影響科爾沁沙地不同林齡樟子松徑向生長的關(guān)鍵氣候因子

中圖分類號：Q948 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號：1000-3142（2025）08-1405-10

Abstract：Toidentifythekeytemperatureandprecipitation factorsinfluencing thegrowthofPinussylvestris var. mongolica at diferent standages，this study focused on middle-aged（33a），mature（48a），and overmature（64 a） （204號 P .sylvestris var.mongolicastandsin Zhanggutai，Liaoning Province，andanalyzed the growth patterns（basal area increment，BAI）of these stands during 2003-2022 .Inaddition，the Climwin model was used to determine the temperatureand precipitationseasonal windows withthe strongest influence on growth at different stand ages.The results were as follows：（1）The BAI of the three stands significantly increased from 2003 to 2022.（2）The current JulyAugust mean temperature had the strongest influence on the growth of middle-aged stand.When the July ^- August mean temperature reached 23.7 ΔC ，the growth-climate relationship of middle-aged stand turned from positive to negative; precipitation from last November to current May showed the strongest impactonboth middle-agedand mature stands.The thresholdsofprecipitationat whichthegrowth-climaterelationshipof midle-agedand maturestandsshifted frompositive to negative were 124 mmand122 mm，respectively.（3）Temperaturehadno significant impacton the growthof mature stand;conversely，thetemperatureinthecurrent Julyandthe precipitation fromlast August tocurrent Maywere significantlypositivelycorelated withthegrowthofovermaturestand.（4）Precipitation was proved tobemore important than temperature across all stand ages，with its relative importance increasing as the stands aged.Therefore，inthe context of a warming and drying climate trend，developing forest management strategies that consider water balance and varywith stand age will facilitate the growth of P .sylvestrisvar.mongolicaplantations.

Key words：Pinus sylvestris var. mongolica，growth-climate relationship，key climate factors，seasonal window，stand age difference

樟子松（Pinus sylvestrisvar.mongolica）是我國北方干旱半干旱沙地主要的防風(fēng)固沙樹種，在生態(tài)恢復(fù)、固碳增匯和改善社會(huì)宜居環(huán)境方面貢獻(xiàn)突出（朱教君和鄭曉，2019）。然而，由于管理措施未能適應(yīng)當(dāng)?shù)貧夂颦h(huán)境，加上氣候變化等多因素的共同作用，樟子松存在生長衰退和成片死亡等退化問題（于豐源等，2022）。研究樟子松生長對氣候變化的響應(yīng)有助于預(yù)測森林生長變化并為森林可持續(xù)經(jīng)營提供科學(xué)依據(jù)，遂成為研究熱點(diǎn)（Sun etal.，2021）。

降水和溫度是決定沙地樟子松生長和分布的關(guān)鍵氣候因子（張世林等，2023），也是建立樹木生長模型最常用的氣候預(yù)測變量。在沙地環(huán)境中，降水可調(diào)節(jié)土壤和空氣水分條件進(jìn)而調(diào)控樟子松水分可利用性，溫度可以通過改變植物生理活動(dòng)和調(diào)節(jié)水文循環(huán)來影響植物的生長（Wenetal.，2023）。已有研究報(bào)道了對沙地樟子松生長具有顯著影響的降水和溫度的月份（Songetal.，2017;Sunetal.，2021），但受研究方法或所選樣本的影響，不同研究最終確定的月份存在差異。同時(shí)，這些研究沒有明確對樹木生長影響最強(qiáng)的溫度和降水因子季節(jié)窗口，導(dǎo)致構(gòu)建包含氣候因子的樹木生長模型時(shí)需考慮的潛在氣候變量較多，而且所用季節(jié)窗口不確定，在一定程度上增加了模型構(gòu)建難度。此外，這些研究大多重點(diǎn)分析樹木生長與氣候間的線性關(guān)系，而自然中樹木生長與氣候關(guān)系不穩(wěn)定，可能為非線性關(guān)系，甚至存在生長-氣候關(guān)系發(fā)生轉(zhuǎn)變的氣候閥值，這一偏差可能影響氣候變化背景下樹木生長的預(yù)測，但相關(guān)研究較少（Anderson-Teixeira etal.，2022；Xue et al.，2023）。

樟子松徑向生長與氣候關(guān)系存在林齡差異。例如，張曉等（2022）發(fā)現(xiàn)遼寧章古臺(tái)地區(qū)樟子松徑向生長對氣候的敏感性隨林齡增大而增強(qiáng)，并且發(fā)現(xiàn)影響低齡樹木生長的降水季節(jié)比高齡樹木晚；周鳳艷等（2025）發(fā)現(xiàn)，影響不同林齡樟子松生長的極端氣候因子發(fā)生在不同季節(jié)；Liu等（2023）發(fā)現(xiàn)中齡樟子松的氣候適應(yīng)性明顯強(qiáng)于幼齡。樟子松生長對氣候因子的響應(yīng)模式（線性vs非線性）是否隨林齡發(fā)生改變，以及影響不同林齡的關(guān)鍵氣候因子是否存在差異，有待進(jìn)一步探索。

