Skip to main content

Growth performance of planted population of Pinus roxburghii in central Nepal



Climate change has altered the various ecosystem processes including forest ecosystem in Himalayan region. Although the high mountain natural forests including treelines in the Himalayan region are mainly reported to be temperature sensitive, the temperature-related water stress in an important growth-limiting factor for middle elevation mountains. And there are very few evidences on growth performance of planted forest in changing climate in the Himalayan region. A dendrochronological study was carried out to verify and record the impact of warming temperature tree growth by using the tree cores of Pinus roxburghii from Batase village of Dhulikhel in Central Nepal with sub-tropical climatic zone. For this total, 29 tree cores from 25 trees of P. roxburghii were measured and analyzed.


A 44-year long tree ring width chronology was constructed from the cores. The result showed that the radial growth of P. roxburghii was positively correlated with pre-monsoon (April) rainfall, although the correlation was not significant and negatively correlated with summer rainfall. The strongest negative correlation was found between radial growth and rainfall of June followed by the rainfall of January. Also, the radial growth showed significant positive correlation with that previous year August mean temperature and maximum temperature, and significant negative correlation between radial growth and maximum temperature (Tmax) of May and of spring season (March-May), indicating moisture as the key factor for radial growth. Despite the overall positive trend in the basal area increment (BAI), we have found the abrupt decline between 1995 and 2005 AD.


The results indicated that chir pine planted population was moisture sensitive, and the negative impact of higher temperature during early growth season (March-May) was clearly seen on the radial growth. We emphasize that the forest would experience further moisture stress if the trend of warming temperatures continues. The unusual decreasing BAI trend might be associated with forest management processes including resin collection and other disturbances. Our results showed that the planted pine forest stand is sub-healthy due to major human intervention at times. Further exploration of growth climate response from different climatic zones and management regimes is important to improve our understanding on the growth performance of mid-hill pine forests in Nepal.


Earth’s climate has never been static and has shown great variability since its origin, the recent change, however, is accelerated by greenhouse effect causing abrupt temperature rise and unpredictable patterns of precipitation (Houghton 2004). Nepal, a Himalayan country, is highly threatened by impact of climate change. The atmospheric temperature in Nepal has been increasing at a rate of 0.04 to 0.06 °C/year, with a higher rate than the global average (Shrestha and Aryal 2011; Shrestha et al. 1999; Shrestha et al. 2000; Shrestha et al. 2017). Environmental changes including rapidly warming temperature and uncertainty in rainfall pattern have widely affected forest ecosystem all over the world (IPCC 2014). These climatic factors cause changes in growth rate of forest trees, which is directly linked with forest economics and ecosystem services provided by the forests. Temperature usually limits tree growth in poleward and high altitude treelines (Körner 2003; Holtmeier and Broll 2007; Harsch et al. 2009); hence, trees are expected to increase growth and shift poleward and upward with warming temperatures and increased CO2 uptake (Körner 2000; Borgaonkar et al. 2011; Gaire et al. 2017). However, tree growth in high mountains and dry regions could also be impaired by increases in temperature, mostly due to warming-induced water deficit (Cook et al. 2003, Sano et al. 2005, Lv and Zhang 2012, Dawadi et al. 2013, Liang et al. 2014, Qi et al. 2015, Gaire et al. 2017).

Nepal comprises nearly 45% of land mass (6.61 million ha) as forest including other woodlands (DFRS 2015). The government of Nepal has handed over about one-third of the total forest area to over 30,000 forest user groups, and the community forestry program in Nepal is considered as one of the highly successful avenues in nature conservation as a result there has been substantial increase in forest area in the middle and high mountains of Nepal (DFRS 2015). However, the technical aspect of forest management in Nepal is poorly understood due to less scientific studies on the growth performance of forest trees including forest health and harvesting modality of forest products in the context of changing climate and human interference.

Nepal is more prone to climate change due to its topographic structure and highly fragile ecosystem, and of course, the absence of long-term instrumental climatic data has become a major problem of studying spatiotemporal trends of climate change in forests ecosystem of Nepal (Cook et al. 2003). To overcome this situation, several alternative tools can be helpful for the precise estimation of past climate, analyze the current status, and speculate the future trends in environmental changes. Dendroclimatology offers one of such alternative tools for reconstruction of forest stand structure and can date the time at which tree rings were formed in many types of wood, to an exact calendar year, as well as in understanding climate dynamics and forest health (Chhetri and Thapa 2010, Cook and Kairiukstis 1990, Fritts 1976, Gaire et al. 2014, Speer 2010, Thapa et al. 2014, Tiwari et al. 2017). Hence, the climatic information preserved in the annual ring of trees is a valuable resource for studying the environmental changes such as temperature and precipitation from the climate-sensitive localities, particularly prior to the interval covered by direct climatic measurements (Cook and Kairiiukstis 1990; Speer 2010).

