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

    栓皮櫟幼苗葉氣溫差隨太陽輻射和風(fēng)速變化的列線研究*

    2016-05-27 02:58:42王佩舒王威紅陳景玲李國懷河南農(nóng)業(yè)大學(xué)林學(xué)院鄭州450002
    中國農(nóng)業(yè)氣象 2016年2期
    關(guān)鍵詞:栓皮櫟太陽輻射風(fēng)速

    王佩舒,王威紅,陳景玲,王 謙,李國懷(河南農(nóng)業(yè)大學(xué)林學(xué)院,鄭州 450002)

    ?

    栓皮櫟幼苗葉氣溫差隨太陽輻射和風(fēng)速變化的列線研究*

    王佩舒,王威紅,陳景玲**,王謙,李國懷
    (河南農(nóng)業(yè)大學(xué)林學(xué)院,鄭州 450002)

    摘要:為了探索太陽輻射、氣溫和風(fēng)速對(duì)栓皮櫟幼苗葉溫(LT)的影響,以及LT與Gates氣候空間列線的一致性,在晴天條件下,控制小氣候風(fēng)速因子,測(cè)定氣溫、LT和太陽輻射,研究不同風(fēng)速條件下栓皮櫟幼苗葉氣溫差(LATD)隨太陽輻射的變化。并依據(jù)Gates氣候空間理論和圖解方法,確定栓皮櫟幼苗的氣候空間列線。結(jié)果表明:靜風(fēng)時(shí)栓皮櫟幼苗LATD隨太陽輻射的增加而增加,二者呈線性關(guān)系;冠層溫度衰減(CTD)隨太陽輻射強(qiáng)度增加而增加,但隨著風(fēng)速的增大,CTD隨太陽輻射而增加的幅度逐漸減小。風(fēng)速為0.8m·s?1時(shí),CTD隨太陽輻射變化的趨勢(shì)線斜率逐漸趨近于0,表明不同風(fēng)速條件下的栓皮櫟幼苗葉片CTD隨著太陽輻射強(qiáng)度的增加呈線性增加趨勢(shì);隨著風(fēng)速的增加,該線性關(guān)系形成斜率逐漸變小的列線,列線與Gates氣候空間理論的上限趨勢(shì)一致,列線上的點(diǎn)與生物學(xué)上限溫度(BMT)相加可以表達(dá)栓皮櫟幼苗葉片氣候空間上限。

    關(guān)鍵詞:栓皮櫟;CTD;太陽輻射;風(fēng)速;氣候空間;列線

    王佩舒,王威紅,陳景玲,等.栓皮櫟幼苗葉氣溫差隨太陽輻射和風(fēng)速變化的列線研究[J].中國農(nóng)業(yè)氣象,2016,37(2):213-221

    干旱和半干旱地區(qū)栓皮櫟實(shí)生苗造林的1-5齡幼苗死苗現(xiàn)象較嚴(yán)重[1-2],存活曲線為DeeveyⅢ型,幼齡期死亡率很高[3]。死苗直接影響造林,但現(xiàn)有研究對(duì)死苗現(xiàn)象關(guān)注較多[1-2],而針對(duì)死苗原因研究較4少。高溫干旱是該氣候區(qū)的主要?dú)夂蛱攸c(diǎn),因而高溫干旱脅迫被認(rèn)為是死苗的主要原因,如華北山區(qū)自然生態(tài)環(huán)境下易發(fā)生長期干旱,栓皮櫟幼苗在這樣的生存空間上遇到小氣候高溫,易導(dǎo)致死苗發(fā)生。另一方面,以往研究針對(duì)單一小氣候因素對(duì)植物的影響較多[4-8],多因素復(fù)合影響的研究相對(duì)較少[9]。而葉溫(Leaf temperature, LT)與小氣候中的太陽輻射、氣溫、風(fēng)速等因素均密切相關(guān)[10-14],受多個(gè)小氣候因子的影響。氣溫和太陽輻射是栓皮櫟幼苗葉片溫度形成的主要因子;土壤濕度影響葉片蒸騰(Transpiration, Tr),潛熱交換作用對(duì)LT有一定影響[15-16]。加之以往研究少有考慮風(fēng)的作用,已有研究報(bào)道關(guān)注的是風(fēng)速影響Tr速率[17]和栓皮櫟植株形態(tài)建成[18]。事實(shí)上,風(fēng)對(duì)栓皮櫟的影響也表現(xiàn)在增加葉片表面的熱量平衡方程中對(duì)流換熱項(xiàng)的比重,可直接影響葉片表面熱量的交換,因而對(duì)LT有較明顯的影響。文獻(xiàn)[9]在栓皮櫟冠層溫度多因子復(fù)合相關(guān)分析中即考慮了風(fēng)速。劉明虎等[19]對(duì)干旱區(qū)植物L(fēng)T模型的研究表明,在給定的條件下,風(fēng)速為10和0.1m·s?1,對(duì)應(yīng)LT為32℃和53.5℃,葉氣溫差LATD 為21.5℃,說明風(fēng)速較小時(shí)LT與氣溫有較大差異,而風(fēng)速較大時(shí),LATD則會(huì)降低。植物受高溫脅迫是在高LT下發(fā)生的,所以定量研究風(fēng)速的影響很有必要。

    研究者認(rèn)為葉片是對(duì)小氣候變化較敏感的器官,最能體現(xiàn)小氣候因子的影響或植物對(duì)小氣候的適應(yīng)性[4]。所以,研究LT變化特征及其與光、溫、濕、風(fēng)等氣候因子的關(guān)系,對(duì)深入揭示高溫脅迫下死苗的機(jī)制具有重要意義。光、溫、濕、風(fēng)小氣候因子對(duì)LT的影響是多因子復(fù)合影響,是個(gè)復(fù)雜的過程。本文從氣候空間理論[20]和圖解法入手,將土壤濕度設(shè)定為保持適宜狀態(tài),并設(shè)定不同的風(fēng)速等級(jí),研究LT隨不同光溫因子的變化,確定不同風(fēng)速的列線形態(tài),以分析光、溫、風(fēng)對(duì)LT的復(fù)合影響。用光、溫、濕、風(fēng)常規(guī)氣象因子監(jiān)測(cè)確定栓皮櫟葉片LT,并與葉片生物學(xué)上限溫度(Biological maximum temperature,BMT)進(jìn)行對(duì)比,據(jù)此判斷栓皮櫟幼苗葉片是否受到高溫脅迫,以及受脅迫的程度和持續(xù)時(shí)間。以期充實(shí)植物氣候空間研究理論,并為栓皮櫟實(shí)生苗造林時(shí)減少高溫干旱脅迫死苗提供依據(jù)。

    1 材料與方法

    1.1試驗(yàn)材料

    試驗(yàn)在河南農(nóng)業(yè)大學(xué)林學(xué)院試驗(yàn)基地(34° 43′N,112°42′E)進(jìn)行。為控制小氣候風(fēng)速因子,試驗(yàn)采用文獻(xiàn)[21]中的開頂箱設(shè)計(jì),用鍍鋅鋼管做骨架,覆塑料薄膜以控制外界風(fēng)對(duì)試驗(yàn)區(qū)的影響。包圍試驗(yàn)區(qū)大小為8m×10m,高3m。試材為一年生盆栽栓皮櫟幼苗,花盆上口直徑14cm,盆底直徑10cm,盆高24cm。2014年10月10日播種,每盆4粒種子,出苗后間苗,保證每盆一株幼苗。共10行,每行10盆共100盆,緊密排列并擺放在試驗(yàn)區(qū)南側(cè)中部。試驗(yàn)區(qū)其余區(qū)域種植植株高度30cm以下的栓皮櫟幼苗作保護(hù)行,使人工通風(fēng)形成的風(fēng)流線平滑通過待測(cè)試栓皮櫟幼苗。

    1.2測(cè)定方法

    2015年3-5月選典型晴天進(jìn)行測(cè)定。用北京產(chǎn)yaxin-0233多點(diǎn)生物溫度儀測(cè)定LT。同時(shí)配置3個(gè)熱電偶傳感器至3個(gè)栓皮櫟葉片,1個(gè)熱電偶探頭放入小型百葉箱內(nèi),并置于LT測(cè)點(diǎn)旁同高度測(cè)量氣溫。該儀器測(cè)溫精度為0.1℃,最小采集時(shí)間間隔0.2s。本試驗(yàn)設(shè)置采集時(shí)間間隔15s。在試驗(yàn)區(qū)中部選生長狀況良好,株高均勻(約15cm),平均成熟葉片為3~4片的花盆,每盆中選取一片上部無其它葉片遮擋的成熟葉片,用特制的鵝頸管支撐裝置將熱電偶探頭固定在葉片背面中部。由于熱電偶端點(diǎn)為測(cè)溫感應(yīng)部分,其直徑僅0.1mm,所以,接觸葉片即可感應(yīng)葉片細(xì)胞實(shí)際溫度。用3個(gè)葉片的平均溫度代表栓皮櫟幼苗LT。該LT是葉片本身的溫度,與紅外線測(cè)溫儀測(cè)得的冠層溫度不同,后者是葉片溫度和葉片間隙(可能是空氣、土壤或下層葉片)溫度的平均。在與測(cè)點(diǎn)同高度處安裝太陽輻射和風(fēng)速傳感器。太陽輻射(W·m?2)采用美國產(chǎn)LI-1400數(shù)據(jù)采集器,連接LI-191傳感器。設(shè)定LI-1400數(shù)據(jù)采集器采集時(shí)間間隔為15s,與熱電偶測(cè)溫儀一致。風(fēng)速測(cè)定采用RR-9200小氣候梯度自動(dòng)監(jiān)測(cè)系統(tǒng)(北京),配置有英國產(chǎn) Gill WindSonic超聲波風(fēng)速風(fēng)向傳感器,精度±2%。由于風(fēng)速的測(cè)定僅用于確定試驗(yàn)環(huán)境平均風(fēng)速等級(jí),所以設(shè)定RR-9200對(duì)Gill風(fēng)速傳感器通道每30s讀取一次通道數(shù)據(jù),每10min記錄一次。全天試驗(yàn)結(jié)束后當(dāng)日試驗(yàn)時(shí)段風(fēng)速的平均值即為當(dāng)日平均風(fēng)速。Gill風(fēng)速傳感器的主要特點(diǎn)是能測(cè)定微小風(fēng)速,避免三杯風(fēng)速傳感器起動(dòng)風(fēng)速過大,無法研究無風(fēng)或微風(fēng)小氣候的問題。試驗(yàn)區(qū)風(fēng)速控制采用2個(gè)不同功率的普通風(fēng)扇,在試驗(yàn)區(qū)北側(cè)由西向東吹風(fēng)。流動(dòng)的空氣沿試驗(yàn)區(qū)邊緣的薄膜流動(dòng)并轉(zhuǎn)向,風(fēng)扇持續(xù)工作,則可在試驗(yàn)區(qū)南部放置試驗(yàn)樣本的花盆處形成穩(wěn)定的風(fēng)速。通過調(diào)節(jié)風(fēng)扇功率可調(diào)節(jié)試驗(yàn)樣本測(cè)點(diǎn)風(fēng)速大小。風(fēng)速控制設(shè)0、0.3、0.6、0.7、0.8、0.9、1.0m·s?1等,每個(gè)風(fēng)速檔測(cè)定3d。實(shí)際測(cè)定時(shí)當(dāng)風(fēng)速達(dá)0.8m·s?1以上時(shí),CTD隨太陽輻射變化已較小,所以不再進(jìn)行更高風(fēng)速的測(cè)定。土壤濕度采用美國產(chǎn)LR5數(shù)據(jù)采集器連接EC-5土壤濕度傳感器監(jiān)測(cè)。使土壤濕度保持在體積比含水率15%~25%,即一直保持土壤適宜水分狀態(tài)。每日試驗(yàn)結(jié)束時(shí)下載土壤濕度數(shù)據(jù),低于15%即澆水,保證次日試驗(yàn)土壤為適宜水分狀態(tài)。

    1.3分析方法

    Gates在分析野生動(dòng)物生存的小氣候條件時(shí)給出了氣候空間圖解[21],如圖1。該圖是由輻射和氣溫形成的一個(gè)二維空間,空間內(nèi)的任意點(diǎn)對(duì)應(yīng)橫坐標(biāo)上不同的太陽輻射值和縱坐標(biāo)上不同的氣溫值而維持植物不同的LT。這個(gè)空間有由輻射因素決定的左右界限,和由溫度決定的上下界限。左右界限即圖中標(biāo)有晴空夜晚和晴空白天的輻射曲線,上下界限受風(fēng)速的影響而有不同的邊界線,即圖中標(biāo)有V=0.1m·s?1和V=2.0m·s?1的直線。位于上下限和左右限之間的區(qū)域?yàn)樯锷娴臍夂蚩臻g。對(duì)于氣候空間上限而言,其上限界線即為植物體溫等于BMT的點(diǎn)的連線,界線以上說明植物體溫超過了BMT而處于氣候空間之外;界線以下說明植物體溫低于BMT而處于氣候空間之內(nèi)。根據(jù)Gates的觀點(diǎn),生物長時(shí)間處于該氣候空間之外,將導(dǎo)致死亡。生物受小氣候因子脅迫均發(fā)生在氣候空間邊緣,超過氣候空間邊緣以外越遠(yuǎn),受脅迫影響的受害程度越大。氣候空間上下限與左右限不同。左右限由物體的輻射的物理學(xué)特性決定。上下限首先由生物能耐受的體溫最大值和最小值(相當(dāng)于生物學(xué)上、下限溫度)決定,在此基礎(chǔ)上,受風(fēng)速影響,上下限界線斜率不同。風(fēng)速大,則斜率?。伙L(fēng)速小,則斜率大。對(duì)于風(fēng)速為0.1m·s?1的氣候空間上限,當(dāng)氣溫為30℃時(shí),如果輻射為700W·m?2,則植物體處于氣候空間范圍以外,不能正常生存;如果環(huán)境溫度仍為30℃,而輻射為500W·m?2,則植物體處于氣候空間范圍內(nèi),能夠正常生存。Gates的氣候空間理論對(duì)生物能否生存的分析不僅僅用氣溫來衡量,同時(shí)考慮了太陽輻射因子和風(fēng)速因子。

    圖1 Gates的氣候空間圖Fig. 1 The climate space graph of Gates

    氣候空間理論為研究干旱半干旱地區(qū)栓皮櫟幼苗高溫脅迫死苗提供了理論基礎(chǔ),祖元?jiǎng)偟萚22]認(rèn)為應(yīng)該在植物生態(tài)領(lǐng)域研究氣候空間問題。這使植物的氣候空間界限指標(biāo)不是簡(jiǎn)單的一個(gè)溫度值,而是一定風(fēng)速下的太陽輻射和氣溫的復(fù)合指標(biāo)值。干旱時(shí)在太陽輻射照射下,LT大于冠層空氣溫度CAT (Canopy Air Temperature)(圖1中為Ta),CAT與LT之差為冠層溫度衰減CTD,即

    CTD=CAT-LT(1)

    CTD在國外文獻(xiàn)中稱為冠層溫度衰減[23-24],CTD絕對(duì)值越大表示冠層溫度衰減越多。它與國內(nèi)文獻(xiàn)中所稱葉氣溫差[25]、LATD[26]或ΔT[27]、Tc-a[28]絕對(duì)值相等,符號(hào)相反。即

    LATD=LT-CAT=-CTD(2)

    本文使用CTD表述冠層溫度變化,是因?yàn)樗m于氣候空間上下限界線分析,因?yàn)镃TD隨太陽輻射增加而增加的趨勢(shì)與Gates的氣候空間上限界線斜率趨勢(shì)一致。以氣候空間上限分析為例,葉片能耐受的LT的最高值是一定的,就是葉片的生物學(xué)上限溫度BMT,它是由植物本身的生物學(xué)特性決定的。把這個(gè)固定的BMT放入氣候空間上限界線上的某一點(diǎn)上,對(duì)BMT需要的氣溫和太陽輻射值進(jìn)行分析。依據(jù)氣候空間圖解(圖1),理論上,氣候空間上限界線延長線在縱軸上的交點(diǎn)即為在太陽輻射為0時(shí)形成BMT所需要的氣溫CAT,即此時(shí)一定風(fēng)速下,CTD與太陽輻射呈線性變化關(guān)系,不同風(fēng)速下二者的關(guān)系即形成一組列線。如果確定了栓皮櫟幼苗葉片的BMT和CTD列線,Gates的氣候空間上限即可表示為氣候空間圖上的具體位置,即在氣溫坐標(biāo)軸上找到等于BMT的點(diǎn),將CTD列線平移至該BMT的位置,則可得到栓皮櫟幼苗生存上限溫度的列線,即得出栓皮櫟幼苗的生存氣候空間上限。陳景玲等[29]基于葉綠素?zé)晒鈪?shù)研究了離體栓皮櫟幼苗葉片的BTM,進(jìn)一步將測(cè)定在體栓皮櫟幼苗葉片的BTM,與本文CTD列線相加,即可確定栓皮櫟幼苗葉片的氣候空間上限界線。

    2 結(jié)果與分析

    2.1靜風(fēng)條件下LATD隨太陽輻射強(qiáng)度的變化趨勢(shì)

    靜風(fēng)條件下典型晴天栓皮櫟葉片蒸騰速率(Tr)的日變化如圖2a。由圖可見,每日上午隨著太陽輻射的增加Tr呈增加趨勢(shì),12:00左右達(dá)到最大值1.6mmol·m?2·s?1,16:00之后隨著太陽輻射的降低Tr值迅速降低[29]。相應(yīng)地,由4月15、16和26日的葉氣溫差(LATD)隨太陽輻射的變化可見(圖2b、2c、2d),總體上LATD隨太陽輻射增加而增大,二者呈較明顯的線性關(guān)系。但從數(shù)據(jù)點(diǎn)分布情況看,圖2b、2c在太陽輻射強(qiáng)度為120W·m?2以上和以下時(shí)回歸直線斜率不同。這主要與葉片的Tr有關(guān),太陽輻射強(qiáng)度在120W·m?2以下時(shí),Tr速率很低,在0.2mmol·m?2·s?1左右,說明Tr耗熱很小。隨著太陽輻射的增加LATD增加很快,回歸直線斜率較大。當(dāng)太陽輻射增至120W·m?2以上時(shí),由于Tr迅速增加,Tr耗熱增加,故LATD增加減緩。由于Tr耗熱的作用低于太陽輻射的增溫作用,LATD仍隨著太陽輻射的增加而增加,但回歸直線斜率降低。圖2d則無這種回歸線斜率的差異,主要是因?yàn)榇杭旧郎乜欤?月26日比15、16日氣溫高,4月26日LT測(cè)定時(shí)沒有觀測(cè)到Tr速率很?。o蒸騰)的時(shí)段。

    對(duì)4月15、16日LATD隨太陽輻射變化曲線進(jìn)行2段分段線性模擬。當(dāng)太陽輻射<120W·m?2時(shí)LATD為0℃左右,此時(shí)的點(diǎn)較為集中,回歸直線斜率較大,15日為0.0317,16日為0.0461,LATD變化迅速;>120W·m?2時(shí)回歸直線斜率降低,15日為0.0122,16日為0.0167。即太陽輻射<120W·m?2時(shí),太陽輻射小,栓皮櫟葉片Tr還沒有開始或很小,為無蒸騰階段,因而LATD隨太陽輻射上升較快;>120W·m?2為植物持續(xù)蒸騰階段,但潛熱消耗對(duì)LATD的降低作用低于太陽輻射對(duì)LATD的增加作用。LATD仍隨著太陽輻射的上升而上升,但上升速率低于太陽輻射<120W·m?2時(shí)段。此階段數(shù)據(jù)點(diǎn)分散是由于持續(xù)蒸騰階段葉片Tr波動(dòng)造成的。

    對(duì)于持續(xù)蒸騰階段,圖2b、2c、2d斜率較接近,分別為0.0122、0.0167和0.0175。可知,LATD隨太陽輻射變化呈較穩(wěn)定的線性變化趨勢(shì)。

    圖2 葉片蒸騰速率(Tr)和典型晴天靜風(fēng)條件下葉氣溫差(LATD)隨太陽輻射強(qiáng)度(R)的變化Fig. 2 The change trend of transpiration rate (Tr) and leaf air temperature difference (LATD) with radiation in aero static conditions

    2.2有風(fēng)條件下CTD隨太陽輻射的變化趨勢(shì)

    圖3 不同風(fēng)速下測(cè)得的冠層溫度衰減(CTD)與太陽輻射(R)的關(guān)系(2015年)Fig. 3 The relationship between canopy temperature depression (CTD) and solar radiation (R) at different wind speeds(2015)

    平均風(fēng)速為0.3、0.6、0.7、0.8m·s?1時(shí)CTD與太陽輻射的關(guān)系如圖3。由圖可以看出,各風(fēng)速下CTD均隨太陽輻射呈線性變化趨勢(shì),但隨著風(fēng)速的增加,回歸直線的斜率越來越小,無風(fēng)時(shí)斜率為0.9353,接近于1,風(fēng)速為0.8m·s-1時(shí)斜率為0.089,接近0。說明隨著風(fēng)速的增加,CTD受太陽輻射的影響越來越小,同時(shí)顯示了風(fēng)速的作用越來越大。風(fēng)速增加加大了葉片表面的對(duì)流熱交換,風(fēng)速達(dá)0.8m·s-1以上時(shí),CTD幾乎不隨太陽輻射增加而變化。試驗(yàn)均在4月進(jìn)行,故Tr速率的日變化仍參見圖2a。由圖3可見,2015年4月23日(圖3a3)和4月13日(圖3b2)均有與圖2無蒸騰階段的特征時(shí)段相類似的階段,即在太陽輻射較低時(shí),CTD隨太陽輻射的增加迅速減少(相當(dāng)于圖2中的LATD迅速增加)階段。但圖3中其它均未觀測(cè)到這一現(xiàn)象。原因之一可能是葉片蒸騰從夜間到白天,從無蒸騰到持續(xù)蒸騰過渡階段發(fā)生時(shí)間不在試驗(yàn)的觀測(cè)時(shí)段內(nèi);原因之二可能是葉片蒸騰從夜間到白天變化不大造成的。如圖3a2即2015年4月9日也有與圖2無蒸騰階段的特征時(shí)段相類似的階段,但前者不太明顯。蒸騰對(duì)CTD的影響有待進(jìn)一步研究。

    由式(2)可知,圖2中靜風(fēng)條件下LATD與太陽輻射關(guān)系的回歸直線方程斜率的負(fù)值即為CTD的斜率。綜合分析風(fēng)速為0、0.3、0.6、0.7和0.8m·s-1時(shí)CTD與太陽輻射的關(guān)系回歸直線的斜率,將同一風(fēng)速下的數(shù)據(jù)回歸方程的斜率取平均值,得到風(fēng)速為0、0.3、0.5、0.6、0.7、0.8m·s-1時(shí)栓皮櫟CTD的列線如圖4。

    圖4 不同風(fēng)速條件下CTD隨R變化列線Fig. 4 The series relationship curves between CTD and R at different wind speeds

    由圖4可以看出,風(fēng)速從0增至0.3m·s-1,CTD與太陽輻射的回歸直線的斜率由-0.0155降至-0.0067,變化較大。而風(fēng)速從0.7m·s-1增至0.8m·s-1,CTD與太陽輻射的回歸直線的斜率由-0.0021 降至-0.0004,變化則較小。由此可見,栓皮櫟LT受風(fēng)速影響較大,從無風(fēng)轉(zhuǎn)至有風(fēng)條件時(shí)尤受影響。當(dāng)風(fēng)速逐漸增大時(shí)回歸直線的斜率的絕對(duì)值逐漸減小,且減少的幅度也逐漸變小,直線分布漸趨密集。

    從圖4還可看出,風(fēng)速為0時(shí),回歸直線斜率絕對(duì)值最大,其斜率為-0.01547。該CTD線在縱軸上截距約為0.8℃,即當(dāng)靜風(fēng)且太陽輻射為0時(shí),CTD 有0.8℃的初始值。這是由于靜風(fēng)夜晚栓皮櫟幼苗葉片長波有效輻射作用,使LT低于CAT,因而CTD為正值是合理的。

    當(dāng)風(fēng)速為0.8m·s-1時(shí),回歸直線幾乎呈水平狀態(tài),即斜率-0.00037,已接近0。該水平線在縱軸上截距約為-0.9℃,即當(dāng)風(fēng)速為0.8m·s-1時(shí),盡管CTD不隨太陽輻射增加而增加,但CTD仍維持平均-0.9℃左右。即LT總比CAT高0.9℃。

    將圖4的V=0和V=0.8m·s-1CTD與太陽輻射的回歸直線與圖1中的氣候空間上限界線V=0.1m·s-1和V=2.0m·s-1進(jìn)行比較,可以看出,圖4中V=0和V=0.8m·s-1有一個(gè)交叉點(diǎn),與圖1氣候空間上限界線形式非常接近。

    3 結(jié)論與討論

    (1)栓皮櫟幼苗葉片葉氣溫差LATD隨著太陽輻射的增強(qiáng)逐漸增大,葉溫LT明顯偏離氣溫,LATD隨著太陽輻射的變化呈直線正相關(guān)關(guān)系。冠層溫度衰減CTD隨著太陽輻射的變化呈直線負(fù)相關(guān)關(guān)系。劉云鵬等[26]對(duì)玉米干旱脅迫日變化特征研究得到了與本文一致的結(jié)果,即LATD先隨太陽輻射增加而上升,后隨著其減弱而降低的規(guī)律。孟平等[27]研究蘋果樹LATD與微氣象因子的關(guān)系,在不考慮土壤水分因子的條件下,LATD與冠層凈輻射呈正相關(guān),這也與本文結(jié)果一致。但在考慮土壤水分因子影響時(shí),太陽輻射的作用可能使植物Tr增加,而導(dǎo)致LATD與冠層凈輻射的關(guān)系為負(fù)相關(guān)關(guān)系[27]。

    (2)受蒸騰作用消耗潛熱的影響,在晴朗無風(fēng)天氣下,栓皮櫟幼苗葉片LATD隨太陽輻射的變化在清晨、傍晚無蒸騰與白天持續(xù)有蒸騰時(shí)段的負(fù)相關(guān)直線的斜率不同,無蒸騰時(shí)段斜率的絕對(duì)值較大,而持續(xù)蒸騰階段較小。張衛(wèi)紅等[31]對(duì)不同土壤濕度下小麥LATD對(duì)光照強(qiáng)度的反應(yīng)的研究結(jié)果認(rèn)為,隨著土壤濕度的增加,LATD隨光照強(qiáng)度增加的回歸系數(shù)逐漸減少。文獻(xiàn)[31]將這種現(xiàn)象解釋為太陽輻射用于光合積累,而不再是影響LATD的有效因素,本文認(rèn)為值得商榷,因?yàn)楦叩戎参锕夂献饔美玫奶栞椛湔继栞椛淇偭康谋壤话憔^低,不足以影響到葉片表面的熱量平衡。這種現(xiàn)象更可能的原因是隨著土壤濕度的增加,葉片蒸騰速率加大,消耗到潛熱的量增加,此時(shí),同樣的太陽輻射能量輸入,LT的增加卻較少。本文的無蒸騰階段LATD隨太陽輻射變化的斜率較高的結(jié)果印證了文獻(xiàn)[31]觀測(cè)到的現(xiàn)象。

    (3)在晴朗有風(fēng)的天氣條件下,栓皮櫟幼苗的CTD隨太陽輻射的增強(qiáng)而逐漸增大。不同風(fēng)速下二者的回歸直線形成列線,其斜率隨著風(fēng)速的增大而逐漸減小,尤以從無風(fēng)到有風(fēng)條件時(shí)的變化明顯;當(dāng)風(fēng)速繼續(xù)增大時(shí)回歸直線斜率減小幅度逐漸變小,直線分布越來越密集,斜率越來越趨近于0。也即隨著風(fēng)速的增大,栓皮櫟LT越來越接近氣溫。前人研究蘋果樹葉片LATD與風(fēng)的關(guān)系[27],不考慮土壤水分時(shí)二者為正相關(guān),考慮土壤水分時(shí)為負(fù)相關(guān),且其偏相關(guān)系數(shù)較低。說明把風(fēng)速作為一個(gè)氣象因子直接參與回歸分析還不能很好地表現(xiàn)風(fēng)速對(duì)LATD 或CTD的影響。因?yàn)轱L(fēng)的對(duì)流換熱作用和促進(jìn)蒸騰水分向大氣擴(kuò)散作用對(duì)LATD或CTD的影響是雙重的,且不呈線性變化。本文研究結(jié)果風(fēng)速從0增至0.8m·s-1的列線之間斜率變化逐漸減少,而非均勻分布證實(shí)了這一點(diǎn)。所以本文用不同風(fēng)速的列線分析風(fēng)速對(duì)LT的影響更合理。

    (4)當(dāng)風(fēng)速為0時(shí),CTD與太陽輻射的回歸直線斜率的絕對(duì)值最大,即斜率-0.01547。該CTD線在縱軸上截距約為0.8℃,即當(dāng)靜風(fēng)且太陽輻射為0時(shí),CTD有0.8℃的初始值。文獻(xiàn)[9]也觀測(cè)到栓皮櫟葉片夜間LT低于CAT。這可以從靜風(fēng)夜晚栓皮櫟幼苗葉片長波有效輻射作用使LT低于CAT,因而CTD為正值,而得到合理的解釋。

    當(dāng)風(fēng)速為0.8m·s-1時(shí),回歸直線幾乎水平,其斜率為-0.00037已接近于0。該水平線在縱軸上截距約為-0.9,即當(dāng)風(fēng)速為0.8m·s-1時(shí),盡管CTD不隨太陽輻射增加而增加,但會(huì)維持平均-0.9℃左右的CTD。即LT總比CAT高0.9℃。許多植物葉溫與氣象條件關(guān)系研究的文獻(xiàn)未考慮風(fēng)的影響[32-35],如果試驗(yàn)環(huán)境風(fēng)速變化,勢(shì)必影響試驗(yàn)結(jié)果。從本文研究結(jié)果可以看出,風(fēng)是重要的CTD影響因子。

    (5)栓皮櫟幼苗葉片的CTD隨風(fēng)速的變化列線形式與能量生態(tài)學(xué)中的植物體溫氣候空間[21]上限界線相似。風(fēng)速從0增至0.8m·s-1時(shí)栓皮櫟幼苗葉片CTD列線與Gates氣候空間理論的上限趨勢(shì)一致。該研究實(shí)踐了Gates的氣候空間理論,給植物氣候空間上限界線賦予了實(shí)際含義,即疊加在栓皮櫟幼苗葉片生物學(xué)上限溫度BMT之上的不同風(fēng)速的CTD列線可以表達(dá)氣候空間上限。

    (6)風(fēng)速為0與0.8m·s-1的CTD回歸直線和氣候空間上限界線風(fēng)速為0.1m·s-1與2.0m·s-1都有一個(gè)交叉點(diǎn),二者形式非常接近。但風(fēng)速為0.3、0.6和0.7m·s-1時(shí)各回歸線并未交叉于該點(diǎn)。按照氣候空間理論,這個(gè)交叉點(diǎn)發(fā)生在氣候空間上限時(shí),是黑體輻射線通過的點(diǎn)[21]。本試驗(yàn)觀測(cè)到了這個(gè)交叉點(diǎn),但不是在氣候空間上限溫度下測(cè)定的CTD,所以該交叉點(diǎn)的含義還有待進(jìn)一步探討。

    (7)土壤濕度對(duì)LATD隨太陽輻射的變化規(guī)律有較大影響,如張衛(wèi)紅等[31]對(duì)小麥葉片的研究。本文是在適宜土壤濕度下進(jìn)行的試驗(yàn),土壤濕度條件主要影響Tr,進(jìn)而影響LT。若設(shè)定不同土壤濕度重復(fù)本文試驗(yàn),則可得到不同土壤濕度條件、對(duì)應(yīng)不同風(fēng)速的列線。此外,對(duì)于空氣濕度、葉片特性等因素對(duì)LT的影響,以及進(jìn)而產(chǎn)生的對(duì)栓皮櫟葉片CTD列線的定量影響也有待于進(jìn)一步研究。

    參考文獻(xiàn)References

    [1]馬莉薇,張文輝,薛瑤芹,等.秦嶺北坡不同生境栓皮櫟實(shí)生苗生長及其影響因素[J].生態(tài)學(xué)報(bào),2010,30(23):6512-6520.

    Ma L W,Zhang W H,Xue Y Q,et al.Growth characteristics and influencing factors of Quercus variabilis seedlings on the north slope of Qinling Mountains[J].Acta Ecologica Sinica,2010,30(23):6512-6520.(in Chinese)

    [2]陳景玲,王靜,王謙,等.基于葉綠素?zé)晒獾那G條灌叢對(duì)栓皮櫟幼苗庇蔭效應(yīng)研究[J].西北林學(xué)院學(xué)報(bào),2014,29(4):46-53.

    Chen J L,Wang J,Wang Q,et al.Patronage effect of verain family to quercus variabilis saplings based on leaf chlorophyll fluorescence[J].Journal of Northwest Forestry University, 2014,29(4):46-53.(in Chinese)

    [3]周健民.秦嶺北坡栓皮櫟種群有性生殖及生存策略研究[D].楊凌:西北農(nóng)林科技大學(xué),2010.

    Zhou J M.Study on the sexual reproduction and survival strategy of the Quercus variablis population in north slope of Qinling mountains[D].Yangling:Northwest A & F University, 2010.(in Chinese)

    [4]趙晶.溫度和光照對(duì)兩種櫟屬植物幼苗葉形態(tài)與結(jié)構(gòu)性狀的影響[D].濟(jì)南:山東大學(xué),2013.

    Zhao J.Influence of temperature and light intensity on morphological and structural trait of leaf seedings of Quercus variabilis and Q. mongolica[D].Jinan:Shandong University, 2013.(in Chinese)

    [5]吳敏.栓皮櫟種群實(shí)生苗定居過程及其對(duì)干旱脅迫的響應(yīng)[D].楊凌:西北農(nóng)林科技大學(xué),2013.

    Wu M.Seedlings settlement of Quercus variabilis BL.and seedlings response to drought stress[D].Yangling:North West Agriculture and Forestry University,2013.(in Chinese)

    [6]王婧,王少波,康宏樟,等.東亞地區(qū)栓皮櫟的地理分布格局及其氣候特征[J].上海交通大學(xué)學(xué)報(bào)(農(nóng)業(yè)科學(xué)版),2009,(3): 235-241.

    Wang J,Wang S B,Kang H Z,et al.Distribution pattern of oriental oak (Quercus variabilis Blume) and the characteristics of climate of distribution area in eastern Asia[J].Journal of Shanghai Jiaotong University(Agricultural Science),2009,(3):235-241.(in Chinese)

    [7]王婧.區(qū)域尺度上栓皮櫟葉片氮、磷含量的地理變異及其與環(huán)境因子的關(guān)系[D].上海:上海交通大學(xué),2009. Wang J.The geographic variation of leaf nitrogen andphosphorus in oriental oak (Quercus variabilis) and its relation to environmental factor[D].Shanghai :Shanghai Jiaotong University,2009.(in Chinese)

    [8]楊文卿,張志銘,侯娜,等.太行山低山丘陵區(qū)3種樹木葉溫對(duì)土壤水分脅迫的響應(yīng)[J].河南農(nóng)業(yè)大學(xué)學(xué)報(bào),2015,49(3): 325-330.

    Yang W Q,Zhang Z M,Hou N,et al.Responses of leaf temperature of three types of trees to soil water stress in hilly areas of Taihang Mountain[J].Journal of Henan Agricultural University,2015,49(3):325-330.(in Chinese)

    [9]魏丹丹,張勁松,孟平,等.栓皮櫟人工林冠層溫度變化特征及其與微氣象因子的關(guān)系[J].應(yīng)用生態(tài)學(xué)報(bào),2012,23(7):1767-1773.

    Wei D D,Zhang J S,Meng P,et al.Variations of canopy temperature in Quercus variabilis plantation and their relations with micrometeorological factors[J].Chinese Journal of Applied Ecology,2012,23(7):1767-1773.(in Chinese)

    [10]李永華,盧琦,吳波,等.干旱區(qū)葉片形態(tài)特征與植物響應(yīng)和適應(yīng)的關(guān)系[J].植物生態(tài)學(xué)報(bào),2012,36(1) :88-98.

    Li Y H,Lu Q,Wu B,et al.A review of leaf morphology plasticity linked to plant response and adaptation characteristics in arid ecosystems[J].Chinese Journal of Plant Ecology,2012,36(1):88-98.(in Chinese)

    [11]Smith W K,Geller G N.Leaf and environmental parameters influencing transpiration:theory and field measurements[J]. Oecologia,1980,46:308-313.

    [12]Jones H G.Plants and microclimate:a quantitative approach to environmental plant physiology[M].Cambridge, UK: Cambridge University Press,1992.

    [13]Nobel P S.Physicochemical and environmental plant physiology[M].San Diego,USA:Academic Press,2005.

    [14]Vogel S.Leaves in the lowest and highest winds:temperature, force and shape[J].New Phytologist,2009,183:13-26.

    [15]陳慧新.北京山區(qū)主要樹種光合蒸騰與耗水特性研究[D].北京:北京林業(yè)大學(xué),2008.

    Chen H X.Study on photosynthetic and transpiration of dominant tree species in Beijing mountainou area[D]. Beijing:Beijing Forestry University,2008.(in Chinese)

    [16]謝會(huì)成.栓皮櫟光合生理生態(tài)的研究[D].南京:南京林業(yè)大學(xué),2002.

    Xie H C.A study on the photosynthetic physioecology of oriental oak[D].Nanjing:Nanjing Forestry University,2002. (in Chinese)

    [17]徐湘婷.風(fēng)與光對(duì)于歐美107楊,栓皮櫟與側(cè)柏生理生態(tài)影響[D].北京:中國林業(yè)科學(xué)研究院,2014.

    Xu X T.Physiological and ecologic traits of hybrid popular(Populus euramericana),Oak(Quercus variabilis) and Chinese arborvitae(Platycladus orientalis)in response to wind and light[D].Beijing:Chinese Academy of Forestry, 2014.(in Chinese)

    [18]徐湘婷,王林,黃平,等.風(fēng)和光照對(duì)側(cè)柏、栓皮櫟生長和形態(tài)建成的影響[J].林業(yè)科學(xué),2014,50(7):164-168.

    Xu X T,Wang L,Huang P,et al.Effects of wind and irradiation on growth and morphogenesis of Platycladus orientalis and Quercus variabilis[J].Scientia Silvae Sinicae, 2014,50(7):164-168.(in Chinese)

    [19]劉明虎,辛智鳴,徐軍,等.干旱區(qū)植物葉片大小對(duì)葉表面蒸騰及葉溫的影響[J].植物生態(tài)學(xué)報(bào),2013,37(5):436-442.

    Liu M H,Xin Z M,Xu J,et al.Influence of leaf size of plant on leaf transpiration and temperature in arid regions[J].Chinese Journal of Plant Ecology,2013,37(5):436-442.(in Chinese)

    [20]Gates D M.Biophysical Ecology[M].New York:Springer Verleg,1980:15-16.

    [21]牛書麗,韓興國,馬克平,等.全球變暖與陸地生態(tài)系統(tǒng)研究中的野外增溫裝置[J].植物生態(tài)學(xué)報(bào),2007,31(2):262- 271.

    Niu S L,Han X G,Ma K P,et al.Field facilities in global warming and terrestrial ecosystem research[J].Journal of Plant Ecology,2007,31(2):262- 271.(in Chinese)

    [22]祖元?jiǎng)?能量生態(tài)學(xué)引論[M].長春:吉林科學(xué)技術(shù)出版社,1990:127-132.

    Zu Y G.Introduction to energy ecology[M].Changchun:Jinlin Science and Technology Publish House,1990:127-132.(in Chinese)

    [23]Amani I,Fischer R A,Reynolds M P.Canopy temperature depression association with yield of irrigated spring wheat cultivars in a hot climate[J].Journal of Agronomy and Crop Science,1996,176(2):119-129.

    [24] Kottmann L,Schittenhelm S,Wittich K P,et al.Suitability of canopy temperature depression in a temperate climate with drought-stressed winter rye,determined with three infrared measurement devices[J].Journal of Agronomy and Crop Science,2013,199(6):385-394.

    [25]程林,江勝國,韓正英,等.安徽省桐城市一季水稻抽穗開花期葉片溫度分析[J].中國農(nóng)業(yè)氣象,2009,30(S2): 218-220.

    Cheng L,Jiang S G,Han Z Y,et al.Analysis on leaf temperature of single-season rice on ear sprouting to flowering period in Tongcheng[J].Chinese Journal of Agrometeorology,2009,30(S2):218-220.(in Chinese)

    [26]劉云鵬,申思,潘余強(qiáng),等.干旱脅迫下玉米葉-氣溫差與葉溫差日變化特征及其品種差異[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào), 2014, 19(5):13-21.

    Liu Y P,Shen S,Pan Y Q,et al.Diurnal variation in leaf air temperature difference and leaf temperature difference and the hybrid difference in maize under different drought stress[J].Journal of China Agricultural University,2014,19(5): 13-21.(in Chinese)

    [27]孟平,張勁松,高峻,等.蘋果樹冠層-空氣溫差變化及其與環(huán)境因子的關(guān)系[J].應(yīng)用生態(tài)學(xué)報(bào),2007,18(9):2030-2034.

    Meng P,Zhang J S,Gao J,et al.Variation of apple tree canopy-air temperature difference and its relations to environment factors[J].Chinese Journal of Applied Ecology, 2007,18(9):2030-2034.(in Chinese)

    [28]劉瑞文,董振國.小麥葉溫對(duì)籽粒灌漿的影響[J].中國農(nóng)業(yè)氣象,1992,(3):1-5.

    Liu R W,Dong Z G.Effect of leaf temperature on grain milking of winter wheat[J].Chinese Journal of Agrometeorology, 1992,(3):1-5.(in Chinese)

    [29]陳景玲,宋曉明,王謙,等.基于葉綠素?zé)晒鈪?shù)的栓皮櫟葉片PSⅡ失活高溫指標(biāo)[J].中國農(nóng)業(yè)氣象,2013,34(5):563-568. Chen J L, Song X M,Wang Q,et al.High temperature index of PSⅡinactivation according to chlorophyll fluorescence of Quercus variabilis leaves[J].Chinese Journal of Agrometeorology, 2013,34(5):563-568.(in Chinese)

    [30]謝會(huì)成,姜志林.栓皮櫟蒸騰速率的日變化、季節(jié)動(dòng)態(tài)及其對(duì)遮蔭的響應(yīng)[J].山東林業(yè)科技,2010,(4):70-71.

    Xie H C,Jiang Z L.The diurnal variation and seasonal dynamics of transpiration rate of Quercus variabilis, and its response to shade[J].Shandong Forestry Science and Technology,2010,(4):70-71.(in Chinese)

    [31]張紅衛(wèi),申雙和,陳懷亮.不同土壤濕度下小麥葉氣溫差對(duì)光照強(qiáng)度的反應(yīng)及對(duì)蒸騰的影響[J].河南農(nóng)業(yè)科學(xué),2008, (7):25-28.

    Zhang H W,Shen S H,Chen H L.Influences of leaf-air temperature on illumination intensity and transpiration-rate with various soil moistures[J].Journal of Henan Agricultural Sciences,2008, (7):25-28.(in Chinese)

    [32]陳金華,岳偉,楊太明.水稻葉溫與氣象條件的關(guān)系研究[J].中國農(nóng)學(xué)通報(bào),2011,27(12):19-23.

    Chen J H,Yue W,Yang T M.A study on the relationship between leaf temperature of rice and meteorological factors[J].Chinese Agricultural Science Bulletin,2011,27(12): 19-23.(in Chinese)

    [33]姚振坤,羅新蘭,李天來,等.日光溫室番茄葉溫的模擬及與環(huán)境因子的關(guān)系[J].江蘇農(nóng)業(yè)學(xué)報(bào),2010,26(3):587-592.

    Yao Z K,Luo X L,Li T L,et al.A Simulation model of the relationship between tomato leaf temperature and ambient factors in solar greenhouse[J].Jiangshu J. of Agr. Sci.,2010, 26(3):587-592.(in Chinese)

    [34]吳強(qiáng),須暉,韓亞東.日光溫室番茄葉溫變化特性研究[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào),2008,38(5):618-620.

    Wu Q,Xu H,Han Y D.Study on temperature change of tomato leaf in facility cultivation[J].Journal of Shenyang Agricultural University,2008,38(5):618-620.(in Chinese)

    [35]高明超.水稻冠層溫度特性及基于冠層溫度的水分脅迫指數(shù)研究[D].沈陽:沈陽農(nóng)業(yè)大學(xué),2013.

    Gao M C.A research on rice canopy temperature properties and in view of the water stress index of canopy temperature[D]. Shenyang:Shenyang Agricultural University, 2013.(in Chinese)

    Series Lines on Change of Leaf Temperature Difference of Quercus variabilis with Solar Radiation and Wind Speed

    WANG Pei-shu, WANG Wei-hong, CHEN Jing-ling, WANG Qian, LI Guo-huai
    (College of Forestry, Henan Agricultural University, Zhengzhou 450002, China)

    Abstract:In order to explore the effect of solar radiation, air temperature and wind speed on leaf temperature (LT) of Quercus variabilis, and to determine the consistency of LT and the series lines of Gates climate space, LT and air temperature (Ta) and solar radiation (R) were measured in clear day with the conditions of controlled wind speeds. The relationship between leaf air temperature difference (LATD) and solar radiation at different wind speeds was studied. The series lines of climate space of Quercus variabilis were determined by using Gates theory climate space and graphic method. The results showed LATD of Quercus seedlings increased linearly with the radiation increasing. Canopy temperature depression (CTD) of Quercus variabilis also increased with radiation increasing, but affected by the wind speed. CTD tended to be zero when wind speed was 0.8m·s?1. CTD of Quercus seedlings changed as R linearly with the radiation increasing and different wind speeds. The linear relationship was series lines and decreased with wind speed increasing. The series lines were in accordance with the upper limits of Gates theory climate space. The upper limit of climate space could be expressed by adding up biological maximum temperature (BMT) and CTD value on series lines with different wind speeds.

    Key words:Quercus variabilis; CTD; Solar radiation; Wind speed; Climate space; Series lines

    doi:10.3969/j.issn.1000-6362.2016.02.011

    * 收稿日期:2015-07-24**通訊作者。E-mail:chenjingling5@163.com

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目“基于氣候空間上限列線的栓皮櫟幼苗高溫疊加干旱脅迫機(jī)制研究”(31370621);國家自然科學(xué)基金項(xiàng)目“番茄葉片低溫脅迫氣象指標(biāo)的氣候空間列線研究”(31071321)

    作者簡(jiǎn)介:王佩舒(1991-),女,碩士生,主要從事農(nóng)林生態(tài)研究。E-mail:15093258404@163.com

    猜你喜歡
    栓皮櫟太陽輻射風(fēng)速
    邯鄲太陽輻射時(shí)空分布特征
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測(cè)
    基于最優(yōu)TS評(píng)分和頻率匹配的江蘇近海風(fēng)速訂正
    栓皮櫟定向培育技術(shù)要點(diǎn)淺析
    栓皮櫟研究進(jìn)展與未來展望
    豫南山區(qū)栓皮櫟生物防火林帶營造技術(shù)
    基于PCA 的太陽輻射觀測(cè)算法研究
    種植 栓皮櫟有前景
    太陽輻射作用下鋼筒倉結(jié)構(gòu)溫度場(chǎng)分析研究
    基于GARCH的短時(shí)風(fēng)速預(yù)測(cè)方法
    午夜老司机福利片| 久久久久国产精品人妻一区二区| 一级片'在线观看视频| e午夜精品久久久久久久| 成人国产麻豆网| 两性夫妻黄色片| 亚洲av中文av极速乱| 大话2 男鬼变身卡| 国产高清不卡午夜福利| 精品国产露脸久久av麻豆| 99re6热这里在线精品视频| 国产成人一区二区在线| 午夜免费男女啪啪视频观看| 少妇精品久久久久久久| 精品人妻一区二区三区麻豆| 天天躁夜夜躁狠狠躁躁| 久久精品久久久久久久性| 国产日韩欧美视频二区| 日韩中文字幕欧美一区二区 | 在线观看国产h片| 亚洲欧洲日产国产| 伦理电影免费视频| 黑人猛操日本美女一级片| 国产97色在线日韩免费| 国产成人免费观看mmmm| 亚洲精品乱久久久久久| 最黄视频免费看| 香蕉国产在线看| 久久性视频一级片| 69精品国产乱码久久久| 视频在线观看一区二区三区| 国产精品 国内视频| 久久免费观看电影| 久久免费观看电影| 成人毛片60女人毛片免费| 久久精品国产亚洲av高清一级| 亚洲av电影在线进入| 国产成人精品在线电影| 夫妻性生交免费视频一级片| 精品国产露脸久久av麻豆| 色精品久久人妻99蜜桃| 大香蕉久久成人网| 亚洲欧美一区二区三区黑人| 亚洲精品aⅴ在线观看| 男女下面插进去视频免费观看| 久久天躁狠狠躁夜夜2o2o | 人人妻人人澡人人爽人人夜夜| 赤兔流量卡办理| 亚洲精品久久成人aⅴ小说| 男人添女人高潮全过程视频| 日韩精品免费视频一区二区三区| 国产精品久久久av美女十八| 天天躁狠狠躁夜夜躁狠狠躁| 国产淫语在线视频| 另类亚洲欧美激情| 一本一本久久a久久精品综合妖精| 国产欧美亚洲国产| 国产精品av久久久久免费| 亚洲欧美一区二区三区黑人| 永久免费av网站大全| 亚洲精品自拍成人| 一级片免费观看大全| 亚洲一区二区三区欧美精品| 中文欧美无线码| 国产极品粉嫩免费观看在线| 性少妇av在线| 日本猛色少妇xxxxx猛交久久| 国产日韩一区二区三区精品不卡| 国产精品三级大全| xxxhd国产人妻xxx| 男男h啪啪无遮挡| 2021少妇久久久久久久久久久| 亚洲一区中文字幕在线| 久久人人爽人人片av| 免费看av在线观看网站| 成年av动漫网址| 国产成人系列免费观看| h视频一区二区三区| 亚洲精华国产精华液的使用体验| 精品酒店卫生间| 午夜福利网站1000一区二区三区| 日韩av在线免费看完整版不卡| 91aial.com中文字幕在线观看| 亚洲国产av新网站| 国产一区亚洲一区在线观看| 久久久久人妻精品一区果冻| 精品少妇久久久久久888优播| 超碰成人久久| 亚洲欧美日韩另类电影网站| 亚洲国产欧美日韩在线播放| 少妇精品久久久久久久| 久久久欧美国产精品| 国产免费现黄频在线看| 交换朋友夫妻互换小说| 1024视频免费在线观看| 91aial.com中文字幕在线观看| 国产成人91sexporn| 国产男女内射视频| 久久99热这里只频精品6学生| 国产免费一区二区三区四区乱码| 女性被躁到高潮视频| 成年av动漫网址| 国产精品香港三级国产av潘金莲 | 1024香蕉在线观看| 美女高潮到喷水免费观看| 成人毛片60女人毛片免费| 91国产中文字幕| 男人操女人黄网站| 久久这里只有精品19| 男女边吃奶边做爰视频| 中文字幕高清在线视频| 亚洲一区中文字幕在线| 国产精品久久久久久精品古装| 涩涩av久久男人的天堂| 亚洲人成网站在线观看播放| 在线天堂中文资源库| 最近最新中文字幕大全免费视频 | av在线观看视频网站免费| 日本一区二区免费在线视频| 欧美中文综合在线视频| 黑人欧美特级aaaaaa片| 99九九在线精品视频| 日韩 亚洲 欧美在线| 欧美 日韩 精品 国产| 国产免费又黄又爽又色| 国产精品二区激情视频| 天天添夜夜摸| 欧美亚洲 丝袜 人妻 在线| 久久这里只有精品19| 免费久久久久久久精品成人欧美视频| 啦啦啦在线免费观看视频4| 亚洲精品国产av成人精品| 亚洲三区欧美一区| 在线亚洲精品国产二区图片欧美| 丝袜喷水一区| 欧美日本中文国产一区发布| 久久精品国产亚洲av高清一级| 国产黄频视频在线观看| 各种免费的搞黄视频| 久久热在线av| 啦啦啦在线免费观看视频4| 男人爽女人下面视频在线观看| 男人舔女人的私密视频| 1024香蕉在线观看| 色网站视频免费| 三上悠亚av全集在线观看| 麻豆乱淫一区二区| 18禁动态无遮挡网站| 午夜影院在线不卡| 99re6热这里在线精品视频| 欧美变态另类bdsm刘玥| 捣出白浆h1v1| 电影成人av| 老汉色av国产亚洲站长工具| 久久热在线av| 另类精品久久| 国产精品久久久久久人妻精品电影 | 免费在线观看完整版高清| 亚洲国产看品久久| 日本黄色日本黄色录像| 精品免费久久久久久久清纯 | 秋霞伦理黄片| 人人妻人人澡人人爽人人夜夜| 电影成人av| 国产成人免费观看mmmm| 精品一区二区三区av网在线观看 | 国精品久久久久久国模美| 国产老妇伦熟女老妇高清| 女人被躁到高潮嗷嗷叫费观| 一本大道久久a久久精品| 国产精品蜜桃在线观看| 午夜免费男女啪啪视频观看| 狠狠婷婷综合久久久久久88av| 91国产中文字幕| 午夜激情久久久久久久| 亚洲精品久久午夜乱码| 9191精品国产免费久久| 亚洲三区欧美一区| 亚洲欧美色中文字幕在线| 亚洲精品美女久久久久99蜜臀 | 99国产精品免费福利视频| 欧美日韩一区二区视频在线观看视频在线| 中文欧美无线码| 一边亲一边摸免费视频| 好男人视频免费观看在线| 婷婷成人精品国产| 人妻 亚洲 视频| 美女脱内裤让男人舔精品视频| 日韩人妻精品一区2区三区| 国产熟女午夜一区二区三区| 男女国产视频网站| 大片免费播放器 马上看| 性色av一级| 精品人妻在线不人妻| 亚洲av日韩精品久久久久久密 | 成年动漫av网址| 国产有黄有色有爽视频| 日本色播在线视频| 亚洲国产精品一区三区| av不卡在线播放| 久久久精品区二区三区| 精品一品国产午夜福利视频| 国产成人精品在线电影| av免费观看日本| 777米奇影视久久| netflix在线观看网站| 国产爽快片一区二区三区| 久久精品国产亚洲av涩爱| 日韩一本色道免费dvd| av片东京热男人的天堂| 啦啦啦视频在线资源免费观看| 亚洲欧美激情在线| 亚洲成人免费av在线播放| videosex国产| 国产精品一国产av| 国产精品三级大全| av女优亚洲男人天堂| 国产精品久久久av美女十八| 久久女婷五月综合色啪小说| 国产精品一国产av| 国产亚洲最大av| 老司机影院成人| 涩涩av久久男人的天堂| 国产精品.久久久| 欧美日韩成人在线一区二区| 啦啦啦在线观看免费高清www| 亚洲色图综合在线观看| 国产精品国产三级国产专区5o| 亚洲第一av免费看| 久久女婷五月综合色啪小说| 亚洲色图 男人天堂 中文字幕| av天堂久久9| 日韩一区二区三区影片| 国产淫语在线视频| 欧美xxⅹ黑人| 国产成人91sexporn| 久久97久久精品| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 1024视频免费在线观看| 777米奇影视久久| 一区二区av电影网| 我的亚洲天堂| 亚洲精品国产区一区二| 熟女av电影| 这个男人来自地球电影免费观看 | 国产精品久久久av美女十八| 国产精品国产av在线观看| 亚洲欧美精品自产自拍| 国产成人欧美在线观看 | av免费观看日本| 免费观看性生交大片5| 性色av一级| kizo精华| 69精品国产乱码久久久| 黄色怎么调成土黄色| 国产无遮挡羞羞视频在线观看| 亚洲国产精品999| 欧美精品一区二区免费开放| av天堂久久9| 妹子高潮喷水视频| 国产精品久久久人人做人人爽| 亚洲第一区二区三区不卡| 女人精品久久久久毛片| 一级毛片我不卡| 黑人巨大精品欧美一区二区蜜桃| 精品酒店卫生间| 看非洲黑人一级黄片| 看免费av毛片| 久久久久久久大尺度免费视频| 亚洲精品美女久久久久99蜜臀 | 国产极品天堂在线| 搡老乐熟女国产| 青春草国产在线视频| 亚洲成人国产一区在线观看 | 中文欧美无线码| 中国国产av一级| 侵犯人妻中文字幕一二三四区| 美女中出高潮动态图| 国产成人精品无人区| 自拍欧美九色日韩亚洲蝌蚪91| 婷婷色综合www| 精品亚洲成国产av| 国产av码专区亚洲av| av片东京热男人的天堂| 99九九在线精品视频| 久久久国产欧美日韩av| 午夜福利免费观看在线| 亚洲熟女精品中文字幕| 丁香六月欧美| 亚洲国产精品国产精品| 国产伦人伦偷精品视频| 色综合欧美亚洲国产小说| 久久精品亚洲熟妇少妇任你| 精品亚洲乱码少妇综合久久| 久久人妻熟女aⅴ| 久久久国产一区二区| 麻豆乱淫一区二区| 久久亚洲国产成人精品v| 国产精品成人在线| 青春草亚洲视频在线观看| 亚洲成人一二三区av| 亚洲中文av在线| 精品国产一区二区久久| 日韩欧美精品免费久久| 国产乱人偷精品视频| 免费女性裸体啪啪无遮挡网站| 国产精品亚洲av一区麻豆 | 黄频高清免费视频| 国产亚洲欧美精品永久| 免费黄色在线免费观看| 国产一卡二卡三卡精品 | 色视频在线一区二区三区| 国产成人欧美| 久久久久久久久久久免费av| 亚洲美女视频黄频| 精品少妇久久久久久888优播| 晚上一个人看的免费电影| 亚洲中文av在线| 成人国产av品久久久| av电影中文网址| 国产成人系列免费观看| 纵有疾风起免费观看全集完整版| 欧美中文综合在线视频| 国产精品国产三级国产专区5o| 十八禁网站网址无遮挡| 国产精品秋霞免费鲁丝片| av电影中文网址| 1024香蕉在线观看| 日本一区二区免费在线视频| 又黄又粗又硬又大视频| h视频一区二区三区| 香蕉国产在线看| 久久精品国产综合久久久| 在线观看一区二区三区激情| 午夜福利网站1000一区二区三区| 国产精品一区二区在线不卡| 久久久久久人人人人人| 国产精品一区二区在线观看99| 777米奇影视久久| 丰满饥渴人妻一区二区三| 日韩熟女老妇一区二区性免费视频| 精品久久久久久电影网| 久久久久网色| 熟女av电影| 国产麻豆69| 黄色视频在线播放观看不卡| 1024香蕉在线观看| 美女主播在线视频| 女性被躁到高潮视频| 高清视频免费观看一区二区| 少妇人妻久久综合中文| 欧美精品亚洲一区二区| 在线精品无人区一区二区三| 少妇人妻久久综合中文| 午夜激情av网站| 中文字幕人妻丝袜一区二区 | 国产精品久久久久久久久免| 亚洲成人av在线免费| 日韩中文字幕视频在线看片| 麻豆乱淫一区二区| 香蕉丝袜av| 成人毛片60女人毛片免费| 国产日韩一区二区三区精品不卡| 亚洲图色成人| 国产熟女欧美一区二区| 亚洲av福利一区| 国产一级毛片在线| 香蕉丝袜av| 丰满迷人的少妇在线观看| 乱人伦中国视频| 免费看不卡的av| 天堂中文最新版在线下载| 日韩大片免费观看网站| 一级a爱视频在线免费观看| 久久精品人人爽人人爽视色| 中文精品一卡2卡3卡4更新| 在线观看www视频免费| 国产免费现黄频在线看| 一级片'在线观看视频| 亚洲精品在线美女| 成人亚洲欧美一区二区av| 亚洲综合色网址| 另类精品久久| 亚洲七黄色美女视频| 制服人妻中文乱码| 日日啪夜夜爽| 久久免费观看电影| 毛片一级片免费看久久久久| 欧美变态另类bdsm刘玥| 精品久久久久久电影网| 色婷婷久久久亚洲欧美| 国产精品二区激情视频| 国产成人欧美在线观看 | svipshipincom国产片| 国产欧美日韩一区二区三区在线| 亚洲第一区二区三区不卡| 亚洲情色 制服丝袜| 乱人伦中国视频| 日韩大码丰满熟妇| 亚洲国产最新在线播放| 宅男免费午夜| 亚洲一级一片aⅴ在线观看| 999精品在线视频| 欧美在线一区亚洲| 91精品国产国语对白视频| 久久婷婷青草| 日本vs欧美在线观看视频| 亚洲熟女毛片儿| 大香蕉久久网| 国产视频首页在线观看| 极品少妇高潮喷水抽搐| 午夜日本视频在线| 无限看片的www在线观看| 国产欧美亚洲国产| av一本久久久久| 老熟女久久久| 久久精品亚洲av国产电影网| 亚洲精品国产av成人精品| 在线精品无人区一区二区三| 精品久久蜜臀av无| 国产一级毛片在线| 纵有疾风起免费观看全集完整版| 丝袜美足系列| 黄片小视频在线播放| 亚洲综合色网址| 熟妇人妻不卡中文字幕| 亚洲av国产av综合av卡| 又黄又粗又硬又大视频| 亚洲国产日韩一区二区| 激情五月婷婷亚洲| 国产亚洲av片在线观看秒播厂| av免费观看日本| 国产黄色视频一区二区在线观看| 丝袜美腿诱惑在线| 亚洲精品国产色婷婷电影| 十八禁高潮呻吟视频| 如日韩欧美国产精品一区二区三区| 免费日韩欧美在线观看| 美女高潮到喷水免费观看| 丰满乱子伦码专区| 午夜福利视频在线观看免费| av福利片在线| 高清欧美精品videossex| 久久精品国产亚洲av涩爱| 夫妻午夜视频| 涩涩av久久男人的天堂| 婷婷色综合大香蕉| 国产精品免费视频内射| 香蕉丝袜av| 亚洲成人手机| 老汉色av国产亚洲站长工具| 久久久久精品国产欧美久久久 | 青草久久国产| 丝袜脚勾引网站| 亚洲成人av在线免费| 日日爽夜夜爽网站| 人妻一区二区av| 国产精品一区二区在线不卡| 欧美97在线视频| 日韩一卡2卡3卡4卡2021年| 如何舔出高潮| 又大又爽又粗| 国产女主播在线喷水免费视频网站| 2018国产大陆天天弄谢| 18禁动态无遮挡网站| 免费高清在线观看视频在线观看| 色婷婷久久久亚洲欧美| 人人妻人人澡人人爽人人夜夜| 国产99久久九九免费精品| 国产午夜精品一二区理论片| 一级片'在线观看视频| 极品少妇高潮喷水抽搐| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品,欧美精品| 下体分泌物呈黄色| 免费观看人在逋| av在线播放精品| 国产不卡av网站在线观看| av又黄又爽大尺度在线免费看| 亚洲国产成人一精品久久久| 久久免费观看电影| 热99国产精品久久久久久7| 波野结衣二区三区在线| 看免费av毛片| 亚洲国产精品999| 中文字幕最新亚洲高清| 成年动漫av网址| 成人国产av品久久久| 嫩草影视91久久| 欧美亚洲 丝袜 人妻 在线| 又黄又粗又硬又大视频| 国产精品一区二区在线观看99| 久久精品久久精品一区二区三区| 欧美黑人欧美精品刺激| 成人影院久久| 国产午夜精品一二区理论片| 亚洲精品国产一区二区精华液| 91成人精品电影| 色婷婷av一区二区三区视频| 亚洲国产精品国产精品| 亚洲欧洲国产日韩| 午夜精品国产一区二区电影| 青草久久国产| 在线免费观看不下载黄p国产| 老鸭窝网址在线观看| 欧美人与性动交α欧美软件| 精品国产一区二区三区久久久樱花| 国产免费一区二区三区四区乱码| 伦理电影大哥的女人| 日本爱情动作片www.在线观看| 黄片小视频在线播放| 国产日韩一区二区三区精品不卡| 精品视频人人做人人爽| 又粗又硬又长又爽又黄的视频| 天天影视国产精品| 亚洲精品乱久久久久久| 在线观看一区二区三区激情| 中国三级夫妇交换| 亚洲精品国产av蜜桃| 国产极品粉嫩免费观看在线| 久久久久久人人人人人| 免费观看a级毛片全部| 丝袜在线中文字幕| 99香蕉大伊视频| 欧美日韩亚洲高清精品| 69精品国产乱码久久久| 黄色 视频免费看| 如日韩欧美国产精品一区二区三区| 亚洲熟女毛片儿| 只有这里有精品99| 亚洲欧美激情在线| 汤姆久久久久久久影院中文字幕| 久久久久国产精品人妻一区二区| 亚洲激情五月婷婷啪啪| 一边摸一边抽搐一进一出视频| 久久精品国产亚洲av高清一级| 亚洲四区av| 亚洲av日韩在线播放| 欧美国产精品va在线观看不卡| 高清黄色对白视频在线免费看| 女性被躁到高潮视频| 国产福利在线免费观看视频| a级毛片黄视频| 男女下面插进去视频免费观看| 午夜激情av网站| 汤姆久久久久久久影院中文字幕| 国产精品一二三区在线看| 日本色播在线视频| 中文字幕亚洲精品专区| 中文字幕人妻熟女乱码| 女的被弄到高潮叫床怎么办| 久久久久精品久久久久真实原创| 狠狠精品人妻久久久久久综合| 精品久久久久久电影网| 波多野结衣av一区二区av| 国产精品免费大片| 80岁老熟妇乱子伦牲交| 国产高清国产精品国产三级| 一本一本久久a久久精品综合妖精| 成人漫画全彩无遮挡| 国产一区二区三区综合在线观看| 黄色毛片三级朝国网站| 免费高清在线观看日韩| 高清欧美精品videossex| 亚洲av欧美aⅴ国产| 日韩一区二区三区影片| 亚洲美女黄色视频免费看| 欧美日韩亚洲高清精品| 亚洲美女视频黄频| 亚洲成人免费av在线播放| 国产精品免费大片| 欧美激情 高清一区二区三区| 两个人看的免费小视频| 久久久久精品人妻al黑| 嫩草影院入口| 99国产综合亚洲精品| 2021少妇久久久久久久久久久| 欧美精品人与动牲交sv欧美| 日韩中文字幕欧美一区二区 | 欧美精品一区二区免费开放| 亚洲精品国产av蜜桃| 最新在线观看一区二区三区 | 在线观看www视频免费| 99热国产这里只有精品6| 新久久久久国产一级毛片| 久久久久视频综合| 精品人妻在线不人妻| 天堂8中文在线网| 波野结衣二区三区在线| 日本欧美国产在线视频| 国产精品 国内视频| 热re99久久精品国产66热6| 美女主播在线视频| 国产亚洲av片在线观看秒播厂| 熟女少妇亚洲综合色aaa.| 亚洲国产精品999| 丝袜美足系列| 国产成人精品久久久久久| 国产精品一区二区精品视频观看| 超色免费av| 久久久久久久大尺度免费视频| e午夜精品久久久久久久| 国产精品香港三级国产av潘金莲 | 91老司机精品| 亚洲人成77777在线视频| 男人爽女人下面视频在线观看| 最新的欧美精品一区二区| 欧美成人精品欧美一级黄| 亚洲伊人久久精品综合| 美女中出高潮动态图| 亚洲av电影在线观看一区二区三区| 亚洲中文av在线| 亚洲四区av| 久久99精品国语久久久|