本文以遼寧章古臺(tái)樟子松中齡林（ ^?33a） 、成熟林（48a）和過熟林（64a）為例，利用樹木年輪學(xué)方法，旨在探討以下3個(gè)問題：（1）樟子松生長速度和趨勢的林齡差異；（2）影響不同林齡樟子松徑向生長的關(guān)鍵溫度和降水發(fā)生的季節(jié)窗口；（3）關(guān)鍵溫度和降水因子的相對重要性隨林齡變化特征。

1材料與方法

1.1研究區(qū)概況

研究區(qū)位于遼寧省阜新市章古臺(tái)鎮(zhèn)L （122^°22^′-123^°22^′E，42^°43^′-43^°20^′N ，平均海拔 226.5m ），該地區(qū)位于科爾沁沙地東南緣，屬于溫帶亞濕潤大陸性季風(fēng)氣候（劉軒等，2022），1961—2022年平均溫度 7% ，平均年降水量518mm 。章古臺(tái)地區(qū)雨熱同季，四季分明，夏季炎熱相對多雨，冬季寒冷干燥，春冬兩季多大風(fēng)天氣。1月是最冷月，溫度范圍為 -18.3～-9.0%，7 月為最熱月，溫度范圍為 22.2～26.1°C （圖1）。降雨主要集中在6一8月（圖1）。土壤主要為風(fēng)沙土，土壤有機(jī)質(zhì)、氮、磷、鉀等養(yǎng)分含量較低。地貌類型有流動(dòng)沙丘、半固定沙丘、固定沙丘和平緩沙地等。主要造林樹種有樟子松、油松（Pinustabuliformis）、楊樹屬（Populusspp.）、榆樹（Ulmuspumila）、五角楓（Acerpictum subsp.mono）等耐性樹種。主要灌木物種有小黃柳（Salixgordejevii）、差巴嘎蒿（Artemisiahalodendron）、小葉錦雞兒（Caraganamicrophylla）等（王景璽等，2022）。林下主要為1年生草本植物，常見有狗尾草（Setariaviridis）、馬唐（Digitariasanguinalis）、中華隱子草（Cleistogeneschinensis）（韓艷剛等，2018；周晏平等，2018）。

圖11961—2022年研究區(qū)氣候特征

Fig.1Climate characteristics of the study area from 1961 to 2022

1.2標(biāo)準(zhǔn)地設(shè)置與年輪樣本采集

于2023年4月，在片狀樟子松人工林中選擇密度和長勢均較典型中齡林（33a）、成熟林（48a）和過熟林（64a）（表1），各設(shè)置1塊 20m×20 m 臨時(shí)標(biāo)準(zhǔn)地。對標(biāo)準(zhǔn)地內(nèi)所有樟子松進(jìn)行每木檢尺后，在胸高（ ）處沿4個(gè)方向鉆取樹芯。使用砂紙將樹芯打磨至所有年輪清晰可見后，在放大鏡下交叉定年。使用LINTAB6年輪儀測量每個(gè)年輪寬度（精度 0.01mm ）。使用COFECHA程序檢驗(yàn)交叉定年結(jié)果，去除無法定年的樣芯后共保留72棵樹的254個(gè)樹芯用于進(jìn)一步分析（表1）。考慮到中齡林和成熟林未經(jīng)過撫育間伐，而過熟林在2001—2003年嚴(yán)重干旱后進(jìn)行了隔行間伐，并且所有樣本最后的完整年輪形成于2022年，因此本研究的分析時(shí)段為2003—2022年。我們的前期研究發(fā)現(xiàn)，本研究所用的年輪樣本符合樹木年輪學(xué)分析相關(guān)要求，適合分析樹木生長對氣候因子的響應(yīng)（周鳳艷等，2025）。考慮到胸高斷面積增量（basalarea increment，BAI）包含更多的生長信息，與生物量密切相關(guān)，因此本研究利用BAI表示生產(chǎn)力，利用林分總BAI表示林分生產(chǎn)力（Labrecque-Foyetal.，2023；李瓊瑩等，2024），計(jì)算公式如下：

表1不同林齡樟子松林分特征 Table1Stand characteristics for Pinus sylvestris var. mongolica with different stand ages

注：A.林齡；D.算術(shù)平均胸徑；MT.平均樹高； BAI_tree .2003—2022年樹木年平均BAI; BAI_stand .2003—2022年林分平均 BAI 數(shù)值為平均值 ?± 標(biāo)準(zhǔn)誤。不同字母表示林齡間差異顯著（ Plt;0.05? 。下同。Note：A. Stand age;D.Arithmetic mean diameter at breast height；MT.Mean tree height; BAI_tree . Average BAI of trees from 2003 to2022 ; BAI_stand .Mean stand BAI from 2003 to 2022.Values are . Different letters indicate significant differences among stand ages （ Plt; 0.05）.The same below.

式中： R_t 和 R_t-1 分別為年份 χ_t 和t-1時(shí)的樹木半徑。

2 結(jié)果與分析

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

