• <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)短波
    綠野短波
    欧美日韩福利视频一区二区| 国产精品 欧美亚洲| 日韩免费av在线播放| 精品久久久久久成人av| 久久天躁狠狠躁夜夜2o2o| 国产精品,欧美在线| 悠悠久久av| √禁漫天堂资源中文www| 18美女黄网站色大片免费观看| 天堂动漫精品| 悠悠久久av| 两个人看的免费小视频| 欧美成人一区二区免费高清观看 | 可以在线观看毛片的网站| 满18在线观看网站| 国产真人三级小视频在线观看| 一本综合久久免费| 很黄的视频免费| 91成年电影在线观看| 免费看a级黄色片| 中文字幕高清在线视频| 女警被强在线播放| 国产黄a三级三级三级人| 999精品在线视频| 国产激情欧美一区二区| 欧美中文日本在线观看视频| 99国产精品一区二区蜜桃av| 美女免费视频网站| 欧美 亚洲 国产 日韩一| 日本撒尿小便嘘嘘汇集6| 国产成年人精品一区二区| 大型黄色视频在线免费观看| 脱女人内裤的视频| 日本 欧美在线| 日韩精品免费视频一区二区三区| 亚洲精品在线美女| 一二三四在线观看免费中文在| 国产日韩一区二区三区精品不卡| 欧美成狂野欧美在线观看| 91精品三级在线观看| 悠悠久久av| 看免费av毛片| 色av中文字幕| 给我免费播放毛片高清在线观看| 久久久久久国产a免费观看| 欧美老熟妇乱子伦牲交| av欧美777| 777久久人妻少妇嫩草av网站| 国产精品一区二区免费欧美| 亚洲欧美激情综合另类| 免费人成视频x8x8入口观看| 亚洲熟妇熟女久久| av视频免费观看在线观看| 嫩草影视91久久| 亚洲精品一卡2卡三卡4卡5卡| 日日夜夜操网爽| 国产激情久久老熟女| 欧美色视频一区免费| 久久精品人人爽人人爽视色| 日韩精品免费视频一区二区三区| 国产亚洲av高清不卡| 国产成+人综合+亚洲专区| 色综合欧美亚洲国产小说| 欧美色视频一区免费| 满18在线观看网站| 美女高潮到喷水免费观看| 丰满的人妻完整版| 日本五十路高清| 国产一区二区三区视频了| av免费在线观看网站| www日本在线高清视频| 9色porny在线观看| 首页视频小说图片口味搜索| 久久性视频一级片| 亚洲无线在线观看| 中文字幕av电影在线播放| 国产区一区二久久| 成人亚洲精品一区在线观看| 人人妻人人爽人人添夜夜欢视频| 欧美激情久久久久久爽电影 | 精品少妇一区二区三区视频日本电影| 99精品在免费线老司机午夜| 国产成人欧美在线观看| 丝袜美腿诱惑在线| 欧美午夜高清在线| 波多野结衣av一区二区av| 老司机在亚洲福利影院| 亚洲自偷自拍图片 自拍| 亚洲第一电影网av| 久久午夜亚洲精品久久| videosex国产| 男女下面进入的视频免费午夜 | 欧美日韩福利视频一区二区| 免费观看人在逋| 国产91精品成人一区二区三区| 国产男靠女视频免费网站| 又紧又爽又黄一区二区| 99国产精品99久久久久| 亚洲av电影不卡..在线观看| 国产成人一区二区三区免费视频网站| 黄片大片在线免费观看| 免费女性裸体啪啪无遮挡网站| 一二三四社区在线视频社区8| 亚洲国产精品sss在线观看| 日本黄色视频三级网站网址| 97人妻天天添夜夜摸| 超碰成人久久| 国产激情欧美一区二区| 国产成人av激情在线播放| 午夜久久久在线观看| 国产成人欧美在线观看| 老司机午夜十八禁免费视频| 18禁裸乳无遮挡免费网站照片 | 少妇熟女aⅴ在线视频| 男男h啪啪无遮挡| 精品国产乱子伦一区二区三区| 国产三级黄色录像| 人人妻人人澡欧美一区二区 | 成人亚洲精品一区在线观看| 欧美老熟妇乱子伦牲交| 一a级毛片在线观看| 国产男靠女视频免费网站| 国产精品久久久久久亚洲av鲁大| 妹子高潮喷水视频| 国产亚洲av高清不卡| 色老头精品视频在线观看| 黄色毛片三级朝国网站| 欧美另类亚洲清纯唯美| 成人亚洲精品av一区二区| 国产激情久久老熟女| 女警被强在线播放| bbb黄色大片| 午夜精品久久久久久毛片777| 亚洲五月婷婷丁香| 51午夜福利影视在线观看| 色综合站精品国产| 免费高清视频大片| 国产不卡一卡二| 国产三级在线视频| 久久欧美精品欧美久久欧美| 成人特级黄色片久久久久久久| 侵犯人妻中文字幕一二三四区| 91精品国产国语对白视频| 黄色a级毛片大全视频| av天堂久久9| 久久伊人香网站| 日本vs欧美在线观看视频| 午夜精品久久久久久毛片777| 麻豆av在线久日| 国内毛片毛片毛片毛片毛片| 国产亚洲欧美在线一区二区| 美女高潮到喷水免费观看| 在线免费观看的www视频| 一级a爱片免费观看的视频| 精品电影一区二区在线| 一区二区三区激情视频| 亚洲精品在线美女| 精品少妇一区二区三区视频日本电影| 久久久久国产精品人妻aⅴ院| 国产熟女xx| 一区二区三区国产精品乱码| 日日爽夜夜爽网站| 亚洲在线自拍视频| 两个人看的免费小视频| 黑人欧美特级aaaaaa片| 欧美大码av| 国产精品永久免费网站| 可以在线观看的亚洲视频| 国产精品久久视频播放| 国产成人精品无人区| 国产亚洲欧美在线一区二区| 国产av在哪里看| 99久久久亚洲精品蜜臀av| 超碰成人久久| 高清在线国产一区| 操出白浆在线播放| 欧美成狂野欧美在线观看| 久久精品成人免费网站| 免费高清在线观看日韩| 亚洲中文av在线| 黄片小视频在线播放| 1024香蕉在线观看| 免费av毛片视频| av在线播放免费不卡| 日韩欧美一区视频在线观看| 精品久久久久久,| 神马国产精品三级电影在线观看 | 91麻豆精品激情在线观看国产| 在线永久观看黄色视频| 国产91精品成人一区二区三区| 国产av一区二区精品久久| 国产精品久久久人人做人人爽| 国产av一区二区精品久久| 精品人妻在线不人妻| 国产av一区二区精品久久| 国产精品久久久久久精品电影 | 午夜免费成人在线视频| av天堂在线播放| 一级a爱片免费观看的视频| 精品一品国产午夜福利视频| 岛国视频午夜一区免费看| 99re在线观看精品视频| 国产不卡一卡二| 丰满的人妻完整版| 国产精品免费一区二区三区在线| 久久婷婷人人爽人人干人人爱 | 香蕉国产在线看| 欧美日韩一级在线毛片| 91字幕亚洲| 日韩欧美国产一区二区入口| 一级a爱片免费观看的视频| 久久亚洲精品不卡| 亚洲,欧美精品.| 99在线人妻在线中文字幕| 天天添夜夜摸| 国产亚洲精品久久久久5区| 亚洲精品中文字幕一二三四区| 1024视频免费在线观看| 国产色视频综合| 一区在线观看完整版| 在线观看日韩欧美| 无人区码免费观看不卡| 国产区一区二久久| 亚洲国产精品久久男人天堂| 欧美久久黑人一区二区| 日日夜夜操网爽| 国产熟女xx| 国产精品99久久99久久久不卡| 国产一级毛片七仙女欲春2 | 999久久久精品免费观看国产| 国产91精品成人一区二区三区| 在线观看免费日韩欧美大片| 精品一区二区三区四区五区乱码| 精品午夜福利视频在线观看一区| 丰满人妻熟妇乱又伦精品不卡| 窝窝影院91人妻| 老司机靠b影院| 精品久久久久久久毛片微露脸| 97超级碰碰碰精品色视频在线观看| 午夜成年电影在线免费观看| 真人做人爱边吃奶动态| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品一区二区在线不卡| 窝窝影院91人妻| 国产精品乱码一区二三区的特点 | 欧美日韩乱码在线| 一本久久中文字幕| 成人国语在线视频| 两个人看的免费小视频| 香蕉久久夜色| 亚洲欧美激情在线| 久久国产乱子伦精品免费另类| 亚洲第一青青草原| 国产亚洲欧美在线一区二区| 免费搜索国产男女视频| 两性夫妻黄色片| 老司机午夜十八禁免费视频| 久久午夜综合久久蜜桃| 成人欧美大片| 午夜激情av网站| 久久久久久人人人人人| 国产不卡一卡二| 欧美日韩中文字幕国产精品一区二区三区 | 精品福利观看| 欧美日韩福利视频一区二区| 亚洲男人天堂网一区| 久久精品亚洲熟妇少妇任你| 国产高清激情床上av| 久久久久久久精品吃奶| 亚洲国产欧美网| 精品久久久久久,| 免费搜索国产男女视频| 亚洲欧美日韩高清在线视频| 午夜a级毛片| 午夜福利视频1000在线观看 | 国产男靠女视频免费网站| 亚洲熟妇中文字幕五十中出| 成人国语在线视频| 天天一区二区日本电影三级 | 欧美日韩福利视频一区二区| 欧美乱码精品一区二区三区| 老汉色av国产亚洲站长工具| 亚洲成国产人片在线观看| 久久久久国产精品人妻aⅴ院| av在线播放免费不卡| 别揉我奶头~嗯~啊~动态视频| 窝窝影院91人妻| 无遮挡黄片免费观看| 人成视频在线观看免费观看| 啦啦啦观看免费观看视频高清 | 97超级碰碰碰精品色视频在线观看| 欧美乱色亚洲激情| 欧美午夜高清在线| 婷婷丁香在线五月| 亚洲中文字幕一区二区三区有码在线看 | 精品免费久久久久久久清纯| 久久久久国产一级毛片高清牌| 国产精品二区激情视频| 人人妻人人澡人人看| 国产精品久久久久久亚洲av鲁大| 九色国产91popny在线| 成人特级黄色片久久久久久久| 亚洲性夜色夜夜综合| 国产精品亚洲一级av第二区| 久久午夜综合久久蜜桃| 亚洲国产精品成人综合色| 搡老岳熟女国产| 超碰成人久久| tocl精华| 国产一区二区在线av高清观看| 日本五十路高清| 国产精品久久久久久精品电影 | a在线观看视频网站| 免费在线观看日本一区| 美女午夜性视频免费| 给我免费播放毛片高清在线观看| 中文亚洲av片在线观看爽| 琪琪午夜伦伦电影理论片6080| www国产在线视频色| 国产99白浆流出| av在线天堂中文字幕| 丁香六月欧美| 久久精品国产亚洲av香蕉五月| 在线观看日韩欧美| 777久久人妻少妇嫩草av网站| 亚洲熟妇熟女久久| 欧美中文综合在线视频| 日本免费一区二区三区高清不卡 | 午夜激情av网站| 国产欧美日韩精品亚洲av| 91麻豆av在线| 国产一区在线观看成人免费| 国产精品美女特级片免费视频播放器 | 亚洲 国产 在线| 高清黄色对白视频在线免费看| 淫妇啪啪啪对白视频| 亚洲欧美精品综合久久99| 亚洲中文字幕日韩| 88av欧美| 日韩欧美一区视频在线观看| 99久久99久久久精品蜜桃| 在线国产一区二区在线| 91九色精品人成在线观看| 精品卡一卡二卡四卡免费| 老司机深夜福利视频在线观看| 欧美日本视频| 亚洲欧美精品综合一区二区三区| 成年版毛片免费区| 777久久人妻少妇嫩草av网站| 又紧又爽又黄一区二区| 韩国精品一区二区三区| 国产私拍福利视频在线观看| 亚洲国产精品成人综合色| 18禁黄网站禁片午夜丰满| 日本在线视频免费播放| 国产成人免费无遮挡视频| 国产精品 欧美亚洲| 最新在线观看一区二区三区| 久久久精品国产亚洲av高清涩受| 亚洲精品美女久久久久99蜜臀| 男女午夜视频在线观看| 国产精品98久久久久久宅男小说| 久久青草综合色| 久久午夜亚洲精品久久| 亚洲 欧美 日韩 在线 免费| 欧美在线一区亚洲| 免费观看精品视频网站| 波多野结衣av一区二区av| 女警被强在线播放| 亚洲aⅴ乱码一区二区在线播放 | 又大又爽又粗| 国产欧美日韩精品亚洲av| 日韩国内少妇激情av| 国产精品野战在线观看| 精品福利观看| 此物有八面人人有两片| 激情视频va一区二区三区| 嫩草影视91久久| 日韩精品免费视频一区二区三区| 在线观看免费视频网站a站| 久久精品国产亚洲av香蕉五月| 国产精品久久久av美女十八| 国产在线观看jvid| 欧美日本亚洲视频在线播放| av在线播放免费不卡| 嫩草影院精品99| 欧美成人一区二区免费高清观看 | 久久精品国产清高在天天线| 精品免费久久久久久久清纯| 后天国语完整版免费观看| www国产在线视频色| 精品国产国语对白av| 女同久久另类99精品国产91| 丝袜在线中文字幕| 岛国在线观看网站| 老司机深夜福利视频在线观看| 色播在线永久视频| www国产在线视频色| √禁漫天堂资源中文www| 人人澡人人妻人| 我的亚洲天堂| 在线观看舔阴道视频| 九色亚洲精品在线播放| 国产99白浆流出| 最好的美女福利视频网| 亚洲精品美女久久av网站| 日本免费a在线| 69精品国产乱码久久久| 久久久国产成人精品二区| 成人欧美大片| 欧美精品亚洲一区二区| 91成人精品电影| 在线天堂中文资源库| 欧美成人性av电影在线观看| 亚洲 欧美一区二区三区| 美女高潮喷水抽搐中文字幕| 欧美中文综合在线视频| 免费少妇av软件| 精品无人区乱码1区二区| 国产亚洲精品久久久久久毛片| 18美女黄网站色大片免费观看| 亚洲自偷自拍图片 自拍| 成人免费观看视频高清| 中文字幕久久专区| 午夜免费观看网址| 国产亚洲欧美精品永久| 亚洲自拍偷在线| 久久这里只有精品19| 97碰自拍视频| 欧美日韩一级在线毛片| 欧美亚洲日本最大视频资源| 叶爱在线成人免费视频播放| 国产日韩一区二区三区精品不卡| 一a级毛片在线观看| 法律面前人人平等表现在哪些方面| 亚洲少妇的诱惑av| 国产欧美日韩综合在线一区二区| 国产男靠女视频免费网站| 日韩有码中文字幕| 亚洲精品中文字幕在线视频| 女性被躁到高潮视频| 午夜激情av网站| av欧美777| 中亚洲国语对白在线视频| 国产单亲对白刺激| 99久久99久久久精品蜜桃| 麻豆久久精品国产亚洲av| 日韩av在线大香蕉| 欧美日韩精品网址| 午夜免费激情av| 久热爱精品视频在线9| 中文字幕另类日韩欧美亚洲嫩草| 咕卡用的链子| 国产精品自产拍在线观看55亚洲| 女人爽到高潮嗷嗷叫在线视频| 久久久久久大精品| 丝袜美足系列| 久久人妻熟女aⅴ| 国产伦一二天堂av在线观看| 黄色视频,在线免费观看| 亚洲最大成人中文| 搡老妇女老女人老熟妇| 国产精品乱码一区二三区的特点 | 国产成人免费无遮挡视频| 在线天堂中文资源库| 我的亚洲天堂| 淫秽高清视频在线观看| 亚洲国产欧美网| 色综合站精品国产| 无限看片的www在线观看| 丝袜美腿诱惑在线| netflix在线观看网站| 国产精品精品国产色婷婷| 欧美久久黑人一区二区| 91精品三级在线观看| 黄片播放在线免费| 岛国视频午夜一区免费看| 9色porny在线观看| 色综合站精品国产| 亚洲中文日韩欧美视频| 国产亚洲欧美在线一区二区| 中文字幕精品免费在线观看视频| а√天堂www在线а√下载| 日韩国内少妇激情av| 俄罗斯特黄特色一大片| 国产精品一区二区精品视频观看| 一级片免费观看大全| 久久精品亚洲熟妇少妇任你| 日韩高清综合在线| 动漫黄色视频在线观看| 在线观看免费午夜福利视频| 日韩免费av在线播放| 一进一出抽搐动态| 亚洲全国av大片| 国产单亲对白刺激| 免费少妇av软件| 成年女人毛片免费观看观看9| 丰满的人妻完整版| 国产片内射在线| 午夜免费观看网址| 99国产精品99久久久久| 亚洲精品国产精品久久久不卡| 一级毛片女人18水好多| 欧美日韩亚洲综合一区二区三区_| 欧美乱码精品一区二区三区| 久久人妻福利社区极品人妻图片| 国语自产精品视频在线第100页| 久久久久久久精品吃奶| 成年版毛片免费区| 久久香蕉激情| 亚洲全国av大片| 一进一出抽搐动态| avwww免费| 亚洲第一电影网av| 1024香蕉在线观看| 精品国产国语对白av| 黄色视频,在线免费观看| 好看av亚洲va欧美ⅴa在| 久久精品91无色码中文字幕| 欧美成人午夜精品| 日韩大码丰满熟妇| 亚洲三区欧美一区| 日本vs欧美在线观看视频| 久久久久亚洲av毛片大全| 国产麻豆成人av免费视频| 精品国产一区二区三区四区第35| 中出人妻视频一区二区| 国产不卡一卡二| 一级毛片女人18水好多| 亚洲av日韩精品久久久久久密| 国产1区2区3区精品| 久久久久久久精品吃奶| 亚洲午夜精品一区,二区,三区| 国产精品久久视频播放| 男人舔女人的私密视频| 非洲黑人性xxxx精品又粗又长| 91成人精品电影| 国产日韩一区二区三区精品不卡| 亚洲成av片中文字幕在线观看| 两个人看的免费小视频| 亚洲熟妇熟女久久| 欧美色欧美亚洲另类二区 | 校园春色视频在线观看| 国产亚洲av高清不卡| 亚洲中文av在线| 国产激情久久老熟女| 精品国产超薄肉色丝袜足j| 长腿黑丝高跟| 久久国产精品影院| 99国产精品免费福利视频| 精品欧美国产一区二区三| 啦啦啦韩国在线观看视频| 日本黄色视频三级网站网址| 久久香蕉国产精品| 国产成人欧美在线观看| 亚洲熟女毛片儿| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频| 激情视频va一区二区三区| 91精品三级在线观看| 国产成人av教育| 人人妻人人澡人人看| 久久久久国产精品人妻aⅴ院| 亚洲av日韩精品久久久久久密| 精品久久久久久久久久免费视频| 国产成人欧美| 18禁美女被吸乳视频| svipshipincom国产片| 日本免费一区二区三区高清不卡 | 美女免费视频网站| 热99re8久久精品国产| 国产高清视频在线播放一区| 久久久久久久久久久久大奶| 精品久久久久久久人妻蜜臀av | 国产免费av片在线观看野外av| 极品人妻少妇av视频| 国产一区二区三区视频了| av超薄肉色丝袜交足视频| 欧美不卡视频在线免费观看 | 激情视频va一区二区三区| 18禁裸乳无遮挡免费网站照片 | 欧美国产日韩亚洲一区| 国产真人三级小视频在线观看| 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 欧美丝袜亚洲另类 | 欧美成狂野欧美在线观看| 午夜福利高清视频| 在线观看免费午夜福利视频| 亚洲少妇的诱惑av| 精品无人区乱码1区二区| 亚洲全国av大片| 欧美黄色淫秽网站| 亚洲精品一区av在线观看| 成熟少妇高潮喷水视频| 午夜福利欧美成人| 亚洲av电影不卡..在线观看| 一区二区三区高清视频在线| 国产欧美日韩一区二区精品| 少妇粗大呻吟视频| 一进一出抽搐动态| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 久久狼人影院| 国产日韩一区二区三区精品不卡| 1024香蕉在线观看| 精品一区二区三区av网在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 法律面前人人平等表现在哪些方面| 日韩欧美一区视频在线观看| 亚洲情色 制服丝袜| 人人妻人人爽人人添夜夜欢视频| 中文字幕高清在线视频| 黄色视频,在线免费观看| 欧美绝顶高潮抽搐喷水|