Dendrochronological studies from Nepal have examined more than 20 tree species including the most common treeline trees Abies spectabilis and Betula utilis from the higher Himalayan region (Tiwari and Jha 2018). Most of these species are conifers, only a few studies are covered for broad-leaved trees like Alnus, Betula, Castanopsis, and Rhodondendron. Tree ring studies in Nepal Himalaya have been more restricted to the lower temperature zones including high mountain forest and subalpine treelines; however, some studies also have been carried out from subtropical region (Speer et al. 2017, Sigdel et al. 2018, Aryal et al. 2018).

Previously, chir pine was examined for its dendrochronological potential in only a few studies in Himalaya (Bhattacharya et al. 1992, Brown et al. 2011, Borgaonkar et al. 1999). However, there have been a number of studies for this species covering both high and middle mountains of Nepal Himalaya (Bhuju and Gaire 2012, Shrestha et al. 2014, Speer et al. 2017, Sigdel et al. 2018, Aryal et al. 2019). Although P. roxburghii is also one of the most important tree species for the local people for its multipurpose uses such as timber, resins, firewood, the growth climate relationship, and forest stand health in panted forest are unknown in Nepal. Studies indicated that chir pine natural forests are highly sensitive to moisture availability in early growing season in relatively drier regions of western Nepal (Sigdel et al. 2018, Aryal et al. 2019); hence, we could expect more moisture-related physiological stress to chir pine forests in mid-hills, as central Nepal is experiencing intensified spring drought (Panthi et al. 2017; Nepal et al. 2020). Generally, the local people in the mid-hills of Nepal blame chir pine for transforming the hills to drier condition with drying of spells and suppression of species turnover, and of course, the blame is more intense for planted population. However, we do not have strong evidence to support the blame and to address the controversial issue regarding the impact of chir pine for drying of the ecosystem. We sought to access the growth-limiting climatic factor for radial growth and trend in aboveground biomass production in the form of BAI (basal area increment) using tree ring width of Pinus roxburghii planted forest from mid-hills of central Nepal, so that we could access the ecological health of these forests.

Materials and methods

Study area

The study was carried out at the Batase community forest (27.5952° N, 85.5477° E) of Kavrepalanchowk district in central Nepal (Fig. 1). The study area is a planted monoculture chir pine community forestin, the sub-tropical lower temperate zone in central Nepal. It occupies an area of 61.50 hectares area extending from 1560 to 1627 m elevation, the hill was bare five decades ago, at which seedlings were planted through Australian funding projects in collaboration with Department of Forests, Government of Nepal (DFRS 2015). The forest is now dominated by Pinus roxburghii with very few individuals of Alnus nepalensis and Ageretina sp., the area looks dry with drying of streams and water resources (local people reports) with white loose soil type with erosion in steep slope. Almost from all pine trees, barks were extracted for resin collection, and frequent cut stumps were observed during field visits.

Fig. 1
figure 1

Map of study area in Nepal showing sampling site and weather stations nearby


Pinus roxburghii Sarg. (chir pine) is one of the most common conifers in the sub-tropical region of Nepal and is distributed in all aspects of Western Himalaya but is generally found in well-exposed southern slopes in Central and Eastern Nepal. It can grow reasonably well in almost all types of soil (Jackson 1994) and has been proved to be a successful pioneer even at most degraded sites due to its high survival rate. Chir pine occupies 8.54% of total forest cover, being the fifth most dense tree species in Nepal with 7.05% of total volume in forest, and has also been an important planted tree species linked with highly successful community forestry program in Nepal (DFRS 2015). It possesses a biomass of 9.09 tons/hectare in Nepal, which accounts for 5.09% of total woody biomass (Aryal et al. 2018).

Sample collection

The tree core samples were collected from relatively less disturbed forest stand of Batase with bigger pine trees (Fig. 2). We collected a total of 85 cores from 40 Pinus roxburghii randomly using increment borer, i.e., two cores and sometimes three cores from the selected trees. The trees were selected with a healthy crown and without visible fire scars. The cores were then immediately kept in the plastic straw pipe with proper labeling. In addition, the GPS location, DBH, and height of each individual trees were measured. The cores were air dried for a few days. The cores were brought to the Dendrolab of NAST (National Academy of Science and Technology) Khumaltar, Kathmandu, for lab analysis. The air-dried cores were mounted in a wooden frame with transverse surface facing up (Speer 2010). The cores then were sanded with sanding papers of grids ranging from 120-400 until the ring boundaries were visible under binocular microscope. The cores with visible annual ring were dated to the calendar year.

Fig. 2
figure 2

Pinus roxburghii forest stand at study area (Batase Sindhupalchok Nepal) (a), individual tree with resin tapping mark (b)