日均溫和降水量來源于中國氣象局氣候研究開放實(shí)驗(yàn)室公開的格點(diǎn)數(shù)據(jù)集（ 0.25^°×0.25^° ）（ Xu etal.，2009；吳佳和高興杰，2013；Wuetal.，2017），計(jì)算后得到各月平均溫（Tmean）和總降水量（Prcp）。針對每個(gè)林齡樟子松，使用R軟件中的Climwin氣候響應(yīng)模型識(shí)別影響樹木生長的關(guān)鍵溫度和降水窗口。參照樹木生長與氣候關(guān)系的散點(diǎn)分布圖及Anderson-teixeira等（2022）方法，最終選擇線性函數(shù)和一元二次函數(shù)分別表示生長與氣候因子間的線性和非線性關(guān)系。考慮到樹木生長可能受當(dāng)年和上一年氣候因子的影響，所選的月份為上一年6月至當(dāng)年9月。Climwin模型在該時(shí)段內(nèi)運(yùn)行所有可能的溫度和降水連續(xù)月份組合的時(shí)間窗口，篩選出△AICc值最小的組合作為最佳模型，進(jìn)而確定對樟子松生長影響最強(qiáng)的溫度和降水的時(shí)間窗口（簡稱關(guān)鍵溫度和降水因子）（van dePol et al.，2016）。

進(jìn)一步使用廣義可加模型（generalizedadditivemodels，GAM）分析BAI與關(guān)鍵降水和溫度因子間關(guān)系，該模型可以擬合生長與各影響因子間的線性和非線性關(guān)系，具有較高的靈活度（式2）。使用R語言mgcv程序包中g(shù)am函數(shù)進(jìn)行模型擬合和參數(shù)估計(jì)，并使用gam.hp函數(shù)確定2個(gè)變量對生長的相對重要性（Laietal.，2022）。

BAI=β₀+f₁（Tmean）+f₂（Prcp）

式中： β₀ 為截距項(xiàng)： f₁ 和 f₂ 分別為因子Tmean和 Prep 對應(yīng)的光滑擬合函數(shù)。

2.1不同林齡樟子松生長

2003—2022年，3個(gè)林齡樟子松樹木和林分生產(chǎn)力均顯著升高（圖2）。就樹木生產(chǎn)力而言，3個(gè)林齡樹木生產(chǎn)力無顯著性差異（表1）。林分生產(chǎn)力方面，未經(jīng)撫育間伐的中齡林和成熟林間無顯著性差異，經(jīng)過間伐的過熟林總生產(chǎn)力顯著低于中齡林和成熟林（表1）。

2.2氣候因子對樟子松生長的影響

影響樟子松生長的關(guān)鍵氣候因子季節(jié)窗口，以及關(guān)鍵氣候因子與生長的函數(shù)關(guān)系，存在明顯林齡差異（圖3、圖4）。關(guān)鍵溫度因子的季節(jié)窗口較短且均出現(xiàn)在當(dāng)年（圖3）。最佳溫度窗口與中齡林、成熟林和過熟林生長回歸模型的 R² 分別為0.41、0.16和0.37。當(dāng)年7—8月平均溫與中齡林生長為非線性關(guān)系且存在生長與氣溫由正轉(zhuǎn)負(fù)的溫度閾值（ 23.7^°C ；圖4：A）。當(dāng)年3月溫度與成熟林生長呈線性關(guān)系，但未達(dá)到顯著水平（ P= 0.077；圖4：B）。當(dāng)年7月平均溫與過熟林生長呈顯著線性正相關(guān)（ P=0.005 ;圖4：C）。

對3個(gè)林齡樟子松生長影響最強(qiáng)的降水因子包括上一年和當(dāng)年部分月份，而且所有林齡的響應(yīng)均顯著（ Plt;0.05 ;圖4）。上一年11月至當(dāng)年5月降水量與中齡林和成熟林生長關(guān)系為非線性，對應(yīng)模型的 R² 分別是0.72和0.55，2個(gè)林齡樹木生長與降水量關(guān)系由正轉(zhuǎn)負(fù)的閾值分別為124mm 和 122mm （ Plt;0.001 和 P=0.001 ;圖4：D、E）。上一年8月至當(dāng)年5月降水量與過熟林生長為線性正相關(guān)，回歸模型 R² 為0.46（圖4：F）。

圖22003—2022年不同林齡樟子松生產(chǎn)力變化Fig.2Variations in productivity of Pinus sylvestris var.mongolica with diferent stand ages from 2OO3 to 2022

A.林分總生產(chǎn)力；B.平均單木生產(chǎn)力。 表示 Plt;0.05 ， ^** 表示 Plt;0.01 ， *** 表示 Plt;0.001 。A. Total stand productivity；B. Average individual tree productivity. * indicates Plt;0.05 ， ^** indicates Plt;0.01 ，and ?=??? indicates Plt; （204號0.001.

圖3Climwin模型的AICc值Fig.3△AICc values for Climwin models

A-C 和D-F分別表示33、48、64a林分的平均溫度和降水結(jié)果。坐標(biāo)軸中字母 p 和c分別表示年輪形成的上一年和當(dāng)年月份。 下同。

A-C and D-F indicate mean temperature and precipitation results for 33，48，64 a stands，respectively；letters p and c on the coordinate axes indicate the previous and current year of tree ring formation. The same below.

圖4不同林齡樹木生長對關(guān)鍵氣候因子的響應(yīng)Fig. 4Response of tree growth to the key climate factors at different stand ages

