Assessing cross-datable distinct annual growth rings in non-native Pinus kesiya Royle ex Gordon in Zambia

The response of non-native forest plantation trees to climate change remains poorly understood.We hypothesized that precipitation and temperature modulate tree-ring width chronology at each site and that higher tree growth is exhibited at remote sites than sites near copper mines. This study investigates if the annual tree-ring boundaries in non-native Pinus kesiya Royle ex Gordon in Zambia are distinct, cross-datable, and coherent with climate signal. We collected increment cores from live trees and climate data near and further away from emission sources and developed site tree-ring width chronologies. Based on cross-dating and chronology building statistics (i.e., ESP > 0.85; Glk > 0.6 and series inter-correlation > 0.4), P. kesiya posses cross-datable distinct annual growth ring boundaries that exhibited a high climate signal at both sites. The tree-ring width chronology was positively modulated by precipitation and negatively by solar radiation and temperature. The dry season precipitation was the limiting factor for the growth of P. kesiya. The predicted decrease in dry season precipitation and increase in temperature and solar radiation may reduce tree growth of P. kesiya, reduce productivity, and extend the rotation age. The mean ring width in P. kesiya was not significantly (p = 0.296) different between sites. However, the mean basal area increment at the site near the emission source (Ichimpe) was significantly (p < 0.001) higher than at the remote site (Chati), suggesting site-specific influences that require investigation. We recommend evaluating the causes and consequences of tree growth variation between sites and their relation to environmental variation, including microclimate, soils, and pollution. In this regard, an assessment of site-specific ring-width chronology and tree growth variation in this study directly contributes to an improved understanding of non-native P. kesiya ecology, and it offers the potential to study trees' responses to edaphic and climatic factors. Knowing these responses deepens our understanding of non-native pine tree growth in the face of climate change, given the significant role of tropical forests in the global carbon cycle.     

Individual and time-varying tree-ring growth to climate sensitivity of Pinus tabuliformis Carr. and Sabina przewalskii Kom. in the eastern Qilian Mountains, China

Individual tree-ring width chronologies and mean chronologies from Pinus tabuliformis Carr. (Chinese pine) and Sabina przewalskii Kom. (Qilian juniper) tree cores were collected and analyzed from two sites in the eastern Qilian Mountains of China. The chronologies were used to analyze individual and time-varying tree-ring growth to climate sensitivity with monthly mean air temperature and total precipitation data for the period 1958–2008. Climate–growth relationships were assessed with correlation functions and their stationarity and consistency over time were measured using moving correlation analysis. Individuals’ growth–climate correlations suggested increased percentages of individuals are correlated with certain variables (e.g., current June temperature at the P. tabuliformis site; previous June, December and current May temperature and May precipitation at the S. przewalskii site). These same climatic variables also correspond to the mean chronology correlations. A decreased percentage of individuals correlated with these climatic variables indicates a reduced sensitivity of the mean chronology. Moving correlation analysis indicated a significant change over time in the sensitivity of trees to climatic variability. Our results suggested: (1) that individual tree analysis might be a worthwhile tool to improve the quality and reliability of the climate signal from tree-ring series for dendroclimatology research; and (2) time-dependent fluctuations of climate growth relationships should be taken into account when assessing the quality and reliability of reconstructed climate signals.     

Development of tree-ring width chronologies and tree-growth response to climate in the mountains surrounding the Issyk-Kul Lake,Central Asia

Three tree-ring width chronologies were developed from 75 Picea schrenkiana trees ranging from low- to high-elevation in the mountains surrounding the Issyk-Kul Lake, Northeast Kyrgyzstan. The reliable chronologies extend back to the mid-18th and late-19th centuries. Spatial correlation analysis indicates that the chronologies for the relatively high-elevation trees contain large-scale climatic signals, while the chronology at relatively low elevation may reflect the local climate variability. The results of the response of tree growth to climate show that these chronologies contain an annual precipitation signal. Furthermore, the influence of temperature indicates mainly moisture stress that is enhanced with rising elevation. The tree-ring records also captured a wetting trend in eastern Central Asia over the past decades. These new tree-ring width chronologies provide reliable proxies of precipitation variability in Central Asia and contribute to the International Tree-Ring Data Bank.     

Establishment of a 4650-year-long eigenvalue chronology based on tree-ring cores from Qilian junipers (Juniperus przewalskii Kom.) in Western China

