Nechitaetal. 2018
Nechitaetal. 2018
net/publication/328205347
A 781-year oak tree-ring chronology for the Middle Ages archaeological dating
in Maramureș (Eastern Europe)
CITATIONS READS
0 112
4 authors:
Some of the authors of this publication are also working on these related projects:
Evaluation and analysis of forest ecosystems from Romania in the pan-European (Level I) and intensive (Level II) grids of forest monitoring View project
All content following this page was uploaded by Nechita Constantin on 15 December 2018.
Dendrochronologia
journal homepage: www.elsevier.com/locate/dendro
A R T I C LE I N FO A B S T R A C T
Keywords: From 1997 to the present, a sustained project in the Maramureș region in Romania was completed with the
Tree-ring chronology construction of a 781-year oak tree-ring chronology. A total of 395 samples from living trees and 429 from
Quercus archaeological wood were analysed with dendrochronological methods. The study aimed to provide the scien-
Sapwood tific community with a new oak chronology that could be applied in dendroarchaeology, dendroclimatology and
Mean segment length
dendroecology studies and also for interpreting past socio-economic events. Furthermore, we studied the number
Average growth rate
of sapwood rings and the growth pattern for different time periods. The chronology covered the continuous
period of A.D. 1236‒2016. From the 824 samples collected, we separated 271 tree-ring series into a new
chronology (A.D. 1406–2016), which fulfilled all the criteria necessary to reconstruct past climate and en-
vironmental changes. The resulting new tree-ring chronology indicated robust signal series intercorrelation
(r = 0.55) and an average mean sensitivity of 0.21. Based on an analysis of the sapwood, we recommend es-
timating a number between 8 and 32 rings to the date of final ring for 95% confidence range of mean on the for
felling in the Maramureș region. Additionally, we observed that the number of sapwood rings was not correlated
with altitude or oak tree species.
⁎
Corresponding author.
E-mail address: [email protected] (I. Popa).
https://doi.org/10.1016/j.dendro.2018.10.006
Received 28 May 2018; Received in revised form 5 October 2018; Accepted 8 October 2018
Available online 10 October 2018
1125-7865/ © 2018 Published by Elsevier GmbH.
C. Nechita et al. Dendrochronologia 52 (2018) 105–112
in this region can be highlighted using dendrochronological dating, countries and N-America, especially the sessile oak from the Rona re-
because information is lacking about the time of construction of many gion, which has high economic value (Horj, 2006).
of them. The material provided can also improve the existing chron-
ology by increasing the replication and therefore the quality of the 3. Material and methods
chronology.
In Maramureș, the constructing of wooden churches is an old tra- In this study, we collected samples from living trees and historical
dition. UNESCO recognizes eight of these churches (in total 43) as constructions. The dendrochronological material from living trees was
World Heritage Sites. The wooden churches in Bârsana, Budești Josani collected from old-growth oak (Quercus robur L.) and sessile oak
and Ieud were dated using an oak chronology constructed at the end of (Quercus petraea (Matt.) Liebl.) from eight sites (Table 1).
the 20th century (Eggertsson and Babos, 2003). The “Church on the From each site, 20 to 50 trees were cored randomly with a Pressler
Hill” in Ieud, is the oldest church, which is dated to 1364. The precise borer following a standard dendrochronological methodology (Cook
dating information established regarding the felling date of timber used and Kairiukstis, 1990; Speer, 2010). Ultimately, a database summar-
for the different parts of the buildings is possible only by the use of the izing 395 samples was obtained representing one core from each tree
dendrochronological dating method, which indicates the importance of (Fig. 1a).
this type of research. Two sites with both oak and sessile oak were investigated (Borlești
Part of this historical wood is sold to various regions of the world in and Coștui), to obtain additional information regarding interspecies
which the law requires that the restoration of monuments must be compatibility for the newly formed chronology, because in terms of
performed with wood from periods close to those of the building con- wood anatomy, the species could not be separated. Most of the sam-
struction. In Maramureș, much of the historical monuments require pling sites were chosen from the central depression of Maramureș, be-
restoration. By highlighting the importance of old wood through den- cause the building timber was most likely provided from nearby forests
drochronological techniques, we demonstrated to the regional autho- (Fig. 1b). Additionally, four sampling sites located on the eastern cliff of
rities and national decision makers the importance of providing funds the volcanic Mountains of Maramureș were included in the chronology
for the conservation and restoration of the wooden monuments. to increase the area of representativity (Nechita et al., 2017).
To extend the tree-ring chronology back in time, we also collected
2. Study area samples from historical buildings located in the region. The timber for
the extension of the Maramureș regional oak tree-ring chronology ori-
The study area is located in the northern part of Romania, i.e in the ginated from old historical buildings (houses, barns and mills), included
Romanian part of the historical Maramureș region, (Fig. 1b). The cli- actually in the patrimony of the Maramureș museum in Sighetul
mate is characterized by Atlantic influence with a mean annual pre- Marmaţiei (Fig. 1c). The sampled buildings represented the oldest and
cipitation of 770 mm/year and mean annual temperature of 5.5 °C well-preserved monuments in the region. The samples from wooden
based on average values extracted from CRU TS3.23 (0.5° x 0.5°) for the constructions were collected in 1997 and 2016. In most cases, samples
region (Harris et al., 2014). Predominant tree species are pedunculate were obtained by coring, but we also collected larger samples using a
oak (Quercus robur L.), sessile oak (Quercus petraea (Matt.) Liebl.), chain saw. At present, a market exists for old wood. Accordingly, we
European beech (Fagus sylvatica L), European hornbeam (Carpinus be- also collected wood from various private old wood stacks that were
tulus L) and sweet chestnut (Castanea sativa Mill.). From the first land open for purchase.
survey and mapping of the region it is known that the natural oak sites The wooden material was processed in the dendrochronological
have been heavily exploited for economical interest (1769–1773 laboratory at the National Research and Development Institute for
Franzisco-Josephinische Landesaufnahme) (Sabau, 2016). Oaks and Silviculture, Marin Drăcea, using standard methods (Baillie, 1982;
beech were clear‐cut in the mid -1800s to provide wood for the Austro- Speer, 2010). After preparation, the wooden materials were scanned at
Hungarian Empire, and later wood was sold to other European 3200 dpi (interpolated from 2400 true dpi) using an Epson Expression
106
C. Nechita et al. Dendrochronologia 52 (2018) 105–112
Table 1
Descriptive statistics of the living raw tree-ring width chronologies.
Site names Species Long. E Lat. N Series intercorrelation Average mean sensitivity No. trees Length (yrs.) RW (mm/yr.)
12000XL flatbed scanner and measured to the nearest 0.001 mm using over 20-year windows lagged by 10 years, using ARSTAN software
image analysis program Cybis™ CooRecorder and CDendro software (Cook and Krusic, 2005). Individual raw measurements were detrended
(Larsson and Larsson, 2017). CooRecorder 9.0 allowed us to import an using ARSTAN software by applying a 67% cubic smoothing spline with
image of the wood sample, identify the ring boundaries and measure a 50% cut-off frequency to remove effects of disturbance events and
the tree-ring parameters. stand competition. The autocorrelation was removed by applying a
Visual and statistical crossdating were both performed and checked filter to the detrended measurements. For removing biases caused by
using TSAP-Win (Rinn, 2010) and COFECHA (Holmes, 1983) software. extreme values, we used bi-weight robust mean to obtain a mean
The regional chronology was compiled from living trees and historical chronology.
wood. The chronology was constructed from oaks samples (combined
both sessile oak and pedunculate oak) with the best visual pattern-
matching and statistical verification results. First, we crossdated his- 4. Results
torical wood internally where multiple samples were available from the
same site/building, since for some cored buildings the date of con- 4.1. The Maramureș oak tree-ring chronology
struction/reconstruction is known. The chronology construction pro-
gresses step by step with caution, using dating checks such as: iterative In the Maramureș (MM) a set of 824 samples were studied, and the
replacement and removal of series to test the stability of the statistical absolutely dated samples resulted in a 781-year chronology covering
results, comparison with others regional tree-ring chronologies. For the time span 1236–2016 suitable for dendrochronological dating
dating quality control and teleconnection assessment, the graphical- (Fig. 1a, d). Summary chronology descriptive statistics were presented
chronology match was statistically verified using the common statistical in Table 2. The length of individual tree-ring series cored both dated
parameters: coefficient of coincidence Gleichläufigkeit (Glk), correla- and undated varied between 23 and 293 years with the mean length of
tion coefficient (ccoeff), t-value, t-value according to Hollstein, Baillie the individual tree-ring series of 125 years. The statistical parameters in
and Pilcher, and cross-date Index (CDI) (Knibbe, 2004; Rinn, 2010; several cases were unsatisfactory for the chronology to be utilized later
Baillie and Pilcher, 1973). as a proxy for climate reconstruction. Altogether, 553 samples re-
The Glk represents the percentage slope equivalence of a given presenting over 67% of the total data set could not be successfully used
sample and reference with the overlapping of both series (Eckstein and in the new tree-ring chronology, because of low correlation. These in-
Bauch, 1969). The correlation coefficient is a percentage of overall si- dividual tree-ring series had to be excluded from the final MM chron-
milarity between two series. When considering the length (N) of over- ology. Therefore, we eliminated individual samples with weak corre-
lapping series, the correlation coefficient becomes more valuable and is lation with the mean chronology (below 0.4).
incorporated in a t-test. We also used the t-values calculated after de- Based on these criteria, we extracted a subset of tree ring width
trending, as proposed by Hollstein (1980), to exclude any age trend that (TRW) data and we created a chronology from 271 wood samples, in-
might significantly affect the resulting statistics. Tree-ring series were cluding 134 living trees and 137 historical wood samples, covering the
graphically matched to one another using statistical parameters; at the 611-year period between 1406 and 2016, including only the samples
end, they were averaged to produce a mean sample chronology. CDI is a that could be strongly correlated (Fig. 2a, b). Several descriptive sta-
coefficient calculated in TSAPWin© that includes Glk and t-values cal- tistics commonly used in dendrochronology were calculated for subset
culated with different detrending options (Rinn, 1996). of tree-ring width data to evaluate the chronology strength. These in-
Haneca et al. (2005) proposed t-test scores ranging between 5.0 and cluded the intercorrelation of the mean series (SI = 0.55) which is a
7.5 resulting in an acceptable t-test threshold of 6.0 for matching tree- measure of the strength of the common signal to sampled trees at the
ring series to a reference chronology. However, Sass-Klassen et al. site.; mean sensitivity (MS), an indicator of the mean relative change
(2008) in their work for establishing new reference chronologies for the between consecutive ring widths showed values of 0.21.
dendroprovenancing of timber indicated the Glk must be at least 50% to The newly established regional chronology fulfils all the constraints
have any statistical significance, in the tree-ring series matching. In required for use in climate reconstruction for the Maramureș region
addition, an empirically tested 99.9% confidence threshold (tBP > 3.7) (Buras, 2017): EPS > 0.85, different age classes occurred, and the
can be also considered as a useful guidance in evaluation of the cross-
dating results and teleconnection pattern (Fowler and Bridge, 2017). Table 2
For our reference chronologies, several overlapping records were Full dataset chronology summary.
relatively short but the correlation was high, which explained why we Sub-sets No of series Start End Period MSL AGR Rbar EPS
carefully analysed the t-test and the other statistical parameters men-
tioned to justify good matches. At the final stage, absolutely dated in- Living 395 1773 2016 244 73 1.60 0.13 0.94
Historical 429 1236 1999 764 52 1.36 0.11 0.81
dividual tree-ring series were included into a regional chronology. To All 824 1236 2016 781 62 1.54 0.29 0.94
assess the quality of the chronology, expressed population signal (EPS)
(Wigley et al., 1984) and inter-series correlation (Rbar) were computed *mean segment length (MSL), average growth rate (AGR), inter-series correla-
tion (Rbar) and the Expressed Population Signal (EPS).
107
C. Nechita et al. Dendrochronologia 52 (2018) 105–112
correlation of living trees with climate was high and stable over time, 1987) and 95% confidence interval of mean falls between 8 and 32
according to Nechita et al. (2017). The Rbar value of the new chron- rings. The average value for all sites in the region was 14 sapwood
ology ranged from 0.16 (1566 y) to 0.36 (1816 y) (Fig. 2d). The signal rings, both for sessile and pedunculate oak. A slightly lower number of
to noise ratio (SNR) varied between 5.91 and 58.5. The chronology was sapwood rings (12 rings) was identified for sessile oak located at lower
replicated with at least 14 samples, with that lowest replication be- elevations when individual sites were analyzed. The total number of
tween 1780 and 1820 and the maximum replication in the 20th century rings counted for living trees from bark to pith varied between 86 and
of 135 samples. t-values calculated using the Baillie-Pilcher (tBP) and 244 rings, and the average site value varied between 118 and 220 rings
Hollstein (tH) algorithm highlight the correlation strength between (Fig. 3b).
chronologies represented by living and historical material in their
overlapping period (Fig. 2c).
4.3. Age related growth trends based on the MM dataset
4.2. Sapwood rings of oaks in NW Romania For the full dataset of living trees, age-related ring width trends
indicated a brief increase in ring width over time (Fig. 4). For the first
We analysed the average number of sapwood rings in the full da- 50 years, the trend of juvenile wood was almost the same for both the
taset of living tree samples. The results indicated that the number of living and the historical wood. Between 50 and 150 years from first tree
sapwood rings ranged from 6 to 40, with an average for each sampling ring was observed a linear trend of growing both in living and historical
site between 10 and 16 rings (Fig. 3a). The distribution of the number wood. The only difference in growth trend observed was an increasing
of sapwood rings conforms to a log-normal distribution (Hillam et al., trend after 150 years in living trees for a small period, followed by a
108
C. Nechita et al. Dendrochronologia 52 (2018) 105–112
Fig. 4. Square symbol: age versus average growth rate (green, historical wood; blue, living trees). Continuous line: mean average growth in report with year starting
from the first tree ring (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.).
downward trend. A similar upward trend was not observed in the his- were included in the new Maramureș (MM) oak TRW chronology sui-
torical material (Fig. 4). table as proxy. A total of 553 samples collected from both living and
historical wood samples were excluded from the proxy chronology.
4.4. Teleconnections of MM chronology These samples were excluded for multiple reasons (low number of tree
rings, low correlation with mean chronology). Some samples of his-
Comparison of MM chronology with regional chronology of SE torical tree-ring data, were excluded because of an insufficient number
Slovenia (Čufar et al., 2008a, b) (site #12) indicate a moderate corre- of tree rings. The suggestion of Cook and Kairiukstis (1990) is not to
lation value (tv = 5.5 / tBP = 4.6 / tH = 4.8). However, using an up- include samples with less than 40 rings per sample. Other studies
dated chronology of Slovenia (Čufar et al., 2014 unpublished), we can propose the minimum of even 50 tree rings per sample (Baillie, 1995;
observe higher t values (tv = 11.2 / tBP = 5.7 / tH = 6.4) (site #13) Haneca et al., 2009). Experience with precision dating of art objects
(Table 3.). attests to the fact that even time series with more than 100 tree rings
Teleconnection statistics with living tree ring chronologies from SE may not be sufficient for successful dating (Horovitz, 2012). In our case,
Europe highlight wide range of variability (Table 3). The area covering only 18 individual series contained less than 40 tree rings, and 35
the sites with tBP / tH > 6 is indicative of the region where newly MM contained less than 50 rings, which was less than 4% of the entire data
updated chronology has the highest potential to be used as a reference set. However, the abundant historical timber and subfossil wood of
chronology. Outside this core region (including sites #1, #4, #8, #9), different ages in the region offer plausible opportunities to extend the
the dating potential of MM is lower. While the areas with tBP < 3.7 chronology further back in time and to improve the replication.
(including sites #5, #6, #10, #11) are already outside the zone were Several historical samples with more than 200 tree rings had low
MM can be used as a reference chronology (Fig. 5). crossdating statistics. The explanation might be the provenance of
timber outside the Maramureș territory or more likely that the samples
were ancient. Some beams might be from older structures and then re-
5. Discussion
used in a younger building. Several holes, carving or joining stem was
observed in different position of the cored beams compared to the
5.1. The Maramureș oak tree-ring chronology
current one which provide structural evidence to support the hypoth-
esis. Historical documents indicate a long presence of civilization in the
In this study, we collected and processed 824 samples resulting in
study area (Drimba, 1987; Filipascu, 1940; Ivanciuc, 1998). The tu-
103,371 rings measured and dated, resulting a 781-year long chron-
multuous historical wars, burns and destructions were the main factor
ology suitable for dendrochronological dating. A total of 271 samples
Table 3
Teleconnection between Maramureș oak chronology and living oak chronologies from SE Europe / regional oak chronologies.
No Location Length (year) Date end Glk tv tBP tH Source
*Glk – signature Gleichläufigkeit; tv – t-value computed without any standardization; t-values calculated using the Baillie-Pilcher algorithm (tBP); t-values calculated
using the Hollstein algorithm (tH). *#12 - #15 represent Regional European Oak Chronology.
109
C. Nechita et al. Dendrochronologia 52 (2018) 105–112
of periodic reconstruction of the buildings. Growth anomalies and Chain, induces a local pattern on the response of trees to climate
narrow rings were also reasons for excluding samples from the proxy (Nechita et al., 2017). At a large scale, the former MM chronology was
chronology. Several other studies highlight similar problems (Haneca strongly correlated (tH = 9.6) with East Austrian Oak Chronology re-
et al., 2009; Kolář et al., 2013; Prokop et al., 2017). ference (Grabner et al., 2014) at a distance of almost 600 km, but
The resulting chronology will have a potential to be used in pa- showed a lower correlation (tBP = 8.61) in the comparison with the
leoclimate reconstruction owing the robust common signal indicated by much closer Czech oak chronology (Prokop et al., 2017). The compar-
the EPS (Buras, 2017) and the well-expressed and time stable correla- ison between updated MM oak chronology and East Austrian Oak
tion between oak radial growth and precipitation found in this region Chronology (site #14) indicate tBP = 7.3 (tH = 6.6). Furthermore, the
(Nechita et al., 2017). The quality of the chronology could be improved North Central Hungary oak chronology (site #15) located between
at the end of 18th and early the 19th century when the EPS reached the Austria and Maramureș indicate almost the same indices tBP = 7.4
lowest thresholds of 0.87 and 0.86, respectively, by adding more his- (tH = 7.3), justifying the spatial stability of teleconnection.
torical wood series. The parameters of the chronology indicated robust According with tBP -value (3.1) was observed insignificant correla-
signal strength as analysed by statistical parameter values: inter-series tion with oak chronology from Bulgaria - Zvezdec (#10), located in
correlation (Rbar) and the expressed population signal (EPS) (Fig. 2b) protected nature reserve natural forest. Contrary to our expectations,
(Wigley et al., 1984). Additionally, we noted the existence of subfossil correlation of oak chronology from Sremska Mitrovica, Serbia, (Čufar
oaks imbedded in the riverbed of the Iza River, which indicate the et al., 2014) with MM chronology show significant values for tBP (4.5)
potential to extend the chronology further back in time. and tH (4.7). We noticed that comparing with South Austria and Hun-
Two widespread species, pedunculate oak and sessile oak, pre- gary the correlation is lower, but even for this area MM chronology has
dominate within the oak forests of Maramureș. However, both species realistic potential for being used as a reference chronology.
could not be distinguished from each other by anatomy Ważny et al. (2014) attest to the necessity for constructing a long
(Schweingruber, 1993; Schoch et al., 2004) in archaeological and his- oak chronology in south-eastern Europe for "bridging gaps" to link north
torical wood, and we can only suppose they existed in mixed popula- central European and east Mediterranean tree-ring networks. The new
tions. Feuillat et al. (1997) and Mosedale et al. (1998) attempted to MM chronology come to fill this gap.
distinguish Quercus robur L. and Quercus petraea (Matt.) Liebl. based on Ważny et al. (2014) and Nechita et al. (2017) noticed the role of
the ratio between the number of earlywood vessels and the latewood Carpathians Chain to induce low correlation between eastern and
ratio; however, this approach was not used in this study. Nechita et al. western side of this mountains, also in the south. MM chronology seems
(2017) studied the common climatic signal for both species. The results to cross-date well with SE Poland oak chronology (Krąpiec, 1998) (tBP /
indicated irrelevant differences between the species responses to pre- tH > 3.7), explain by extended influence of North Atlantic air mases in
cipitation in April and June. This information highlights the importance our study area. Even so, tree ring signal can be modified by local
of the MM oak tree-ring chronology, which can be used as a proxy for condition, genetic factors, provenance of the regeneration material,
the reconstruction of past climate and environmental changes for the disturbances or anthropogenic influences, common in lower altitude.
last 611 years.
Several studies, analysing the teleconnection between oak chron-
ologies in Europe, have referred on the former Maramureș oak tree-ring 5.2. Sapwood rings of oaks in NW Romania
chronology in an unpublished form. Information outlined in Kolář et al.
(2012) and Ważny et al. (2014) for the past 600 years corroborates with Commonly, when dendrochronological samples are collected for
that of a previous study on the Maramureș living tree-ring data (Nechita dating from constructions, they do not contain the complete sapwood,
et al., 2017), indicating the necessity of building a long oak chronology because most builders remove a part of the outermost sapwood rings.
in eastern Europe. The primary reason is because of the Carpathians The estimation of missing sapwood tree rings is of great importance in
dendroarchaeological dating, which differ in number from west to east
110
C. Nechita et al. Dendrochronologia 52 (2018) 105–112
Europe (Baillie et al., 1985; Prokop et al., 2017). The number of sap- earliest/latest felling date. For timbers with no sapwood the pos-
wood tree rings for the Maramureș region varied between 6 and 40. sible earliest felling date could be estimated by the addition of 8
Mean sapwood ring number ranged from 10.6 to 16.8. According to our years to the date of the final ring.
results, sapwood tree-ring variability did not depend on altitude or
interspecific characteristics. Special thanks also go to anonymous reviewers which help us to
The number of mean sapwood rings found in our study (n = 14.5) is improve quality of this manuscript, as well as Katarina Čufar, Zoltán
consistent with the results of Prokop et al. (2017) in the Czech Republic Kern and Tomasz Ważny for share with us their chronologies.
(n = 14.6). Compared with other studies in different parts of the Eur-
opean continent such as: Finland and western Estonia (4–21), eastern Acknowledgements
Estonia, Latvia, and Lithuania (6–18) (Sohar et al., 2012), south Finland
(7–24) (Baillie et al., 1985), Ireland (14–62) (Baillie, 1982), the British This work was supported by a grant from the Romanian National
Isles (14–66) (Hillam et al., 1987), north west England and north Wales Authority for Scientific Research and Innovation, CNCS/
(10–55) (Hughes et al., 1981), north England (10–60), south England CCCDI–UEFISCDI, project number PN-III-P2-2.1-PED-2016-1058,
(4–57) (Miles, 1997), north France (12–49) (Pilcher, 1987), Belgium within PNCDI III, “A new technique regarding dendrochronological
(13–39) (Haneca et al., 2009), Germany (7–66) (Hollstein, 1965,1980), dating. Statistical, biological and chemical approach” (DendroTECH).
Hungary (7–36) (Grynaeus, 1998), and Poland (9–36) (Ważny, 1990), EEA Financial Mechanism 2009–2014, contract no. 18SEE “CLIMFOR”
our results confirmed a decreasing trend of the number of sapwood also supported this work. Thanks to Roibu Catalin Constantin and
rings, expressed in absolute range, from west to east. Mursa Andrei for their support and for assistance in fieldwork. Special
thanks to Alexandru Dumitru Babos who initiated the dating work of
5.3. Age related growth trends based on the MM dataset the Maramureș churches.
111
C. Nechita et al. Dendrochronologia 52 (2018) 105–112
112