圖5基于廣義可加模型的關(guān)鍵溫度和降水因子對徑向生長的相對重要性

Fig.5Relative importance of the key temperature andprecipitation factorsonradial growth based on generalized additive model

2.3溫度和降雨對樟子松生長的相對重要性

使用上述確定的關(guān)鍵溫度和降水因子為預(yù)測變量時(shí)，GAM模型解釋了3個(gè)林齡樟子松2003—2022年間BAI年際變化的 0.55～0.72 （表2）。中齡林和成熟林生長與降水因子呈顯著非線性關(guān)系（估計(jì)自由度 gt;1 Plt;0.05 ）。過熟林生長與關(guān)鍵溫度和降水因子均呈顯著線性關(guān)系（估計(jì)自由度 Σ=Σ 1，Plt;0.05） （表2）。

關(guān)鍵降水因子對3個(gè)林齡樟子松生長的影響大于溫度因子且隨著林齡增加，降水的相對重要性也逐漸增加（圖5）。具體而言，降水因子對中齡林、成熟林和過熟林生長的重要性分別是溫度因子的2.2倍、4.7倍和6.0倍（圖5）。

3討論

3.1不同林齡樟子松生長特征

本研究發(fā)現(xiàn)，2003—2022年期間3個(gè)林齡樟

表2廣義可加模型擬合結(jié)果

Table 2Fitting results for generalized additive model

注：A.林齡；Tmean.平均溫；Prcp.總降水量。 * 表示 Plt;0.05 ， ^** 表示 Plt;0.01 ， *** 表示 Plt;0.001 。 f₁ 和 f₂ 為對應(yīng)變量擬合函數(shù)。Note：A. Stand age；Tmean. Mean temperature；Prcp. Total precipitation. * indicates Plt;0.05 ， ** indicates Plt;0.01 ，and ***indicates Plt;0.001 f₁ and f₂ are fitting functions for responding variables.

子松（對應(yīng)生長年齡分別為14～33a、29～48 a、45～64a 林分生產(chǎn)力和樹木生產(chǎn)力均顯著上升，顯示不同林齡樟子松都具有較大生長潛力，與雷澤勇等（2018）認(rèn)為遼寧章古臺(tái)樟子松14年后胸徑連年生長量下降的結(jié)果不同，還與以往研究中科爾沁樟子松隨林齡增大普遍生長衰退存在一定差異（宋立寧等，2009，2017;Song etal.，2019）。這一結(jié)果的原因主要包括以下3個(gè)方面。（1）氣候條件改善：研究階段內(nèi)（2003—2022年），章古臺(tái)地區(qū)氣溫變暖減緩，降水量和相對濕度增加，風(fēng)速降低（韓輝等，2023），同時(shí)本研究所判定的關(guān)鍵季節(jié)降水量也呈上升趨勢（ Prcp_p11-?5 和 Prcp_p8-c5 年均變化率分別為 1.5mm 和 8.7mm ），減輕了樟子松受到的干旱脅迫，促進(jìn)了樟子松生長。（2）經(jīng)營管理措施優(yōu)化：撫育間伐和病蟲害防治等多種措施有效提升了現(xiàn)有樹木生長速率和生長穩(wěn)定性，延遲了樟子松生長速率下降的年齡（Duanetal.，2023）。（3）選擇指標(biāo)差異：BAI反映了樹木實(shí)際生產(chǎn)力，其對氣候和非氣候等多因子響應(yīng)敏感，無量綱化年輪指數(shù)主要反映生長受氣候因子的影響，難以反映其他因素的影響。例如，我們使用年輪指數(shù)分析這些樟子松生長時(shí)，發(fā)現(xiàn)樟子松生長并未發(fā)生顯著變化（周鳳艷等，2025）。

3.2不同林齡樟子松生長對關(guān)鍵溫度和降水因子的響應(yīng)

盡管3個(gè)林齡樟子松生長受降水和溫度的共同影響，但降水因子對樟子松生長的重要性強(qiáng)于溫度因子且其重要性隨著林齡增加而逐漸增強(qiáng)。章古臺(tái)地區(qū)植物生長所利用的土壤水分主要靠降水補(bǔ)給。隨著林齡增加，樟子松人工林的生物量和冠幅均增大，因而加大蒸騰量（韓艷剛等，2018；童郁強(qiáng)，2023），同時(shí)木質(zhì)部導(dǎo)水性能降低（Liuetal.，2018），這些變化可能使樹木生長對降水響應(yīng)敏感，該結(jié)果與半干旱區(qū)側(cè)柏（Platycladusorientalis）、油松和云杉（Piceacrassifolia）等結(jié)果一致（趙瑩等，2021；黨宏忠等，2022；SuJetal.，2024）。本研究結(jié)果驗(yàn)證了“以水定綠”沙地生態(tài)修復(fù)理論的正確性（宋立寧等，2017；黨宏忠等，2022）。