This paper describes a tree-ring width chronology that spans the past 4650 years, established using the recently developed eigenanalysis technique. The aim is to show whether this eigenanalysis method enables the extraction of long-term tree-growth variations that are due to climatic changes, from a large dataset comprising 1263 tree-ring width records sampled from the highlands of Western China. In order to exclude the so-called growth rate/life span association effect, tree-ring width records were sorted into six subsets, based on the life spans of the trees sampled: 200–400, 400–600, 600–800, 800–1000, 1000–1500 years old, and trees older than 1500 years. Some partial chronologies were created, by pairing the set of all tree samples (living, dead, archaeological remains) with the living trees belonging to each subset. We computed the contours of tree-growth variations (on both 100-year and longer time scales) for each subset, ending with six pairs of these partial subset chronologies. Two sums of all these partial chronologies thus yielded a record of precipitation variations over a period ranging from 2627 BCE up to 2012 CE. It was found that this record shows a high degree of similarity to the existing chronology produced using the regional curve standardization (RCS) method applied to the same dataset, indicating that the eigenvalue chronology is capable of faithfully extracting long-term climatic variations. This also confirms that the first eigenvector represents the growth pattern that is characteristic of each biologically unique tree as well as the micro-environment of each tree stand. The variations observed over the last millennium seem to be connected to a cycle of solar activity with a period of ∼200 years. However, a clear lack of correspondence between solar activity and tree growth prior to 1000 CE indicates that the recent consistency may be coincidental. We believe that the eigenanalysis technique is readily applicable to other kinds of tree-ring datasets from different parts of the world.     

长白山北坡不同年龄红松年表及其对气候的响应

运用树木年轮气候学方法,研究了长白山北坡红松(Pinus koraiensis)不同年龄年表特征及其与气候因子间的关系,以期揭示年龄因素对年表的潜在影响。结果表明,平均年龄为63a的红松低龄年表与平均年龄为184a的高龄年表对气候的响应明显不同:低龄红松径向生长与当年1、2月月平均温度负相关(P<0.05),同时也受到上年及当年多个月份的月平均最高温度或最低温度的影响,但与降水的相关性未达到显著水平;高龄红松径向生长则与月平均温度间的关系不明显,而与当年1、2、4、6、7、9月的月平均最高温度正相关,与当年4月、9月的月平均最低温度负相关,同时受到上年5月及当年5月月总降水量的影响。因此,年龄因素对红松年表的气候响应方面存在一定影响,且高龄年表对气候响应的敏感性更高,包含有更多的气候信息。     

用西伯利亚落叶松年轮最大密度重建和布克赛尔5—8月份平均最高温度

利用采自和布克塞尔铁布克山二号沟的西伯利亚落叶松树轮样本,研制出树轮最大密度年表(MXD)和年轮宽度年表(TRW),分析了其年表特征和气候响应特点。结果表明,该样点的落叶松树轮最大密度年表与和布克塞尔气象站5-8月份平均温度和平均最高温度度均具有很好的正相关关系,最高单相关系数为0.660。用铁布克山二号沟的最大密度差值年表序列,可较好地重建和布克塞尔地区1795-2007年来春夏季平均最高温度距平,47 a重建值对实测值的解释方差达43.5%,且方程稳定。重建结果揭示,在和布克赛尔地区,20世纪平均最高温度距平要高于20世纪以前的平均最高温度距平,20世纪前中期平均最高温度距平出现了明显的上升,并且在重建的时段的末期,5-8月份平均最高温度距平表现出上升趋势。     

川西米亚罗林区不同海拔岷江冷杉生长对气候变化的响应

为分析青藏高原东缘半湿润区树木生长与气候关系随海拔变化的规律,于川西米亚罗林区分别在高、中、低海拔选取3个采样点,共采集132棵岷江冷杉年轮样芯,建立了最长达170a(1842-2011年)3个海拔高度的差值年表.岷江冷杉年轮指数与气候因子的相关分析表明:随海拔高度降低,温度与生长的负相关呈增加趋势.高海拔岷江冷杉径向生长与前—年冬季最低温呈显著正相关,中低海拔与当年春季均温、最高温和年均最高温呈显著负相关.低海拔岷江冷杉与当年4月降水呈显著正相关,随海拔升高降水与岷江冷杉生长的相关性降低.中低海拔岷江冷杉年表与4、5月帕尔默干旱指数(P DSI)呈显著正相关,表明在中低海拔存在春季干旱胁迫,抑制了岷江冷杉的生长.另外,大龄树木比小龄树木对气候变化的响应更敏感.     

