• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2017-11-06 09:54:45白昕欣黃蕭霖秦孟晟張?jiān)婪?/span>陳留根
    關(guān)鍵詞:反照率短波晴天

    白昕欣, 黃蕭霖, 秦孟晟, 張?jiān)婪? 陳留根, 郭 智, 郝 璐**

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

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

    白昕欣1, 黃蕭霖1, 秦孟晟1, 張?jiān)婪?, 陳留根2, 郭 智2, 郝 璐1**

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

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

    水稻田; 地表反照率; 生長(zhǎng)季; 不對(duì)稱; 氣象因子

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

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

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

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

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

    1 研究方法

    1.1 研究區(qū)概況

    研究站點(diǎn)位于江蘇省農(nóng)業(yè)科學(xué)院溧水試驗(yàn)基地,該基地位于溧水縣白馬鎮(zhèn)(119.2°E, 31.6°N, 海拔高度38 m), 總面積為 0.8 km2, 站點(diǎn)四周均為水稻田, 水源來(lái)自于白馬湖水庫(kù)?;镜孛差愋蜑榍鹆? 土質(zhì)類型為白漿土。基地屬于秦淮河流域, 為亞熱帶季風(fēng)氣候, 雨量充沛, 年降水量為 1 048 mm, 四季分明, 年平均溫度15.4 ℃。根據(jù)溧水站多年氣象觀測(cè)數(shù)據(jù)(1985—2016年, 包括氣溫、降水量以及相對(duì)濕度)分析, 溧水2016年5月和6月的平均氣溫較歷年偏低, 而其他月份的平均氣溫較歷年偏高。2016年水稻生長(zhǎng)季平均氣溫略高于歷年同期值。除8月份外,2016年水稻生長(zhǎng)季內(nèi)各月的相對(duì)濕度和降水量均大于歷年同期值??傮w上, 2016年水稻生長(zhǎng)季屬于偏暖濕年份(表1)。觀測(cè)站點(diǎn)種植的水稻品種為‘南粳9108’, 2016年的插秧時(shí)間為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 研究站點(diǎn)(溧水)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 觀測(cè)儀器與數(shù)據(jù)預(yù)處理

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

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

    2 結(jié)果與分析

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

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

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

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

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

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

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

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

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

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

    圖2 稻田不同時(shí)期典型晴天地表反照率的不對(duì)稱性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)表明, 稻田地表反照率整體上晴天>陰天>雨天。晴天和陰天的地表反照率不對(duì)稱性較為明顯, 均為上午較高, 下午較低; 而雨天的地表反照率較為對(duì)稱, 這可能與晴天和陰天上午地表附著較多的露水有關(guān)。晴天地表接收的太陽(yáng)短波輻射較大,土壤濕度較小, 因此地表反照率較高。受太陽(yáng)高度角的影響, 地表反照率的日變化較為明顯。陰天到達(dá)地表的太陽(yáng)短波輻射較少土壤濕度較大, 地表反照率較小, 反照率的日變化幅度最小。雨天條件下, 到達(dá)地表的太陽(yáng)輻短波射最小, 土壤濕度較大, 反照率也相應(yīng)減小, 地表反照率的日變化較幅度較小。

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

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

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

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

    圖3 稻田不同時(shí)期典型晴天的平均露點(diǎn)溫度與實(shí)際溫度(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 稻田不同時(shí)期典型晴天平均風(fēng)速日變化Fig.4 Daily dynamics of average wind speed of paddy fields in typical sunny days in different growth periods of rice

    圖5 稻田不同時(shí)期典型晴天的風(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)。陰天地表反照率與相對(duì)濕度之間的相關(guān)程度最高。雨天地表反照率與三者之間的相關(guān)程度均比較低。

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

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

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

    圖7 水稻生長(zhǎng)季內(nèi)的稻田地表反照率(Ⅰ為水稻插秧到分蘗期, Ⅱ?yàn)樗痉痔Y期, Ⅲ為水稻拔節(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 水稻生長(zhǎng)季內(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 觀測(cè)期內(nèi)的稻田地表反照率與溫度的變化Fig.9 Changes of surface albedo and temperature in paddy field during the observation period

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

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

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

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

    3 討論

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

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

    3.2 影響地表反照率的因子

    Liu等[25]發(fā)現(xiàn)玉米(Zea mays)農(nóng)田地表反照率與太陽(yáng)高度角和土壤濕度呈負(fù)相關(guān)關(guān)系, 與植被指數(shù)呈正相關(guān)關(guān)系。鄭志遠(yuǎn)等[29]發(fā)現(xiàn)裸土地表反照率隨太陽(yáng)高度角及土壤濕度的增加而減小。張樂(lè)樂(lè)等[30]研究發(fā)現(xiàn), 土壤水分增加, 地表反照率明顯減小,而且這種關(guān)系還受到植被生長(zhǎng)周期的影響。Usowicz等[31]發(fā)現(xiàn)隨著生物炭用量的增加, 地表反照率明顯減小, 并且降雨過(guò)后, 休耕地的地表反照率降低,而草地則相反。李慧芳等[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較小時(shí), 地表反照率的差異較大, 而隨著 LAI的增加, 地表反照率趨于穩(wěn)定。由此可見(jiàn), 地表反照率的大小受到多種因素的共同影響, 如下墊面的顏色、植被覆蓋率、土壤溫度、土壤濕度、太陽(yáng)高度角和大氣透過(guò)率等。

    我們的結(jié)果表明, 典型晴天狀況下, 稻田地表反照率與太陽(yáng)高度角有關(guān), 白天中午地表反照率最低, 上午和下午地表反照率較大, 因此地表反照率隨太陽(yáng)高度角的增加而減小。同時(shí), 晴天地表反照率每日內(nèi)的變化特征還受到當(dāng)天晝夜溫差及風(fēng)的影響, 太陽(yáng)高度角小于60°時(shí), 晝夜溫差導(dǎo)致上午植被表面附著較多露水, 同一太陽(yáng)高度角下上午地表反照率較下午大; 而太陽(yáng)高度角小于60°時(shí), 風(fēng)使得植被葉面傾斜, 同一太陽(yáng)高度角下下午地表反照率較上午大。地表反照率也與天氣狀況有關(guān), 晴天大氣透過(guò)率較大, 地表反照率也較大, 而陰雨天大氣透過(guò)率較小, 地表反照率也較小。地表反照率還與水稻的生長(zhǎng)狀況有關(guān)。水稻生長(zhǎng)初期, 地表反照率隨水稻的生長(zhǎng)而增加, 在水稻接近成熟時(shí), 稻田地表反照率變化較為平緩。當(dāng)水稻成熟時(shí), 地表反照率開(kāi)始下降, 水稻收獲之后, 地表反照率進(jìn)一步降低。對(duì)于土壤濕度及植被覆蓋度等因子對(duì)地表反照率的影響研究, 將是未來(lái)需要開(kāi)展的重點(diǎ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”型分布, 中午較低, 下午和上午較高。隨著水稻的生長(zhǎng), 地表反照率的日變化特征變得更為明顯。晴天稻田地表反照率的不對(duì)稱性主要是由露水和風(fēng)速風(fēng)向引起, 太陽(yáng)高度角較小時(shí), 露水的散射作用使得上午的地表反照率值較下午高, 而太陽(yáng)高度角較大時(shí), 西南風(fēng)促使作物葉面傾斜, 從而使得下午的地表反照率值較上午高。

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

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

    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種水稻趨勢(shì)產(chǎn)量擬合方法的比較分析[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 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年氣候變化對(duì)西南冬小麥潛在和雨養(yǎng)產(chǎn)量影響的模擬分析[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),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]楊娟, 陳洪濱, 王開(kāi)存, 等.利用 MODIS衛(wèi)星資料分析北京地區(qū)地表反照率時(shí)空分布及變化特征[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]張超, 宋海清, 吳國(guó)周, 等.NCAR/CLM系列陸面模式對(duì)內(nèi)蒙古地表溫度的模擬評(píng)估[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 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]陳云浩, 李曉兵, 謝鋒.我國(guó)西北地區(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]蔡福, 祝青林, 何洪林, 等.中國(guó)月平均地表反照率的估算及其時(shí)空分布[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ā)育期反照率變化及其不對(duì)稱性[J].農(nóng)業(yè)工程學(xué)報(bào), 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)云, 童應(yīng)祥, 等.壽縣地區(qū)小麥和水稻田地表反照率觀測(cè)分析[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]王勝, 張強(qiáng), 張鴻.稀疏植被地表反照率及土壤熱傳導(dǎo)率特征研究[J].中國(guó)沙漠, 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]文麗娟, 呂世華, 陳世強(qiáng), 等.干旱區(qū)綠洲地表反照率不對(duì)稱觀測(cè)研究[J].太陽(yáng)能學(xué)報(bào), 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]張?chǎng)? 謝輝, 張平, 等.扁桃||冬小麥間作模式下樹(shù)冠結(jié)構(gòu)對(duì)間作區(qū)域光環(huán)境的影響[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 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]張樂(lè)樂(lè), 趙林, 李韌, 等.青藏高原唐古拉地區(qū)暖季土壤水分對(duì)地表反照率及其土壤熱參數(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

    白昕欣, 黃蕭霖, 秦孟晟, 張?jiān)婪? 陳留根, 郭智, 郝璐.亞熱帶地區(qū)水稻田地表反照率變化特征[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 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

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

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

    白昕欣, 主要研究方向?yàn)閼?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

    猜你喜歡
    反照率短波晴天
    基于藍(lán)天空反照率的氣溶膠輻射強(qiáng)迫模擬
    薩吾爾山木斯島冰川反照率時(shí)空變化特征研究
    冰川凍土(2022年6期)2022-02-12 08:31:06
    它若安好,便是晴天
    長(zhǎng)江三角洲地區(qū)大氣氣溶膠柱單次散射反照率特性研究
    一起來(lái)寫(xiě)“雪”
    晴天有時(shí)下豬
    樂(lè)海短波
    工運(yùn)短波
    工運(yùn)短波
    綠野短波
    久久人妻熟女aⅴ| 精品国产一区二区三区四区第35| 亚洲午夜精品一区,二区,三区| 亚洲欧美精品综合一区二区三区| 国产精品综合久久久久久久免费 | 国产精品久久电影中文字幕| 性欧美人与动物交配| 亚洲avbb在线观看| 亚洲情色 制服丝袜| 黄片播放在线免费| 免费在线观看日本一区| 久久久久久人人人人人| 久久精品影院6| 中亚洲国语对白在线视频| 成年女人毛片免费观看观看9| 亚洲国产欧美网| 亚洲国产欧美日韩在线播放| 久久中文看片网| 久久天躁狠狠躁夜夜2o2o| 人妻久久中文字幕网| 无限看片的www在线观看| 国产深夜福利视频在线观看| 国产精品av久久久久免费| bbb黄色大片| 免费观看精品视频网站| 伊人久久大香线蕉亚洲五| 国产午夜精品久久久久久| 麻豆国产av国片精品| 黑人欧美特级aaaaaa片| 免费在线观看影片大全网站| 五月开心婷婷网| 中文字幕人妻熟女乱码| 一级毛片精品| 人人妻,人人澡人人爽秒播| 国产精品日韩av在线免费观看 | 黄色怎么调成土黄色| 亚洲欧美激情在线| 亚洲熟妇中文字幕五十中出 | 少妇粗大呻吟视频| 亚洲久久久国产精品| 午夜免费激情av| 国产精品 欧美亚洲| 手机成人av网站| 国产不卡一卡二| 最近最新中文字幕大全免费视频| 亚洲一区高清亚洲精品| 老熟妇乱子伦视频在线观看| 欧美老熟妇乱子伦牲交| 桃红色精品国产亚洲av| 脱女人内裤的视频| 精品人妻1区二区| 成人亚洲精品一区在线观看| 亚洲人成伊人成综合网2020| 国产伦一二天堂av在线观看| 日韩精品免费视频一区二区三区| 国产黄色免费在线视频| 国产激情欧美一区二区| av中文乱码字幕在线| 大陆偷拍与自拍| 日本vs欧美在线观看视频| 国产亚洲精品综合一区在线观看 | 色尼玛亚洲综合影院| 777久久人妻少妇嫩草av网站| 国产主播在线观看一区二区| 大香蕉久久成人网| 18禁黄网站禁片午夜丰满| 午夜老司机福利片| 精品欧美一区二区三区在线| 一区二区三区国产精品乱码| 一级黄色大片毛片| 午夜免费成人在线视频| 一区二区日韩欧美中文字幕| 91成年电影在线观看| 黄频高清免费视频| 99国产综合亚洲精品| 欧美日本中文国产一区发布| 成年人免费黄色播放视频| а√天堂www在线а√下载| 最近最新免费中文字幕在线| 成人黄色视频免费在线看| 国产亚洲欧美98| 国产日韩一区二区三区精品不卡| 亚洲国产欧美网| 亚洲 欧美 日韩 在线 免费| 美女大奶头视频| 国产黄a三级三级三级人| 在线看a的网站| 18禁裸乳无遮挡免费网站照片 | 亚洲精品中文字幕在线视频| 国产真人三级小视频在线观看| 黄色视频不卡| 在线av久久热| 国产亚洲av高清不卡| 欧美黑人欧美精品刺激| 极品人妻少妇av视频| 亚洲专区字幕在线| av欧美777| 亚洲人成电影观看| 色播在线永久视频| 午夜福利一区二区在线看| a在线观看视频网站| 亚洲情色 制服丝袜| 身体一侧抽搐| 欧美日韩乱码在线| 成人影院久久| 亚洲精品在线美女| 亚洲精品粉嫩美女一区| 国产成人一区二区三区免费视频网站| 成年版毛片免费区| 岛国视频午夜一区免费看| 国产欧美日韩一区二区精品| 午夜福利影视在线免费观看| 午夜免费观看网址| 99久久精品国产亚洲精品| 欧美日本中文国产一区发布| 性欧美人与动物交配| 手机成人av网站| 亚洲成人国产一区在线观看| 亚洲欧美激情综合另类| av免费在线观看网站| 欧美黄色片欧美黄色片| 最新在线观看一区二区三区| a级毛片黄视频| 久久亚洲精品不卡| 国产欧美日韩综合在线一区二区| 美女国产高潮福利片在线看| 亚洲成av片中文字幕在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 99国产精品一区二区三区| av视频免费观看在线观看| 中文欧美无线码| 一级毛片精品| 国产高清激情床上av| 高清毛片免费观看视频网站 | 免费av毛片视频| 国产精品成人在线| 成人精品一区二区免费| 午夜福利一区二区在线看| 国产av又大| 丝袜美足系列| 午夜福利一区二区在线看| 波多野结衣一区麻豆| 看黄色毛片网站| 亚洲国产精品999在线| а√天堂www在线а√下载| 免费av中文字幕在线| 身体一侧抽搐| 欧美在线一区亚洲| 久久国产亚洲av麻豆专区| 亚洲精品国产区一区二| 国产午夜精品久久久久久| 一级作爱视频免费观看| 极品人妻少妇av视频| 国产欧美日韩一区二区三区在线| 亚洲国产精品sss在线观看 | 一进一出抽搐gif免费好疼 | 啪啪无遮挡十八禁网站| 操美女的视频在线观看| 亚洲欧美日韩高清在线视频| xxx96com| 国产乱人伦免费视频| 99久久久亚洲精品蜜臀av| 欧美激情 高清一区二区三区| 超色免费av| 精品国产乱子伦一区二区三区| 精品国产亚洲在线| 一级a爱片免费观看的视频| 嫩草影院精品99| 麻豆国产av国片精品| 国产成人啪精品午夜网站| 亚洲av成人一区二区三| 欧美日韩亚洲高清精品| 国产91精品成人一区二区三区| 电影成人av| 精品久久久久久电影网| 99久久精品国产亚洲精品| 两个人看的免费小视频| 亚洲精品av麻豆狂野| 欧美精品啪啪一区二区三区| 999久久久精品免费观看国产| 国产野战对白在线观看| 亚洲精品中文字幕一二三四区| 久久国产精品人妻蜜桃| 啦啦啦免费观看视频1| 欧美乱妇无乱码| 日韩免费高清中文字幕av| 亚洲成人国产一区在线观看| 97超级碰碰碰精品色视频在线观看| 一级a爱视频在线免费观看| 久久香蕉国产精品| 国产精品一区二区在线不卡| 精品久久久久久,| 男人的好看免费观看在线视频 | 韩国精品一区二区三区| 午夜免费鲁丝| 午夜日韩欧美国产| 欧美一区二区精品小视频在线| 日本欧美视频一区| 中文字幕人妻丝袜制服| 女人被狂操c到高潮| 免费一级毛片在线播放高清视频 | 动漫黄色视频在线观看| 亚洲,欧美精品.| 国产黄色免费在线视频| 国产精品一区二区在线不卡| 亚洲精品成人av观看孕妇| 18禁国产床啪视频网站| 在线观看免费日韩欧美大片| 高清在线国产一区| 欧美日韩av久久| 色哟哟哟哟哟哟| 亚洲男人天堂网一区| 嫩草影院精品99| 午夜老司机福利片| 欧美亚洲日本最大视频资源| 视频在线观看一区二区三区| 亚洲专区字幕在线| 成人亚洲精品一区在线观看| 人人妻人人添人人爽欧美一区卜| 国产蜜桃级精品一区二区三区| 很黄的视频免费| 亚洲五月色婷婷综合| 成在线人永久免费视频| 久久青草综合色| 一区福利在线观看| 十分钟在线观看高清视频www| 999久久久精品免费观看国产| 精品免费久久久久久久清纯| 亚洲av片天天在线观看| 亚洲精品美女久久久久99蜜臀| 中文字幕色久视频| 人人妻,人人澡人人爽秒播| 久久久久九九精品影院| 欧美日韩精品网址| 免费观看人在逋| xxxhd国产人妻xxx| 久久中文看片网| 久久草成人影院| 少妇粗大呻吟视频| 一边摸一边抽搐一进一小说| 欧美激情久久久久久爽电影 | 亚洲欧美激情综合另类| 欧美一区二区精品小视频在线| 日韩 欧美 亚洲 中文字幕| 操出白浆在线播放| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 看免费av毛片| 日韩欧美三级三区| 伦理电影免费视频| 亚洲精品久久成人aⅴ小说| 国产成人av教育| 侵犯人妻中文字幕一二三四区| 亚洲性夜色夜夜综合| 国产日韩一区二区三区精品不卡| 一级a爱片免费观看的视频| 免费在线观看黄色视频的| 精品少妇一区二区三区视频日本电影| 狠狠狠狠99中文字幕| 91精品国产国语对白视频| av视频免费观看在线观看| 成人国语在线视频| 国产av在哪里看| 麻豆国产av国片精品| 黄色成人免费大全| 亚洲一码二码三码区别大吗| 亚洲午夜理论影院| 日本黄色日本黄色录像| 久久久久久久久免费视频了| 国产精品久久久人人做人人爽| 国产区一区二久久| 乱人伦中国视频| 欧美日韩视频精品一区| 久久精品91蜜桃| 亚洲黑人精品在线| 不卡av一区二区三区| 久久精品人人爽人人爽视色| 国产精品野战在线观看 | 国产亚洲精品综合一区在线观看 | 色综合欧美亚洲国产小说| 好看av亚洲va欧美ⅴa在| 一级毛片女人18水好多| 久久欧美精品欧美久久欧美| 久久人妻福利社区极品人妻图片| 精品久久久久久成人av| 美国免费a级毛片| 99国产极品粉嫩在线观看| 国产精品日韩av在线免费观看 | 久久久久国内视频| 黄色怎么调成土黄色| 久久久久亚洲av毛片大全| 亚洲七黄色美女视频| 成人18禁在线播放| 亚洲国产欧美日韩在线播放| 99精国产麻豆久久婷婷| 亚洲avbb在线观看| 一区二区三区激情视频| 国产一区二区在线av高清观看| 夜夜夜夜夜久久久久| 变态另类成人亚洲欧美熟女 | 黄色视频,在线免费观看| 午夜久久久在线观看| 99久久久亚洲精品蜜臀av| 99国产综合亚洲精品| www.精华液| 亚洲国产欧美日韩在线播放| 欧美色视频一区免费| 午夜两性在线视频| 在线观看免费日韩欧美大片| 丰满饥渴人妻一区二区三| av超薄肉色丝袜交足视频| 国产免费男女视频| 18禁美女被吸乳视频| 咕卡用的链子| 精品国产乱码久久久久久男人| 看黄色毛片网站| 黄色视频不卡| 真人做人爱边吃奶动态| 久久精品人人爽人人爽视色| 亚洲国产欧美网| 琪琪午夜伦伦电影理论片6080| 夜夜躁狠狠躁天天躁| 免费av毛片视频| 日韩欧美在线二视频| 国产精品久久久av美女十八| 亚洲欧美精品综合久久99| 亚洲美女黄片视频| 最新在线观看一区二区三区| 色在线成人网| 女人被狂操c到高潮| 国产av在哪里看| 热99re8久久精品国产| 12—13女人毛片做爰片一| 如日韩欧美国产精品一区二区三区| 变态另类成人亚洲欧美熟女 | 成人影院久久| 亚洲第一青青草原| 国产成人av激情在线播放| a级片在线免费高清观看视频| 日本三级黄在线观看| 男女下面进入的视频免费午夜 | 女性生殖器流出的白浆| 免费av中文字幕在线| 亚洲av第一区精品v没综合| 日本撒尿小便嘘嘘汇集6| 欧美乱码精品一区二区三区| 亚洲国产看品久久| 欧美色视频一区免费| 日韩免费高清中文字幕av| 一进一出好大好爽视频| 久久香蕉激情| 久久久久国产精品人妻aⅴ院| 欧美老熟妇乱子伦牲交| 亚洲免费av在线视频| 无限看片的www在线观看| 日日摸夜夜添夜夜添小说| 男女午夜视频在线观看| 久久人妻福利社区极品人妻图片| 亚洲精品国产色婷婷电影| 自线自在国产av| 18美女黄网站色大片免费观看| 久久 成人 亚洲| 久久亚洲精品不卡| 精品国产一区二区久久| 国产亚洲av高清不卡| 国产精品自产拍在线观看55亚洲| 成人18禁在线播放| 国产男靠女视频免费网站| 丝袜人妻中文字幕| 日韩免费高清中文字幕av| 99国产精品免费福利视频| 婷婷丁香在线五月| 欧美亚洲日本最大视频资源| 99国产精品一区二区蜜桃av| xxxhd国产人妻xxx| 久久狼人影院| 制服诱惑二区| 水蜜桃什么品种好| 欧美日韩亚洲综合一区二区三区_| 黑人操中国人逼视频| 99精国产麻豆久久婷婷| av电影中文网址| 免费女性裸体啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av高清一级| 色老头精品视频在线观看| 久久久久久久精品吃奶| 亚洲五月色婷婷综合| 高清毛片免费观看视频网站 | 国产精品一区二区三区四区久久 | 无人区码免费观看不卡| 亚洲美女黄片视频| 精品久久久久久成人av| 宅男免费午夜| 日日夜夜操网爽| 国产精品日韩av在线免费观看 | 国产亚洲精品第一综合不卡| 十分钟在线观看高清视频www| 人人妻,人人澡人人爽秒播| 精品久久久久久电影网| 日韩成人在线观看一区二区三区| 校园春色视频在线观看| 欧美午夜高清在线| 国产99久久九九免费精品| 国产午夜精品久久久久久| 亚洲中文日韩欧美视频| 香蕉丝袜av| 色婷婷久久久亚洲欧美| 在线观看午夜福利视频| 亚洲欧美激情在线| 在线观看免费午夜福利视频| 国产男靠女视频免费网站| 国产单亲对白刺激| 精品人妻在线不人妻| 在线播放国产精品三级| 91精品三级在线观看| 久久国产亚洲av麻豆专区| 成在线人永久免费视频| 美女国产高潮福利片在线看| 亚洲五月婷婷丁香| 好男人电影高清在线观看| 两人在一起打扑克的视频| 一级a爱片免费观看的视频| 欧美中文综合在线视频| 亚洲男人天堂网一区| 一级毛片精品| 久久人妻熟女aⅴ| 色老头精品视频在线观看| 亚洲一码二码三码区别大吗| 亚洲熟女毛片儿| 一级作爱视频免费观看| 色在线成人网| 成人18禁在线播放| 在线天堂中文资源库| 99久久99久久久精品蜜桃| 国产亚洲精品久久久久5区| 曰老女人黄片| 青草久久国产| 嫩草影视91久久| 国产精品亚洲av一区麻豆| av欧美777| 欧美久久黑人一区二区| 淫秽高清视频在线观看| 女性被躁到高潮视频| 日韩大码丰满熟妇| 大型黄色视频在线免费观看| 大码成人一级视频| 日韩欧美三级三区| 看黄色毛片网站| 日本a在线网址| 又紧又爽又黄一区二区| 操出白浆在线播放| 亚洲第一av免费看| 欧美久久黑人一区二区| 久久久精品欧美日韩精品| 又黄又爽又免费观看的视频| 两性夫妻黄色片| 精品电影一区二区在线| www国产在线视频色| 少妇被粗大的猛进出69影院| 中文字幕av电影在线播放| 久久人妻福利社区极品人妻图片| 久久99一区二区三区| 国产亚洲精品一区二区www| 男女下面插进去视频免费观看| 在线观看免费视频日本深夜| 日韩有码中文字幕| 欧美在线黄色| 午夜免费鲁丝| 欧美日韩乱码在线| 黄网站色视频无遮挡免费观看| 一进一出抽搐动态| 国产不卡一卡二| 久久久久久久久免费视频了| 极品教师在线免费播放| 一级黄色大片毛片| 视频区图区小说| 精品一区二区三卡| 高潮久久久久久久久久久不卡| 国产精品自产拍在线观看55亚洲| 黄色成人免费大全| 99热国产这里只有精品6| 丁香欧美五月| 人人澡人人妻人| 极品教师在线免费播放| 亚洲欧美日韩另类电影网站| 操出白浆在线播放| 一a级毛片在线观看| 淫妇啪啪啪对白视频| 极品教师在线免费播放| 亚洲精品成人av观看孕妇| 国产国语露脸激情在线看| 国产真人三级小视频在线观看| xxx96com| 精品国产乱码久久久久久男人| 又紧又爽又黄一区二区| 午夜老司机福利片| 久久精品国产亚洲av香蕉五月| 久久午夜综合久久蜜桃| 91成年电影在线观看| 美女国产高潮福利片在线看| 69精品国产乱码久久久| 美女午夜性视频免费| 夜夜躁狠狠躁天天躁| 国产一区二区三区综合在线观看| 脱女人内裤的视频| 女生性感内裤真人,穿戴方法视频| 免费一级毛片在线播放高清视频 | 视频在线观看一区二区三区| 夜夜看夜夜爽夜夜摸 | av电影中文网址| 男女之事视频高清在线观看| 又黄又爽又免费观看的视频| 麻豆av在线久日| 精品一区二区三区视频在线观看免费 | 欧美精品一区二区免费开放| 亚洲成人免费av在线播放| 制服人妻中文乱码| 亚洲狠狠婷婷综合久久图片| 国产精品免费视频内射| 18禁黄网站禁片午夜丰满| 午夜a级毛片| 可以在线观看毛片的网站| 精品免费久久久久久久清纯| 999精品在线视频| 91麻豆av在线| 国产午夜精品久久久久久| 国产精品国产高清国产av| 别揉我奶头~嗯~啊~动态视频| 欧美在线黄色| 亚洲精品av麻豆狂野| 淫秽高清视频在线观看| 欧美日韩国产mv在线观看视频| 日韩欧美免费精品| 一区二区三区国产精品乱码| 国产亚洲精品久久久久久毛片| 最近最新中文字幕大全免费视频| 高清欧美精品videossex| 亚洲人成77777在线视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲性夜色夜夜综合| 又黄又粗又硬又大视频| 少妇裸体淫交视频免费看高清 | 色播在线永久视频| 免费女性裸体啪啪无遮挡网站| 岛国在线观看网站| 亚洲五月天丁香| 在线观看免费视频日本深夜| 久久九九热精品免费| 每晚都被弄得嗷嗷叫到高潮| 国产成人av教育| 亚洲成人精品中文字幕电影 | 黄网站色视频无遮挡免费观看| 如日韩欧美国产精品一区二区三区| 一进一出抽搐gif免费好疼 | 久久久国产欧美日韩av| 国产成人av激情在线播放| 亚洲七黄色美女视频| 国产精品综合久久久久久久免费 | 免费一级毛片在线播放高清视频 | 9色porny在线观看| 欧美日韩乱码在线| 午夜视频精品福利| 99久久精品国产亚洲精品| 免费少妇av软件| 国产人伦9x9x在线观看| 日韩三级视频一区二区三区| 久久国产亚洲av麻豆专区| 少妇裸体淫交视频免费看高清 | 精品福利观看| 黄片大片在线免费观看| 老熟妇仑乱视频hdxx| aaaaa片日本免费| av在线播放免费不卡| 久久国产精品人妻蜜桃| 女生性感内裤真人,穿戴方法视频| 香蕉久久夜色| 电影成人av| 欧美日韩黄片免| 日韩一卡2卡3卡4卡2021年| 搡老岳熟女国产| 成人影院久久| 国产一区二区激情短视频| 又黄又粗又硬又大视频| 日韩人妻精品一区2区三区| 夜夜看夜夜爽夜夜摸 | 欧美乱码精品一区二区三区| 国产成人精品无人区| 日韩中文字幕欧美一区二区| 国产单亲对白刺激| 国产精品日韩av在线免费观看 | 天堂动漫精品| 777久久人妻少妇嫩草av网站| 两人在一起打扑克的视频| 长腿黑丝高跟| 久久国产精品影院| e午夜精品久久久久久久| 欧美黑人欧美精品刺激| 激情在线观看视频在线高清| 日韩一卡2卡3卡4卡2021年| 亚洲久久久国产精品| 水蜜桃什么品种好| 在线av久久热| 五月开心婷婷网| 日韩中文字幕欧美一区二区| 国产亚洲欧美精品永久| 国产精品自产拍在线观看55亚洲| 亚洲精品中文字幕在线视频| 欧美中文综合在线视频| 国产有黄有色有爽视频| 三上悠亚av全集在线观看| 欧美大码av| 久久久久久久精品吃奶| 免费看a级黄色片| 日日夜夜操网爽|