不同林齡樟子松生長對氣候因子的響應(yīng)模式也存在差異，中齡林和成熟林樟子松生長與關(guān)鍵溫度和降水因子的關(guān)系主要為非線性，而過熟林樟子松生長與溫度和降水為線性關(guān)系，表明生長與氣候關(guān)系隨樹木年齡存在不穩(wěn)定性，這與劉亞靜等（2022）研究結(jié)果一致。中齡林和成熟林樟子松生長與上一年11月至當(dāng)年5月降水量呈顯著非線性關(guān)系，并且存在生長與降水關(guān)系由正轉(zhuǎn)負(fù)的降水閾值。這是因?yàn)檎鹿排_(tái)地區(qū)樟子松生長季開始于5月上旬（劉軒等，2022），上一年11月至當(dāng)年5月包括了樟子松的休眠季和生長季早期，休眠季充足的降水能夠增加積雪覆蓋，保護(hù)樹木免受凍害，減少凍融產(chǎn)生的木質(zhì)部導(dǎo)管氣穴化栓塞（Charrieretal.，2021），同時(shí)可殺死蟲害促進(jìn)樟子松生長。進(jìn)一步，該時(shí)段內(nèi)降水增加可有效供給樟子松生長季早期所需的土壤水分，提前達(dá)到形成層細(xì)胞開始活動(dòng)的水分閥值（Renetal.，2015），并加快細(xì)胞生長和分裂，進(jìn)而容易形成寬輪（Heldetal.，2021）。然而，過量降水可能加劇早春融雪和降雨造成的土壤侵蝕（SuLetal.，2024），導(dǎo)致養(yǎng)分流失，同時(shí)過多的降雨意味陰天時(shí)間增加，減少光合有效輻射，降低光合速率，不利于樹木生長。對于過熟林，由于林齡增大伴隨的樹冠蒸騰耗水量和耗水速率增強(qiáng)（童郁強(qiáng)，2023），樹木對水分依賴性也增強(qiáng)（Songetal.，2016），因此上一年8月至當(dāng)年5月與樟子松過熟林生長呈線性正相關(guān)。此外，可能現(xiàn)階段降水量并未達(dá)到過熟林樟子松生長與降水負(fù)相關(guān)的氣候閾值，導(dǎo)致生長與降水關(guān)系為顯著正相關(guān)。

溫度對不同林齡樹木生長的影響存在明顯差異。其中，成熟林樟子松生長與溫度無顯著相關(guān)性。中齡林樟子松生長與當(dāng)年7一8月平均溫為非線性關(guān)系，生長與關(guān)鍵溫度因子的關(guān)系存在由正轉(zhuǎn)負(fù)的閾值（ 23.7°C ）。一定范圍內(nèi)升溫有利于細(xì)胞分裂和光合作用，加速樹木生長。但溫度較高會(huì)加劇林地蒸散發(fā)，造成土壤水分虧缺，同時(shí)可能造成樹木高溫脅迫，進(jìn)而限制樹木生長（Lietal.，2020）。此外，我們發(fā)現(xiàn)過熟林BAI與7月溫度顯著正相關(guān)，相同年輪樣本所建立的年輪指數(shù)也與7月溫度正相關(guān)（周鳳艷等，2025），這不同于以往多數(shù)研究結(jié)果中樟子松生長與7月溫度負(fù)相關(guān)（Lietal.，2024；李雅婧等，2025），然而由于未找到可參考的研究，我們暫時(shí)無法給出準(zhǔn)確的解釋，這有待于在后續(xù)結(jié)合同位素和樹干液流等生理指標(biāo)解析相關(guān)機(jī)制。

3.3本研究的局限性

本研究成功確定了影響不同林齡樟子松生長的關(guān)鍵氣候因子—溫度和降水的最強(qiáng)季節(jié)窗口，這一發(fā)現(xiàn)對于構(gòu)建樟子松生長模型和預(yù)測未來氣候變化下的樟子松生產(chǎn)力具有重要的實(shí)際應(yīng)用價(jià)值（李明陽等，2023）。然而，影響樹木生長及其對氣候因子響應(yīng)關(guān)系的因子較多，本研究未涉及會(huì)影響生長的密度、撫育措施等林分因子，以及其他氣候環(huán)境因子中的水分因子，如地下水、地表水及空氣中的水分，也未考慮降水形式（如積雪與降雨）、降水強(qiáng)度及頻率等變量（Gebaueretal.，2023）。此外，極端氣候因子也未在研究中考慮，這些極端氣候因子可能對樟子松的生長有潛在影響。因此，未來研究應(yīng)納入這些多樣化的水分和溫度變量，以更精準(zhǔn)地識(shí)別影響樟子松生長的關(guān)鍵氣候因子。這將有助于全面理解樟子松對氣候變化的適應(yīng)性與響應(yīng)機(jī)制，為林業(yè)管理和應(yīng)對氣候變化提供更有效的理論支撐。

4結(jié)論

遼寧章古臺(tái)地區(qū)樟子松人工林2003—2022 年生產(chǎn)力顯著升高。樟子松生長對氣候的響應(yīng)受林齡和氣候因子的共同影響。中齡林與影響其生長的關(guān)鍵溫度和降水因子之間，以及成熟林與關(guān)鍵降水因子之間，均為顯著非線性關(guān)系，生長-氣候關(guān)系存在由正轉(zhuǎn)負(fù)的氣候閾值。過熟林生長與關(guān)鍵溫度和降水因子呈顯著線性正相關(guān)。降水對不同林齡樟子松生長的重要性強(qiáng)于溫度。隨著林齡增加，樹木生長對降水的依賴性逐漸增強(qiáng)。在未來氣候暖干化趨勢下，應(yīng)側(cè)重于制定針對林齡變化的林木水分平衡管理措施，以促進(jìn)樟子松人工林的生長和持續(xù)發(fā)展。