昌岭山两个优势树种径向生长对气候变化的响应

在气候变暖背景下,树木径向生长对气候变化的响应存在不稳定性。利用采自祁连山东部余脉昌岭山两个优势树种油松和青海云杉的树轮样芯,建立树轮宽度标准年表,通过分析树轮宽度年表与气候要素的相关关系,探讨两个树种径向生长对气候变化的响应。结果表明：（1）油松年表比青海云杉年表包含更多的气候信息,其平均敏感度、标准差、信噪比和样本对总体的代表性等统计量均高于青海云杉标准年表。（2）气候要素对不同树种径向生长限制程度不同,油松径向生长主要与降水（前一年9月和当年3-8月）和气温（前一年9月）有关,但对降水的响应更为敏感,而青海云杉径向生长则受到气温（当年9月）和降水（前一年9月、当年3月和7月）的共同作用。（3）气温突变后,油松和青海云杉年表与各气温要素的相关性显著增强,而青海云杉年表与气温要素的相关性变化更明显,指示了青海云杉径向生长对气温的响应更不稳定。（4）生长季平均最低气温的升高诱导的干旱胁迫是油松和青海云杉树木径向生长-气温响应变化的主要原因。     

江西大岗山地区7-9月降水量的重建与分析

根据江西大岗山地区4个采样点的马尾松年轮样本,建立了本地区的综合年轮年表,分析了马尾松径向生长与气候因子变化的相关及响应关系,发现大岗山地区树木径向生长受生长季7—9月降水量影响较显著,呈负相关关系。在响应分析的基础上,首次重建了大岗山地区1892年以来7—9月的降水量,交叉检验表明重建序列是可靠的。重建结果表明,大岗山地区7—9月份降水变化在过去的117年中总体经历了3个偏干时期：1895—1902年,1908—1926年和1944—1985年,和3个偏湿阶段：1903—1907年,1927—1943年及1986—2008年。重建的降水量序列在1921年,1937年及1977年发生明显的方差突变,表明百年以来该地区降水量变化趋势存在阶段性变化。     

祁连山中部地区树轮宽度年表特征随海拔高度的变化

利用采自祁连山中部地区不同海拔高度的四个采样点的青海云杉树轮样芯 ,分别建立了树木年轮宽度年表。发现随海拔高度的上升 ,树轮宽度指数的振幅减小 ,年表的平均敏感性降低 ,样本间的一致性也逐步减小 ,上限年表与气候因子的相关性最低 ,这与目前大家普遍认同的上限树木的生长受温度控制的概念并不一致。进一步的分析表明 ,年表的敏感性随海拔高度降低主要是由于该区域树木生长的限制因子是春季降水 ,而降水随海拔高度的升高而增加 ,从而使得春季降水对树木生长的限制作用随海拔升高而逐步减弱 ;生物学指标的测定结果表明 ,生长在高海拔的树木对环境的生态适应策略发生变化 ,其生理代谢维持在较低水平 ,以避免环境变化带来的影响 ,因此生长在高海拔的青海云杉对环境变化的敏感性较差。     

Variability of climate–growth relationships along an elevation gradient in the Changbai Mountain, northeastern China

In order to explore climate–growth relationships at different elevations, tree-ring width chronologies of larch (Larix olgensis) were developed from three sampling sites on the northern slope of the Changbai Mountain, northeastern China. There were no consistent trends in statistical characteristics of the tree-ring chronologies along the elevation gradient, since trees in the forest interior had a complacent growth pattern. Monthly mean temperature and monthly total precipitation were used for the analysis. Correlation analysis indicated that temperatures in winter had negative correlations with tree growth (previous November, December and current March for the low-, mid- and high-elevation sites, respectively). The correlations between tree growth and June temperature varied from weakly negative at low elevations to significantly positive at high elevations. Precipitation in June of the growth year had negative relationship with the high-elevation chronology. However, high precipitation was associated with low temperature in early growing season, further supporting that temperature is a growth-limiting factor at high elevations. Our results suggest that along the elevation gradient, L. olgensis may respond in different ways to local climate change.     