Then, the tree ring measurement was done using a hardware called LINTAB attached having TSAP-Win program. Every single ring in each series was counted from bark toward the pith under the microscope adjusting the resolution for clear visualization. After the completion of ring width measurement, the individual tree-ring series were cross dated using the alignment technique, looking at the math graph and cross dating statistics as explained by Rinn (2003). The errors in the cross dating was rechecked and confirmed by using the computer program COFECHA (Grissino-Mayer 2001; Holmes 1983); standardization was carried out applying a computer program ARSTAN (Cook 1985) in order to remove age and other long-term variability that are considered as noise in dendrochronological study (Fritts 1976). The detrending of each sample was done using a negative exponential curve in order to remove the non-climatic age trends (Cook and Peters 1981). After detrending each series, ring-width chronology of both stands was developed by averaging the ensemble of detrended tree ring indices across the series for each year using arithmetic mean (Fritts 1976) which produces a mean value function that concentrates on signal and averages out the noise (Cook and Kairiiukstis 1990). Thus, standard, residual, and Arstan chronologies were developed. The residual chronology statistics of the overall period and standard chronology statistics of the common period were analyzed. Arstan chronology includes the residual chronology incorporated with the pooled autoregression. The major chronological statistics like, mean sensitivity, standard deviation, chronology length, mean index, serial correlation, and expresses population signal (EPS) were also computed in order to verify the dendroclimatic analysis. Growth-climate relationships were determined by correlating site standard chronology with monthly climatic variables (total precipitation, mean air temperature) from June of the previous growth year until October of the current growth year (Fritts 1976). Significance was tested using 1000 bootstrap replicates and with 95% confidence intervals in R with the R package bootRes (Zang and Biondi 2013; R Core Team 2019).

Basal area index (BAI) chronology

BAI is commonly used to assessing tree and stand growth as it allows precise quantification of tree productivity (Rubino and McMarthy 2000; Peñuelas et al. 2011). The unstandarized BAI sigmoidal growth model avoids detrending and standardizing employed in traditional calculation of ring-width indices (RWI) (Esper et al. 2002; Salzer et al. 2009). We employed individual tree BAI to produce mean unstandardized BAI series of all sampled trees for each year.

Ring widths were converted into tree BAI according to the following standard formula:

BAI = π (Rn2 − Rn − 12),

where “R” is the radius of the tree and “n” is the year of tree ring formation. The BAI chronology was produced using the bai.out function in the dplR package in R as some tree cores had missed pith and almost every core had intact bark.

Climatic data

The temperature and rainfall data of Panchkhal and Khopasi Meteorological station were taken from the Department of Hydrology and Meteorology (DHM), Kathmandu, respectively, as they are the nearest climate station from the study sites. Panchkhal station (27.6527° N, 85.6169° E) is situated in north-east direction (16 km), whereas Khopasi station (27.5663° N, 85.5566° E) in south-east direction (6 km) from study area. The rainfall data was recorded from Khopasi (1971-2016 AD) and the temperature record of Panchkhal (1978-2016 AD). These instrumental data of monthly mean, minimum, maximum temperature, and rainfall were used for growth climate response analysis. The paucity of long meteorological records is one of the major difficulties in undertaking dendroclimatic research in Nepal, which is associated with statistical calibrating the tree rings (Bhattacharyya et al. 1992; Shrestha et al. 1999, 2000; Aryal et al. 2018).


Climatic trend

The climate of this area is sub-tropical to temperate as indicated by the climate station data. Climatic diagram of Panchkhal and Khopasi is shown in Fig. 3, where average temperature (1978-2016) and total monthly rainfall (1971-2016) were taken into consideration, which showed that the mean maximum temperature in the Panchkhal area is 32.49 °C and the mean minimum temperature is 6.53 °C, with the highest rainfall in the month of July. There was a gradual increase in rainfall and temperature from May to July, and declining from August to December. However, there was an increasing trend in rainfall and temperature from January to May.

Fig. 3
figure 3

Monthly patterns of average temperature (Panchkhal) and monthly total rainfall (Khopasi) (DHM, Nepal)

The trends of annual average minimum, maximum, and mean temperature are shown in (Fig. 4). The annual average maximum and mean temperature showed increasing trend with the rate of 0.032 °C/year and 0.014 °C/year respectively, and the trend is significant only for maximum temperature. The rainfall showed non-significant declining trend.

Fig. 4
figure 4

Trend of annual temperature (maximum, mean, and minimum) and annual rainfall in climate stations nearby study site

Tree ring chronology

We produced a 44-year long tree ring width chronology from 25 trees (29 cores) of Pinus roxburghii planted population from Batase (Kavrepalanchowk) in the central Nepal (Fig. 5). The chronology has fulfilled all the statistical parameters used in standard dendrochronological studies such as mean sensitivity, series correlation, and expressed population signal (EPS) values (Table 1). The main difficulty of cross dating the tree ring width series was the most frequently occurring false rings, due to which many tree cores could not be cross dated. Similarly, the stand was relatively younger as the oldest tree was of 44 years; hence, it was obvious that the ring growth showed more fluctuation through time. The chronology statistics are presented in Table 1.

Fig. 5
figure 5

Tree ring-width chronology (standard and residual) of Pinus roxburghii (using ARSTAN), the left vertical panel showed ring-width indices and the right vertical panel showed number of tree cores measured. Both standard and residual chronologies were smoothed for a 5-year period

Table 1 Chronology statistics

Growth-climate relationship

We found a negative correlation between summer rainfall and tree radial growth. The strongest negative correlation was found between radial growth and rainfall of June (r = −0.47) followed by the rainfall of January (r = −0.31) (Fig. 6). It was also seen that the rainfall of November also showed negative significant correlation with radial growth; however, this is not more influential for radial growth as to be the ending phase of growth of the given year. It was found that previous year August mean temperature and maximum temperature showed a significant positive correlation with radial growth (r = 0.33, 0.20 respectively). However, significant negative correlation was observed between radial growth and maximum temperature (Tmax) of May (r = −0.32) and of spring season (March-May) (r = −0.29) (Fig. 7).