參考文獻(xiàn)：

ANDERSON-TEIXEIRA KJ，HERRMANN V，ROLLINSON CR，et al.，2O22. Joint effects ofclimate，tree size，and year on annual tree growth derived from tree-ring records of ten globally distributed forests [J].Global Change Biology， 28（1）：245-266.

BABST F，BOURIAUD O，POULTER B，etal.，2019. Twentieth centuryredistribution inclimatic driversof global tree growth[J].Science Advances，5（1）：eaat4313.

CHARRIERG，MARTIN-STPAUL N，DAMESIN C， et al.， 2021.Interaction of drought and frost in tree ecophysiology ： Rethinking the timing of risks[J].Annals of Forest Science，78（2）：40.

DANG HZ，ZHANG XL，HAN H，et al.，2O22. Research advances on forest-water relationships in Pinus sylvestris var.mongolica plantations for sand dune immobilization and guidance to forest management practices [J]. Chinese Journal ofPlantEcology，46（9）：971-983.[黨宏忠，張學(xué)利，韓輝， 等，2022.樟子松固沙林林水關(guān)系研究進(jìn)展及對營林實(shí)踐 的指導(dǎo)[J].植物生態(tài)學(xué)報(bào)，46（9）：971-983.

DUAN CY，LI MY，ZHANG C，et al.，2023.Thinning effectively mitigates the decline of aging Mongolian pine plantationsby alleviatingdrought stressand enhancing plant carbon balance[J].Environmental and Experimental Botany，213：105434.

GEBAUER R，VOLARIK D，HOUSKOVAK， etal.， 2023.Sensitivity of physiological traits to different short-term drought events and subsequent recovery at the sapling stage inEuropean whiteelm （UlmuslaevisPall）[J]. Environmental and Experimental Botany，214：105469.

HAN H，LANG MH，ZHANG YW，et al.，2023.Preliminary study on meteorological factor change of Zhang Gutai from

1954 to 2O21[J].Protection Forest Science and Technology （1）：39-41，68.[韓輝，郎明翰，張亞偉，等，2023.章 古臺(tái) 1954—2021年氣象因子變化初步研究［J]．防護(hù)林 科技（1）：39-41，68.] HAN YG，LEI ZY，ZHAO GJ，et al.，2018.Canopy-DBH models for sand-fixing plantation of Pinus sylvestris var. mongolica[J]. Arid Zone Research，35（5）：1129-

1137.［韓艷剛，雷澤勇，趙國軍，等，2018.樟子松人工 固沙林冠幅一 —胸徑模型［J].干旱區(qū)研究，35（5）：

1129-1137.] HELD M，GANTHALER A，LINTUNEN A， et al.，2021. Tracheid and pit dimensions hardly vary in the xylem of Pinussylvestrisundercontrastinggrowing conditions [J].Frontiers in Plant Science，12：786593. LABRECQUE-FOY JP，ANGERS-BLONDIN S，ROPARS P， et al.，2O23.The use of basal area increment to preserve the multi-decadalclimaticsignalinshrubgrowthring chronologies：A case study of Betula glandulosa in a rapidly warming environment[J].Atmosphere，14（2）：319. LAI JS，ZOU Y， ZHANG S，et al.，2022. glmm. hp：an R package for computing individual effect of predictors in generalized linear mixed models [J].Journal ofPlant Ecology，15（6）：1302-1307. LEI ZY，HAN YG，ZHAO GJ，et al.，2018.Growth of Mongolian pine in Zhanggutai of Liaoning Province[J]. Journal of Zhejiang Aamp; F University，35（2）：324-

330.[雷澤勇，韓艷剛，趙國軍，等，2018.遼寧章古臺(tái)樟 子松生長過程分析［J]．浙江農(nóng)林大學(xué)學(xué)報(bào)，35（2）：

324-330.] LI MY，YANG LL，CAO Y，et al.，2024.Aging Mongolian pine plantations face high risks of drought-induced growth decline：Evidence from both individual tree and forest stand measurements[J].JournalofForestryResearch，

35（1） ： 38. LI MY，F(xiàn)ANG LD，DUAN CY，et al.，2O2O.Greater risk of hydraulic failure due to increased drought threatens pine plantations in Horqin Sandy Land of northern China [J].Forest Ecology and Management，461：117980. LI MY，DANG HZ，CHEN S，et al.，2O23. Sap flow variation in Pinus sylvestrisvar. mongolica plantation at the southern margin of Mu Us sandy land and its response to environmental factors [J]. Journal of Arid Land Resources and Environment，

37（4）：153-161.［李明陽，黨宏忠，陳帥，等，2023.毛烏 素沙地南緣樟子松樹干液流特征及其對環(huán)境因子的響應(yīng) [J].干旱區(qū)資源與環(huán)境，37（4）：153-161.] LI QY，CHEN RJ，LI ZW，et al.，2024.Eects ofPinus massoniana and Castanopsis hystrix mixture on the productivity and temporal stability of Pinus massoniana based on tree ring analysis[J].Chinese Journal of Ecology，43（12）：3605-

