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

    紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng)

    2017-10-13 12:32:53董伊晨劉艷紅
    生態(tài)學(xué)報 2017年17期
    關(guān)鍵詞:紅松林苗齡紅松

    董伊晨, 劉艷紅

    北京林業(yè)大學(xué)森林資源與生態(tài)系統(tǒng)過程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083

    紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng)

    董伊晨, 劉艷紅*

    北京林業(yè)大學(xué)森林資源與生態(tài)系統(tǒng)過程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083

    通過分析不同苗齡紅松幼苗在控制溫度和光照條件下葉性狀的變化差異,以及各控制因素交互作用對紅松幼苗葉性狀的影響程度和紅松葉性狀的線性相關(guān)關(guān)系,探討紅松幼苗葉性狀響應(yīng)氣候環(huán)境因素變化的特征和所采取的物質(zhì)分配上的變化策略。實(shí)驗(yàn)采取人工氣候箱培養(yǎng)紅松幼苗,設(shè)置3組溫度為10、20、30℃,每組溫度下設(shè)置有效光照輻射4組:45.3、90.5、181.1、271.6(μmol m-2s-1),每組觀測苗齡1—4a。觀測性狀為:比葉重(LMA)、單位面積光合能力(Aarea)、單位質(zhì)量葉氮含量(Nmass)和單位質(zhì)量葉磷含量(Pmass)。結(jié)果顯示:(1)紅松幼苗的LMA、Aarea、Nmass、Pmass受到溫度、光照和苗齡的影響,各控制組葉性狀間差異顯著(P<0.05);(2)溫度、光照和苗齡對紅松幼苗葉性狀的交互作用效應(yīng)顯著(P<0.05),影響程度依次為:苗齡>溫度>溫度×光照×苗齡>光照>光照×苗齡>溫度×光照>溫度×苗齡;(3)LMA與溫度和苗齡呈正相關(guān)性,Aarea與苗齡呈正相關(guān)性,與溫度呈負(fù)相關(guān)性,Nmass與溫度、光照、苗齡、LMA和Aarea均呈負(fù)相關(guān)性,Pmass與溫度、Nmass呈正相關(guān)性而與苗齡和LMA呈負(fù)相關(guān)性。通過結(jié)果分析可以得出以下結(jié)論:(1)紅松幼苗葉性狀對溫度、光照和苗齡3個因素均有明顯響應(yīng)但響應(yīng)程度不一,影響性狀差異的主要因素是溫度和苗齡,光照在其溫度和苗齡二者的作用基礎(chǔ)上對性狀產(chǎn)生影響;(2)紅松幼苗葉性狀之間的線性關(guān)系實(shí)際反映的是對溫度、光照和苗齡響應(yīng)的結(jié)果,葉性狀間的密切聯(lián)系反映出葉片在響應(yīng)環(huán)境因素時所采取的分配策略的調(diào)整,其目的是維持植物體內(nèi)同化作用和異化作用的平衡,確保紅松幼苗穩(wěn)定生長。

    紅松幼苗; 控制實(shí)驗(yàn); 比葉重; 光合能力; 葉氮磷含量; 因素交互作用

    Abstract: In the present study, changes in leaf traits ofPinuskoraiensis(Korean pine) seedlings in response to seedling age, temperature, and light conditions were investigated, as well as the correlation of these factors with Korean pine seedling leaf traits. Korean pine seedlings were cultivated in artificial climate chambers with three levels of temperature settings (10, 20, and 30 ℃) and four levels of effective light radiation applied to each temperature level (45.3, 90.5, 181.1, and 271.6 μmol m-2s-1); the observed seedlings in each group were aged between 1 and 4 years. The observed leaf mass traits were per area (LMA), area-based photosynthetic capacity (Aarea), mass-based nitrogen concentration (Nmass), and mass-based phosphorus concentration (Pmass). Our results showed that (1) LMA,Aarea,Nmass, andPmassof Korean pine seedlings were affected by temperature, light, and seedling age. The differences in leaf traits among different controlled groups were statistically significant (P<0.05). (2) Temperature, light, and seedling age in combination affected Korean pine seedling leaf traits. The order of significance of effects was seedling age > temperature > temperature × light × seedling age > light > light × seedling age > temperature × light > temperature × seedling age. (3) LMA was positively correlated with temperature and seedling age;Aareawas positively correlated with seedling age but negatively correlated with temperature;Nmasswas negatively correlated with temperature, light, seedling age, LMA, andAarea; andPmasswas positively correlated with temperature andNmassbut was negatively correlated with seedling age and LMA. Our results suggest the following: (1) Korean pine seedling leaf traits responded to all three studied factors (temperature, light, and seedling age) to differing degrees. The main factors that caused differences among traits were temperature and seedling age, while the effect of light was influenced by these two factors. (2) The correlation between Korean pine seedling leaf traits reflected their response to temperature, light, and seedling age. The close relationship between leaf traits is consistent with the allocation strategies that the leaves use in response to environmental factors, which maintain the balance between assimilation and respiration to ensure stable growth in Korean pine seedlings.

    KeyWords: Korean pine seedlings; controlled experiments; leaf mass per area; photosynthetic capacity; leaf nitrogen and phosphorus concentrations; factor interaction

    紅松針闊混交林是中國東北東部山區(qū)的典型森林群落類型,紅松作為主要優(yōu)勢樹種在該地區(qū)森林演替更新、群落結(jié)構(gòu)動態(tài)和生物多樣性組成等方面起到關(guān)鍵作用,因此其生態(tài)特性一直是人們關(guān)注和研究的對象[1-4]。此外,紅松林群落及其物種多樣性對氣候動態(tài)響應(yīng)關(guān)系等方面的研究又使其與全球氣候變化、多樣性保護(hù)等熱點(diǎn)問題[5-6]緊密聯(lián)系在一起。近年來,在對前輩研究不斷深入和補(bǔ)充的基礎(chǔ)上,利用植物功能性狀(plant functional traits)和葉功能性狀(leaf functional traits)評價紅松群落生態(tài)學(xué)特征和氣候適應(yīng)性的方法,對紅松群落數(shù)量化理論和模型研究起到推動作用[7-8]。通過對紅松葉性狀的長期監(jiān)測不僅可以掌握紅松林群落動態(tài)變化,而且可以對紅松生長的環(huán)境影響程度和物質(zhì)分配策略進(jìn)行評價[9- 10],為紅松響應(yīng)環(huán)境變化方面的研究提供理論依據(jù)和方法。

    植物功能性狀和葉功能性狀研究的主要目的是闡述植物在進(jìn)化過程中所形成的環(huán)境適應(yīng)機(jī)制[11]。植物性狀在環(huán)境干擾評價、量化解釋生態(tài)系統(tǒng)功能和多樣性以及易于調(diào)查獲取等方面的優(yōu)勢廣泛應(yīng)用于生態(tài)學(xué)不同尺度的研究[9,12-17]。其中葉經(jīng)濟(jì)型譜(leaf economics spectrum)已形成較完整的理論體系[18- 19],在其基礎(chǔ)上建立的全球植物性狀數(shù)據(jù)庫[20- 21]也廣受關(guān)注。此外,植物葉性狀與植物代謝和植物相關(guān)生長等生態(tài)學(xué)理論聯(lián)系緊密,通過模型擬合的方法可以進(jìn)一步深化植物生長與氣候環(huán)境適應(yīng)關(guān)系的研究[22-28]。

    目前氣候環(huán)境因素影響植物葉性狀的研究多以全球、區(qū)域或群落尺度為主,而對特定種植物葉性狀的研究卻較鮮見[21],盡管如此,個體物種研究依然可以通過其性狀的生理特征闡述植物的環(huán)境適應(yīng)性和生態(tài)系統(tǒng)過程和功能[11,17]。因此研究紅松葉性狀無論是對紅松樹種還是紅松林群落意義都很明顯:首先通過葉性狀描述紅松對氣候環(huán)境的適應(yīng)性過程符合目前植物葉性狀的主要研究方向;其次通過葉性狀間的線性關(guān)系論述紅松內(nèi)部組織的物質(zhì)分配策略符合葉經(jīng)濟(jì)型譜的主要研究內(nèi)容[18];紅松及其群落作為東北亞地區(qū)典型森林群落類型,是世界葉經(jīng)濟(jì)型譜中重要的組成部分[20,29]。

    氣候環(huán)境因素對植物生長的影響是綜合性的[30-32],溫度和光照是其中重要的影響因素。溫度限制植物個體生長、生物產(chǎn)量和物種分布范圍[33-36],植物葉性狀對溫度變化的響應(yīng)機(jī)制反映了其生長適應(yīng)性,其形態(tài)、物候及生長型都是長期適應(yīng)進(jìn)化的結(jié)果[27,37]。光照強(qiáng)度影響植物光合作用中光響應(yīng)過程,進(jìn)而影響植物的生長和生物量積累[38-40]。溫度和光照對紅松生長分布及更新的意義明顯,紅松葉性狀對溫度和光照變化的響應(yīng)特征在紅松林群落結(jié)構(gòu)和物種多樣性研究中作用突出,尤其是光照對天然紅松林更新的影響,作為紅松林群落演替過程的重點(diǎn)問題一直備受關(guān)注[2,41-47]。

    基于上述原因,本文旨在通過分析紅松幼苗葉性狀在氣候環(huán)境因素影響下的響應(yīng)特征來判斷:(1)能否通過葉性狀的變化差異判斷紅松受到環(huán)境因素影響;(2)能否通過葉性狀的響應(yīng)差異判斷環(huán)境因素間影響程度是否相同;(3)能否通過葉性狀相關(guān)性判斷環(huán)境因素的影響是否改變紅松物質(zhì)分配的策略。此外,我國東北地區(qū)的紅松闊葉林所面臨現(xiàn)實(shí)困境依然嚴(yán)峻:如紅松林分布減少、紅松林?jǐn)?shù)量萎縮等風(fēng)險[3]。因此進(jìn)行紅松個體性狀研究除可以為紅松群落功能性狀提供數(shù)據(jù)支持,也為紅松的保護(hù)和栽培提供實(shí)驗(yàn)依據(jù)。

    1 實(shí)驗(yàn)材料與方法

    1.1 實(shí)驗(yàn)材料

    紅松幼苗產(chǎn)自黑龍江小興安嶺南坡。先將紅松幼苗在北京林業(yè)大學(xué)氣象站實(shí)驗(yàn)地盆栽培養(yǎng),其栽培用土為產(chǎn)苗當(dāng)?shù)赝寥馈_x擇苗齡在1—4a的健壯紅松幼苗作為實(shí)驗(yàn)樣本。

    1.2 實(shí)驗(yàn)控制條件

    實(shí)驗(yàn)時間為2014年3—9月間,將紅松幼苗放置在人工氣候箱(RXZ- 280C型)中進(jìn)行溫度和光照處理。溫度、光照、樹齡三者對紅松幼苗是否存在交互作用是本實(shí)驗(yàn)關(guān)注的問題之一,為有效判斷這些因素的功效作用,實(shí)驗(yàn)采取的控制條件設(shè)置為:(1)根據(jù)東北地區(qū)紅松針闊混交林在生長季時的溫度范圍[3,33]設(shè)置10、20、30℃,3個溫度處理組,視10—20℃為紅松生長季溫度,將30℃視為極端溫度的臨界值,低于10℃可能因生理低溫的影響不作為考慮溫度[48-49]。(2)根據(jù)天然紅松林下更新的遮蔭影響[45- 47],將光照范圍設(shè)定在實(shí)際光有效輻射40—300μmol m-2s-1之間,以此測定紅松幼苗在微光環(huán)境下的葉性狀變化。根據(jù)每個溫度組分別設(shè)置4個光照控制組:2200、4400、8800、13200lx(氣候箱實(shí)際光照度最大值22000lx),文中將光強(qiáng)換算為光有效輻射(換算公式:1μE-1·s-1=48.6lx)[50],后文均以光輻射值45.3、90.5、181.1、271.6 μmol m-2s-1表示光照條件。(3)按照實(shí)驗(yàn)條件的設(shè)置,每個溫度下均有4個不同光照的氣候箱,每個氣候箱中放置1—4a苗齡的紅松幼苗,每苗齡5株重復(fù),全部實(shí)驗(yàn)組共計(jì)48組(即48水平處理組)。每日光照時長為14h,箱內(nèi)空氣濕度(RH%)保持在(60±5)%。

    1.3 葉性狀測量方法

    在紅松幼苗生長旺盛的8—9月進(jìn)行性狀測定??紤]到針葉壽命長和往年生葉片的差異,故選擇生長情況較一致的當(dāng)年生葉片進(jìn)行性狀測量,以排除選擇干擾。葉性狀測定方法:(1)比葉重LMA(g/m2),取針葉3簇,視紅松葉為5等分圓柱體,視等分體外側(cè)弧狀面為光輻射面,葉面積即為光輻射面積總和。利用尺量葉長,掃描葉橫截面(Image-Pro Plus 6.0計(jì)算弧長),葉面積=弧長×葉長×5(未計(jì)楔狀面);烘干葉片稱取葉干重, LMA=葉干重/葉面積,LMA計(jì)算曾通過經(jīng)驗(yàn)?zāi)P陀?jì)算法比對[16,51]。(2)單位面積葉光合能力Aarea(μmolC m-2s-1),利用LI6400便攜式光合測定儀(LiCor, Lincoln, NE,USA)和針狀葉室(needle chamber, IRGA)測定紅松幼苗凈光合速率,每株取葉3簇。測定環(huán)境與處理環(huán)境相同,測定前記錄測定葉寬,測定時手動輸入:測量面積=葉寬×葉室長(若葉不及葉室長則輸入葉面積),選擇最大瞬時光合速率(10:00—11:00)作為單位面積光合能力Aarea[16,52]。(3)單位質(zhì)量葉氮含量Nmass(%)利用凱式定氮法測定。(4)單位質(zhì)量葉磷含量Pmass(%)利用鉬銻抗火焰光度計(jì)法測定。

    1.4 數(shù)據(jù)分析方法

    紅松幼苗葉性狀數(shù)據(jù)利用SPSS 18.0進(jìn)行多因素方差分析(GLM Multivariate),得到各組均值(mean±SE)、Duncan一致性子集檢驗(yàn)各組差異顯著值以及多因素功效檢驗(yàn)解釋值。利用相關(guān)性分析(Bivariate, pearson correlation)、一元線性回歸(Linear Regression, least-squares)分析性狀與因素關(guān)系、性狀間關(guān)系等,其中葉性狀指標(biāo)均進(jìn)行對數(shù)轉(zhuǎn)化(log10-transformed)[18-19],以降低方差非齊性可能造成的數(shù)據(jù)離散和判定概率系數(shù)P值偏小的干擾。利用Origin 8.0和Excel 2007制作圖表。

    2 結(jié)果分析

    2.1 不同苗齡紅松幼苗葉性狀對溫度、光照條件變化的響應(yīng)

    紅松幼苗葉性狀指標(biāo):LMA、Aarea、Nmass和Pmass對溫度、光照、苗齡3因素48水平(溫度3水平、光照4水平、苗齡4水平)的響應(yīng)呈現(xiàn)不同的變化趨勢(表1—表4)。實(shí)驗(yàn)結(jié)果表明溫度、光照、苗齡的變化對紅松幼苗葉性狀均產(chǎn)生明顯影響。

    表1 紅松幼苗比葉重(LMA)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    表2 紅松幼苗單位面積光合能力(Aarea)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    (1)紅松幼苗LMA的變化(表1):在不同溫度或光照組條件下,紅松幼苗LMA隨著苗齡的增長而上升,各苗齡組間差異顯著(P<0.05)。在固定苗齡組比較不同光照條件中LMA均值隨溫度的升高而呈現(xiàn)上升的趨勢,各溫度組間差異顯著(P<0.05)。相同溫度和苗齡組,紅松幼苗LMA受光輻射強(qiáng)度影響的組間差異顯著(P<0.05)。

    (2)紅松幼苗Aarea變化趨勢(表2):在不同的光照和苗齡組,溫度由10℃上升到20℃時Aarea呈現(xiàn)隨之升高的變化,而溫度持續(xù)升高到30℃時Aarea出現(xiàn)明顯下降,溫度組間差異顯著(P<0.05)。在相同的溫度和光照組中,隨苗齡的增長Aarea呈現(xiàn)上升的變化趨勢,各苗齡組差異顯著(P<0.05)。而Aarea隨光照有效輻射的升高,同一苗齡組有上升的趨勢,各光照組間差異顯著(P<0.05)。

    (3)紅松幼苗Nmass的變化趨勢(表3):在不同的溫度和光照組,Nmass隨苗齡的增長而下降,各苗齡間差異顯著(P<0.05)。在相同的溫度和苗齡組Nmass隨光有效輻射的增高而下降,各光照組間差異顯著(P<0.05)。Nmass在同一苗齡組隨溫度上升而減少,各溫度組間差異顯著(P<0.05)。

    (4)紅松幼苗Pmass的變化趨勢(表4):在不同的溫度和光照組中,Pmass隨苗齡的增長呈現(xiàn)小幅降低,各苗齡組間差異顯著(P<0.05)。在相同的光照和苗齡組中,Pmass隨溫度的升高而出現(xiàn)波動變化,但仍呈現(xiàn)下降的趨勢,各溫度組間差異顯著(P<0.05)。在同一的溫度組中Pmass隨光有效輻射的增強(qiáng)出現(xiàn)不同的變化,但在同一的苗齡組中隨光照變化的趨勢不明顯,各光照組差異顯著(P<0.05)。

    表3 紅松幼苗單位質(zhì)量葉氮含量(Nmass)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    表4 紅松幼苗單位質(zhì)量葉磷含量(Pmass)在不同溫度和光照條件下的變化

    表中數(shù)值為 GLM Multivariate方差檢驗(yàn)均值和標(biāo)準(zhǔn)誤差,n=5,F檢驗(yàn)概率P<0.05,字母為Duncan一致性子集檢驗(yàn),ABCD為苗齡組,abc為光照組,xyz為溫度組

    2.2 溫度、光照和苗齡對紅松幼苗葉性狀影響的交互作用功效分析

    紅松幼苗葉性狀對溫度、光照和苗齡3個因素的響應(yīng)變化趨勢不同,通過交互作用分析和功效檢驗(yàn)對三者的影響程度做進(jìn)一步分析。表5給出的是溫度、光照和苗齡3個因素對紅松幼苗葉性狀的交互作用和功效分析,3個單項(xiàng)(溫度、光照、苗齡)表示因素交互作用的主效應(yīng),3個二維項(xiàng)(溫度×光照、溫度×苗齡、光照×苗齡)表示的是兩兩因素的一階交互效應(yīng),1個三維項(xiàng)(溫度×光照×苗齡)則表示的是3因素的二階交互效應(yīng)。Wilks′ Lambda值為基于組間無差別的零級假設(shè)判別系數(shù),系數(shù)越接近于0則說明因素功效越有差異,檢驗(yàn)概率P<0.05表明各因素的影響差異顯著。

    主效應(yīng)判別中溫度和苗齡的系數(shù)最小(0.003、0.035),低于光照的系數(shù)(0.163),說明溫度和苗齡的獨(dú)立作用效率高于光照。二維交互項(xiàng)判別中溫度×光照系數(shù)最小(0.115),其次是光照×苗齡(0.136),溫度×苗齡系數(shù)最大(0.223),說明溫度×光照和光照×苗齡的交互作用效率高于溫度×苗齡,也即是說溫度和苗齡二者獨(dú)立作用效率高于二者的交互作用,光照×溫度、光照×苗齡的交互作用效率高于光照的獨(dú)立作用。三維交互項(xiàng)溫度×光照×苗齡系數(shù)(0.087)高于其他所有二維項(xiàng),說明溫度×光照×苗齡三者交互作用效率高于任意兩兩交互作用。通過交互作用判別可以看出光照的獨(dú)立作用程度受到溫度和苗齡的制約,同時光照也使溫度和苗齡的交互作用效率有所降低。

    效應(yīng)歸因解釋(Partial Eta-Squared)給出的因素解釋概率η2值指出溫度、光照和苗齡3因素對紅松幼苗葉性狀影響程度的強(qiáng)弱關(guān)系。各因素交互作用對葉性狀的解釋率分別為:苗齡對LMA的主效應(yīng)解釋概率最高(0.861)。溫度對Aarea的主效應(yīng)解釋率最高(0.990),其次為苗齡(0.885)、溫度×光照(0.703)、光照(0.629),說明紅松幼苗光合能力對環(huán)境因素的響應(yīng)過程比較復(fù)雜。光照和苗齡對Nmass的主效應(yīng)解釋率較高(0.630、0.794),其次是光照×苗齡(0.604),說明其二者的交互效應(yīng)對紅松幼苗葉性狀的影響既有獨(dú)立性又有相互促進(jìn)性。苗齡對Pmass的主效應(yīng)解釋率最高(0.792),其次為溫度(0.645)和溫度×光照×苗齡(0.613)??偨Y(jié)表5功效解釋率依次為:苗齡>溫度>溫度×光照×苗齡>光照>光照×苗齡>溫度×光照>溫度×苗齡。

    表5 各影響因素功效檢驗(yàn)和對葉性狀的交互作用解釋

    Wilks′ Lambda值為各因素差異性檢驗(yàn),F檢驗(yàn)概率P<0.05說明差異顯著,效應(yīng)估計(jì)η2值表示各因素功效解釋度;*因素為溫度,光照,苗齡,溫度×光照,溫度×苗齡,光照×苗齡,溫度×光照×苗齡

    通過分析比較紅松幼苗各葉性狀與溫度、光照和苗齡的交互作用和功效解釋率,發(fā)現(xiàn)各因素間的主效應(yīng)、二維交互效應(yīng)和三維交互效應(yīng)都較顯著(P<0.05),且對葉性狀的影響過程較為復(fù)雜。說明溫度、光照和苗齡3因素對紅松幼苗葉性狀的影響程度不同,因素間的獨(dú)立性和依賴性在交互過程中由于分組水平的變化而產(chǎn)生差異,因此在描述葉性狀的差異變化上顯得龐雜,然而這也是客觀現(xiàn)實(shí)的體現(xiàn)。利用因素的解釋度來分析溫度、光照和苗齡對紅松幼苗葉性狀影響程度,可以使結(jié)果更加清晰。

    2.3 紅松幼苗葉性狀間及與溫度、光照及苗齡等因素的相關(guān)關(guān)系

    在肯定溫度、光照和苗齡對紅松幼苗葉性狀影響程度存在差別的基礎(chǔ)上,通過相關(guān)性分析對各葉性狀與因素之間的密切程度,以及葉性狀間的相互關(guān)聯(lián)進(jìn)行分析。表6為紅松幼苗葉性狀與溫度、光照、苗齡的相關(guān)關(guān)系,以及葉性狀之間的相關(guān)關(guān)系。采用Pearson相關(guān)性檢驗(yàn)和雙尾顯著性t-檢驗(yàn),概率分為P<0.01、P<0.05表示顯著性。其中LMA與溫度和苗齡存在顯著正相關(guān)關(guān)系;Aarea與苗齡存在顯著正相關(guān)關(guān)系,與溫度成顯著負(fù)相關(guān)關(guān)系;Nmass與溫度、光照、苗齡、LMA和Aarea均成顯著負(fù)相關(guān)關(guān)系;Pmass與溫度、Nmass成顯著正相關(guān)關(guān)系而與苗齡和LMA則呈顯著負(fù)相關(guān)關(guān)系。

    表6 紅松幼苗葉性狀與溫度、光照間的相關(guān)性

    各葉性狀對數(shù)值采用Pearson Correlation 相關(guān)性檢驗(yàn)(2- tailed),**表示P<0.01顯著,*表示P<0.05顯著

    3 討論

    3.1 紅松幼苗葉性狀與溫度、光照和苗齡的響應(yīng)關(guān)系

    通過實(shí)驗(yàn)結(jié)果分析可知,紅松幼苗葉性狀對溫度、光照和苗齡響應(yīng)差異顯著,溫度、光照和苗齡既有獨(dú)立作用也有交互作用,且影響過程復(fù)雜,光照的作用是在溫度和苗齡的作用基礎(chǔ)上體現(xiàn)的。紅松幼苗葉性狀與溫度、光照和苗齡相關(guān)性表現(xiàn)為:溫度和苗齡對4種葉性狀都具有顯著相關(guān)性,而光照除與Nmass相關(guān)外與其他性狀沒有顯著相關(guān)性,說明紅松葉性狀更多表現(xiàn)為響應(yīng)溫度和苗齡的變化,響應(yīng)關(guān)系明顯。這些關(guān)系反映出:一方面苗齡對Aarea變化功效顯著,兩者顯著正相關(guān),即苗齡越高紅松的生物積累越多,另一方面溫度與LMA、Pmass呈正顯著相關(guān)性,說明紅松幼苗因?yàn)闇囟鹊纳仙訌?qiáng)了呼吸和代謝的作用,LMA上升是葉片呼吸造成水分散失加劇,導(dǎo)致葉面積收縮相對減少的結(jié)果。而LMA和苗齡的顯著正相關(guān)性又說明紅松幼苗生長和抵御環(huán)境干擾的能力也隨苗齡增強(qiáng)。磷素在植物細(xì)胞核質(zhì)和線粒體為補(bǔ)充轉(zhuǎn)錄RNA提供構(gòu)件[53-54],在群落研究中屬于生境限制因素[19,29],部分研究也利用磷素在膜質(zhì)結(jié)構(gòu)的多少來判斷植物速生與否[21,26,55]。對于紅松樹種而言,其本身非速生樹種,Pmass更多反映的是線粒體的代謝水平,它與LMA都反映紅松幼苗對溫度升高的響應(yīng)措施:隨著溫度逐步升高代謝增強(qiáng)導(dǎo)致Pmass下降,水分物質(zhì)流失加劇,對紅松生長帶來不利影響。而Pmass與苗齡呈負(fù)顯著相關(guān)性,即在溫度持續(xù)升高條件下紅松幼苗降低磷素代謝的水平隨年齡遞增而降低,這樣的措施抑制了部分水分物質(zhì)散失量,光合作用所積累的物質(zhì)消耗相對減少,可以避免紅松生長時因環(huán)境干擾所帶來不利影響。

    紅松幼苗的生長需光量的研究國內(nèi)已有不少報道[2,45- 47,56],涉及到紅松林下更新和幼苗生長適應(yīng)性等問題,對紅松幼苗需光或需庇蔭的討論較多[46-47,56]。從本實(shí)驗(yàn)的數(shù)據(jù)結(jié)果上看:首先功效檢驗(yàn)溫度×光照、光照×苗齡、溫度×光照×苗齡交互作用效率要比光照單獨(dú)作用顯著,其次紅松幼苗葉性狀中的除Nmass與光照顯著負(fù)相關(guān)外,LMA、Aarea、Pmass與光照均無顯著相關(guān)性,說明光照對紅松幼苗的生長獨(dú)立影響不明顯,而是與溫度和苗齡共同影響。此外在本實(shí)驗(yàn)中對紅松幼苗Aarea的功效解釋較高的因素有溫度、光照、苗齡、溫度×光照,也可以看出紅松光合作用受到多種因素限制,在已知的研究中也體現(xiàn)出這一特點(diǎn)[57-58]。不僅如此,本實(shí)驗(yàn)中Aarea在有效光輻射45.3,90.5,181.1,271.6(μmol m-2s-1)4種條件下數(shù)值差異顯著且沒有低于補(bǔ)償值范圍,同化依然高于呼吸,說明即使在這種光照偏低的條件下紅松幼苗依然可以進(jìn)行光合作用。在前人的研究中通過紅松幼苗的光響應(yīng)曲線判斷其在林下光適應(yīng)的大致范圍從20到900不等,全光可接近2000 μmol m-2s-1,說明紅松幼苗的光適應(yīng)范圍較寬[56,58-59]。

    綜上所述,溫度和光照對不同苗齡的紅松幼苗葉性狀影響差異顯著,葉性狀表現(xiàn)是溫度、光照和苗齡的獨(dú)立作用和交互作用共同影響的結(jié)果。溫度和苗齡影響功效更顯著,光照在其溫度和苗齡二者的影響基礎(chǔ)上產(chǎn)生作用。

    3.2 紅松幼苗葉性狀間線性相關(guān)關(guān)系

    通過相關(guān)性分析得知紅松幼苗葉性狀間存在顯著相關(guān)性,說明在響應(yīng)環(huán)境因素影響時葉性狀之間存在密切聯(lián)系。表7為紅松幼苗葉性狀logLMA、logAarea、logNmass和logPmass的一元線性擬合結(jié)果。其中l(wèi)ogLMA與logNmass、logPmass為負(fù)線型相關(guān),解釋概率分別為35%和30%,即紅松幼苗LMA與Nmass和Pmass變化趨勢相反。紅松幼苗比葉重在溫度、光照影響下隨苗齡上升而上升,而葉片中氮磷元素含量卻隨之減少,這說明紅松幼苗在維持光合呼吸過程中消耗的元素增多,而對氮磷元素吸收的能力因?yàn)榄h(huán)境因素的影響而降低[45]。logNmass和logPmass的正線性相關(guān)說明二者的變化趨勢是類似的(解釋概率12%),原因可能是在響應(yīng)外界干擾時,紅松幼苗需要同時調(diào)節(jié)光合和呼吸水平[58- 59]。除此以外,logNmass與logAarea之間呈顯著負(fù)線型關(guān)系,氮素含量的下降并沒有使紅松幼苗光合能力降低,說明說明光合作用在氮含量降低的情況下仍能維持。同時logNmass與logAarea之間的關(guān)系也并非簡單的數(shù)量上相關(guān)(解釋概率2%),還與氮參與光合作用的形式、機(jī)制以及利用效率有關(guān)[60],氮磷元素含量在光合作用和葉形態(tài)上相互關(guān)聯(lián)[61- 62],這其中氮素提供光合作用所需原件,磷素分別為膜系統(tǒng)和分解代謝提供所需原件[63-64],它們含量的降低主要是溫度上升導(dǎo)致LMA上升,水分和元素散失加劇的影響。紅松幼苗比葉重和光合能力在溫度和光照條件影響下對氮磷含量降低的表現(xiàn)反映的是其維持光合、呼吸之間的平衡的結(jié)果。

    表7 紅松幼苗葉性狀線性回歸分析

    根據(jù)性狀相關(guān)關(guān)系的顯著性進(jìn)行線性回歸分析,表中各項(xiàng)數(shù)值基于最小二乘法(least-squares regression)的結(jié)果

    紅松維持光合和呼吸平衡的過程中LMA和Aarea間雖無顯著線性相關(guān)性,但二者均與溫度、苗齡顯著相關(guān),存在間接關(guān)聯(lián)。LMA和Aarea都反映植物干物質(zhì)積累與葉生長關(guān)系[18-19,55],光合作用又受到胞間CO2濃度、水分利用能力和氣孔等限制[60,65]。紅松幼苗LMA的變化是由溫度的升高引起的,葉片水分在高溫下散失加劇,導(dǎo)致葉表面積收縮干物質(zhì)比重上升,而葉片通過降低氣體流量和水分散失來應(yīng)對這種干擾,導(dǎo)致光合能力下降,當(dāng)生長年限增加時紅松對這種干擾的抵抗能力都有所加強(qiáng)。所以盡管LMA與Aarea之間不存在線性相關(guān)性,但對相同影響因素所做出的響應(yīng)都是降低物質(zhì)水分消耗,維持紅松幼苗正常生長的過程,具有一定的相互影響。

    綜上所述,紅松幼苗葉性狀間的線性關(guān)系是對溫度、光照和苗齡響應(yīng)的結(jié)果。葉性狀間緊密聯(lián)系反應(yīng)的是紅松為響應(yīng)環(huán)境因素影響所采取的物質(zhì)分配策略的變化。在這個過程中同一因素對紅松葉性狀的影響方式和程度都有差別,同一性狀所響應(yīng)的因素也不相同,其目的都是維持植物體內(nèi)同化和異化作用的平衡,即保障紅松幼苗穩(wěn)定生長。

    4 結(jié)論

    通過上述討論可以明確紅松幼苗葉性狀對溫度、光照和苗齡的響應(yīng)本質(zhì)是其在不同條件下維持自身生長所采取的措施。葉性狀與3種因素間的關(guān)系,以及葉性狀之間的變化關(guān)系是其在調(diào)節(jié)同化和異化作用平衡所采取的相應(yīng)策略。因此總結(jié)數(shù)據(jù)分析和討論的結(jié)果可以得出的結(jié)論為:

    (1)紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng)差異顯著,葉性狀差異是溫度、光照和苗齡的獨(dú)立作用和交互作用共同影響的結(jié)果。溫度和苗齡影響功效更顯著,光照在其溫度和苗齡二者的影響基礎(chǔ)上產(chǎn)生作用。

    (2)紅松幼苗葉性狀間的線性關(guān)系是對溫度、光照和苗齡響應(yīng)的結(jié)果。葉性狀間緊密聯(lián)系反應(yīng)的是紅松為響應(yīng)環(huán)境因素影響所采取的物質(zhì)分配策略的變化,其目的是維持植物體內(nèi)同化和異化作用的平衡保障紅松幼苗穩(wěn)定生長。

    除此以外,從實(shí)驗(yàn)中不同苗齡紅松幼苗葉性狀做出的響應(yīng)過程看氣候環(huán)境因素的作用是綜合的。如何對氣候環(huán)境因素分別進(jìn)行定量劃分范圍和歸因解釋是一個難點(diǎn),而進(jìn)行人為控制實(shí)驗(yàn)正是要將不同影響因素進(jìn)行逐一驗(yàn)證,尋找準(zhǔn)確的解釋方法。本實(shí)驗(yàn)選取的溫度、光照是重要的氣候因素,對紅松生長發(fā)育意義明顯,而苗齡(或樹齡)既是性狀的限定因素,又可以指示紅松林群落生物積累,是聯(lián)系個體紅松生長和種群分布、更新的橋梁。對幾種因素間的獨(dú)立、交互作用效果分別進(jìn)比較,可以較為準(zhǔn)確地判斷引起性狀變化的原因,為今后更有針對性地開展紅松的性狀研究提供了依據(jù)。

    [1] 朱濟(jì)凡, 劉慎諤, 王戰(zhàn), 馮宗煒, 劉同生. 小興安嶺紅松針闊葉混交林. 林業(yè)科學(xué), 1958, 4(4): 355- 369.

    [2] 李俊清. 闊葉紅松林中紅松的分布格局及其動態(tài). 東北林業(yè)大學(xué)學(xué)報, 1986, 14(1): 33- 38.

    [3] 王業(yè)蘧. 闊葉紅松林. 哈爾濱: 東北林業(yè)大學(xué)出版社, 1995.

    [4] 于振良, 趙士洞, 董立榮, 張?jiān)? 闊葉紅松林的研究綜述. 吉林林學(xué)院學(xué)報, 1996, 12(4): 235- 238.

    [5] 鄧慧平, 吳正方, 周道瑋. 全球氣候變化對小興安嶺闊葉紅松林影響的動態(tài)模擬研究. 應(yīng)用生態(tài)學(xué)報, 2000, 11(1): 43- 46.

    [6] 吳曉莆, 朱彪, 趙淑清, 樸世龍, 方精云. 東北地區(qū)闊葉紅松林的群落結(jié)構(gòu)及其物種多樣性比較. 生物多樣性, 2004, 12(1): 174- 181.

    [7] 胡耀升, 么旭陽, 劉艷紅. 長白山不同演替階段森林植物功能性狀及其與地形因子間的關(guān)系. 生態(tài)學(xué)報, 2014, 34(20): 5915- 5924.

    [8] 么旭陽, 胡耀升, 劉艷紅. 長白山闊葉紅松林不同群落類型的植物功能性狀與功能多樣性. 西北農(nóng)林科技大學(xué)學(xué)報: 自然科學(xué)版, 2014, 42(3): 77- 84.

    [9] McIntyre S, Lavorel S, Landsberg J, Forbes T D A. Disturbance response in vegetation-towards a global perspective on functional traits. Journal of Vegetation Science, 1999, 10(5): 621- 630.

    [10] Wright I J, Reich P B, Cornelissen J H C, Falster D S, Groom P K, Hikosaka K, Lee W, Lusk C H, Niinemets ü, Oleksyn J, Osada N, Poorter H, Warton D I, Westoby M. Modulation of leaf economic traits and trait relationships by climate. Global Ecology and Biogeography, 2005, 14(5): 411- 421.

    [11] 孟婷婷, 倪健, 王國宏. 植物功能性狀與環(huán)境和生態(tài)系統(tǒng)功能. 植物生態(tài)學(xué)報, 2007, 31(1): 150- 165.

    [12] 毛偉, 李玉霖, 張銅會, 趙學(xué)勇, 黃迎新, 宋琳琳. 不同尺度生態(tài)學(xué)中植物葉性狀研究概述. 中國沙漠, 2012, 32(1): 33- 41.

    [13] Díaz S, Cabido M C, Casanoves F. Plant functional traits and environmental filters at a regional scale. Journal of Vegetation Science, 1998, 9(1): 113- 122.

    [14] Díaz S, Cabido Marcelo C, Zak M, Carretero E M, Araníbar J. Plant functional traits, ecosystem structure and land-use history along a climatic gradient in central-western Argentina. Journal of Vegetation Science, 1999, 10(5): 651- 660.

    [15] Díaz S, Lavorel S, Mcintyre S, Falczuk V, Casanoves F, Milchunas D G., Skarpe C, Rusch G, Sternberg M, Noy-Meir I, Landsberg J, Zhang W, Clark H, Campbell B. Plant trait responses to grazing-a global synthesis. Global Change Biology, 2007, 13(2): 313- 341.

    [16] Cornelissen J H C, Lavorel S, Garnier E, Díaz S, Buchmann N, Gurvich D E, Reich P B, Steege H T, Morgan H D, van der Heijden M G A, Pausas J G, Poorter H. A Handbook of Protocols for Standardised and Easy Measurement of Plant Functional Traits Worldwide. Australian Journal of Botany, 2003, 51(4): 335- 380.

    [17] 肖玉, 謝高地, 安凱, 魯春霞. 基于功能性狀的生態(tài)系統(tǒng)服務(wù)研究框架. 植物生態(tài)學(xué)報, 2012, 36(4): 353- 362.

    [18] Wright I J, Reich P B, Westoby M, Ackerly D D, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen J H C, Diemer M, Flexas J, Garnier E, Groom P K, Gulias J, Hikosaka K, Lamont B B, Lee T, Lee W, Lusk C, Midgley J J, Marie-Laure N, Niinemets ü, Oleksyn J, Osada N, Poorter H, Poot P, Prior L, Pyankov V I, Roumet C, Thomas S C, Tjoelker M G, Veneklaas E J, Villar R. The worldwide leaf economics spectrum. Nature, 2004, 428(6985): 821- 827.

    [19] Wright I J, Reich P B, Cornelissen J H C, Falster D S, Garnier E, Hikosaka K, Lamont B B, Lee W, Oleksyn J, Osada N, Poorter H, Villar R, Warton D I, Westoby M. Assessing the generality of global leaf trait relationships. New Phytologist, 2005, 166(2): 485- 496.

    [20] Kattge J, Díaz S, Lavorel S, Prentice I C, Leadley P, B?nisch G, Garnier E, Westoby M, Reich P B, Wright I J, Cornelissen J H C, Violle C, Harrison S P, van Bodegom P M, Reichstein M, Enquist B J, Soudzilovskaia N A, Ackerly D D, Anand M, Atkin O, Bahn M, Baker T R, Baldocchi D, Bekker R, Blanco C C, Blonder B, Bond W J, Bradstock R, Bunker D E, Casanoves F, Cavender-Bares J, Chambers J Q, ChapinF S, Chave J, Coomes D, Cornwell W K, Craine J M, Dobrin B H, Duarte L, Durka W, Elser J, Esser G, Estiarte M, Fagan W F, Fang J, Fernández-Méndez F, Fidelis A, Finegan B, Flores O, Ford H, Frank D, Freschet G T, Fyllas N M, Gallagher R V, Green W A, Gutierrez A G, Hickler T, Higgins S I, Hodgson J G, Jalili A, Jansen S, Joly C A, Kerkhoff A J, Kirkup D, Kitajima K, Kleyer M, Klotz S, Knops J M H, Kramer K, Kühn I, Kurokawa H, Laughlin D, Lee T D, Leishman M, Lens F, Lenz T, Lewis S L, Lloyd J, Llusià J, Louault F, Ma S, Mahecha M D, Manning P, Massad T, Medlyn B E, Messier J, Moles A T, Müller S C, Nadrowski K, Naeem S, Niinemets ü, N?llert S, Nüske A, Ogaya R, Oleksyn J, Onipchenko V G, Onoda Y, Ordoez J, Overbeck G, Ozinga W A, Patio S, Paula S, Pausas J G, Peuelas J, Phillips O L, Pillar V, Poorter H, Poorter L, Poschlod P, Prinzing A, Proulx R, Rammig A, Reinsch S, Reu B, Sack L, Salgado-Negret B, Sardans J, Shiodera S, Shipley B, Siefert A, Sosinski E, Soussana F J, Swaine E, Swenson N, Thompson K, Thornton P, Waldram M, Weiher E, White M, White S, Wright S J, Yguel B, Zaehle S, Zanne A E, Wirth C. TRY-a global database of plant traits. Global Change Biology, 2011, 17(9): 2905- 2935.

    [21] 陳瑩婷, 許振柱. 植物葉經(jīng)濟(jì)譜的研究進(jìn)展. 植物生態(tài)學(xué)報, 2014, 38 (10): 1135- 1153.

    [22] Swenson N G, Enquist B J. Ecological and evolutionary determinants of a key plant functional trait: wood density and its community-wide variation across latitude and elevation. American Journal of Botany, 2007, 94(3): 451- 459.

    [23] van Kleunen M, Weber E, Fischer M. A meta-analysis of trait differences between invasive and non-invasive plant species. Ecology Letters, 2010, 13(2): 235- 245.

    [24] Sack L, Scoffoni C, Grace P J, Poorter H, Mason C M, Mendez-Alonzo R, Donovan L A. How do leaf veins influence the worldwide leaf economic spectrum? Review and synthesis. Journal of Experimental Botany, 2013, 64(13): 4053- 4080.

    [25] Mason C M, McGaughey S E., Donovan L A. Ontogeny strongly and differentially alters leaf economic and other key traits in three diverse Helianthus species. Journal of Experimental Botany, 2013, 64(13): 4089- 4099.

    [26] Price C A, Wright I J, Ackerly D D, Niinemets ü, Reich P B, Veneklaas E J. Are leaf functional traits ‘invariant’ with plant size and what is ‘invariance’ anyway? Functional Ecology, 2014, 28(6): 1330- 1343.

    [27] Atkin O K, Bloomfield K J, Reich P B, Tjoelker M G, Asner G P, Bonal D, B?nisch G, Bradford M G, Cernusak L A, Cosio E G, Creek D, Crous K Y, Domingues T F, Dukes J S, Egerton J J G, Evans J R, Farquhar G D, Fyllas N M, Gauthier P P G, Gloor E, Gimeno T E, Griffin K L, Guerrieri R, Heskel M A, Huntingford C, Ishida F Y, Kattge J, Lambers H, Liddell M J, Lloyd J, Lusk C H, Martin R E, Maksimov A P, Maximov T C, Malhi Y, Medlyn B E, Meir P, Mercado L M, Mirotchnick N, Ng D, Niinemets ü, O′Sullivan O S, Phillips O L, Poorter L, Poot P, Prentice I C, Salinas N, Rowland L M, Ryan MG, Sitch S, Slot M, Smith N G, Turnbull M H, VanderWel M C, Valladares F, Veneklaas EJ, Weerasinghe L K, Wirth C, Wright I J, Wythers K R, Xiang J, Xiang S, Zaragoza-Castells J. Global variability in leaf respiration in relation to climate, plant functional types and leaf traits. New Phytologist, 2015, 206(2): 614- 636.

    [28] 祝介東, 孟婷婷, 倪健, 蘇宏新, 謝宗強(qiáng), 張守仁, 鄭元潤, 肖春旺. 不同氣候帶間成熟林植物葉性狀間異速生長關(guān)系隨功能型的變異. 植物生態(tài)學(xué)報, 2011, 35(7): 687- 698.

    [29] Chen Y H, Han W X, Tang L Y, Tang Z Y, Fang J Y. Leaf nitrogen and phosphorus concentrations of woody plants differ in responses to climate, soil and plant growth form. Ecography, 2013, 36(2): 178- 184.

    [30] 方精云. 我國森林植被帶的生態(tài)氣候?qū)W分析. 生態(tài)學(xué)報, 1991, 11(4): 377- 387.

    [31] 周廣勝, 張新時, 高素華, 白克智, 延曉東, 鄭元潤. 中國植被對全球變化反應(yīng)的研究. 植物學(xué)報, 1997, 39(9): 879- 888.

    [32] 黎磊, 陳家寬. 氣候變化對野生植物的影響及保護(hù)對策. 生物多樣性, 2014, 22(5): 549- 563.

    [33] 衛(wèi)民, 王輝民, 王其冬, 劉允芬, 賀慶棠, 袁嘉祖, 邵海榮, 宋從和. 氣候變化對我國紅松林的影響. 地理研究, 1995, 14(1): 17- 26.

    [34] 郭泉水, 閻洪, 徐德應(yīng), 王兵. 氣候變化對我國紅松林地理分布影響的研究. 生態(tài)學(xué)報, 1998, 18(5): 484- 488.

    [35] 郝占慶, 代力民, 賀紅士, David J.Malandnoff, 邵國凡. 氣候變暖對長白山主要樹種的潛在影響. 應(yīng)用生態(tài)學(xué)報.2001, 12(5):653- 658.

    [36] 冷文芳, 賀紅士, 布仁倉, 胡遠(yuǎn)滿. 氣候變化條件下東北森林主要建群種的空間分布. 生態(tài)學(xué)報, 2006, 26(12): 4257- 4266.

    [37] 馮秋紅, 史作民, 董莉莉. 植物功能性狀對環(huán)境的響應(yīng)及其應(yīng)用. 林業(yè)科學(xué), 2008, 44(4): 125- 131.

    [38] 陳悅, 許大全. 植物光合速率對光強(qiáng)從飽和到有限轉(zhuǎn)變響應(yīng)方式的物種依賴性. 植物生理與分子生物學(xué)學(xué)報, 2007, 33(6): 538- 546.

    [39] 王凱, 朱教君, 于立忠, 孫一榮, 張金鑫. 光環(huán)境對胡桃楸幼苗生長與光合作用的影響. 應(yīng)用生態(tài)學(xué)報, 2010, 21(4): 821- 826.

    [40] 馮玉龍, 曹坤芳, 馮志立, 馬玲. 四種熱帶雨林樹種幼苗比葉重, 光合特性和暗呼吸對生長光環(huán)境的適應(yīng). 生態(tài)學(xué)報, 2002, 22(6): 901- 910.

    [41] 劉慎諤. 關(guān)于大小興安嶺的森林更新問題. 林業(yè)科學(xué), 1957, 3(3): 263- 280.

    [42] 王戰(zhàn). 對于小興安嶺紅松林更新和主伐方式的意見. 林業(yè)科學(xué), 1957, 3(3): 281- 288.

    [43] 李俊清, 朱春全, 柴一新, 張東力, 劉傳照, 姚成濱. 闊葉紅松林的營養(yǎng)結(jié)構(gòu)與動態(tài)特性. 吉林林學(xué)院學(xué)報, 1989, 5(2): 1- 16.

    [44] 李俊清, 柴一新, 張東力. 人工闊葉紅松林的結(jié)構(gòu)與生產(chǎn)力. 林業(yè)科學(xué), 1990, 26(1): 1- 8.

    [45] 李俊清, 柴一新, 朱春全, 王業(yè)蘧. 不同光強(qiáng)下紅松幼樹光合作用和營養(yǎng)物質(zhì)含量的季節(jié)模式. 生態(tài)學(xué)雜志, 1991, 10(5): 1- 5.

    [46] 劉傳照, 李俊清, 金奎剛. 林下光照條件與紅松幼樹生長的相關(guān)性研究. 東北林業(yè)大學(xué)學(xué)報, 1991, 19(3): 103- 108.

    [47] 詹鴻振, 劉傳照, 徐桂林, 程文俊, 李茹秀, 李景文. 不同庇蔭的光照特點(diǎn)及對紅松苗的影響. 東北林業(yè)大學(xué)學(xué)報, 1986, 14(2): 25- 33.

    [48] 李晶, 閻秀峰, 祖元剛. 低溫脅迫下紅松幼苗活性氧的產(chǎn)生及保護(hù)酶的變化. 植物學(xué)報, 2000, 42(2): 148- 152.

    [49] 閻秀峰, 李晶, 祖元剛. 干旱脅迫對紅松幼苗保護(hù)酶活性及脂質(zhì)過氧化作用的影響. 生態(tài)學(xué)報, 1999, 19(6): 850- 854.

    [50] 聶修和, 聶宜茂, 聶俊華, 姜黎霞. 光合有效輻射測量原理及其單位換算. 山東農(nóng)業(yè)大學(xué)學(xué)報, 1992, 23(3): 247- 258.

    [51] Vile D, Garnier é, Shipley B, Laurent G R, Marie-Laure N, Roumet C, Lavorel S, DíAz S, Hodgson J G, Lloret F, Midgley G F, Poorter H, Rutherford M C, Wilson P J, Wright I J. Specific leaf area and dry matter content estimate thickness in laminar leaves. Annals of Botany, 2005, 96(6): 1129- 1136.

    [52] 許大全. 光合作用測定及研究中一些值得注意的問題. 植物生理學(xué)通訊, 2006, 42(6): 1163- 1167.

    [53] López-Bucio J, Cruz-Ramírez A, Herrera-Estrella L. The role of nutrient availability in regulating root architecture. Current Opinion in Plant Biology, 2003, 6(3): 280- 287

    [54] Reid K E, Olsson N, Schlosser J, Peng F, Lund S T. An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biology, 2006, 6: 27.

    [55] Santiago L S, Wright S J. Leaf functional traits of tropical forest plants in relation to growth form. Functional Ecology, 2007, 21(1): 19- 27.

    [56] 孫一榮, 朱教君, 于立忠, 閆巧玲, 王凱. 不同光環(huán)境對紅松幼苗光合生理特征的影響. 生態(tài)學(xué)雜志, 2009, 28 (5): 850- 857.

    [57] 張雷明, 曹沛雨, 朱亞平, 李慶康, 張軍輝, 王曉凌, 戴冠華, 李金功. 長白山闊葉紅松林生態(tài)系統(tǒng)光能利用率的動態(tài)變化及其主控因子. 植物生態(tài)學(xué)報, 2015, 39(12): 1156- 1165.

    [58] 韓麗冬, 牟長城, 張軍輝. 透光撫育對長白山闊葉紅松林冠下紅松光合作用的影響. 東北林業(yè)大學(xué)學(xué)報, 2016, 44(4): 38- 40.

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

    [60] 鄭淑霞, 上官周平. 不同功能型植物光合特性及其與葉氮含量、比葉重的關(guān)系. 生態(tài)學(xué)報, 2007, 27(l): 171- 181.

    [61] Garnier E, Gobin O, Poorter H. Nitrogen productivity depends on photosynthetic nitrogen use efficiency and on nitrogen allocation within the plant. Annals of Botany, 1995, 76(6): 667- 672.

    [62] Wright I J, Reich P B, Westoby M. Strategy shifts in leaf physiology, structure and nutrient content between species of high- and low-rainfall and high- and low-nutrient habitats. Functional Ecology, 2001, 15(4): 423- 434.

    [63] Elser J J, Urabe J. The stoichiometry of consumer-driven nutrient recycling: theory, observations, and consequences. Ecology, 1999, 80(3): 735- 751.

    [64] Elser J J, Sterner R W, Gorokhova E, Fagan W F, Markow T A, Cotner J B, Harrison J F, Hobbie S E, Odell G M, Weider L J. Biological stoichiometry from genes to ecosystems. Ecology Letters, 2000, 3(6): 540- 550.

    [65] 許大全. 光合作用氣孔限制分析中的一些問題. 植物生理學(xué)通訊, 1997, 33(4): 241- 244.

    ChangesintheresponseofleaftraitsinPinuskoraiensis(Koreanpine)seedlingsofdifferentagestocontrolledtemperaturesandlightconditions

    DONG Yichen, LIU Yanhong*

    TheKeyLaboratoryforForestResourcesandEcosystemProcessesofBeijing,BeijingForestryUniversity,Beijing100083,China

    國家“十二五”科技支撐計(jì)劃研究(2012BAC01B03)

    2015- 11- 12;

    2016- 08- 24

    10.5846/stxb201511122293

    *通訊作者Corresponding author.E-mail: liuyh@ bjfu.edu.cn

    董伊晨, 劉艷紅.紅松不同苗齡幼苗葉性狀對溫度和光照變化的響應(yīng).生態(tài)學(xué)報,2017,37(17):5662- 5672.

    Dong Y C, Liu Y H.Changes in the response of leaf traits inPinuskoraiensis(Korean pine) seedlings of different ages to controlled temperatures and light conditions.Acta Ecologica Sinica,2017,37(17):5662- 5672.

    猜你喜歡
    紅松林苗齡紅松
    繞口令
    沒有紅松的紅松林
    不同苗齡大花序桉造林效果分析
    不同接穗苗齡對小果型西瓜嫁接育苗效果的影響
    紅松林采伐更新的確定方式
    東折棱河楓樺紅松林與柞樹紅松林結(jié)構(gòu)特征對比分析
    苗齡、光照強(qiáng)度和施氮量對水稻營養(yǎng)物質(zhì)的影響及其與抗褐飛虱的關(guān)系
    選擇紅松寶就是選擇財(cái)富
    選擇紅松寶就是選擇財(cái)富
    高見的詩
    av视频在线观看入口| 亚洲中文字幕一区二区三区有码在线看 | 欧美精品亚洲一区二区| www.自偷自拍.com| 男女做爰动态图高潮gif福利片| 十八禁网站免费在线| 夜夜爽天天搞| 国内精品久久久久久久电影| netflix在线观看网站| 亚洲三区欧美一区| 宅男免费午夜| 高清毛片免费观看视频网站| 中文在线观看免费www的网站 | 午夜免费成人在线视频| 嫩草影院精品99| 在线观看日韩欧美| 欧美日韩亚洲国产一区二区在线观看| 国产1区2区3区精品| 久久精品91蜜桃| 两个人免费观看高清视频| 免费观看人在逋| 欧美色视频一区免费| 久久午夜综合久久蜜桃| 亚洲电影在线观看av| 草草在线视频免费看| 日本免费a在线| 非洲黑人性xxxx精品又粗又长| 91老司机精品| 成人av一区二区三区在线看| 女性生殖器流出的白浆| 精品无人区乱码1区二区| 国产精品电影一区二区三区| 成人国产一区最新在线观看| 国产一区二区三区在线臀色熟女| 国产日本99.免费观看| 亚洲精品av麻豆狂野| 国产一级毛片七仙女欲春2 | 搞女人的毛片| 99re在线观看精品视频| 18美女黄网站色大片免费观看| 色av中文字幕| 黄色成人免费大全| 国产一卡二卡三卡精品| 亚洲色图 男人天堂 中文字幕| 国产片内射在线| 老司机在亚洲福利影院| 丰满的人妻完整版| 又大又爽又粗| 男女那种视频在线观看| 熟女电影av网| 亚洲精品色激情综合| 久久久久国内视频| 精品一区二区三区av网在线观看| 熟女电影av网| 久久久久久久久久黄片| 真人一进一出gif抽搐免费| 久久人妻av系列| av免费在线观看网站| 男女下面进入的视频免费午夜 | 色播在线永久视频| 精品久久久久久,| 久久久国产精品麻豆| 亚洲人成电影免费在线| 看片在线看免费视频| 国产高清视频在线播放一区| 国产精品精品国产色婷婷| 精品一区二区三区四区五区乱码| 国产成+人综合+亚洲专区| 神马国产精品三级电影在线观看 | 久9热在线精品视频| 亚洲人成网站高清观看| www.自偷自拍.com| 免费高清在线观看日韩| 久久久久国产一级毛片高清牌| 免费在线观看视频国产中文字幕亚洲| 中文字幕人成人乱码亚洲影| 中亚洲国语对白在线视频| 色综合欧美亚洲国产小说| 亚洲国产精品合色在线| 国产欧美日韩一区二区三| 国产精品久久视频播放| www.熟女人妻精品国产| 久久久久久免费高清国产稀缺| 久久中文看片网| 在线观看午夜福利视频| 在线播放国产精品三级| 国产一区二区三区视频了| 99精品欧美一区二区三区四区| 91国产中文字幕| 国产成人av激情在线播放| 麻豆久久精品国产亚洲av| 十八禁网站免费在线| 欧美人与性动交α欧美精品济南到| 中文字幕av电影在线播放| 精品国产超薄肉色丝袜足j| 美女免费视频网站| 免费高清视频大片| 久久久久久国产a免费观看| 国产成年人精品一区二区| 淫秽高清视频在线观看| 日韩av在线大香蕉| 狠狠狠狠99中文字幕| av福利片在线| 中文字幕精品亚洲无线码一区 | 国产精品亚洲一级av第二区| 韩国精品一区二区三区| 日韩欧美一区视频在线观看| 深夜精品福利| 搡老妇女老女人老熟妇| 久久久久国产一级毛片高清牌| 婷婷亚洲欧美| 99热只有精品国产| 香蕉国产在线看| 亚洲 欧美一区二区三区| 伊人久久大香线蕉亚洲五| 国产精品av久久久久免费| 在线av久久热| 国产亚洲精品久久久久久毛片| 亚洲男人天堂网一区| 1024手机看黄色片| 成人永久免费在线观看视频| 日韩精品免费视频一区二区三区| tocl精华| 欧美乱码精品一区二区三区| 久热这里只有精品99| 又大又爽又粗| 在线国产一区二区在线| 日韩欧美一区视频在线观看| 国产av又大| 亚洲欧美精品综合一区二区三区| 男人舔女人下体高潮全视频| 精品欧美国产一区二区三| 国产激情欧美一区二区| 久久国产精品男人的天堂亚洲| 一个人观看的视频www高清免费观看 | 欧美黄色片欧美黄色片| 三级毛片av免费| 国产伦一二天堂av在线观看| 亚洲国产看品久久| 一本一本综合久久| 国产亚洲精品一区二区www| 美女午夜性视频免费| 99久久无色码亚洲精品果冻| 啦啦啦韩国在线观看视频| 97碰自拍视频| 免费人成视频x8x8入口观看| 免费一级毛片在线播放高清视频| 日韩精品免费视频一区二区三区| 成人手机av| 热re99久久国产66热| 国产精品久久久av美女十八| 久久精品成人免费网站| 午夜日韩欧美国产| 欧美午夜高清在线| 亚洲人成伊人成综合网2020| 一级黄色大片毛片| 90打野战视频偷拍视频| 黄色丝袜av网址大全| 亚洲国产欧美日韩在线播放| 免费看美女性在线毛片视频| 亚洲在线自拍视频| 97人妻精品一区二区三区麻豆 | 中文字幕精品免费在线观看视频| 久久久久久久久久黄片| 久久久久国产一级毛片高清牌| 日本成人三级电影网站| 国产一区二区三区在线臀色熟女| 国产人伦9x9x在线观看| av中文乱码字幕在线| 一个人观看的视频www高清免费观看 | 日本免费a在线| 伦理电影免费视频| 91成人精品电影| 久久久久久久精品吃奶| 丁香六月欧美| 草草在线视频免费看| 国产熟女午夜一区二区三区| 亚洲国产精品合色在线| www.999成人在线观看| 久久这里只有精品19| 超碰成人久久| 长腿黑丝高跟| 成人亚洲精品av一区二区| 久久精品91蜜桃| 嫩草影院精品99| 日本三级黄在线观看| 麻豆久久精品国产亚洲av| 免费看美女性在线毛片视频| 两个人看的免费小视频| 国产又爽黄色视频| 别揉我奶头~嗯~啊~动态视频| av电影中文网址| 草草在线视频免费看| 麻豆成人午夜福利视频| 国产精品99久久99久久久不卡| 黄色视频不卡| 国产v大片淫在线免费观看| 亚洲中文字幕日韩| 免费女性裸体啪啪无遮挡网站| 中文在线观看免费www的网站 | 国产精品av久久久久免费| 在线播放国产精品三级| 在线视频色国产色| 国产麻豆成人av免费视频| 亚洲午夜精品一区,二区,三区| 亚洲国产欧洲综合997久久, | 91麻豆精品激情在线观看国产| 午夜精品久久久久久毛片777| 韩国av一区二区三区四区| 精品国产美女av久久久久小说| bbb黄色大片| 色老头精品视频在线观看| 亚洲精品美女久久久久99蜜臀| 国产97色在线日韩免费| 熟女电影av网| 国产成人精品久久二区二区91| 国产成人av激情在线播放| 久久人人精品亚洲av| 免费观看人在逋| 狂野欧美激情性xxxx| 亚洲天堂国产精品一区在线| 中亚洲国语对白在线视频| 久久九九热精品免费| 国产高清videossex| 色婷婷久久久亚洲欧美| 免费在线观看影片大全网站| 久久久久久久午夜电影| 国产亚洲欧美98| 国产亚洲精品一区二区www| 日韩精品免费视频一区二区三区| 国产精品av久久久久免费| 久久久久久亚洲精品国产蜜桃av| 91九色精品人成在线观看| 天天躁夜夜躁狠狠躁躁| 天堂动漫精品| av超薄肉色丝袜交足视频| 99久久精品国产亚洲精品| 国产成+人综合+亚洲专区| 色在线成人网| a级毛片在线看网站| 亚洲在线自拍视频| 不卡一级毛片| 日韩三级视频一区二区三区| 欧美色欧美亚洲另类二区| 亚洲一区中文字幕在线| 巨乳人妻的诱惑在线观看| 国产亚洲精品久久久久5区| 桃色一区二区三区在线观看| 757午夜福利合集在线观看| 欧美在线黄色| 久久精品成人免费网站| 50天的宝宝边吃奶边哭怎么回事| 久久久久久免费高清国产稀缺| 午夜福利免费观看在线| 中出人妻视频一区二区| 侵犯人妻中文字幕一二三四区| 欧美激情 高清一区二区三区| 丝袜人妻中文字幕| 91老司机精品| 人人妻,人人澡人人爽秒播| 国产精品av久久久久免费| 特大巨黑吊av在线直播 | 中文字幕人妻丝袜一区二区| 三级毛片av免费| 韩国av一区二区三区四区| 人人妻,人人澡人人爽秒播| 满18在线观看网站| 久久久久精品国产欧美久久久| 欧美zozozo另类| 日韩一卡2卡3卡4卡2021年| 亚洲自偷自拍图片 自拍| 在线观看一区二区三区| 91国产中文字幕| 麻豆一二三区av精品| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 精品不卡国产一区二区三区| 露出奶头的视频| 国产精品久久电影中文字幕| 亚洲国产欧美日韩在线播放| 国产黄色小视频在线观看| 国产aⅴ精品一区二区三区波| 国语自产精品视频在线第100页| 欧美日本亚洲视频在线播放| 日本a在线网址| 男女那种视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 男女床上黄色一级片免费看| 欧美日韩福利视频一区二区| 国产精品乱码一区二三区的特点| 精品人妻1区二区| 18美女黄网站色大片免费观看| 成人精品一区二区免费| 亚洲片人在线观看| 嫁个100分男人电影在线观看| 中文字幕av电影在线播放| 亚洲av成人不卡在线观看播放网| 久久精品国产亚洲av高清一级| 精品久久久久久成人av| 国产精品亚洲一级av第二区| 精品久久久久久久毛片微露脸| 午夜福利成人在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲性夜色夜夜综合| 欧美一级毛片孕妇| 在线观看www视频免费| 黑丝袜美女国产一区| www.精华液| 不卡一级毛片| 一夜夜www| 淫妇啪啪啪对白视频| 午夜两性在线视频| 99精品欧美一区二区三区四区| 美女 人体艺术 gogo| 欧美精品啪啪一区二区三区| 久久精品国产清高在天天线| 最近在线观看免费完整版| 久久国产亚洲av麻豆专区| 日本五十路高清| 一夜夜www| 男女那种视频在线观看| 精品人妻1区二区| 成人午夜高清在线视频 | 好看av亚洲va欧美ⅴa在| 啦啦啦免费观看视频1| 青草久久国产| 国产激情偷乱视频一区二区| 久久精品夜夜夜夜夜久久蜜豆 | 一个人免费在线观看的高清视频| 成熟少妇高潮喷水视频| 一区二区三区高清视频在线| 免费女性裸体啪啪无遮挡网站| 久久国产精品人妻蜜桃| 999久久久国产精品视频| 亚洲av成人av| 精品一区二区三区四区五区乱码| 757午夜福利合集在线观看| 97碰自拍视频| 亚洲久久久国产精品| aaaaa片日本免费| 香蕉久久夜色| 女人高潮潮喷娇喘18禁视频| 亚洲国产精品成人综合色| 日韩大尺度精品在线看网址| 午夜福利在线在线| 一夜夜www| 色尼玛亚洲综合影院| 国产一级毛片七仙女欲春2 | 久久热在线av| 男人舔女人的私密视频| 日韩欧美免费精品| 欧美日韩亚洲国产一区二区在线观看| 亚洲 欧美 日韩 在线 免费| 色综合站精品国产| 国产国语露脸激情在线看| 精品卡一卡二卡四卡免费| 欧美+亚洲+日韩+国产| 黑丝袜美女国产一区| 桃色一区二区三区在线观看| 超碰成人久久| 婷婷六月久久综合丁香| 免费看十八禁软件| 香蕉国产在线看| 久久青草综合色| 国产亚洲精品久久久久久毛片| 免费电影在线观看免费观看| 俄罗斯特黄特色一大片| 日韩一卡2卡3卡4卡2021年| 男男h啪啪无遮挡| 美女国产高潮福利片在线看| 18禁黄网站禁片免费观看直播| 这个男人来自地球电影免费观看| 亚洲国产精品久久男人天堂| 国产精品国产高清国产av| 免费电影在线观看免费观看| 一本一本综合久久| 黑人操中国人逼视频| 色在线成人网| 香蕉国产在线看| 免费一级毛片在线播放高清视频| 2021天堂中文幕一二区在线观 | 50天的宝宝边吃奶边哭怎么回事| 18禁观看日本| 满18在线观看网站| 亚洲天堂国产精品一区在线| 亚洲精品国产区一区二| 久久久久久国产a免费观看| 成人国语在线视频| 老司机在亚洲福利影院| 国产精品,欧美在线| 女警被强在线播放| 亚洲va日本ⅴa欧美va伊人久久| 啦啦啦免费观看视频1| 俺也久久电影网| 日日夜夜操网爽| 国产三级在线视频| 欧美在线一区亚洲| 亚洲人成网站高清观看| 18禁观看日本| 久久中文字幕一级| 丰满人妻熟妇乱又伦精品不卡| 精品福利观看| 国内揄拍国产精品人妻在线 | 一区二区三区精品91| 中文字幕最新亚洲高清| 亚洲成人免费电影在线观看| 午夜福利免费观看在线| 午夜福利成人在线免费观看| 老司机深夜福利视频在线观看| 国产91精品成人一区二区三区| 国产欧美日韩一区二区精品| 国产精品野战在线观看| 亚洲国产精品999在线| 久久久久久久精品吃奶| 国产99久久九九免费精品| 此物有八面人人有两片| 满18在线观看网站| 国产又黄又爽又无遮挡在线| ponron亚洲| 亚洲一码二码三码区别大吗| 日本一本二区三区精品| 精品人妻1区二区| cao死你这个sao货| 我的亚洲天堂| 亚洲精品一区av在线观看| 色精品久久人妻99蜜桃| 日本一本二区三区精品| 午夜福利在线在线| av视频在线观看入口| xxx96com| 99国产极品粉嫩在线观看| av超薄肉色丝袜交足视频| 成人欧美大片| 一a级毛片在线观看| 这个男人来自地球电影免费观看| 看片在线看免费视频| 一级a爱片免费观看的视频| 国产精品久久视频播放| 亚洲 欧美一区二区三区| 18禁黄网站禁片免费观看直播| 免费人成视频x8x8入口观看| 黄色a级毛片大全视频| 成人欧美大片| 9191精品国产免费久久| 日韩欧美一区视频在线观看| 精品一区二区三区四区五区乱码| 亚洲五月色婷婷综合| 天天添夜夜摸| 最新美女视频免费是黄的| 丁香欧美五月| 国产真实乱freesex| 免费在线观看视频国产中文字幕亚洲| 国产精品久久久久久精品电影 | 老熟妇乱子伦视频在线观看| 久久中文看片网| 日韩一卡2卡3卡4卡2021年| 国语自产精品视频在线第100页| 国产av不卡久久| 美女扒开内裤让男人捅视频| 日韩欧美 国产精品| 69av精品久久久久久| 一级毛片精品| 国产主播在线观看一区二区| 日韩有码中文字幕| 白带黄色成豆腐渣| 亚洲九九香蕉| 淫妇啪啪啪对白视频| av欧美777| 免费在线观看成人毛片| 精品少妇一区二区三区视频日本电影| 女性生殖器流出的白浆| 午夜成年电影在线免费观看| 久久久久亚洲av毛片大全| 久久亚洲真实| 亚洲国产欧美网| 草草在线视频免费看| 欧美日韩福利视频一区二区| 婷婷亚洲欧美| 天堂√8在线中文| 法律面前人人平等表现在哪些方面| 老熟妇仑乱视频hdxx| 两个人视频免费观看高清| 免费观看精品视频网站| 精品少妇一区二区三区视频日本电影| 国产精品一区二区三区四区久久 | 亚洲人成网站高清观看| 中文字幕人成人乱码亚洲影| 三级毛片av免费| 欧美一区二区精品小视频在线| 久久精品91无色码中文字幕| 可以在线观看毛片的网站| 一区福利在线观看| 狂野欧美激情性xxxx| 精品人妻1区二区| 久久久国产成人免费| 女同久久另类99精品国产91| 日日干狠狠操夜夜爽| 国产成人精品久久二区二区免费| 一进一出好大好爽视频| 老司机午夜福利在线观看视频| 好男人电影高清在线观看| 日韩高清综合在线| 好看av亚洲va欧美ⅴa在| 亚洲熟妇熟女久久| 观看免费一级毛片| 一区福利在线观看| 在线av久久热| 大型黄色视频在线免费观看| 欧美成人一区二区免费高清观看 | 精品久久久久久久毛片微露脸| 最新美女视频免费是黄的| 一区二区日韩欧美中文字幕| 免费在线观看亚洲国产| 午夜视频精品福利| 一个人观看的视频www高清免费观看 | 91麻豆av在线| 国产爱豆传媒在线观看 | 免费观看人在逋| 免费在线观看日本一区| 欧美在线黄色| 黄色片一级片一级黄色片| 欧美日韩精品网址| or卡值多少钱| 欧美久久黑人一区二区| 久久国产精品男人的天堂亚洲| 99久久国产精品久久久| 黄色成人免费大全| a在线观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 国产伦一二天堂av在线观看| 亚洲精品色激情综合| 国产成人啪精品午夜网站| 好男人电影高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看66精品国产| 两个人视频免费观看高清| 亚洲国产日韩欧美精品在线观看 | 久久精品夜夜夜夜夜久久蜜豆 | 欧美日韩亚洲综合一区二区三区_| 夜夜躁狠狠躁天天躁| 黄色女人牲交| 亚洲欧美精品综合久久99| 午夜福利在线在线| 成人手机av| 12—13女人毛片做爰片一| 亚洲五月天丁香| 成人国产综合亚洲| 嫩草影院精品99| 黄色毛片三级朝国网站| 天堂√8在线中文| 国产精品一区二区免费欧美| 黄色女人牲交| 精品久久久久久成人av| 无遮挡黄片免费观看| 国产成人欧美在线观看| 无遮挡黄片免费观看| 99久久精品国产亚洲精品| 久久久久久九九精品二区国产 | 男人的好看免费观看在线视频 | 国产伦一二天堂av在线观看| 黄色视频,在线免费观看| 中文字幕高清在线视频| 日韩三级视频一区二区三区| 淫妇啪啪啪对白视频| 国产精品免费一区二区三区在线| 国产片内射在线| 久久天堂一区二区三区四区| √禁漫天堂资源中文www| 美女免费视频网站| aaaaa片日本免费| www.自偷自拍.com| 一进一出抽搐动态| 大型黄色视频在线免费观看| 欧美日韩一级在线毛片| 麻豆久久精品国产亚洲av| tocl精华| 亚洲中文字幕一区二区三区有码在线看 | 免费在线观看成人毛片| 老司机靠b影院| 一级作爱视频免费观看| 亚洲中文字幕一区二区三区有码在线看 | 婷婷丁香在线五月| 国产在线观看jvid| 成熟少妇高潮喷水视频| 国产又色又爽无遮挡免费看| 婷婷丁香在线五月| 久久狼人影院| 美女扒开内裤让男人捅视频| 男人舔女人下体高潮全视频| 国产成人啪精品午夜网站| 亚洲激情在线av| 99精品在免费线老司机午夜| 91成年电影在线观看| 中文资源天堂在线| 男女之事视频高清在线观看| 曰老女人黄片| 亚洲五月天丁香| 亚洲三区欧美一区| 成人手机av| 黄色 视频免费看| 久久香蕉激情| 神马国产精品三级电影在线观看 | 一边摸一边抽搐一进一小说| 亚洲成人免费电影在线观看| 在线观看日韩欧美| 免费看a级黄色片| 久久精品国产99精品国产亚洲性色| 巨乳人妻的诱惑在线观看| 国产片内射在线| 欧美黄色片欧美黄色片| 50天的宝宝边吃奶边哭怎么回事| 久久香蕉国产精品| 亚洲国产欧美网| 久久中文字幕一级|