Fig. 6
figure 6

Correlation coefficients between radial growth and total monthly rainfall of June-December in the previous year (pJ…pD) to June-October of the current year (J…O), spring season, summer, and annual climate, dashed horizontal lines indicate significant correlation at 95% confidence limit for a two-tailed test

Fig. 7
figure 7

Correlation coefficients between radial growth and total monthly maximum (Tmax), mean (Tmean), and minimum (Tmin) temperature of June-December in the previous year (pJ…pD) to June-October (J…O) of the current year spring, annual, and summer months, the dashed horizontal lines indicate significant correlation at 95% confidence level for a two-tailed test

Basal area increment (BAI)

The BAI pattern in the studied chir pine forests appeared abnormal (Fig. 8) due to the clear growth suppression phase during the years of 1990-2000 AD, and due to the absence of juvenile growth increment trend (Jiao et al. 2017; Tiwari et al. 2017) However, the growth release phase after 2000 AD showed that the poor forest health was recovered. In mature forest stands, age-related trends of BAI are generally positive. BAI may continue to increase in healthy stands (LeBlanc et al. 1992; Duchesne et al. 2003), or stabilize (LeBlanc et al. 1992), but it does not show a decreasing trend until trees begin to senesce or unless trees are subject to significant growth stress (Weiner and Thomas 2001; Duchesne et al. 2003; Jump and Hunt 2006).

Fig. 8
figure 8

BAI (basal area increment) through time, the smoothing line operated for a 5-year period


Climatic trend

The warming trend in the study area is associated with higher rate of increase in maximum temperature (Tmax), and the trend is typical for central Himalaya (DHM 2017, Shrestha et al. 2017, Karki et al. 2019). This temperature trend is associated with higher rate of evapotranspiration to increase moisture stress for the forest trees (Dawadi et al. 2013; Liang et al. 2014; Tiwari et al. 2017), given the uncertainty of future precipitation patterns as indicated by IPCC (2014). The climatic trend of the study area representing the sub-tropical region of central Himalaya showed the great similarity with climate of whole Nepal despite some regional differences.

The entire central Himalayan zone has shown remarkable warming of maximum temperature although the rate of warming is site specific (DHM 2017, Shrestha et al. 2017). But the minimum temperature and rainfall were decreasing (Fig. 4) with the rate of 0.016 °C/year and 0.54 mm/year respectively, although the trend is not significant. The warming trend of the study area is a little low in comparison to whole Nepal (0.04 °C/year) for the period of 1980-2012 (Karki et al. 2019) and the global-surface global temperature trend (0.85 °C/year) for the period of 1880-2012 as mentioned by IPCC 2013. Similarly, the cooling trend of the study area is also lower than that of whole Nepal (0.02 °C/year) as mentioned by Karki et al. (2019), indicating that the area is getting warmer and drier with the decreasing trend of rainfall. The higher trend of day temperature (Tmax), more or less stable trend in annual rainfall, and the occurrence and severity of droughts especially in the winter and pre-monsoon seasons have contributed warmer and drier climate in central Himalaya (Panthi et al. 2017; Nepal et al. 2020). Consequently, the warming trend would affect tree growth patterns and trends in aboveground biomass production as it has already been observed in the Trans-Himalayan drier valley of Nepal (Tiwari et al. 2017).

Growth-climate relationship

We found positive correlation between radial growth and rainfall of March and April, although the correlation was not significant. The similar studies from western Nepal (Panchase) showed significant positive relationship between March, April rainfall, and radial growth (Aryal et al. 2018). Many studies in Nepal Himalaya showed spring season moisture sensitivity in tree growth (Dawadi et al. 2013; Tiwari et al. 2017; Aryal et al. 2018; Sigdel et al. 2018). Chir pine is a subtropical conifer most commonly found in dry, exposed, and well-drained southern aspects of sub-tropical hills in Nepal Himalaya. Irrespective of the strong moisture sensitivity of high mountain forest including treelines and timberlines in Nepal Himalaya (Dawadi et al. 2013; Liang et al. 2014; Gaire et al. 2017), pine radial growth in the study area showed weaker correlation with spring season (March-May) rainfall, and negative relationship with summer rainfall indicating that higher rainfall and lower day temperature during the early growing season (March-May) is preferred by chir pine forest.