3614.［李瓊瑩，陳仁杰，李澤偉，等，2024.基于樹輪研究 馬尾松-紅錐同齡混交對馬尾松生產(chǎn)力及穩(wěn)定性的影響 [J].生態(tài)學(xué)雜志，43（12）：3605-3614.]

LI YJ，LI ZJ，SUN SJ，et al.，2025.Diferences in growth and wood anatomy in relation to meteorological variables between healthy and declining Mongolian pine trees [J].Chinese Journal of Ecology，44（6）：1837-1846.[李雅婧，李子靜， 孫守家，等，2025.健康與衰退樟子松生長和解剖特征差 異及其與氣象因子的關(guān)系［J].生態(tài)學(xué)雜志，44（6）： 1837-1846.]

LIU P，HU SY，WEI HX，et al.，2023.Response of radial growth of Pinus sylvestris var. mongolica of different stand ages to climate and extreme drought events in the semi-arid region of western Liaoning，Northeast China [J].Frontiers in Forests and Global Change，6：1272477

LIU X，ZHAO PS，GAO GL，et al.，2O22. Variations in phenology of Pinus sylvestris var. mongolica and the response to climate factors [J]. Journal of Desert Research，42（2）： 25-35.［劉軒，趙珮杉，高廣磊，等，2022.沙地樟子松 （Pinus sylvestris var.mongolica）物候特征及其對氣候的響 應(yīng)［J].中國沙漠，42（2）：25-35.]

LIU YJ，ZHOUL，ZHANGB，etal.，2O21.Radial variation of Cunninghamia lanceolata in different aged forests and its response to meteorological factors[J].Journal of Nanjing Forestry University（Natural Sciences Edition），45（2）： 135-144.［劉亞靜，周來，張博，等，2021.不同林齡杉木 徑向變化及其對氣象因子的響應(yīng)［J].南京林業(yè)大學(xué)學(xué) 報(bào)（自然科學(xué)版），45（2）：135-144.]

LIU YY，WANG AY，AN YN， et al.， 2018. Hydraulics play an important role in causing low growth rate and dieback of aging Pinus sylvestris var. mongolica trees in plantations of Northeast China [J].Plant Cell and Environment，41（7）： 1500-1511.

REN P，ROSSI S，GRICAR J，et al.， 2015. Is precipitation a trigger for the onset of xylogenesis in Juniperus przewalski on the north-eastern Tibetan Plateau？[J].Annalsof Botany， 115（4） ： 629-639.

SONG LN， ZHU JJ， LI MC， et al.， 2016. Sources of water used by Pinus sylvestris var. mongolica trees based on stable isotope measurements in a semiarid sandy region of Northeast China[J].Agricultural Water Management，164：281-290.

SONG LN，LI MC，ZHU JJ，et al.，2017. Comparisons of radial growth and tree-ring cellulose 8¹³C for Pinus sylvestris var. mongolica in natural and plantation forests on sandy lands[J]. Journal of Forest Research，22（3）：160-168.

SONG LN，ZHU JJ，YAN QL，20O9.Review on the shelter forest decline [J]. Chinese Journal of Ecology，28（9）： 1684-1690.［宋立寧，朱教君，閆巧玲，2009.防護(hù)林衰 退研究進(jìn)展［J].生態(tài)學(xué)雜志，28（9）：1684-1690.]

SONG LN，ZHU JJ，ZHANG JX，et al.，2019.Divergent growth responses to warming and drying climates between native and non-native tree species in Northeast China [J].Trees-Structure and Function，33（4）：1143-1155.

SONG LN，ZHU JJ，ZHENG X，2017. Forestation and management scheme ofPinussylvestrisvar. mongolica plantations in sandy lands based on their decline mechanisms [J].Chinese Journal of Ecology，36（11）：3249-3256.[宋立 寧，朱教君，鄭曉，2017.基于沙地樟子松人工林衰退機(jī)制 的營林方案［J].生態(tài)學(xué)雜志，36（11）：3249-3256.]

SUL，HUANG DH，ZHOU LL，2O24.Differences in sediment provenance from rainfall and snowmelt erosion in the Mollisol region of Northeast China [J]. Earth Surface Processes and Landforms，49（8） ： 2442-2457.

SU JR，XIAO SC， PENG XM， et al.，2024. Responses of radial growth to climate change for two dominant artificial coniferous trees [J]. Dendrochronologia，83：126163.

SUN SJ， ZHANG JS， ZHOU J，et al.，2021. Long-term effects of climate and competition on radial growth，recovery，and resistance in Mongolianpines [J].Frontiers in Plant Science，12： 729935.

TONG YQ，2023. Estimation of the transpiration based on the radial patern of sap flow and its response to meteorological factors in sandy Pinus sylvestrisvar. mongolica plantations [D].Shenyang：Shenyang Agricultural University.[童郁 強(qiáng)，2023.基于樹干液流徑向變化的沙地樟子松蒸騰耗 水估算及影響因素研究［D].沈陽：沈陽農(nóng)業(yè)大學(xué).]