芦芽山华北落叶松(Larix principis-rupprechtii)树轮宽度年表对气候因子的响应

在芦芽山地区采集3个不同海拔的华北落叶松(Larix principis-rupprechtii),在传统去趋势的基础上,采用"signal-free"方法对拟合曲线进行修正,避免了中等频率的气候信息引起的拟合偏差,最终建立3个不同海拔树轮宽度标准年表(STD)。同时以10a为界对上述年表进行滤波处理,得到3个低频年表。年表特征值表明,随着海拔升高,年轮平均轮宽变窄,敏感性和高频信息增强,低频信息减弱,这可能与逐渐恶劣的生境有关。中、低海拔年表的低频信息更一致,中、高海拔的高频信息更接近,而高、低海拔无论是标准年表还是高频、低频年表相似性均较差。树轮气候响应分析显示,低海拔STD年表与5月最低温负相关最为显著,STD和低频年表均与5、6月份土壤温度显著负相关,说明生境暖干,树木主要受生长季的干旱胁迫;中海拔STD年表与当年5月最高温正相关最为显著,STD和低频年表与土壤温度相关均不显著,说明生境逐渐变得冷湿,生长季的低温成为树木生长的限制因子;高海拔STD年表与气象要素相关不显著,低频年表与当年4月土壤温度正相关,说明高海拔最为冷湿,并有季节性冻土分布,生长季的土壤低温成为树木生长的限制因子。因此,全球变暖的趋势将更有利于高海拔树木的生长,而低海拔树木的干旱胁迫进一步加剧。     

Climatic response of three tree species growing at different elevations in the Lüliang Mountains of Northern China

In this study, we present the results of a dendroclimatological investigation of three coniferous tree species, Larix principis-rupprechtii, Picea meyeri and Pinus tabulaeformis, growing along an altitudinal gradient at the Lüliang Mountains in Northern China. Totally five tree-ring width chronologies were developed to explore the climate-growth responses of these tree species. No obviously regular trend associated with the increase of elevation was found by comparing the statistical characteristics of the chronologies. Correlation analysis indicated that the chronologies from lowerest to middle-high sites (SZ, BWD, BDGL and BDGP, respectively) were highly correlated, and different species from the same site showed the highest correlation. Growth–climate analysis indicated that the chronology of Larix principis-rupprechti at the uppermost site near the tree line (XWS) did not exhibit a significant response to the seasonal climatic factors, whereas the other four lower chronologies were consistently and significantly influenced by both the mean temperature from May to July and the total precipitation from March to June, regardless of tree species and elevation. The similarity of the tree growth–climate relationships of different species growing at different elevations (except that from the tree line) suggests that the trees in this region can provide common regional climate information, and combinations of multiple species (RC) are more successful in reconstructing the climate data than single species. The results of this research are very crucial for the future forest management and dendroclimatological sampling strategy in the arid to semi-arid area of northern China.     

Climate–growth relationships of the dominant tree species from semi-arid savanna woodland in Ethiopia

Long-term climate–growth relationships, were examined in tree rings of four co-occurring tree species from semi-arid Acacia savanna woodlands in Ethiopia. The main purpose of the study was to prove the presence of annual tree rings, evaluate the relationship between radial growth and climate parameters, and evaluate the association of El Niño and drought years in Ethiopia. The results showed that all species studied form distinct growth boundaries, though differences in distinctiveness were revealed among the species. Tree rings of the evergreen Balanites aegyptiaca were separated by vessels surrounding a thin parenchyma band and the growth boundary of the deciduous acacias was characterized by thin parenchyma bands. The mean annual diameter increment ranged from 3.6 to 5.0 mm. Acacia senegal and Acacia seyal showed more enhanced growth than Acacia tortilis and B. aegyptiaca. High positive correlations were found between the tree-ring width chronologies and precipitation data, and all species showed similar response to external climate forcing, which supports the formation of one tree-ring per year. Strong declines in tree-ring width correlated remarkably well with past El Niño Southern Oscillation (ENSO) events and drought/famine periods in Ethiopia. Spectral analysis of the master tree-ring chronology indicated occurrences of periodic drought events, which fall within the spectral peak equivalent to 2–8 years. Our results proved the strong linkage between tree-ring chronologies and climate, which sheds light on the potential of dendrochronological studies developing in Ethiopia. The outcome of this study has important implications for paleoclimatic reconstructions and in restoration of degraded lands.     