The negative relationship between radial growth and temperatures is in accordance with various studies in Nepal Himalaya both from high mountain forests including treeline and timberline (Dawadi et al. 2013; Tiwari et al. 2017), and from sub-tropical regions at lower elevations (Aryal et al. 2018, Sigdel et al. 2018). However, irrespective of the negative response of maximum temperature (Tmax) on tree growth, the rainfall in the same period did not appear to have a significant positive effect. This could be due to less amount of rainfall to recharge water in the deep soil to be absorbed by the tall trees, failure of steep slope for retaining water, and a higher rate of evaporation during the season. Interestingly, the correlation strength between radial growth and spring season (March-May) rainfall is found decreasing from western to eastern region through mid-hills of Nepal showing that the climatic response to radial growth is site specific (Aryal et al. 2018, Sigdel et al. 2018). The low correlation strength in our study could be due to topographic structure with specific regional environmental conditions, and this corresponds to the fact that the total rainfall trend in Nepal is stronger in the eastern region that becomes weaker while going from eastern to western Nepal (Pokharel et al. 2019). We present that chir pine planted populations at Batase showed the similar response to temperature and rainfall at natural pine forests (Aryal et al. 2018, Sigdel et al. 2018). Although the strength of the correlation is low for spring season rainfall, more rainfall during the month of June showed negative impact on radial growth. Hence, we emphasize that the intensification of spring season drought in the future may have a negative impact on radial growth highlighting the moisture sensitivity of chir pine (Dawadi et al. 2013; Liang et al. 2014; Sigdel et al. 2018).


We found a sharp decline in BAI during 1994-2008, and the lack of juvenile growth release phase in chir pine trees. This may be due to the very slow growth for a few early years as the trees were planted. Further, our tree cores were extracted from the diameter of breast height; hence, we might have missed some early years’ wood formation. The relatively unhealthy BAI trend with decline of BAI over times indicated that the forest was not healthy, which could be associated with premature resin extraction and higher human disturbance pressure in the forest stand (Fig. 2). The unhealthy trend of BAI could be associated with younger age of forest stand (50 years >) and prolonged pre-monsoon drought (common in central Himalayan region), as many studies indicated the vulnerability of younger trees under worsening long-term drought (Vose and Swank 1994, Ruiz-Benito et al. 2015, Pompa-García and Hadad 2016). In comparison to the stronger competitive advantage of old trees for resources, the middle and young trees obtained less moisture and nutrition with lower trunk height, shallower root systems, and less trunk runoff, were more susceptible to long-term drought, and showed growth decline (Kloeppel et al. 1993, Pichler and Oberhuber 2007). However, our results indicated that the decline in BAI for the investigated chir pine forest did not show a strong impact of environmental factors (no remarkable drought events during the growth suppression years), rather it could be the results of human disturbance such as fire, grazing, and resin collection. Climate-induced growth decline has already been observed for drier parts of the Trans-Himalayan zone of Nepal, where tree growth was found to be positively correlated with spring season rainfall (March-May) (Tiwari et al. 2017). The negative impact of temperature-induced drought to pine forest growth has also been reported in different regions of Asia, Pinus taiwanensis, Pinus massoniana in the Dabie Mountains in subtropical China (Cai et al. 2020), Pinus koraiensis (Korean pine) of Changbai Mountain in northeast China (Yu et al. 2013) and Pinus latteri of Northeastern Thailand (Rakthai et al. 2020) all showed decreased growth at lower elevations under climate change characterized by warming and drought.

The BAI trend in the study area showed that the chir pine forest was not in good health and is showing poor performance in terms of carbon assimilation. And climatic factors did not show any strong impact on radial growth as indicated by our growth-climate relationship results. The sharp decline of BAI in times, and growth release afterward indicated that the growth is highly sensitive to human disturbance factors and resin extraction. The impact of forest disturbance in growth suppression was also indicated for P. roxburghii from similar climatic conditions in Kathmandu, Nepal (Speer et al. 2017).

Conclusions and management implications

Chir pine forest in Batase (Kavrepalanchowk) district region was 44-50 years old as evidenced from our ring count, and the stand age matches the plantation history of the forest. The radial growth of P. roxburghii showed significant negative correlation with both rainfall and temperature of different seasons. And the limiting factor for tree radial growth is not a single one. The BAI trend showed deviation from normal sigmoidal pattern (increasing) indicating that the forest is sub-healthy, although there was a substantial increment of BAI following the great decline of BAI after 1995. However, a sharp decline in BAI after 1995 until 2005 was very much clear, which could be the influence of resin extraction and other potential disturbance in the forests. The community users also explained the past trend of resin collection and higher disturbance in the forests. This study indicated the planted pine forests in the mid-hills of Nepal are much affected by human activities than changing climate. The magnitude of the impact of anthropogenic factors should be more closely observed for managing the forests in a sustainable basis. We suggest that forest managers should choose only one option as their desired product from pine forests. If they would wish to get more quantity of logs they should stop resin extraction, and if they would wish to get more resin they should not expect higher timber yield in the forests. Further replication of sites and taking cases of different climatic conditions would refine our knowledge on impact of modern climatic changes and human influence in the mid-hill forests. A detailed study on tree physiological performance and climate sensitivity of tree species is highly recommended before selecting the species for mass scale plantation.

Availability of data and materials

All data involved in this study are provided by the authors upon request.



Basal aArea iIncrement


Diameter at breast height


Department of Forest Resource and Survey


Department of Hydrology and Meteorology


Expresses pPopulation sSignal


Global Positioning System


Intergovernmental Panel on Climate Change






Nepal Academy to Science and Technology




Mmaximum temperature


Mean temperature