VAN DE POL M，BAILEY LD，MCLEAN N，et al.，2016. Identifying the best climatic predictors in ecology and evolution[J]. Methods in Ecology and Evolution，7（10）： 1246-1257.

WANG JX，ZHAO HC，SONG QM，et al.， 2022. Efects of stand factors on herbaceous species diversity in Pinus sylvestris var. mongolica plantations in Zhanggutai area， Liaoning Province[J].Protection Forest Science and Technology（4）：17-20.[王景璽，趙浩唱，宋啟民，等， 2022.林分因子對遼寧章古臺(tái)樟子松人工林草本植物多 樣性的影響［J].防護(hù)林科技（4）：17-20.]

WEN S，SHI ZJ，ZHANGX，et al.，2O23.Effect of climate and competition onradial growth ofPinussylvestrisvar. mongolica forest in Hulunbuir sandy land of Inner Mongolia， China[J].Plants-Basel，12（13）.

WU J，GAO XJ，GIORGI F，et al.， 2017. Changes of effctive temperature and cold/hot days in late decades over China based on a high resolution gridded observation dataset [J]. International Journal of Climatology，37（S1）： 788-800.

WU J，GAO XJ， 2013.A gridded daily observation dataset over China region and comparison with the other datasets [J].Chinese Journal of Geophysics，56（4）：1102- 1111.［吳佳，高學(xué)杰，2013.一套格點(diǎn)化的中國區(qū)域逐日 觀測資料及與其它資料的對比［J]．地球物理學(xué)報(bào)， 56（4）：1102-1111.]

XU Y，GAO XJ， SHEN Y， et al.， 2009. A daily temperature dataset over China and its application in validating a RCM simulation [J].Advances in Atmospheric Sciences，26（4）： 763-773 sensitivity processes of conifers radial growth to climate factors based on tree ring width variations [J].Global Ecology and Conservation，48：e02743. YU FY，ZHANG，JX，SUN YR，et al.，2022. Leaf δ¹³C of the mainafforestation tree species in Horgin Sandy Land [J].Forest Research，35（4）：179-187.[于豐源，張金 鑫，孫一榮，等，2022.科爾沁沙地主要造林樹種葉片 S¹³（ C比較研究［J].林業(yè)科學(xué)研究，35（4）：179-187.] ZHANG SL，GAO RH，GAO ML，et al.，2023. Prediction of the potential distribution pattern of Pinus sylvestris var. mongolica in China under climate change[J].Journal of ZhejiangAamp;FUniversity，40（3）：560-568.[張世林，高 潤紅，高明龍，等，2023.氣候變化背景下中國樟子松潛 在分布預(yù)測［J].浙江農(nóng)林大學(xué)學(xué)報(bào)，40（3）：560-568.] ZHANG X，WU MW，KWON S，et al.，2022.Radial growth responsesofMongolianpine（Pinussylvestrisvar. mongolica） plantations at different ages to climate and groundwater level changes[J].Acta Ecologica Sinica，42 （16）：6827-6837.［張曉，吳夢婉，KWONS，等，

2022.不同林齡樟子松人工林徑向生長對氣候及地下水 位變化的響應(yīng)［J].生態(tài)學(xué)報(bào)，42（16）：6827-6837.] ZHAO Y，CAI LX，JIN YT，et al.，2021.Warming-drying climate intensifies the restriction of moisture on radial growth ofPinus tabuliformisplantationinsemi-arid areaofNortheast China[J]. Chinese Journal of Applied Ecology，32（10） ：

3459-3467.［趙瑩，蔡立新，靳雨婷，等，2021.暖干化加 劇東北半干旱地區(qū)油松人工林徑向生長的水分限制 [J].應(yīng)用生態(tài)學(xué)報(bào)，32（10）：3459-3467.] ZHOU FY，WANG SN，ZHANG F，et al.，2025. Response of radial growth of Pinus sylvestris var.mongolica to climate change at different ages on sandy land［J]. Chinese Journal of Ecology，44（6）：1847-1856.[周鳳艷，王賽男，張帆， 等，2025.不同林齡沙地樟子松徑向生長對氣候變化的 響應(yīng)［J].生態(tài)學(xué)雜志，44（6）：1847-1856.] ZHOU YP，LEI ZY， ZHAO GJ，et al.，2018.Impacts of plant diseases and insect pests on tree height growth of Pinus sylvestris var.mongolica plantations on sandy land[J].Arid ZoneResearch，35（3）：653-661.［周晏平，雷澤勇，趙國 軍，等，2018.病蟲害對沙地樟子松樹高生長的影響 [J].干旱區(qū)研究，35（3）：653-661.] ZHU JJ， ZHENG X， 2019. The prospects of development of the Three-North Afforestation Program（TNAP）： On the basis of the results of the 4O-year construction general assessment of the TNAP[J]. Chinese Journal of Ecology，38（5）：1600-

1610.［朱教君，鄭曉，2019.關(guān)于三北防護(hù)林體系建設(shè)的 思考與展望——基于40年建設(shè)綜合評估結(jié)果［J].生態(tài) 學(xué)雜志，38（5）：1600-1610.]

（責(zé)任編輯周翠鳴）