Recent growth increase in endemic Juglans boliviana from the tropical Andes

The spatial coverage of tree-ring chronologies in tropical South America is low compared to the extratropics, particularly in remote regions. Tree-ring dating from such tropical sites is limited by the generally weak temperature seasonality, complex coloration, and indistinct anatomical morphology in some tree species. As a result, there is a need to complement traditional methods of dendrochronology with innovative and independent approaches. Here, we supplement traditional tree-ring methods via the use of radiocarbon analyses to detect partial missing rings and/or false rings, and wood anatomical techniques to precisely delineate tree-ring boundaries. In so doing we present and confirm the annual periodicity of the first tree-ring width (TRW) chronology spanning from 1814 to 2017 for Juglans boliviana (‘nogal’), a tree species growing in a mid-elevation tropical moist forest in northern Bolivia. We collected 25 core samples and 4 cross-sections from living and recently harvested canopy-dominant trees, respectively. The sampled trees were growing in the Madidi National Park and had a mean age of 115 years old, with certain trees growing for over 200 years. Comparison of (residual and standard) TRW chronologies to monthly climate variables shows significant negative relationships to prior year May-August maximum temperatures (r = −0.54, p < 0.05) and positive relationships to dry season May-October precipitation (r = 0.60, p < 0.05) before the current year growing season. Additionally, the radial growth of Juglans boliviana shows a significant positive trend since 1979. Our findings describe a new and promising tree species for dendrochronology due to its longevity and highlight interdisciplinary techniques that can be used to expand the current tree-ring network in Bolivia and the greater South American tropics.     

Dendrochronology and dendroclimatology of yew in Poland

This study concerns dendrochronology and dendroclimatology of yew growing in Poland. The yew (Taxus baccata L.), a long-lived, slow-growing tree, is regarded as a species threatened with extinction. The eastern boundary of its range transects Poland. The analyses were performed on 34 yew populations which are protected as parts of nature reserves, as nature monuments, or which are planned to be protected. Samples were collected using Pressler borer from a total of 774 trees yielding 1307 profiles. Classical dating techniques (including cross-dating method) were subsequently used to reconstruct 34 local chronologies. Average tree-ring width was 0.84 mm and ranged from 0.27 to 1.47 mm. For most stands, the age of the yew trees is overestimated. The studied populations most frequently are 100–200 years old, and the oldest yew trees in Poland are growing in the KS population (age estimated at 500–600 years). Average tree-ring width was found to be strongly dependent on the chronology length/tree age. The 674 pointer years calculated for the local chronologies served as a basis for determining pointer years for the entire study area. Negative years include: 1862, 1865, 1917, 1927, 1940, 1947, 1956, 1963, 1969, 1976, 1979, 1993, 1996, and 2003. Positive years are: 1884, 1914, 1916, 1965, 1977, 1988, 1997, and 2007. Analysis of weather conditions in those years demonstrated a strong relationship between tree-ring width and air temperature in winter, pre-spring, and early spring. Higher-than-average temperatures during those seasons correlate positively with yew tree-ring width. Response function analysis performed for local and regional chronologies point to a dominant role of air temperature in December of the year preceding growth and in January, February, and March of the current year (linear relationships). June precipitation is an additional factor affecting tree-ring width in some areas of northern and northwestern Poland: the higher the precipitation, the wider the tree-rings. The results obtained, particularly those concerning growth-climate relationships and dendroclimatic regionalisation, can be used in the on-going programme of yew restitution in Poland.     

坡向对海拔梯度上祁连圆柏树木生长的影响

选择青海省同德县南部河北林场的一个连续坡面,根据不同海拔和坡向设置4个采样点,采集祁连圆柏(Sabina przewalskii)树轮数据,分析不同海拔和坡向对树木生长的影响。结果表明:坡面上部3个采样点的树轮年表特征值均呈一定的变化规律——平均敏感值(MS)和标准差(SD)随海拔升高而增大,一阶自相关(AC)随海拔升高而递减,下限年表特征值均表现出与其他3点的不同,都是最值(MS和SD均最大,AC最小);年表间相关和主成分分析结果都显示出海拔梯度上的变化规律,但下限差异显著;树轮指数与当年6–8月平均气温的相关系数呈增强趋势,森林上限受当年7、8月平均气温影响较大,下限树轮指数不仅与当年6月和前一年11月的气温显著负相关,而且受前一年8月和当年5月的月降水量影响显著。与通常情况"下限树木生长受降水制约"比较,这里的温度作用增强而降水限制减弱。显然,坡向扭转是海拔梯度上影响祁连圆柏生长变化的重要因子。     