Minimum temperature


  • Aryal S, Bhuju DR, Kharal DK, Gaire NP, Dyola N. Climatic upshot using growth patternof Pinus roxburghii from western Nepal. Pak J Bot. 2018;50(2):579–88.

    Google Scholar 

  • Bhattacharyya A, La Marche VC, Hughes MK. Tree-ring chronologies from Nepal. Tree Ring Bulletin. 1992;52:59–66.

    Google Scholar 

  • Bhuju DR, Gaire NP. Plantation history and growth of old Pine stands in Kathmandu valley: a dendrochronological approach. FUUAST J. Biol. 2012;2(2):13–7.

    Google Scholar 

  • Borgaonkar HP, Pant GB, Kumar KR. Tree-ring chronologies from Western Himalaya and their dendroclimatic potential. IAWA Journal. 1999;20(3):295–309.

    Article  Google Scholar 

  • Borgaonkar HP, Sikder AB, Ram S. High altitude forest sensitivity to the recent warming: a tree-ring analysis of conifers from Western Himalaya, India. Quat Int. 2011;236:158–66.

    Article  Google Scholar 

  • Brown PM, Bhattacharyya A, Shah SK. Potential for developing fire histories in Chir Pine (Pinus roxburghii) forests in the Himalayan Foothills. Tree-ring Res. 2011;67(1):57–62.

    Article  Google Scholar 

  • Cai Q, Liu Y, Qian H, Liu R. Inverse effects of recent warming on trees growing at the low and high altitudes of the Dabie Mountains, subtropical China. Dendrochronologia. 2020;59:125649.

    Article  Google Scholar 

  • Chhetri PK, Thapa S. Tree ring and climate change in Langtang National Park, central Nepal. Our Nature. 2010;8(1):139–43.

    Article  Google Scholar 

  • Cook ER. A time-series analysis approach to tree-ring standardization. Ph.D. Dissertation. The University of Arizona Press, Tucson; 1985.

  • Cook ER, Kairiiukstis LA. Methods of dendrochronology: Applications in the Environmental Sciences. Dordrecht, The Netherlands: Kluwer Academic Publisher and International Institute for Applied System Analysis; 1990.

    Book  Google Scholar 

  • Cook ER, Krusic PJ, Jones PD. Dendroclimatic signals in long tree - ring chronologies from the Himalayas of Nepal. Int J Climatol. 2003;23(7):707–32.

    Article  Google Scholar 

  • Cook ER, Peters K. The smoothing spline: a new approach to standardizing forest interior tree-ring width series for dendroclimatic studies. Tree-Ring Bulletin. 1981;41:45–53.

    Google Scholar 

  • Core Team R. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2019.

    Google Scholar 

  • Dawadi B, Liang E, Tian L, Devkota LP, Yao T. Premoonsoon precipitation signal in tree ring of timberline Betula utilis in Central Himalayas. Quat Int. 2013;283:72–7.

    Article  Google Scholar 

  • DFRS. State of Nepal’s Forests. Department of Forest Resource and Survey.: Kathmandu; 2015.

    Google Scholar 

  • DHM. Observed Climate Trend Analysis in the Districts and Physiographic Regions of Nepal (1971–2014). Department of Hydrology and Meteorology, Kathmandu; 2017.

  • Duchesne L, Ouimet R, Morneau C. Assessment of sugar maple health based on basal area growth pattern. Can J For Res. 2003,33:2074–2080.

  • Esper J, Cook E, Schweingruber F. Low-frequency signals in long tree-ring chronologies for reconstructing past temperature variability. Science. 2002;295:2250–3.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  • Fritts HC. Tree rings and climate. London: Academic Press; 1976.

    Google Scholar 

  • Gaire NP, Koirala M, Bhuju DR, Borgaonkar HP. Treeline dynamics with climate change at the central Nepal Himalaya. Clim Past. 2014;10(4):1277–90.

    Article  Google Scholar 

  • Gaire, NP, Koirala, M, Bhuju, DR. Site- and species-specific treeline responses to climatic variability in eastern Nepal Himalaya. Dendrochronologia. 2017;41:44–56.

  • Grissino-Mayer HD. Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHA. Tree-Ring Res. 2001;57(2):17.

    Google Scholar 

  • Harsch MA, Hulme PE, McGlone MS, Duncan RP. Are treelines advancing? A global meta-analysis of treeline response to climate warming. Ecol Lett. 2009;12:1040–9.

    Article  PubMed  Google Scholar 

  • Holmes RL. Computer-assisted quality control in tree-ring dating and measurement. Tree -Ring Bulletin. 1983;43:69–78.

    Google Scholar 

  • Holtmeier FK, Broll G. Treeline advance - driving processes and adverse factors. Landscape Online. 2007;1:1–33.

    Article  Google Scholar 

  • Houghton J. Global warming: the complete briefing. Cambridge: Cambridge University Press; 2004.

    Book  Google Scholar 

  • IPCC. Climate change 2013: the physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. New York: Cambridge University Press; 2013. p. 1535.

    Book  Google Scholar 

  • IPCC. Summary for policymakers. In Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects.Contribution of Working Group II to Fifth Assessment Report of the Intergovernmental Panel on Climate Change (eds) FieldCB, BarrosVR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissal ES, Levy AN, MacCracken S, Mastrandrea PR, White LL. Cambridge University Press, Cambridge, and New York. 2014: 1-32.

  • Jackson JK. Manual of afforestation in Kathmandu, Nepal: Forest Research and Survey; 1994. p. 2.

    Google Scholar 

  • Jiao L, Jiang Y, Wang M, Zhang W, Zhang Y. Age-effect radial growth responses of Picea schrenkiana to climate change in the eastern Tianshan Mountains, Northwest China. Forests. 2017;8:294.

    Article  Google Scholar 

  • Jump AS, Hunt JM, Pen̈uelas J. Rapid climate change-related growth decline at the southern range edge of Fagus sylvatica. Glob Chang Biol. 2006, 12:2163–2174. doi:

  • Karki R, Hasson S, Gerlitz L, Talchabhadel R, Schickhoff U, Scholten T, Böhner J. Rising mean and extreme near-surface air temperature across Nepal. Int J Climatol. 2019:1–19.

  • Kloeppel BD, Abrams MD, Kubiske ML. Seasonal ecophysiology and leaf morphology of four successional Pennsylvania barrens species in open versus understory environments. Can J For Res. 1993;23:181–9.

    Article  Google Scholar 

  • Körner C. Biosphere response to CO2 enrichment. Ecol Appl. 2000;10:1590–619.

    Google Scholar 

  • Körner C. Alpine plant life: functional plant ecology of high mountain ecosystems. Berlin: Springer; 2003.

    Book  Google Scholar 

  • LeBlanc DC, Nicholas NS, Zedaker SM. Prevalence of individual-tree growth decline in red spruce populations of the southern Appalachian Mountains. Can J For Res. 1992;22:905–14.

    Article  Google Scholar 

  • Liang E, Dawadi B, Pederson N, Eckstein D. Is the growth of birch at the upper timberline in the Himalayas limited by moisture or by temperature? Ecology. 2014;95:140307191613003.

    Article  Google Scholar 

  • Lv LX, Zhang QB. Asynchronous recruitment history of Abies spectabilis along an altitudinal gradient in the Mt. Everest region. Journal of Plant Ecology 2012,5:147–156.

  • Nepal S, Pradhananga S, Shrestha NK, Kralisch S, Shrestha J, Fink M. Space-time variability of soil moisture droughts in the Himalayan region. Hydrol Earth Syst Sci Discuss. 2020., in review.

  • Panthi S, Bräuning A, Zhou ZK, Fan ZX. Tree rings reveal recent intensified spring drought in the Central Himalaya, Nepal. Glob Planet Chang. 2017;157:26–34.

    Article  Google Scholar 

  • Peñuelas J, Canadell JG, Ogaya R. Increased water-use efficiency during the 20th century did not translate into enhanced tree growth. Glob Ecol Biogeogr. 2011;20:597–608.

    Article  Google Scholar 

  • Pichler P, Oberhuber W. Radial growth response of coniferous forest trees in an inner Alpine environment to heat-wave in 2003. For Ecol Manag. 2007;242:688–99.

    Article  Google Scholar 

  • Pokharel B, Wang S-Y S, Meyer J, Marahatta S, Nepal B, Chikamoto Y, Gillies R. The east-west division of changing precipitation in Nepal. Int J Clim. 2019.

  • Pompa-García M, Hadad MA. Sensitivity of pines in Mexico to temperature varies with age. Atmósfera. 2016;29:209–19.

    Article  Google Scholar 

  • Qi Z, Liu H, Wu X, Hao Q. Climate-driven speedup of alpine treeline forest growth in the Tianshan Mountains, Northwestern China. Glob Chang Biol. 2015;21:816–26.

    Article  PubMed  Google Scholar 

  • Rakthai S, Fu P-L, Fan ZX, Gaire NP, Pumijumnong N, Eiadthong W, Tangmitcharoen S. Increased drought sensitivity results in a declining tree growth of Pinus latteri in Northeastern Thailand. Forests. 2020;11:361.

    Article  Google Scholar 

  • Rinn F. TSAP-Win: Time Series Analysis and Presentation for Dendrochronology and Related Applications. Version 0.55 User reference. Heidelberg, Germany. 2003.

  • Rubino DL, McMarthy BC. Dendroclimatological analysis of white oak (Quercus alba L., Fagaceae) from an old-growth forest of southeastern Ohio, USA. J Torrey Bot Soc. 2000;127:240–50.

    Article  Google Scholar 

  • Ruiz-Benito P, Madrigal-González J, Young S, Mercatoris P, Cavin L, Huang T-J. Climatic stress during stand development alters the sign and magnitude of age-related growth responses in a subtropical mountain pine. 2015;10(5):e0126581.

  • Salzer MW, Hughes MK, Bunnb AG, Kipfmueller KF. Recent unprecedented tree-ring growth in bristlecone pine at the highest elevations and possible causes. PNAS. 2009;106:20348–53.

    CAS  PubMed  Article  Google Scholar 

  • Sano M, Furuta F, Kobayashi O, Sweda T. Temperature variations since the mid-18th century for western Nepal, as reconstructed from tree-ring width and density of Abies spectabilis. Dendrochronologia. 2005;23:83–92.

    Article  Google Scholar 

  • Shrestha AB, Aryal R. Climate change in Nepal and its impact on Himalayan glaciers. Environ Change. 2011;11(1):65–77.

    Article  Google Scholar 

  • Shrestha AB, Bajracharya SR, Sharma AR, Duo C, Kulkarni A. Observed trends and changes in daily temperature and precipitation extremes over the Koshi river basin 1975–2010. Int J Clim. 2017;37:1066–83.

    Article  Google Scholar 

  • Shrestha AB, Wake CP, Dibb JE, Mayewski PA. Precipitation fluctuations in the Nepal Himalaya and its vicinity and relationship with some large scale climatological parameters. Int J Clim. 2000;20:317–27.

    Article  Google Scholar 

  • Shrestha AB, Wake CP, Mayewsk PA, Dibb JE. Maximum temperature trends in the Himalaya and its vicinity: an analysis based on temperature records from Nepal for the period 1971-94. J Clim. 1999;12(9):12.

    Article  Google Scholar 

  • Shrestha KB, Hofgaard A, Vandvik V. Recent treeline dynamics are similar between dry and mesic areas of Nepal, central Himalaya. Journal of Plant Ecology 2014,8: 347–358.

    Article  Google Scholar 

  • Sigdel SR, Dawadi B, Camarero JJ, Liang E, Leavitt SW. Moisture-limited tree growth for a subtropical Himalayan conifer forest in Western Nepal. Forests. 2018;9(6):340.

    Article  Google Scholar 

  • Speer JH. Fundamentals of tree-ring research. Tucson: The University of Arizona Press; 2010.

    Google Scholar 

  • Speer JH, et al. Pinus roxburghii stand dynamics at a heavily impacted site in Nepal: research through an educational fieldweek. Dendrochronologia. 2017;2017(1):2–9.

    Article  Google Scholar 

  • Thapa UK, Shah SK, Gaire NP, Bhuju DR. Spring temperatures in the far-western Nepal Himalaya since AD 1640 reconstructed from Picea smithiana treering widths. Clim Dyn. 2014.

  • Tiwari A, Fan ZX, Jump AS, Zhou ZK. Warming induced growth decline of Himalayan birch at its lower range edge in a semi-arid region of Trans-Himalaya, Central Nepal. Plant Ecol. 2017;218(5):621–33.

    Article  Google Scholar 

  • Tiwari A, Jha PK. An overview of treeline response to environmental changes in Nepal Himalaya. Tropical Ecology 2018,59(2): 273–285.

    Google Scholar 

  • Vose JM, Swank WT. Effects of long-term drought on the hydrology and growth of a white pine plantation in the southern Appalachians, For. Ecol. Manag. 1994,64:(1)25–39.

  • Weiner J, Thomas SC. The nature of tree growth and the age-related decline in forest productivity. Oikos. 2001;94:374–6.

    Article  Google Scholar 

  • Yu D, Liu J, Lewis BJ, Li Z, Wangming Z, Xiangmin F, Yawei W, Shengwei J, Limin D. Spatial variation and temporal instability in the climate–growth relationship of Korean pine in the Changbai Mountain region of Northeast China. For Ecol Manag. 2013;300:96–105.

    Article  Google Scholar 

  • Zang C, Biondi F. Dendroclimatic calibration in R: the bootRes package for response and correlation function analysis. Dendrochronologia. 2013;31:68–74.

    Article  Google Scholar 

