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

    散射輻射對(duì)亞熱帶人工針葉林光能利用率的影響

    2017-06-22 13:47:48王輝民張雷明溫學(xué)發(fā)劉壽東
    生態(tài)學(xué)報(bào) 2017年10期
    關(guān)鍵詞:針葉林亞熱帶光合作用

    衛(wèi) 楠,張 彌,*,王輝民,張雷明,溫學(xué)發(fā),劉壽東

    1 南京信息工程大學(xué)氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室,氣候與環(huán)境變化國(guó)際合作聯(lián)合實(shí)驗(yàn)室,氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,耶魯-南京信息工程大學(xué)大氣環(huán)境研究中心,南京 210044

    ?

    散射輻射對(duì)亞熱帶人工針葉林光能利用率的影響

    衛(wèi) 楠1,張 彌1,*,王輝民2,張雷明2,溫學(xué)發(fā)2,劉壽東1

    1 南京信息工程大學(xué)氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室,氣候與環(huán)境變化國(guó)際合作聯(lián)合實(shí)驗(yàn)室,氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,耶魯-南京信息工程大學(xué)大氣環(huán)境研究中心,南京 210044

    2 中國(guó)科學(xué)院地理科學(xué)與資源研究所,生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測(cè)與模擬重點(diǎn)實(shí)驗(yàn)室,北京 100101

    地表接受的太陽(yáng)輻射中散射輻射的改變是影響森林生態(tài)系統(tǒng)光能利用率(Light Use Efficiency,LUE)的重要因素。以千煙洲亞熱帶人工針葉林為研究對(duì)象,利用30 min通量和常規(guī)氣象觀測(cè)數(shù)據(jù),以晴空指數(shù)(kt)和地表接受的散射輻射(Sf)占太陽(yáng)總輻射(S0)的比值(Sf/S0)為指標(biāo),分析了2003—2012年生長(zhǎng)旺季(6—8月)散射輻射變化對(duì)千煙洲亞熱帶人工針葉林光能利用率的影響,并利用改進(jìn)的光響應(yīng)曲線模型分析了散射輻射變化對(duì)植被光合特性的影響。研究結(jié)果表明:2003—2012年生長(zhǎng)旺季中,kt在0.6—0.7范圍內(nèi)的LUE比kt在0.4—0.5范圍內(nèi)的LUE平均減少了44.66%;Sf/S0在70%—85%之間的LUE比Sf/S0在55%—70%之間的LUE平均提高了22.24%。由此可以看出,與晴朗天空相比,多云及氣溶膠增加導(dǎo)致的散射輻射增加可使該生態(tài)系統(tǒng)的LUE提高。并且,未受到高溫干旱影響的2005、2006、2008、2009、2010及2012年散射輻射下該生態(tài)系統(tǒng)的初始量子效率αf比直接輻射下的αr平均增加了0.63 gCO2/mol;而10年間所有年份散射輻射下的光飽和時(shí)的潛在最大光合作用速率比直接輻射下平均提高了0.81 mgCO2m-2s-1,說(shuō)明散射輻射可使該生態(tài)系統(tǒng)植被光合能力提高。

    光能利用率;散射輻射;晴空指數(shù);千煙洲亞熱帶人工針葉林

    太陽(yáng)輻射是植被光合作用的能量來(lái)源,也是氣候形成及變化的重要驅(qū)動(dòng)因子[1]。然而,當(dāng)天空云量及大氣氣溶膠改變時(shí),會(huì)影響到達(dá)地表的總太陽(yáng)輻射強(qiáng)度及散射輻射與直接輻射之間的比例,從而影響陸地生態(tài)系統(tǒng)的光合生產(chǎn)力[2]。研究發(fā)現(xiàn),近55年來(lái)中國(guó)地面太陽(yáng)輻射的變化特征與全球變化特征相一致,即經(jīng)歷了先逐漸降低變暗,而后逐漸升高變亮的變化過(guò)程[3- 6]。太陽(yáng)輻射的變化尤其是散射輻射比例的改變勢(shì)必會(huì)影響陸地生態(tài)系統(tǒng)對(duì)碳的固定及其光能利用率。

    光能利用率(Light Use Efficiency, LUE)是表征植物利用太陽(yáng)能通過(guò)光合作用轉(zhuǎn)化有機(jī)物的效率[7]。一個(gè)地區(qū)生態(tài)系統(tǒng)生產(chǎn)力的大小與該地區(qū)植被光能利用效率有直接關(guān)系[8-9]。理論上植被葉片中的葉綠體對(duì)光的轉(zhuǎn)化率最高可達(dá)20%—25%,但是自然條件下生長(zhǎng)的植被對(duì)光能的利用率卻不到1%[10-11]。由此可以看出,找到提高光能利用率的途徑成為提高植被生產(chǎn)力的關(guān)鍵問(wèn)題。森林是重要的陸地生態(tài)系統(tǒng),其光能利用率的變化得到很多學(xué)者的關(guān)注。大量研究表明,當(dāng)天空云量及大氣氣溶膠含量增加導(dǎo)致地表接受的太陽(yáng)輻射中散射輻射的比例增加時(shí),會(huì)使得森林生態(tài)系統(tǒng)的光合作用較晴朗天氣條件下高[2,12- 14],從而提高了森林生態(tài)系統(tǒng)的光能利用率及固碳能力。近年來(lái),我國(guó)中東部地區(qū)常出現(xiàn)嚴(yán)重的霧霾天氣[15-17]。霧霾天氣的出現(xiàn)會(huì)導(dǎo)致地表接受的太陽(yáng)總輻射減少[6],而散射輻射的比例增加。因此,在當(dāng)前霧霾天氣經(jīng)常出現(xiàn)的條件下研究散射輻射變化對(duì)森林生態(tài)系統(tǒng)光能利用率的影響變得尤為重要。

    人工林是我國(guó)森林生態(tài)系統(tǒng)的重要組成部分,面積比居世界首位,因此其對(duì)我國(guó)陸地生態(tài)系統(tǒng)的碳收支會(huì)產(chǎn)生巨大影響[18-20]。我國(guó)南方林區(qū)以人工林為主,面積占全國(guó)人工林的54.3%,故在當(dāng)前人為氣溶膠排放增加及氣候變化條件下,研究散射輻射變化對(duì)南方人工林光能利用率的影響能夠?yàn)楦玫卦u(píng)價(jià)人工林生態(tài)系統(tǒng)在區(qū)域碳平衡中的作用提供科學(xué)依據(jù)。

    地處江西省千煙洲的亞熱帶人工針葉林是我國(guó)南方亞熱帶地區(qū)典型的人工林生態(tài)系統(tǒng),大量研究表明,該生態(tài)系統(tǒng)凈CO2交換量(Net Ecosystem Exchange of CO2,NEE)在5—6月出現(xiàn)最高值,而太陽(yáng)輻射強(qiáng)烈少雨的7月份易受高溫干旱影響[18,20- 21]。但是同期在散射輻射提高、氣溫及飽和水汽壓差下降的條件下NEE會(huì)提高[22]。由此可見(jiàn),該生態(tài)系統(tǒng)生產(chǎn)力受光以及溫度、水分等環(huán)境因子的影響。在不受干旱脅迫、溫度適宜的條件下,植被光能利用率決定著該生態(tài)系統(tǒng)的生產(chǎn)力,但是目前該生態(tài)系統(tǒng)光能利用率對(duì)散射輻射變化的響應(yīng),特別是植被在散射輻射及直接輻射條件下光合特性的差異還不明確。本研究以千煙洲亞熱帶人工針葉林為研究對(duì)象,利用2003—2012年的30 minCO2通量及常規(guī)氣象觀測(cè)資料,分析太陽(yáng)輻射中散射輻射的變化對(duì)千煙洲亞熱帶人工針葉林光能利用率的影響,以及光合有效輻射中散射部分和直接輻射部分對(duì)植被光合特性的影響,從而為準(zhǔn)確地評(píng)價(jià)同區(qū)域人工林碳匯能力對(duì)氣候變化的響應(yīng)特征提供科學(xué)依據(jù)。

    1 研究地概況及研究方法

    1.1 研究地概況

    千煙洲亞熱帶人工針葉林位于江西省泰和縣,其通量觀測(cè)塔(26°44′29.1″N,115°03′29.2″E)位于中國(guó)生態(tài)研究網(wǎng)絡(luò)的千煙洲紅壤丘陵農(nóng)業(yè)綜合開(kāi)發(fā)實(shí)驗(yàn)站內(nèi)。千煙洲地處亞熱帶大陸性季風(fēng)氣候區(qū),海拔高度為60—120 m,年平均氣溫約為17.8 ℃,年平均降水量為1360 mm,年平均太陽(yáng)總輻射量為4349 MJ/m2[23]。林內(nèi)主要植被為馬尾松(PinusMassoniana)、濕地松(P.elliottii)、杉木(Cunninghamialanceolata)等,冠層平均高度為12 m,覆蓋率在90%以上。

    1.2 研究方法

    1.2.1 數(shù)據(jù)獲取

    本研究獲取2003—2012年千煙洲站實(shí)時(shí)的30 minCO2通量數(shù)據(jù)及常規(guī)氣象數(shù)據(jù),包括太陽(yáng)總輻射數(shù)據(jù)(S0)、光合有效輻射數(shù)據(jù)(PAR)、溫度(氣溫Ta及土壤5 cm處的溫度Ts)、飽和水汽壓差(VPD)、降水量(P)及土壤表層5 cm處的含水量(SWC)。CO2通量數(shù)據(jù)是由安裝在觀測(cè)塔39 m高度處的渦度相關(guān)系統(tǒng)觀測(cè)得到,該觀測(cè)系統(tǒng)包括三維超聲風(fēng)速儀(Model CSAT- 3, Campbell Scientific)和紅外氣體分析儀(Model Li- 7500, Licor Inc)。儀器原始采樣頻率為10 Hz,數(shù)據(jù)采集器(CR10XTD, CR23XTD, CR5000) 在采集實(shí)時(shí)數(shù)據(jù)的同時(shí)計(jì)算通量數(shù)據(jù)30 min的平均值并進(jìn)行存儲(chǔ)。輻射數(shù)據(jù)來(lái)自42 m,氣溫、水汽壓來(lái)自23.5 m的觀測(cè),詳細(xì)的站點(diǎn)及數(shù)據(jù)材料描述見(jiàn)宋霞等[21]的研究。

    1.2.2 數(shù)據(jù)處理

    在CO2通量觀測(cè)過(guò)程中,由于儀器響應(yīng)誤差、下墊面起伏條件、天氣狀況等因素的影響,會(huì)造成觀測(cè)數(shù)據(jù)出現(xiàn)缺失或異常,因此需要對(duì)獲取的30 min的CO2通量數(shù)據(jù)進(jìn)行處理。具體操作步驟包括:(1)CO2氣體濃度的密度校正;(2)傾斜校正,即坐標(biāo)軸旋轉(zhuǎn);(3)儲(chǔ)存項(xiàng)校正;(4)數(shù)據(jù)篩選和插補(bǔ),剔除由于儀器故障、天氣以及不符合通量觀測(cè)條件等導(dǎo)致的無(wú)效通量數(shù)據(jù),CO2通量數(shù)據(jù)采用查表法進(jìn)行插補(bǔ),常規(guī)氣象數(shù)據(jù)采用滑動(dòng)平均法進(jìn)行缺失數(shù)據(jù)插補(bǔ)[14,18-19,24],最終得到連續(xù)半小時(shí)的凈CO2交換量,即NEE的數(shù)據(jù)。

    1.2.3 NEE組分拆分

    由于渦度相關(guān)法只能直接獲取凈生態(tài)系統(tǒng)CO2交換量(NEE,mgCO2m-2s-1),該量值為正代表生態(tài)系統(tǒng)為碳源,即釋放CO2;該量值為負(fù)代表生態(tài)系統(tǒng)為碳匯,即吸收CO2。生態(tài)系統(tǒng)水平上總初級(jí)生產(chǎn)力(GPP,mgCO2m-2s-1)無(wú)法直接觀測(cè),因此,需對(duì)NEE進(jìn)行組分拆分[25],計(jì)算得到GPP,即使用下式:

    GPP=-NEE+Re

    (1)

    式中,Re(Ecosystem Respiration)表示生態(tài)系統(tǒng)呼吸(mgCO2m-2s-1);Re是利用Lloyd-Taylor方程[26],使用篩選后的全年夜間有效NEE 數(shù)據(jù)與對(duì)應(yīng)有效的土壤5 cm處溫度數(shù)據(jù)對(duì)式2進(jìn)行擬合,并將擬合的方程外推至白天,從而計(jì)算白天的生態(tài)系統(tǒng)呼吸,最終利用式1計(jì)算出GPP。

    Re=RrefeE0(1/(Tref-T0)-1/(Ts-T0))

    (2)

    式中,Ts為土壤5 cm處溫度(K);Rref為生態(tài)系統(tǒng)在參考溫度Tref(283.1 K)下的呼吸值(mgCO2m-2s-1);E0為活化能;T0為常數(shù)(-46.02 ℃,即227.13 K)。

    1.2.4 光能利用率(LUE)

    生態(tài)系統(tǒng)水平的光能利用率為生態(tài)系統(tǒng)總初級(jí)生產(chǎn)力與植被上方接受的光合有效輻射的比值[27-29]:

    (3)

    式中,PAR(photosynthetic active radiation)表示相應(yīng)時(shí)間段生態(tài)系統(tǒng)接受的光合有效輻射量(μmol m-2s-1),因此LUE的單位為(mgCO2/μmol),本文進(jìn)行單位轉(zhuǎn)換,轉(zhuǎn)換為(gCO2/mol)。

    1.2.5 晴空指數(shù)(kt)

    晴空指數(shù)(Clearness index,kt)是指一定太陽(yáng)高度角下地表接受的總太陽(yáng)輻射S0(W/m2) 與地球大氣層上方平行于地表面上接受的總太陽(yáng)輻射Se(W/m2) 的比值[12]:

    (4)

    Se=Ssc[1+0.033cos(360td/365)]sinh

    (5)

    式中,Ssc為太陽(yáng)常數(shù)(1370W/m2),td為日序數(shù),h為太陽(yáng)高度角。晴空指數(shù)是天空晴朗程度的判定指標(biāo),當(dāng)kt接近于0時(shí),天空完全被云層覆蓋,到達(dá)地面的太陽(yáng)輻射非常??;當(dāng)kt=0.5左右時(shí),為多云天空中等輻射條件,太陽(yáng)總輻射下降而散射輻射增加;當(dāng)kt接近于1時(shí),天空晴朗無(wú)云[14]。

    1.2.6 散射輻射

    本文使用散射輻射分割模型簡(jiǎn)化后的形式對(duì)千煙洲人工針葉林接受的散射輻射進(jìn)行估算[2,30-31],方法如下:

    0 ≤kt≤ 0.3;Sf/Se≤kt

    (6)

    Sf/Se=kt(1.020-0.254kt+ 0.0123sinh)

    (7)

    0.3

    (8)

    Sf/Se=kt(1.400-1.749kt+ 0.177sinh)

    (9)

    kt≥0.78;Sf/Se≥ 0.1kt

    (10)

    Sf/Se=kt(0.486kt-0.182sinh)

    (11)

    式中,Sf為地表水平面上接受的散射輻射(W/m2)。

    光合有效輻射的散射輻射部分利用Spitters等[32-33]提出的經(jīng)驗(yàn)公式計(jì)算。

    1.2.7 NEE對(duì)光的響應(yīng)

    本文用Michaelis-Menten方程(式12)對(duì)多云及晴朗天空條件下千煙洲人工林生態(tài)系統(tǒng)的光響應(yīng)曲線進(jìn)行擬合[13-14,34-35]:

    (12)

    式中,α為生態(tài)系統(tǒng)的初始量子效率,表示弱光條件下植被對(duì)光的利用能力(mgCO2/μmol),本文進(jìn)行單位轉(zhuǎn)換,轉(zhuǎn)換為(gCO2/mol);β表示生態(tài)系統(tǒng)潛在最大光合作用速率(mgCO2m-2s-1);由于NEE為負(fù)值表示CO2的吸收,故α和β都為負(fù)值,所以α值的絕對(duì)值越大表明植被弱光下對(duì)光的利用能力越強(qiáng),β值的絕對(duì)值越大表明光飽和時(shí)的潛在最大光合作用速率越大[2]。

    為區(qū)分該森林生態(tài)系統(tǒng)光合作用對(duì)直接輻射、散射輻射的響應(yīng)特征,本文使用Gu等[2]提出的植被光響應(yīng)模型,分析千煙洲人工林生態(tài)系統(tǒng)NEE對(duì)散射輻射與直接輻射的光響應(yīng)特征:

    (13)

    式中,If和Ir分別表示光合有效輻射中散射輻射部分與直接輻射部分,單位:μmol m-2s-1;αf與αr分別表示生態(tài)系統(tǒng)在If和Ir下的光量子初始量子效率(mgCO2/μmol),本文進(jìn)行單位轉(zhuǎn)換,轉(zhuǎn)換為(gCO2/mol);βf與βr分別表示生態(tài)系統(tǒng)在If和Ir下的潛在最大光合作用速率(mgCO2m-2s-1),由于本模型中采用NEE為正值表示CO2吸收,故α、β均為正值。

    本研究采用11天滑動(dòng)法[2]對(duì)生長(zhǎng)旺盛季節(jié)不同時(shí)段的光響應(yīng)曲線進(jìn)行擬合,以得到生長(zhǎng)旺季內(nèi)不同時(shí)段的光合參數(shù),即從6月1日開(kāi)始到8月30日,每11天為一組數(shù)據(jù)進(jìn)行擬合,后一個(gè)11天的開(kāi)頭從前一個(gè)11天的中間開(kāi)始,共擬合17次得出17組不同參數(shù)。

    2 結(jié)果與討論

    2.1 千煙洲亞熱帶人工針葉林環(huán)境因子的季節(jié)變化

    千煙洲亞熱帶人工針葉林氣溫和降水季節(jié)變化如圖1所示,氣溫變化呈單峰型,多年月平均氣溫在7月份達(dá)到最高值,為29.1 ℃,冬季平均最低氣溫出現(xiàn)在1月份,為4.9 ℃。該生態(tài)系統(tǒng)總體夏季降水多,冬季降水少,總體可以看出4—6月氣候濕潤(rùn),其降水量約占全年降水的33.6%,而7—9月相對(duì)干旱,這也與王萌萌等[22]對(duì)2003—2010年以及孫曉敏等[20]對(duì)2003、2004年千煙洲人工林氣溫和降水季節(jié)動(dòng)態(tài)的研究相一致。飽和水汽壓差(VPD)表征環(huán)境大氣的水分虧缺狀況,從圖1可以看出:全年中VPD在7、8月份相對(duì)較高,7月達(dá)到峰值,說(shuō)明在7、8月份,空氣中水分含量相對(duì)較少,比較干旱;土壤表層5 cm處的含水量在1—6月份較多,平均值為0.21 m3/m3,7—12月份較少,平均值為0.14 m3/m3,這與7—12月份降水量較前半年降水量少相一致。

    2003—2012年太陽(yáng)總輻射以及散射輻射季節(jié)變化特征如圖1所示:二者都是夏季大于冬季,在1月份出現(xiàn)最低值,在7月份達(dá)到最大值,太陽(yáng)總輻射最大值為633.8 MJ/m2,散射輻射最大值為318.7 MJ/m2。6月份太陽(yáng)總輻射相比5月份減少而散射輻射卻上升,這是由于當(dāng)月降水增多,多云天空會(huì)導(dǎo)致散射輻射增加。

    綜上所述,該生態(tài)系統(tǒng)光、水、熱資源比較豐富,但是在生長(zhǎng)旺盛季節(jié),即6—8月,受副熱帶高壓控制,7月會(huì)出現(xiàn)高溫少雨的現(xiàn)象,易出現(xiàn)季節(jié)性干旱,對(duì)植被生長(zhǎng)會(huì)造成一定影響。

    圖1 千煙洲亞熱帶人工針葉林2003—2012年氣溫(Ta)與降水(P)、飽和水汽壓差(VPD)與土壤表層5 cm處的含水量(SWC)、太陽(yáng)總輻射(S0)與散射輻射(Sf)的平均季節(jié)變化Fig.1 Averaged seasonal variation of air temperature (Ta) and precipitation (P), vapor pressure deficit (VPD) and soil water content at 5 cm (SWC), global solar radiation (S0) and diffuse radiation (Sf) in Qianyanzhou subtropical plantation coniferous forests from 2003 to 2012

    2.2 千煙洲亞熱帶人工針葉林LUE的季節(jié)變化

    從圖2可以看出,2003—2012年生長(zhǎng)旺盛季節(jié)內(nèi),GPP和LUE都有略微上升的趨勢(shì),但是從P值得出二者都沒(méi)有通過(guò)F檢驗(yàn),這說(shuō)明10年間植被本身生長(zhǎng)變化對(duì)LUE的影響較小,可以忽略。

    光能利用率的季節(jié)變化特征呈雙峰型(圖3),生長(zhǎng)較旺盛季節(jié),兩個(gè)高值分別出現(xiàn)在6月和9月,谷值出現(xiàn)在7月份,這是由于6月份平均GPP為739.4 gCO2/m2,平均PAR為874.3 mol/m2,7月份平均GPP為855.4 gCO2/m2,平均PAR為1257.8 mol/m2,雖然7月份的GPP較高,但7月份PAR相比6月份卻猛增,二者比值減小,故光能利用率低。

    由此可以看出,PAR最高的月份并非LUE最高的季節(jié),而在云雨相對(duì)較多的6月、9月,LUE較高,因此,太陽(yáng)輻射變化尤其是散射輻射比例的改變會(huì)對(duì)LUE產(chǎn)生重要影響。

    圖2 千煙洲亞熱帶人工針葉林2003—2012 年6—8月生態(tài)系統(tǒng)總初級(jí)生產(chǎn)力(GPP)與光能利用率(LUE)年動(dòng)態(tài)Fig.2 Inter-annual variation of gross primary productivity (GPP) and light use efficiency (LUE) in Qianyanzhou subtropical plantation coniferous forests from 2003 to 2012 in June to August

    圖3 千煙洲亞熱帶人工針葉林2003—2012年光能利用率(LUE)的平均季節(jié)變化Fig.3 Averaged seasonal variation of light use efficiency (LUE) in Qianyanzhou subtropical plantation coniferous forests from 2003 to 2012

    2.3 散射輻射變化對(duì)千煙洲亞熱帶人工針葉林LUE的影響

    為去除葉面積指數(shù)變化與太陽(yáng)高度角的改變對(duì)林冠接受太陽(yáng)輻射的影響,本文選擇在該區(qū)域生長(zhǎng)旺盛季節(jié)(6—8月)的最大太陽(yáng)高度角范圍內(nèi)(76.7°—86.7°)分析散射輻射的改變對(duì)LUE的影響,該太陽(yáng)高度角范圍對(duì)應(yīng)的時(shí)段也是一天中植被生理活動(dòng)較為旺盛的時(shí)段。由于各年份變化趨勢(shì)一致,圖4—圖8皆以2011年、2012年為例。

    2.3.1 生長(zhǎng)旺盛季節(jié)LUE隨kt的變化

    如圖4所示,光能利用率隨晴空指數(shù)的變化呈遞減趨勢(shì),kt值越高,天空越晴朗,太陽(yáng)輻射越強(qiáng),LUE越小。2011、2012年kt在0.6—0.7范圍內(nèi)的晴朗天氣條件下平均光能利用率比kt在0.4—0.5范圍內(nèi)的多云條件下光能利用率減小了46.55%、42.95%,其它年份如表1所示,10年中晴朗天空下的LUE較多云天空下平均減小了44.66%。以往研究發(fā)現(xiàn),長(zhǎng)白山闊葉紅松林、鼎湖山亞熱帶常綠闊葉林、西雙版納熱帶雨林森林生態(tài)系統(tǒng)中,kt在0.4—0.6范圍內(nèi)的中等輻射條件下,各生態(tài)系統(tǒng)的光能利用率都會(huì)相應(yīng)的提高[36-37],本文研究結(jié)果與其一致。

    圖4 千煙洲亞熱帶人工針葉林2011、2012年6—8月太陽(yáng)高度角在76.7°—86.7°范圍內(nèi)LUE與kt的關(guān)系Fig.4 Relationship between LUE and kt for solar elevation angle 76.7°—86.7° in Qianyanzhou subtropical plantation coniferous forests from June to August in 2011 and 2012

    2.3.2 生長(zhǎng)旺盛季節(jié)(6—8月)LUE隨散射輻射占總輻射比例的變化

    千煙洲亞熱帶人工針葉林在生長(zhǎng)旺季太陽(yáng)高度角76.7°—86.7°范圍內(nèi),到達(dá)地面的散射輻射(Sf)在太陽(yáng)總輻射中(S0)所占的比例Sf/S0隨kt的增加而減小(圖5)。該結(jié)果表明天空越晴朗,散射輻射占總輻射的比例越小。

    圖5 千煙洲人工針葉林2011、2012年6—8月太陽(yáng)高度角在76.7°—86.7°范圍內(nèi)Sf/S0與kt的關(guān)系Fig.5 Relationship between Sf/S0 and kt for solar elevation angle 76.7°—86.7° in Qianyanzhou subtropical plantation coniferous forests from June to August in 2011 and 2012

    而千煙洲亞熱帶人工針葉林,LUE隨Sf/S0的增大而增大(圖6)。10年間,Sf/S0在70%—85%之間的LUE比Sf/S0在55%—70%之間的平均提高了22.24%(表1)。由此可以看出,在多云及氣溶膠一定程度增多導(dǎo)致總太陽(yáng)輻射減小而散射輻射增加的天氣條件下,LUE增加。

    圖6 千煙洲亞熱帶人工針葉林2011、2012年6—8月太陽(yáng)高度角在76.7°—86.7°范圍內(nèi)LUE與Sf/S0的關(guān)系Fig.6 Relationship between LUE and Sf/S0 (solar elevation 76.7°—86.7°) in Qianyanzhou subtropical plantation coniferous forests from June to August in 2011 and 2012

    光能利用率的增加與散射輻射促使植被光合作用的增加有關(guān)。從圖7中可以看出,2003—2012年GPP隨Sf/S0的增加先增加后減小,最大值出現(xiàn)在60%—70%附近,這是由于散射輻射在林冠中具有更強(qiáng)的透射性[38-39],易入射至林冠下方,因此,散射輻射占太陽(yáng)總輻射增加時(shí),更多的輻射入射到冠層底部的葉片上,使林冠中占比例較大只能接受散射輻射進(jìn)行光合作用的陰生葉的光合作用增加,從而使整個(gè)冠層的光合作用提高[40]。因此,多云天空或氣溶膠一定程度的增加導(dǎo)致散射輻射增加時(shí),會(huì)使植被的光合作用增強(qiáng)。

    表1 2003—2012年千煙洲亞熱帶人工針葉林6—8月太陽(yáng)高度角在76.7°—86.7°范圍內(nèi)LUE與kt以及Sf/S0的定量關(guān)系

    Table 1 Quantitative relationship between LUE andkt, LUE andSf/S0(solar elevation 76.7°—86.7°) in Qianyanzhou subtropical plantation coniferous forests from June to August from 2003 to 2012

    年份Year2003200420052006200720082009201020112012a46.27%45.98%48.48%46.68%30.24%46.28%46.60%46.57%46.55%42.95%b21.49%23.56%21.61%19.11%23.27%38.02%25.71%16.80%22.10%10.76%

    a:kt在0.6—0.7范圍內(nèi)的LUE比kt在0.4—0.5范圍內(nèi)的LUE平均減少的百分比; b:Sf/S0在70%—85%之間的LUE比Sf/S0在55%—70%之間的LUE平均提高的百分比

    圖7 千煙洲亞熱帶人工針葉林2011、2012年6—8月太陽(yáng)高度角在76.7°—86.7°范圍內(nèi)GPP與Sf/S0的關(guān)系Fig.7 Relationship between GPP and Sf/S0 (solar elevation 76.7°—86.7°) in Qianyanzhou subtropical plantation coniferous forests from June to August in 2011 and 2012

    2.4 千煙洲亞熱帶人工針葉林光合作用對(duì)光的響應(yīng)

    2.4.1 不同天空條件下光合作用的光響應(yīng)特征

    利用Michaelis-Menten方程對(duì)2003—2012年生長(zhǎng)旺季(6—8月)多云天氣和晴朗天氣下該生態(tài)系統(tǒng)的光響應(yīng)曲線進(jìn)行擬合發(fā)現(xiàn)(圖8):10年中除2008、2009、2010年之外其他年份多云天氣條件下的初始量子效率α比晴朗天氣條件下平均提高了29.3%(表2)。這說(shuō)明在弱光下多云天氣條件更有利于千煙洲人工針葉林對(duì)光的利用。對(duì)2003—2012年多云天氣條件下和晴朗天氣條件下擬合潛在光合速率參數(shù)進(jìn)行比較發(fā)現(xiàn):除了2008、2009年,其他年份潛在最大光合速率在多云天氣條件下低于晴朗天氣條件下的擬合值(表2)。

    張彌等[14]及Law等[13]的研究發(fā)現(xiàn),光飽和時(shí)潛在最大光合作用速率在多云天空條件下大于晴朗天空條件下的值,與本研究的結(jié)果不一致。這與千煙洲亞熱帶人工針葉林的植被類(lèi)型對(duì)環(huán)境的適應(yīng)性有關(guān)。千煙洲人工林植被類(lèi)型為馬尾松、濕地松、杉木等,梁春等[41]對(duì)馬尾松光響應(yīng)曲線研究表明:1500 μmol m-2s-1光合有效輻射條件下,馬尾松的光合速率達(dá)到最高值,比光合有效輻射在2000 μmol m-2s-1條件下的光合速率高,即強(qiáng)光下,該樹(shù)種會(huì)產(chǎn)生光抑制現(xiàn)象。本文對(duì)光響應(yīng)曲線比較發(fā)現(xiàn),無(wú)論是多云天氣條件還是晴朗天氣條件下,其PAR都會(huì)超過(guò)1500 μmol m-2s-1,即都會(huì)對(duì)植被的光合作用產(chǎn)生抑制作用,故與晴朗天空相比較,千煙洲亞熱帶人工林光飽和時(shí)的潛在光合速率在多云天空條件下并沒(méi)有明顯的增加。

    雖然將多云天氣條件與晴朗天氣條件進(jìn)行區(qū)分,可以得出兩種輻射條件下植被的光合特性,但是該方法并不能完全將散射輻射與直接輻射分開(kāi),為了進(jìn)一步探索散射輻射與直接輻射對(duì)植被光合的影響,將光合有效輻射分為散射有效輻射和直接有效輻射進(jìn)行研究。

    圖8 千煙洲亞熱帶人工針葉林2011、2012年不同天空條件下的光響應(yīng)曲線Fig.8 Light response curves in Qianyanzhou subtropical plantation coniferous forests under clear sky condition and cloudy sky condition in 2011 and 2012

    Table 2 Parameters value of light response curve in Qianyanzhou subtropical plantation coniferous forests under clear sky condition and cloudy sky condition from June to August in the years from 2003 to 2012

    年份Yearα/(gCO2/mol)β/(mgCO2m-2s-1)Re/(mgCO2m-2s-1)相關(guān)性系數(shù)R2晴空條件Clearskycondition多云條件Cloudyskyconditionc/%晴空條件Clearskycondition多云條件Cloudyskyconditiond/%晴空條件Clearskycondition多云條件Cloudyskycondition晴空條件Clearskycondition多云條件Cloudyskycondition2003-1.07-1.3118.32-1.00-0.85-17.830.180.140.660.632004-0.95-1.3630.30-1.74-1.20-44.570.150.150.610.72005-1.04-1.2013.33-1.88-1.18-59.590.250.120.70.742006-0.77-1.4747.70-1.48-1.33-11.030.100.170.670.782007-0.89-1.3131.94-1.31-1.22-7.220.170.170.770.682008-1.59-1.54-3.25-0.88-1.1423.000.110.170.630.732009-1.47-1.35-8.89-1.13-1.162.850.110.140.640.692010-1.37-1.30-5.38-1.30-1.18-9.880.210.150.720.722011-1.16-1.4419.44-1.45-1.23-17.690.240.190.720.762012-0.74-1.3444.42-2.55-1.37-86.190.190.100.660.71

    α: 生態(tài)系統(tǒng)的初始量子效率 The initial quantum efficiency of ecosystem;β:生態(tài)系統(tǒng)潛在最大光合作用速率 The capacity of a canopy to resist photosynthetic saturation at high levels of PAR;Re:生態(tài)系統(tǒng)呼吸 Ecosystem respiration; c表示多云天空條件α比晴朗條件下的α增加的百分比; d表示多云天空條件β比晴朗條件下的β增加的百分比

    2.4.2 植被的光合作用對(duì)散射輻射與直接輻射變化的響應(yīng)特征

    利用公式(13)給出的模型擬合分析了2003—2012年間生長(zhǎng)旺盛季節(jié)光合作用對(duì)散射輻射與直接輻射的響應(yīng)特征。從散射輻射與直接輻射條件下的初始量子效率αf、αr的變化看出(圖9),2005、2006、2008、2009、2010、2012年相類(lèi)似,以2012年為例說(shuō)明,這些年份中70%以上的數(shù)據(jù)點(diǎn)αf>αr,6年中生長(zhǎng)旺盛季節(jié)αf比αr平均高0.63 gCO2/mol,由此說(shuō)明,氣溫降水適宜時(shí),弱光條件下散射輻射比直接輻射更有利于植被對(duì)光的利用。而2003、2004、2007、2011年相類(lèi)似,以2003年為例可以看出6月份以及8月份αf>αr,而在7月份αf<αr,該現(xiàn)象表明,在降水較少的7月,即易發(fā)生高溫干旱的時(shí)期,散射輻射條件下的初始量子效率并未表現(xiàn)出顯著的增加。

    圖9 千煙洲亞熱帶人工針葉林2003、2012年6—8月散射輻射與直接輻射條件下對(duì)應(yīng)的初始量子效率(αf與αr)的對(duì)比Fig.9 The initial quantum efficiency under diffuse radiation condition (αf) and under direct radiation condition (αr) in Qianyanzhou subtropical plantation coniferous forests from June to August in 2003 and 2012

    2003—2012年間(圖10,以2011、2012年為例)擬合所得散射輻射與直接輻射條件下最大潛在光合速率βf和βr相比,10年間生長(zhǎng)旺季散射輻射下的βf比直接輻射下的βr平均提高了0.81 mgCO2m-2s-1,這說(shuō)明散射輻射條件下該人工林生態(tài)系統(tǒng)光合潛力比直接輻射條件下大,即散射輻射條件下,植被的光合潛力增加,這將有利于提高植被的光能利用率。

    本研究中αf和αr、βf和βr的結(jié)果與Gu 等在2002年對(duì)歐洲赤松林、白楊林、混合落葉林、草地和冬小麥作物的研究結(jié)果基本一致。但本研究中初始量子效率在一些年份中出現(xiàn)散射輻射條件下小于直接輻射條件下的現(xiàn)象,而出現(xiàn)較多αf<αr的年份在生長(zhǎng)旺季都出現(xiàn)不同程度的干旱天氣。以2003年為例,如圖11所示,與2012年相比,2003年的7月份土壤含水量偏少,且氣溫高,導(dǎo)致較嚴(yán)重的季節(jié)性伏旱,植被在受到干旱脅迫時(shí),散射輻射條件下,初始量子效率并未表現(xiàn)出顯著的提高。

    10年間生長(zhǎng)旺季βf的平均值比βr的平均值均有所提高(圖10),這是由于散射輻射在林冠中具有更強(qiáng)的透射性[37-38],易入射至林冠下方。因此,在散射輻射增加時(shí),冠層中占比例較大的陰生葉的光合作用增加,同時(shí)也可避免直射輻射過(guò)多導(dǎo)致林冠上部葉片出現(xiàn)光飽和的現(xiàn)象,因此,散射輻射的增加會(huì)使林冠整體上具有更高的光合潛力[42]。

    圖10 千煙洲亞熱帶人工針葉林2011、2012年6—8月散射輻射與直接輻射條件下對(duì)應(yīng)的最大潛在光合速率(βf與βr)的對(duì)比Fig.10 The capacity of a canopy to resist photosynthetic saturation at high levels of PAR under diffuse radiation condition (βf) and under direct radiation condition (βr) in Qianyanzhou subtropical plantation coniferous forests from June to August in 2011 and 2012

    圖11 千煙洲亞熱帶人工針葉林2003、2012年6—8月氣溫與土壤表層5 cm處的含水量配置圖Fig.11 The dynamic of air temperature and soil water content at 5 cm in Qianyanzhou subtropical plantation coniferous forests from June to August in 2003 and 2012

    3 結(jié)論

    通過(guò)對(duì)千煙洲亞熱帶人工針葉林光能利用率隨散射輻射的變化,以及散射輻射與直接輻射條件下植被光合作用光響應(yīng)參數(shù)的特征進(jìn)行分析,并與同類(lèi)型研究的比較討論發(fā)現(xiàn):

    散射輻射更有利于提高該生態(tài)系統(tǒng)的光能利用率。與多云及氣溶膠一定程度增加的中等輻射條件下(0.4

    不受干旱脅迫影響的條件下,植被光合作用對(duì)光的響應(yīng)表現(xiàn)為散射輻射條件下的初始量子效率αf比直接輻射條件下的αr平均提高0.63 gCO2/mol,即弱光條件下散射輻射比直接輻射更有利于植被對(duì)光的利用。在強(qiáng)光下,散射輻射增加更有利于植被光合潛力的增加,從而提高植被的光能利用率。

    太陽(yáng)散射輻射比例改變時(shí),溫度、水分等環(huán)境因子也會(huì)隨之改變,因此要量化太陽(yáng)輻射改變對(duì)植被光能利用率的影響還需借助過(guò)程模型來(lái)加以分析。這也將有助于分析環(huán)境因子脅迫發(fā)生時(shí)直接輻射與散射輻射對(duì)植被光能利用率的影響及其機(jī)制。

    [1] 楊勝朋, 王可麗, 呂世華. 近40年來(lái)中國(guó)大陸總輻射的演變特征. 太陽(yáng)能學(xué)報(bào), 2007, 28(3): 227- 232.

    [2] Gu L H, Baldocchi D, Verma S B, Black T A, Vesala T, Falge E M, Dowty P R. Advantages of diffuse radiation for terrestrial ecosystem productivity. Journal of Geophysical Research, 2002, 107(D6): 4050.

    [3] 馬金玉, 梁宏, 羅勇, 李世奎. 中國(guó)近50年太陽(yáng)直接輻射和散射輻射變化趨勢(shì)特征. 物理學(xué)報(bào), 2011, 60(6): 853- 866.

    [4] Wild M, Gilgen H, Roesch A, Ohmura A, Long C N, Dutton E G, Forgan B, Kallis A, Russak V, Tsvetkov A. From dimming to brightening: decadal changes in solar radiation at earth′s surface. Science, 2005, 308(5723): 847- 850.

    [5] Mercado L M, Bellouin N, Sitch S, Boucher O, Huntingford C, Wild M, Cox P M. Impact of changes in diffuse radiation on the global land carbon sink. Nature, 2009, 458(7241): 1014- 1071.

    [6] 齊月, 房世波, 周文佐. 近50年來(lái)中國(guó)東、西部地面太陽(yáng)輻射變化及其與大氣環(huán)境變化的關(guān)系. 物理學(xué)報(bào), 2015, 64(8): 089201- 1-089201- 10.

    [7] 衛(wèi)亞星, 王莉雯. 青海省植被光能利用率模擬研究. 生態(tài)學(xué)報(bào), 2010, 30(19): 5209- 5216.

    [8] 李澤暉, 王云龍, 魏遠(yuǎn), 馬悅, 陳潤(rùn)芝. 湖南岳陽(yáng)楊樹(shù)人工林光能利用率動(dòng)態(tài)特征分析. 資源科學(xué), 2012, 34(10): 1832- 1838.

    [9] 穆少杰, 李建龍, 周偉, 楊紅飛, 章超斌, 居為民. 2001—2010年內(nèi)蒙古植被凈初級(jí)生產(chǎn)力的時(shí)空格局及其與氣候的關(guān)系. 生態(tài)學(xué)報(bào), 2013, 33(12): 3752- 3764.

    [10] 解鋒. 論提高光能利用率的途徑. 陜西農(nóng)業(yè)科學(xué), 2008, 54(1): 144- 146.

    [11] 徐恒永, 趙君實(shí), 徐茂臻. 高產(chǎn)冬小麥光能利用率的研究. 萊陽(yáng)農(nóng)學(xué)院學(xué)報(bào), 1994, 11(1): 6- 10.

    [12] Gu L H, Fuentes J D, Shugart H H, Staebler R M, Black T A. Responses of net ecosystem exchanges of carbon dioxide to changes in cloudiness: results from two North American deciduous forests. Journal of Geophysical Research: Atmospheres, 1999, 104(D24): 31421- 31434.

    [13] Law B E, Falge E, Gu L, Baldocchi D D, Bakwin P, Berbigier P, Davis K, Dolman A J, Falk M, Fuentes J D, Goldstein A, Granier A, Grelle A, Hollinger D, Janssens I A, Jarvis P, Jensen N O, Katul G, Mahli Y, Matteucci G, Meyers T, Monson R, Munger W, Oechel W, Olson R, Pilegaard K, Paw U K T, Thorgeirsson H, Valentini R, Verma S, Vesala T, Wilson K, Wofsy S. Environmental controls over carbon dioxide and water vapor exchange of terrestrial vegetation. Agricultural and Forest Meteorology, 2002, 113(1/4): 97- 120.

    [14] 張彌, 于貴瑞, 張雷明, 孫曉敏, 溫學(xué)發(fā), 韓士杰. 太陽(yáng)輻射對(duì)長(zhǎng)白山闊葉紅松林凈生態(tài)系統(tǒng)碳交換的影響. 植物生態(tài)學(xué)報(bào), 2009, 33(2): 270- 282.

    [15] 丁一匯, 柳艷菊. 近50年我國(guó)霧和霾的長(zhǎng)期變化特征及其與大氣濕度的關(guān)系. 中國(guó)科學(xué): 地球科學(xué), 2014, 44(1): 37- 48.

    [16] 張人禾, 李強(qiáng), 張若楠. 2013年1月中國(guó)東部持續(xù)性強(qiáng)霧霾天氣產(chǎn)生的氣象條件分析. 中國(guó)科學(xué): 地球科學(xué), 2014, 44(1): 27- 36.

    [17] 吳兌. 近十年中國(guó)灰霾天氣研究綜述. 環(huán)境科學(xué)學(xué)報(bào), 2012, 32(2):257-269.

    [18] 劉允芬, 于貴瑞, 溫學(xué)發(fā), 王迎紅, 宋霞, 李菊, 孫曉敏, 楊風(fēng)亭, 陳永瑞, 劉琪珠. 千煙洲中亞熱帶人工林生態(tài)系統(tǒng)CO2通量的季節(jié)變異特征. 中國(guó)科學(xué) D輯 地球科學(xué), 2006, 34(S1): 91- 102.

    [19] 劉允芬, 宋霞, 孫曉敏, 溫學(xué)發(fā), 陳永瑞. 千煙洲人工針葉林CO2通量季節(jié)變化及其環(huán)境因子的影響. 中國(guó)科學(xué) D輯 地球科學(xué), 2004, 34(S2): 109- 117.

    [20] 孫曉敏, 溫學(xué)發(fā), 于貴瑞, 劉允芬, 劉琪璟. 中亞熱帶季節(jié)性干旱對(duì)千煙洲人工林生態(tài)系統(tǒng)碳吸收的影響. 中國(guó)科學(xué) D輯 地球科學(xué), 2006, 36(S1): 103- 110.

    [21] 宋霞, 劉允芬, 徐小鋒, 于貴瑞, 溫學(xué)發(fā). 紅壤丘陵區(qū)人工林冬春時(shí)段碳、水、熱通量的觀測(cè)與分析. 資源科學(xué), 2004, 26(3): 96- 104.

    [22] 王萌萌, 張彌, 王輝民, 張雷明, 楊沈斌. 太陽(yáng)輻射變化對(duì)千煙洲亞熱帶人工針葉林凈CO2交換量的影響. 生態(tài)學(xué)雜志, 2015, 34(2): 303- 311.

    [23] 韓佳音, 李勝功, 張雷明, 溫學(xué)發(fā), 李慶康, 王輝民. 江西千煙洲地區(qū)散射輻射的模擬與驗(yàn)證. 應(yīng)用生態(tài)學(xué)報(bào), 2015, 26(10): 2991- 2999.

    [24] 張雷明, 于貴瑞, 孫曉敏, 溫學(xué)發(fā), 任傳友, 宋霞, 劉允芬, 關(guān)德新, 閏俊華, 張一平. 中國(guó)東部森林樣帶典型生態(tài)系統(tǒng)碳收支的季節(jié)變化. 中國(guó)科學(xué) D輯 地球科學(xué), 2006, 36(S1): 45- 59.

    [25] 黃昆, 王紹強(qiáng), 王輝民, 儀垂祥, 周蕾, 劉允芬, 石浩. 中亞熱帶人工針葉林生態(tài)系統(tǒng)碳通量拆分差異分析. 生態(tài)學(xué)報(bào), 2013, 33(17): 5252- 5265.

    [26] Lloyd J, Taylor J A. On the temperature dependence of soil respiration. Functional Ecology, 1994, 8(3): 315- 323.

    [27] Gilmanov T G, Tieszen L L, Wylie B K, Flanagan L B, Frank A B, Haferkamp M R, Meyers T P, Morgan J A. Integration of CO2flux and remotely-sensed data for primary production and ecosystem respiration analyses in the Northern Great Plains: potential for quantitative spatial extrapolation. Global Ecology and Biogeography, 2005, 14(3): 271- 292.

    [28] 趙育民, 牛樹(shù)奎, 王軍邦, 李海濤, 李貴才. 植被光能利用率研究進(jìn)展. 生態(tài)學(xué)雜志, 2007, 26(9): 1471- 1477.

    [29] Gilmanov T G, Soussana J F, Aires L, Allard V, Ammann C, Balzarolo M, Barcza Z, Bernhofer C, Campbell C L, Cernusca A, Cescatti A, Clifton-Brown J, Dirks B O M, Dore S, Eugster W, Fuhrer J, Gimeno C, Gruenwald T, Haszpra L, Hensen A, Ibrom A, Jacobs A F G, Jones M B, Lanigan G, Laurila T, Lohila A, Manca G, Marcolla B, Nagy Z, Pilegaard K, Pinter K, Pio C, Raschi A, Rogiers N, Sanz M J, Stefani P, Sutton M, Tuba Z, Valentini R, Williams M L, Wohlfahrt G. Partitioning European grassland net ecosystem CO2exchange into gross primary productivity and ecosystem respiration using light response function analysis. Agriculture, Ecosystems & Environment, 2007, 121(1/2): 93- 120.

    [30] 孫敬松, 周廣勝. 散射輻射測(cè)量及其對(duì)陸地生態(tài)系統(tǒng)生產(chǎn)力影響的研究進(jìn)展. 植物生態(tài)學(xué)報(bào), 2010, 34(4): 452- 461.

    [31] Reindl D T, Beckman W A, Duffie J A. Diffuse fraction correlations. Solar Energy, 1990, 45(1): 1- 7.

    [32] Spitters C J T. Separating the diffuse and direct component of global radiation and its implications for modeling canopy photosynthesis part II. Calculation of canopy photosynthesis. Agricultural and Forest Meteorology, 1986, 38(1/3): 231- 242.

    [33] Spitters C J T, Toussaint H A J M, Goudriaan J. Separating the diffuse and direct component of global radiation and its implications for modeling canopy photosynthesis part I. Components of incoming radiation. Agricultural and Forest Meteorology, 1986, 38(1/3): 217- 229.

    [34] Goulden M L, Daube B C, Fan S M, Sutton D J, Bazzaz A, Munger J W, Wofsy S C. Physiological responses of a black spruce forest to weather. Journal of Geophysical Research: Atmospheres, 1997, 102(D24): 28987- 28996.

    [35] Aubinet M, Chermanne B, Vandenhaute M, Longdoz B, Yernaux M, Laitat E. Long term carbon dioxide exchange above a mixed forest in the Belgian Ardennes. Agricultural and Forest Meteorology, 2001, 108(4): 293- 315.

    [36] 張彌, 吳家兵, 關(guān)德新, 施婷婷, 陳鵬獅, 紀(jì)瑞鵬. 長(zhǎng)白山闊葉紅松林主要樹(shù)種光合作用的光響應(yīng)曲線. 應(yīng)用生態(tài)學(xué)報(bào), 2006, 17(9): 1575- 1578.

    [37] Zhang M, Yu G R, Zhuang J, Gentry R, Fu Y L, Sun X M, Zhang L M, Wen X F, Wang Q F, Han S J, Yan J H, Zhang Y P, Wang Y F, Li Y N. Effects of cloudiness change on net ecosystem exchange, light use efficiency, and water use efficiency in typical ecosystems of China. Agricultural and Forest Meteorology, 2011, 151(7): 803- 816.

    [38] Roderick M L, Farquhar G D, Berry S L, Noble I R. On the direct effect of clouds and atmospheric particles on the productivity and structure of vegetation. Oecologia, 2001, 129(1): 21- 30.

    [39] Yamasoe M A, von Randow C, Manzi A O, Schafer J S, Eck T F, Holben B N. Effect of smoke and clouds on the transmissivity of photosynthetically active radiation inside the canopy. Atmospheric Chemistry and Physics, 2006, (6): 1645- 1656.

    [40] Williams M, Rastetter E B, var der Pol L, Shaver G R. Arctic canopy photosynthetic efficiency enhanced under diffuse light, linked to a reduction in the fraction of the canopy in deep shade. New Phytologist, 2014, 202(4): 1267- 1276.

    [41] 梁春, 林植芳, 孔國(guó)輝. 不同光強(qiáng)下生長(zhǎng)的亞熱帶樹(shù)苗的光合-光響應(yīng)特性的比較. 應(yīng)用生態(tài)學(xué)報(bào), 1997, 8(1): 7- 11.

    [42] 何學(xué)兆, 周濤, 賈根鎖, 張自銀, 李秀娟, 趙超, 馮勝輝. 光合有效輻射總量及其散射輻射比例變化對(duì)森林GPP影響的模擬. 自然資源學(xué)報(bào), 2011, 26(4): 619- 634.

    The impacts of changes in diffuse radiation on light use efficiency in a subtropical plantation coniferous forest

    WEI Nan1, ZHANG Mi1,*, WANG Huimin2, ZHANG Leiming2, WEN Xuefa2, LIU Shoudong1

    1KeyLaboratoryofMeteorologicalDisaster,MinistryofEducation,InternationalJointLaboratoryonClimateandEnvironmentChange,CollaborativeInnovationCenteronForecastandEvaluationofMeteorologicalDisasters,Yale-NUISTCenteronAtmosphericEnvironment,NanjingUniversityofInformationScienceandTechnology,Nanjing210044,China2KeyLaboratoryofEcosystemNetworkObservationandModeling,SynthesisResearchCenterofChineseEcosystemResearchNetwork,InstituteofGeographicSciencesandNaturalResourcesResearch,ChineseAcademyofSciences,Beijing100101,China

    Changes in diffuse radiation

    at the ground surface impact light use efficiency (LUE) of the forest ecosystem. In this study, we focused on a subtropical plantation coniferous forest at Qianyanzhou. Based on 30 min flux data, routine meteorological data, and two indexes, clearness index (kt) and the ratio (Sf/S0) of diffuse radiation (Sf) to global solar radiation (S0), we analyzed the impacts of changes in diffuse radiation on light use efficiency in this ecosystem during the mid-growing season (June to August) from 2003 to 2012. Moreover, applying the modified light response model, we investigated effects of changes in diffuse radiation on vegetation photosynthetic characteristics. The results indicated that the LUE decreased by a mean of 44.66% inktrange of 0.6 to 0.7 during the mid-growing season from 2003 to 2012,compared that inktrange of 0.4 to 0.5. In the same period, the LUE increased by a mean of 22.24% inSf/S0range of 70% to 85%, compared to that inSf/S0range of 55% to 70%. These results indicated that LUE was higher under cloudy skies and aerosol loading increase conditions than that under clear sky conditions. Moreover, the initial quantum efficiency (αf) increased by a mean of 0.63 gCO2/mol under diffuse radiation, comparing with that (αr) under direct radiation in 2005, 2006, 2008, 2009, 2010, and 2012 during the mid-growing season without drought stress. Regardless of where the water stress occurred, the capacity of a canopy to resist photosynthetic saturation at high levels of PAR (βf) was greater under diffuse radiation thanβrunder direct radiation by a mean of 0.81 mgCO2m-2s-1. The results verified that diffuse radiation improved the photosynthetic capacity of the forest ecosystem.

    light use efficiency; diffuse radiation; clearness index; subtropical plantation coniferous forests at Qianyanzhou

    國(guó)家自然科學(xué)基金青年基金(31200377);中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專(zhuān)項(xiàng)項(xiàng)目(XDA05050208);高等學(xué)校博士學(xué)科點(diǎn)專(zhuān)項(xiàng)科研基金新教師類(lèi)課題(20123228120003)

    2016- 03- 08; 網(wǎng)絡(luò)出版日期:2017- 02- 17

    10.5846/stxb201603080411

    *通訊作者Corresponding author.E-mail: zhangm.80@nuist.edu.cn

    衛(wèi)楠,張彌,王輝民,張雷明,溫學(xué)發(fā),劉壽東.散射輻射對(duì)亞熱帶人工針葉林光能利用率的影響.生態(tài)學(xué)報(bào),2017,37(10):3403- 3414.

    Wei N, Zhang M, Wang H M, Zhang L M, Wen X F, Liu S D.The impacts of changes in diffuse radiation on light use efficiency in a subtropical plantation coniferous forest.Acta Ecologica Sinica,2017,37(10):3403- 3414.

    猜你喜歡
    針葉林亞熱帶光合作用
    金沙江上游與其他流域典型植被碳含量差異性探究
    光合作用研究動(dòng)向
    科學(xué)(2022年4期)2022-10-25 02:43:00
    《亞熱帶農(nóng)業(yè)研究》2019年征稿簡(jiǎn)則
    大興安嶺主要針葉樹(shù)種移植容器苗的培育技術(shù)與造林效果的研究
    施肥對(duì)油松人工林針葉養(yǎng)分變化的影響
    結(jié)合Sentinel- 1B和Landsat8數(shù)據(jù)的針葉林葉片含水量反演研究
    森林工程(2018年4期)2018-08-04 03:23:16
    戶(hù)撒刀
    家居布藝的光合作用
    Coco薇(2016年5期)2016-06-03 09:17:41
    例析凈光合作用相關(guān)題型解題中的信息提取錯(cuò)誤與矯正
    考試周刊(2016年6期)2016-03-11 08:14:32
    愛(ài)上光合作用
    老司机午夜福利在线观看视频| 日日干狠狠操夜夜爽| 18禁黄网站禁片免费观看直播| 国产亚洲精品第一综合不卡| 每晚都被弄得嗷嗷叫到高潮| 悠悠久久av| 在线观看66精品国产| 天天一区二区日本电影三级| 日韩精品青青久久久久久| 国产精品一及| 桃红色精品国产亚洲av| 久久久精品大字幕| 99久久久亚洲精品蜜臀av| 丰满人妻一区二区三区视频av | 久久久久免费精品人妻一区二区| 久久欧美精品欧美久久欧美| videosex国产| xxxwww97欧美| 亚洲成av人片在线播放无| 亚洲国产高清在线一区二区三| 久久婷婷人人爽人人干人人爱| 99国产精品一区二区蜜桃av| √禁漫天堂资源中文www| 在线十欧美十亚洲十日本专区| 精品电影一区二区在线| 中文字幕熟女人妻在线| 亚洲免费av在线视频| 亚洲国产欧美网| 91字幕亚洲| 亚洲av成人av| 好男人在线观看高清免费视频| 九色国产91popny在线| 精品不卡国产一区二区三区| 国产精品一区二区三区四区久久| av福利片在线| 久久久久免费精品人妻一区二区| 国产真人三级小视频在线观看| 午夜亚洲福利在线播放| 成人三级做爰电影| 久热爱精品视频在线9| 欧美黑人精品巨大| 一边摸一边抽搐一进一小说| 正在播放国产对白刺激| 身体一侧抽搐| 91成年电影在线观看| 免费在线观看黄色视频的| 99精品在免费线老司机午夜| 日韩免费av在线播放| 变态另类丝袜制服| 国内精品久久久久久久电影| 欧美日韩国产亚洲二区| 男人舔女人下体高潮全视频| a在线观看视频网站| 正在播放国产对白刺激| 一区二区三区高清视频在线| 变态另类丝袜制服| 首页视频小说图片口味搜索| 欧美乱色亚洲激情| 亚洲 欧美 日韩 在线 免费| 18禁观看日本| 少妇人妻一区二区三区视频| av国产免费在线观看| 91麻豆精品激情在线观看国产| 国产欧美日韩一区二区三| 哪里可以看免费的av片| 俄罗斯特黄特色一大片| 国产成人影院久久av| 可以免费在线观看a视频的电影网站| 精品国产乱子伦一区二区三区| 亚洲午夜理论影院| 国产一区二区三区视频了| 国产精品久久久av美女十八| 欧美黄色片欧美黄色片| 国产高清视频在线观看网站| 久久香蕉国产精品| 久久精品国产亚洲av高清一级| 日本一二三区视频观看| 桃色一区二区三区在线观看| svipshipincom国产片| 国产精品野战在线观看| 国产精品一区二区三区四区免费观看 | a在线观看视频网站| 亚洲aⅴ乱码一区二区在线播放 | 在线观看舔阴道视频| 韩国av一区二区三区四区| 亚洲电影在线观看av| 国产精品影院久久| 一进一出好大好爽视频| 国产精品爽爽va在线观看网站| 18禁美女被吸乳视频| 亚洲精品久久国产高清桃花| av片东京热男人的天堂| 91大片在线观看| 亚洲人成网站在线播放欧美日韩| 国产精品九九99| 精品无人区乱码1区二区| 女人爽到高潮嗷嗷叫在线视频| 此物有八面人人有两片| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 人人妻人人澡欧美一区二区| 中国美女看黄片| 黑人欧美特级aaaaaa片| 婷婷精品国产亚洲av在线| 91av网站免费观看| 一本精品99久久精品77| 波多野结衣高清无吗| 精品国产乱码久久久久久男人| 成人三级做爰电影| 日本黄色视频三级网站网址| 极品教师在线免费播放| 亚洲欧洲精品一区二区精品久久久| 久久久久久免费高清国产稀缺| 久久久水蜜桃国产精品网| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品一区二区www| 美女午夜性视频免费| a级毛片在线看网站| 最近最新中文字幕大全电影3| 久久精品91蜜桃| 女生性感内裤真人,穿戴方法视频| 久久久精品国产亚洲av高清涩受| 午夜激情av网站| 国产精品,欧美在线| 欧美人与性动交α欧美精品济南到| 欧美高清成人免费视频www| 亚洲国产精品999在线| 色播亚洲综合网| 久久久久精品国产欧美久久久| 999久久久精品免费观看国产| 女生性感内裤真人,穿戴方法视频| 成年免费大片在线观看| 国产精品久久久久久久电影 | 搡老岳熟女国产| 最近最新免费中文字幕在线| 麻豆一二三区av精品| 亚洲成人精品中文字幕电影| 1024视频免费在线观看| 亚洲成av人片在线播放无| 99re在线观看精品视频| 欧美日韩一级在线毛片| 成年版毛片免费区| 午夜福利免费观看在线| 国产片内射在线| 亚洲成人久久性| 俄罗斯特黄特色一大片| 久久久久久久久免费视频了| 村上凉子中文字幕在线| 嫁个100分男人电影在线观看| 一级毛片精品| 国产主播在线观看一区二区| 三级毛片av免费| 又黄又爽又免费观看的视频| 黄色视频不卡| 日本三级黄在线观看| 国产99白浆流出| 国产精品自产拍在线观看55亚洲| 日韩欧美国产在线观看| 99国产极品粉嫩在线观看| 日本免费a在线| 在线播放国产精品三级| 欧美日本亚洲视频在线播放| 成年人黄色毛片网站| 国产亚洲精品综合一区在线观看 | 黄色a级毛片大全视频| 国产三级中文精品| 一本久久中文字幕| 老司机深夜福利视频在线观看| 国产免费av片在线观看野外av| 99在线视频只有这里精品首页| 韩国av一区二区三区四区| 久久婷婷人人爽人人干人人爱| 午夜精品久久久久久毛片777| 国产成+人综合+亚洲专区| 国产精品久久久久久精品电影| 最近在线观看免费完整版| 香蕉久久夜色| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 色哟哟哟哟哟哟| 香蕉丝袜av| 久久久久精品国产欧美久久久| 18禁观看日本| 男人的好看免费观看在线视频 | 男女视频在线观看网站免费 | 黄频高清免费视频| 99精品在免费线老司机午夜| 少妇被粗大的猛进出69影院| 日本撒尿小便嘘嘘汇集6| 欧美日韩精品网址| 中文字幕熟女人妻在线| 50天的宝宝边吃奶边哭怎么回事| 久久精品夜夜夜夜夜久久蜜豆 | 国产视频一区二区在线看| 久久精品综合一区二区三区| 欧美大码av| 亚洲精品一区av在线观看| 久久精品成人免费网站| 欧美性长视频在线观看| 亚洲无线在线观看| 淫妇啪啪啪对白视频| 亚洲美女黄片视频| 日韩高清综合在线| 亚洲国产精品合色在线| 最近视频中文字幕2019在线8| 在线播放国产精品三级| av国产免费在线观看| 最近在线观看免费完整版| 日本熟妇午夜| 欧美中文综合在线视频| 人人妻,人人澡人人爽秒播| 久久精品91无色码中文字幕| 天天一区二区日本电影三级| 精品久久久久久久毛片微露脸| 高清在线国产一区| 五月伊人婷婷丁香| 欧美色视频一区免费| 搞女人的毛片| 久久久久久大精品| 夜夜躁狠狠躁天天躁| 国产在线观看jvid| 免费在线观看成人毛片| 999久久久精品免费观看国产| 一个人免费在线观看电影 | 少妇粗大呻吟视频| 这个男人来自地球电影免费观看| 国产精品乱码一区二三区的特点| 国产一区二区在线av高清观看| 精品久久久久久成人av| 免费看a级黄色片| 国产精华一区二区三区| 国语自产精品视频在线第100页| 久久久久久国产a免费观看| 亚洲第一电影网av| 国产亚洲精品第一综合不卡| bbb黄色大片| 一二三四在线观看免费中文在| 成在线人永久免费视频| 日韩大尺度精品在线看网址| 香蕉av资源在线| xxxwww97欧美| 久久精品人妻少妇| 1024手机看黄色片| 熟女电影av网| 日本a在线网址| 国产激情欧美一区二区| 国产私拍福利视频在线观看| 国产精品 欧美亚洲| 午夜免费成人在线视频| 久99久视频精品免费| 亚洲欧美精品综合一区二区三区| 成人午夜高清在线视频| 国产精品亚洲美女久久久| 亚洲精品美女久久久久99蜜臀| 91麻豆av在线| 亚洲欧美日韩无卡精品| 午夜激情福利司机影院| 99久久久亚洲精品蜜臀av| 中亚洲国语对白在线视频| 亚洲成av人片在线播放无| 亚洲男人天堂网一区| 老司机福利观看| 午夜福利在线在线| bbb黄色大片| 日本黄色视频三级网站网址| 国产成人精品久久二区二区91| 国产1区2区3区精品| 在线永久观看黄色视频| 国产亚洲欧美98| 国产精品综合久久久久久久免费| 成人手机av| 欧美激情久久久久久爽电影| 国产亚洲精品久久久久5区| 免费一级毛片在线播放高清视频| 十八禁网站免费在线| 欧美中文综合在线视频| av在线播放免费不卡| 免费在线观看亚洲国产| 999精品在线视频| 日本 欧美在线| 欧美日韩国产亚洲二区| 天天添夜夜摸| 男男h啪啪无遮挡| 日本精品一区二区三区蜜桃| 成人18禁高潮啪啪吃奶动态图| 97人妻精品一区二区三区麻豆| 精品一区二区三区四区五区乱码| 亚洲电影在线观看av| 亚洲国产精品成人综合色| netflix在线观看网站| 19禁男女啪啪无遮挡网站| 中文在线观看免费www的网站 | 怎么达到女性高潮| 国产高清视频在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 丁香欧美五月| 久久久国产成人精品二区| 久久精品人妻少妇| 99久久国产精品久久久| 大型黄色视频在线免费观看| 女警被强在线播放| ponron亚洲| 青草久久国产| 精品不卡国产一区二区三区| 一边摸一边抽搐一进一小说| 精品不卡国产一区二区三区| 一级片免费观看大全| 在线十欧美十亚洲十日本专区| 精品久久久久久久毛片微露脸| 99精品在免费线老司机午夜| 亚洲熟妇熟女久久| 国产激情欧美一区二区| 天堂动漫精品| 日韩大尺度精品在线看网址| 老汉色av国产亚洲站长工具| 一级毛片女人18水好多| АⅤ资源中文在线天堂| 国产精品电影一区二区三区| 久久香蕉国产精品| 中文字幕久久专区| 黄色视频不卡| 亚洲av中文字字幕乱码综合| 国产三级黄色录像| 真人做人爱边吃奶动态| 国产一区二区激情短视频| 亚洲av成人精品一区久久| 午夜免费观看网址| 国产又黄又爽又无遮挡在线| 国产私拍福利视频在线观看| 黑人欧美特级aaaaaa片| 国产精品一区二区三区四区久久| 高潮久久久久久久久久久不卡| 1024香蕉在线观看| 黑人欧美特级aaaaaa片| 亚洲国产欧洲综合997久久,| 中文字幕久久专区| 免费高清视频大片| 国产探花在线观看一区二区| 国产一区二区激情短视频| 国产99久久九九免费精品| 深夜精品福利| 免费在线观看日本一区| 九色成人免费人妻av| 国产精品一区二区免费欧美| 十八禁网站免费在线| 亚洲av成人精品一区久久| 18禁国产床啪视频网站| 亚洲在线自拍视频| 九色成人免费人妻av| 国产探花在线观看一区二区| 国产亚洲欧美在线一区二区| 99久久99久久久精品蜜桃| 午夜福利免费观看在线| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久毛片微露脸| 一级毛片高清免费大全| 舔av片在线| 亚洲成av人片在线播放无| 精品久久久久久久毛片微露脸| 中文字幕高清在线视频| 欧美日韩中文字幕国产精品一区二区三区| 人人妻人人看人人澡| 亚洲国产欧美人成| 国产亚洲精品av在线| 中文字幕熟女人妻在线| 欧美绝顶高潮抽搐喷水| 亚洲欧美一区二区三区黑人| 手机成人av网站| 日韩av在线大香蕉| 国产亚洲av高清不卡| 在线观看午夜福利视频| 国产99白浆流出| 久久午夜亚洲精品久久| 久久久久久人人人人人| 老司机午夜十八禁免费视频| 亚洲美女黄片视频| 校园春色视频在线观看| 黄色女人牲交| 最近最新中文字幕大全电影3| 国产精品,欧美在线| 男女那种视频在线观看| 免费高清视频大片| 久久久久久久久中文| 听说在线观看完整版免费高清| 亚洲午夜理论影院| 亚洲七黄色美女视频| 又黄又粗又硬又大视频| cao死你这个sao货| 中文亚洲av片在线观看爽| 一边摸一边抽搐一进一小说| 999精品在线视频| 久久久国产成人免费| 午夜福利高清视频| 桃色一区二区三区在线观看| 99re在线观看精品视频| 中文字幕最新亚洲高清| 老汉色av国产亚洲站长工具| 人人妻人人看人人澡| 热99re8久久精品国产| 此物有八面人人有两片| 可以免费在线观看a视频的电影网站| 五月伊人婷婷丁香| 亚洲欧洲精品一区二区精品久久久| 91麻豆精品激情在线观看国产| 精品第一国产精品| 亚洲精品国产精品久久久不卡| 亚洲中文日韩欧美视频| 黄频高清免费视频| 一二三四在线观看免费中文在| 18禁裸乳无遮挡免费网站照片| 在线观看www视频免费| 真人一进一出gif抽搐免费| 老司机靠b影院| 国产亚洲av嫩草精品影院| 免费在线观看影片大全网站| x7x7x7水蜜桃| www日本黄色视频网| 国产成人啪精品午夜网站| 1024香蕉在线观看| 日日摸夜夜添夜夜添小说| 嫩草影院精品99| 熟妇人妻久久中文字幕3abv| 亚洲精华国产精华精| 成人特级黄色片久久久久久久| 亚洲色图 男人天堂 中文字幕| 亚洲专区中文字幕在线| 国产区一区二久久| 久久这里只有精品中国| 久久国产精品人妻蜜桃| 国产成年人精品一区二区| 一级毛片女人18水好多| av在线天堂中文字幕| 757午夜福利合集在线观看| 久久天躁狠狠躁夜夜2o2o| 日韩中文字幕欧美一区二区| bbb黄色大片| 色在线成人网| 国内久久婷婷六月综合欲色啪| 黄片小视频在线播放| 精品国内亚洲2022精品成人| 90打野战视频偷拍视频| 欧美激情久久久久久爽电影| 日韩欧美国产在线观看| 久久精品国产清高在天天线| 亚洲精品av麻豆狂野| 国产精品av久久久久免费| 在线十欧美十亚洲十日本专区| 日本 欧美在线| 亚洲中文字幕一区二区三区有码在线看 | 热99re8久久精品国产| 久久久国产欧美日韩av| 国内久久婷婷六月综合欲色啪| 亚洲18禁久久av| 看黄色毛片网站| 日本免费a在线| 久热爱精品视频在线9| 欧美色视频一区免费| 午夜免费成人在线视频| 午夜成年电影在线免费观看| 波多野结衣高清作品| 亚洲成av人片在线播放无| 国产成人精品无人区| 夜夜看夜夜爽夜夜摸| 欧美人与性动交α欧美精品济南到| 国产av在哪里看| 成人亚洲精品av一区二区| 亚洲成人国产一区在线观看| 精品久久久久久久人妻蜜臀av| 国产黄色小视频在线观看| 亚洲专区国产一区二区| 久99久视频精品免费| 成人三级黄色视频| 中文字幕最新亚洲高清| 久久久国产精品麻豆| 亚洲欧美日韩东京热| 桃红色精品国产亚洲av| 婷婷精品国产亚洲av在线| 99热这里只有精品一区 | 亚洲狠狠婷婷综合久久图片| 老司机在亚洲福利影院| 国产真实乱freesex| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清videossex| 亚洲精品美女久久av网站| 国产伦一二天堂av在线观看| 99riav亚洲国产免费| 国产爱豆传媒在线观看 | 五月玫瑰六月丁香| 久久久久免费精品人妻一区二区| 中国美女看黄片| 18禁观看日本| 美女扒开内裤让男人捅视频| 亚洲熟妇中文字幕五十中出| 老司机深夜福利视频在线观看| 丝袜人妻中文字幕| 亚洲av成人不卡在线观看播放网| 国产成人啪精品午夜网站| 99在线人妻在线中文字幕| 国产免费男女视频| 两性夫妻黄色片| 人人妻人人澡欧美一区二区| 国产人伦9x9x在线观看| 亚洲人成77777在线视频| 老熟妇乱子伦视频在线观看| 久久香蕉激情| 美女午夜性视频免费| 50天的宝宝边吃奶边哭怎么回事| 日本免费a在线| 一区二区三区激情视频| 桃色一区二区三区在线观看| 18禁国产床啪视频网站| 韩国av一区二区三区四区| 国模一区二区三区四区视频 | 国产精品自产拍在线观看55亚洲| 久久国产精品影院| 非洲黑人性xxxx精品又粗又长| 国产又色又爽无遮挡免费看| 99精品在免费线老司机午夜| 久久精品国产亚洲av高清一级| 精品国产超薄肉色丝袜足j| 亚洲免费av在线视频| 免费电影在线观看免费观看| 好看av亚洲va欧美ⅴa在| 又大又爽又粗| 国产男靠女视频免费网站| 日韩欧美国产在线观看| 国产精品av久久久久免费| 久久香蕉精品热| 桃色一区二区三区在线观看| 99热只有精品国产| 欧美绝顶高潮抽搐喷水| tocl精华| 午夜福利欧美成人| 国产激情欧美一区二区| 夜夜爽天天搞| 此物有八面人人有两片| 中文字幕最新亚洲高清| 亚洲九九香蕉| 亚洲中文字幕一区二区三区有码在线看 | 国产精品 欧美亚洲| 亚洲无线在线观看| 精品少妇一区二区三区视频日本电影| 黄色片一级片一级黄色片| 男女做爰动态图高潮gif福利片| 一区二区三区国产精品乱码| 久久精品91无色码中文字幕| 精品人妻1区二区| 最近最新中文字幕大全电影3| 国内毛片毛片毛片毛片毛片| 一区二区三区激情视频| 中出人妻视频一区二区| 黄色视频不卡| 国产精品98久久久久久宅男小说| 免费在线观看亚洲国产| 免费在线观看完整版高清| 又黄又粗又硬又大视频| 国产欧美日韩一区二区精品| 午夜a级毛片| 中文字幕高清在线视频| 日韩大码丰满熟妇| 亚洲欧美日韩无卡精品| 欧美最黄视频在线播放免费| 一个人观看的视频www高清免费观看 | 黑人巨大精品欧美一区二区mp4| 日本撒尿小便嘘嘘汇集6| 日韩欧美三级三区| 757午夜福利合集在线观看| 免费在线观看黄色视频的| 国产免费av片在线观看野外av| 亚洲国产日韩欧美精品在线观看 | 老鸭窝网址在线观看| 大型黄色视频在线免费观看| 天堂影院成人在线观看| 日本免费a在线| 99久久久亚洲精品蜜臀av| 法律面前人人平等表现在哪些方面| a级毛片a级免费在线| 一区二区三区激情视频| 色综合欧美亚洲国产小说| 男女床上黄色一级片免费看| 又黄又爽又免费观看的视频| 午夜老司机福利片| 宅男免费午夜| 不卡av一区二区三区| 在线观看免费午夜福利视频| 精品国产美女av久久久久小说| 亚洲专区字幕在线| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三| 老司机深夜福利视频在线观看| 亚洲人成电影免费在线| 香蕉av资源在线| 亚洲欧美日韩东京热| 午夜老司机福利片| 久久久久久久久久黄片| 色老头精品视频在线观看| 亚洲人成伊人成综合网2020| 99精品欧美一区二区三区四区| 日韩精品中文字幕看吧| 老司机午夜十八禁免费视频| 国产高清视频在线观看网站| 午夜激情av网站| 欧美3d第一页| 国产精品九九99| 国产精品久久久久久久电影 | 亚洲一区高清亚洲精品| 99在线人妻在线中文字幕| 麻豆成人av在线观看| 成年女人毛片免费观看观看9| 国产精品日韩av在线免费观看| 亚洲中文字幕日韩| 中出人妻视频一区二区| 午夜福利免费观看在线| 757午夜福利合集在线观看| 美女 人体艺术 gogo|