Tree growth and its association with climate between individual tree-ring series at three mountain ranges in north central China

Individual tree-ring series may show changed growth trends and divergent climate–growth associations even within a site, highlighting the need to examine tree growth and its climate association before building a chronology. We provided a case study for the stratification and temporal variability of tree growth and its climate associations of individual cores for three mountain ranges in north central China. Tree growth is mainly limited by moisture conditions in previous July–September and current June–August. Repeated sampling and field investigations of Picea wilsonii at Xinglong Mountain over a growth year of 2004 suggested that the growing season is from about the end of April to the end of September. It appears that the moisture conditions in previous and current growing seasons are crucial for tree growth in this region. However, a decrease in drought limitation was observed for a few tree-ring series. We thereby built the pooled chronology and sub-site chronologies with only drought-sensitive tree rings similar climate–growth relationships from the three mountain slopes. Growth disturbances of tree-ring series are detected by checking the occurrence of successively low values of the biweight series, which are treated by fitting a flexible curve.     

Biases in RCS tree ring chronologies due to sampling heights of trees

The Regional Curve Standardization (RCS) is one of the most employed standardization methods to remove biological signals in long tree ring chronologies. The approach assumes that an overall age-related growth trend typify all tree ring series to be included in a standardized tree ring chronology. Although several potential problems of the method have been examined, the influence of varying the sampling height along tree stems has not been evaluated. Considering that age-related growth trends may vary with stem height, biases may arise when combining samples from unknown or variable sampling heights, a frequent situation with subfossil logs. In this study we perform a detailed stem analysis of 15 lakeshore black spruce (Picea mariana Mill. B.S.P.) trees in the taiga of eastern Canada to describe how the age-related growth trend varies with stem height and evaluate associated biases in RCS chronologies built from living and subfossil trees. Results show that the age-related growth trends vary markedly and systematically along stems, potentially generating large methodological biases in RCS chronologies, especially near the recent chronology end. These biases may lead to erroneous reconstructions of recent climatic trends and cause false divergence between tree ring and climate series. We have developed a correction procedure that appears efficient in removing these biases from chronologies built with the lakeshore trees and associated subfossil logs.     

An investigation of the snowpack signal in moisture-sensitive trees from the Southern Canadian Cordillera

Variations in mountain snowpack in the western Canadian Cordillera have widespread and important impacts on ecosystems, environmental processes and socio-economic activities (e.g. water availability downstream). Historical records of snowpack generally span only the latter half of the 20th century offering a limited perspective on the causes and uniqueness of recently observed changes across the region. This paper explores the potential utility of a network of low elevation Douglas-fir (Pseudotsuga menziesii) and ponderosa pine (Pinus ponderosa) tree ring-width chronologies to reconstruct past snowpack variations. Correlation coefficients between the tree-ring chronologies and a set of snow water equivalent (SWE) records are calculated and mapped. Separate analyses were carried out for total ring- width (TRW) and partial-ring measurements (earlywood and latewood; EW and LW). A set of Adjusted LW chronologies was also developed; in these, the relationship between LW and the preceding EW width has been removed. The ring-width chronologies exhibit moderately strong relationships with SWE records from the western Canadian Cordillera and these relationships vary in sign across the region. Distinctive regional groups are identified where chronologies exhibit same-sign correlations with SWE, in possible accordance with the elevation and characteristics of the tree-ring chronology sample sites. The EW chronologies correlate more strongly and consistently with SWE records in regions where the growth relationship with SWE is negative. The LW chronologies, and particularly the Adjusted LW chronologies, exhibit a greater number of positive correlations with the set of SWE records. Collectively these results offer valuable insights for developing a targeted sampling and/or reconstruction strategy that can exploit these different relationships with SWE to generate more robust estimates of pre-instrumental snowpack for the region.