Download references


We would like to acknowledge the Ministry of Forests and Environment, Government of Nepal for providing financial support, and members of the Batase Community Forest Users Group of Dhulikhel Central Nepal for sample collection permission for this research work. We are grateful to Nepal Academy of Science and Technology (NAST) Nepal for providing lab work permission in the Dendrolab. Further, Saroj Basnet of NAST and Mr. Basu Dev Paudel are acknowledged for their kind cooperation on sample measurement and preparation of study area map.


This research was supported by Ministry of Forests and Environment, Government of Nepal.

Author information




AT and PR conceptualized and designed the research, NT collected sample and performed the laboratory measurement. AT, SA1, and SA2 performed data analysis. The authors read and approved the final manuscript.

Authors’ information

Achyut Tiwari, PhD, is an Assistant Professor at the Central Department of Botany, Tribhuvan University, Kathmandu, Nepal, with research interests in forest ecology (dendroecology). Neeta Thapa is a MSc student at Central Department of Botany, Tribhuvan University; Sugam Aryal is a MSc student at Institute of Geography, University of Erlangen Nürnberg, Dr. Prabina Rana is a senior scientific officer at Nepal Academy of Science and Technology Nepal and Mr. Shankar Adhikari is an undersecretary at Department of Forests, Government of Nepal.

Corresponding author

Correspondence to Achyut Tiwari.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable.

Competing interests

All authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tiwari, A., Thapa, N., Aryal, S. et al. Growth performance of planted population of Pinus roxburghii in central Nepal. j ecology environ 44, 31 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Radial growth
  • Tree-rings
  • Basal area increment (BAI)