MA Thesis Marije Vlaar 2007
MA Thesis Marije Vlaar 2007
MA Thesis Marije Vlaar 2007
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Reconstruction of the Palaeoecology<br />
of the Eem Polder by<br />
means of dendrochronology, pollen<br />
and macrofossil analysis<br />
<strong>Marije</strong> <strong>Vlaar</strong><br />
Master thesis Utrecht University, <strong>2007</strong><br />
Supervision:<br />
- Prof. dr. E. Jansma<br />
The National Service for Archaeology, Cultural<br />
Landscape and Built Heritage (RACM)<br />
Dutch Centre for Dendrochronology (Ring<br />
Foundation)<br />
- Prof. dr. G.J. van der Zwaan<br />
Utrecht University
Index ...........................................................................................................2<br />
Preface.......................................................................................................4<br />
Abstract.....................................................................................................6<br />
1. Introduction........................................................................................7<br />
2. Material and Methods.....................................................................9<br />
2.1 Location and site description...................................................9<br />
Geomorphological development of the Eem Valley<br />
2.2 Material: Sampling sites (Location 1 and 2).......................... 15<br />
2.2.1 Location 1, Groeneweg............................................... 15<br />
2.2.2 Location 2, Lodijk....................................................... 17<br />
2.3 Methods..................................................................................... 17<br />
2.3.1 Dendrochronology....................................................... 17<br />
General<br />
Calendars<br />
Sample preparation and the dating process<br />
Pointer years<br />
2.3.2 Pollen.............................................................................20<br />
General<br />
Sample preparation<br />
Counting and diagrams<br />
Comparison with known data<br />
2.3.3 Macrofossils...................................................................22<br />
General<br />
Sample preparation<br />
Analysis<br />
2.3.4 Radiocarbon dating........................................................23<br />
General<br />
Sample preparation<br />
3. Results and Interpretation.............................................................26<br />
3.1 Dendrochronological results.....................................................26<br />
Pointer years<br />
3.2 Pollen analysis............................................................................36<br />
Description of pollen diagram GP01<br />
Description of pollen diagram GP02<br />
3.3 Macrofossils................................................................................40<br />
3.4 Radiocarbon dating....................................................................41<br />
4. Discussion..............................................................................................43<br />
Dendrochronology<br />
Pointer Years<br />
Pollen and Macrofossil analysis<br />
2
5. Conclusions......................................................................................... 49<br />
Groeneweg<br />
Lodijk<br />
6. Further research............................................................................... 50<br />
7. References...........................................................................................51<br />
Appendix...................................................................................................55<br />
Appendix 1.......................................................................................55<br />
Appendix 2.......................................................................................56<br />
Appendix 3.......................................................................................59<br />
Appendix 4.......................................................................................61<br />
Appendix 5.......................................................................................67<br />
3
Preface<br />
After almost five years of studying Earth Sciences it was very pleasant to apply my<br />
knowledge during a final internship and to combine this with learning about a scientific<br />
discipline that was new to me. Starting my internship at the National Service for Archaeology,<br />
Cultural Landscape and Built Heritage (RACM), I discovered the field of dendrochronology,<br />
a method I had never worked with before. Coming from a background where time scales<br />
stretch out over millions of years, it was interesting to work with a dating method that is exact<br />
to the year. Besides the accuracy of dating, the dendrochronological analysis is also very<br />
understandable. Already during my first day I was able to date an old oak trunk. The<br />
combination of the study of dendrochronology with a pollen and macrofossil analysis<br />
supported my interest for palaeo-ecology and palaeo-botany.<br />
Prior to the different analyses that mostly occurred by using a microscope at the laboratory we<br />
accomplished some fieldwork in the Eem Polder. During these days we searched for tree<br />
remains and we sampled the soil for the pollen analysis. Large trunks were discovered in the<br />
field or were already collected by farmers. The practical activities of collecting, sawing and<br />
sampling made the fieldwork very pleasant as well.<br />
The last six months were very informative and interesting. This report is the result of all<br />
different analyses, discussions and interpretations. I never could have written this report<br />
without the help of many different people and I’m very grateful for their help and enthusiasm<br />
during this project.<br />
First of all I would like to thank my supervisor Esther Jansma (RACM, Ring Foundation).<br />
She introduced me to the field of dendrochronology and shared her enthusiasm about this<br />
science. She also performed many practical activities that were very useful during this<br />
research: the development of a new reference calendar was essential to date some of the trees<br />
by means of dendrochronology. It was very nice that she introduced different persons to me<br />
who were willing to participate in this research. Also, she had spent a lot of time on the<br />
revision of my thesis and she gave me useful comments to improve my report. I also want to<br />
congratulate her with her recent installation at Utrecht University as holder of a special chair<br />
dedicated to dendrochronology and palaeo-ecology of the quaternary.<br />
I could not have performed the pollen analysis without the help of Otto Brinkkemper<br />
(RACM). With great patience he taught me to determine all different vegetation species. He<br />
also scanned all samples on rare types and showed them to me to increase my knowledge.<br />
Besides the pollen analysis he dedicated a lot of his time to help me with many other<br />
problems and it was very pleasant to discuss some of the results together.<br />
Also, great regards go out to Marta Domínguez Delmás (Ring Foundation) for her very<br />
pleasant company, her help with dendrochronological problems and her comforting words<br />
especially during the last phase of this project. I want to thank Tamara Vernimmen (formerly<br />
of Ring Foundation) for her help with the macrofossil analysis and Elsemieke Hanraets<br />
(formerly of Ring Foundation) for her help with some dendrochronological measurements. I<br />
want to thank Bert van der Zwaan (Utrecht University) for his supervision, for bringing me in<br />
contact with Esther Jansma and for showing me the option of doing my final internship at the<br />
RACM.<br />
This whole project would not have existed without the enthusiasm of the Eem Polder<br />
volunteer group with Jantine Strating, Carla Vintges and Gerard Monkhorst. They started the<br />
project ‘Leefbaar Landschap’ in collaboration with the IVN (Association for natural and<br />
4
environmental education) to improve the public interest for the Eem Polder and were<br />
especially interested in the natural history of this area. Also the involved farmers, Kees van<br />
Valkengoed en Jan van Halteren were enthusiastic about this project. They showed us some<br />
tree remains that they found and tolerated that we excavated their grasslands. We couldn’t<br />
have performed the fieldwork without the help of Wim Jong (RACM). Due to his field<br />
activities I obtained some nice wood samples at the laboratory. Thanks to Hans van der Plicht<br />
(Centre for Isotope Research; University of Groningen) some wood and soil samples were<br />
radiocarbon dated. Despite the very busy schedule of his research centre he was willing to<br />
date my samples rapidly. He also welcomed me at his Centre of Isotope Research at the<br />
University of Groningen, showed me the laboratory and explained me the dating methods.<br />
The radiocarbon dating could not have been performed without Jos Deeben (RACM) who<br />
took care for the finances on behalf of the RACM. I want to thank the Laboratory of<br />
Palaeobotany en Palynology (Utrecht University) for giving me the opportunity to prepare the<br />
pollen and macrofossil samples, with a special thanks for Jan van Tongeren (Utrecht<br />
University) who performed most of the chemical treatments in their laboratory. Frans Bunnik<br />
(TNO-NITG) helped me with some problems during the pollen analysis and it was nice to<br />
discuss the vegetation history with him. I want to thank all other persons of the RACM who<br />
were interested in my project; Menne Kosian for giving me nice maps, Hans Peeters for<br />
literature and Ruud de Man for helping me to find my way in the beautiful library.<br />
Finally, I want to thank the National Service for Archaeology, Cultural Landscape and Built<br />
Heritage (RACM) and the Ring Foundation for giving me the opportunity to perform my final<br />
internship in a nice atmosphere and for using their facilities.<br />
5
Abstract<br />
The aim of this research is to reconstruct the palaeo-ecology of the Eem Polder by means of<br />
dendrochronology, pollen and macrofossils. Dendrochronology involves the study of annual<br />
ring-width variations through time. The tree-ring patterns in wood reflect the annual growth<br />
conditions, depending on different factors such as precipitation, temperature, hydrology, soil<br />
conditions and others. Synchronization of tree-ring patterns of sub-fossil wood with<br />
absolutely dated reference calendars makes it possible to date the trees. If trees can not be<br />
dated by means of dendrochronology, it is possible to date wood samples with the<br />
radiocarbon dating method. Pollen and macrofossil analyses are used to study the surrounding<br />
vegetation development. Together with the tree-ring pattern they contain information about<br />
the development of the soil on which the vegetation was growing.<br />
The Eem Polder is an area at the Eem Valley. This valley is formed during the Saalien period<br />
by ice sheets which covered large parts of Europe, including the northern part of the<br />
Netherlands. Tree remains are sampled at two different locations within the Eem Polder. At<br />
location Groeneweg two oaks (GFE01 and GFE02) are found in situ and they are dated by<br />
means of dendrochronology and radiocarbon. One of the studied oaks lived from 2719 until<br />
2570 BC and the other one from 2882 until 2581 BC. Synchronization with oaks from<br />
Ypenburg indicates that they were growing on an unsaturated substrate. The radiocarbon<br />
dates of GFE02 differ from the dendrochronological dates and it might be necessary to<br />
reconsider these radiocarbon dates. The tree-ring patterns of GFE01 and GFE02 show some<br />
simultaneous pointer years as the eight oaks from Ypenburg, caused by similar response to<br />
changing growth conditions.<br />
The pollen analysis that is performed on samples from the Groeneweg show a dominance of<br />
trees over herbs, indicating that the vegetation history of this area is characterized by the<br />
growth of a closed forest. Some macrofossil remains and also the appearance of clumps of<br />
pollen show the local appearance of species. Anthropogenic indicators are present in very low<br />
quantities indicating a minor influence of people in this area. The absence of water plants or<br />
peat forming plants confirm the unsaturated characteristic of the soil as shown by the studied<br />
oaks. From the results of the radiocarbon dating of three soil samples it turns out that the<br />
vegetation development took place during the Subboreal. The dates show that the studied oaks<br />
lived during the same period as the development of the studied profile occurred. The tree-ring<br />
pattern of GFE01 and GFE02 show that the trees died because of a sudden event and the<br />
preservation of the trees show that they were buried soon after their death. The pollen<br />
diagrams did not show any remarkable decrease of tree pollen during the studied period<br />
indicating that not the whole forest was influenced.<br />
Two oaks from location Lodijk are dated with a reference calendar build up with oaks that<br />
were growing in the province of Zuid-Holland. Oak LFE01 lived from 1678 until 1498 BC<br />
and LFE02 from 1688 until 1476 BC. The tree-ring patterns of the oaks show similarities with<br />
pointer years of 32 trees from the province of Zuid-Holland. Especially the rings formed<br />
during 1629 and 1628 BC were remarkable smaller. The trunks of LFE01 and LFE02 indicate<br />
that they were growing in a closed forest and other discovered root systems confirm this. The<br />
relation with undated pines from the same location is unknown, the stem form and shallow<br />
roots of these pines show that they were growing on a saturated, unstable substrate. In<br />
summary it can be stated that 5000 to 3000 years ago a closed-canopy forest grew in what<br />
now is called the Eem Polder.<br />
6
1. Introduction<br />
The surface on which we live is a treasure trove of earth’s own history. Sediments hold<br />
information about activities that once took place, a million years ago but also during the<br />
recent past. There are different ways to search for this information; from large glaciers in the<br />
mountains to the oceans floor, but also on smaller scales in our surrounding environment.<br />
The aim of this research is to reconstruct the palaeo-ecology of a part of the Eem Valley by<br />
performing a multi-proxy analysis on vegetation remains. Palaeo-ecological research<br />
concentrates on fossils and subfossils to reconstruct ecosystems of the past. This includes all<br />
organisms, their life cycle and natural environment. The study presented here only<br />
concentrates on subfossil vegetation remains to reconstruct the vegetation history.<br />
Knowledge about the palaeo-ecology of this valley is interesting from several points of view:<br />
From a biological point of view it is interesting to know which species were growing at a<br />
certain time and how their succession occurred. Looking at larger timescales, migration<br />
patterns of species through time are also interesting.<br />
From an archaeological point of view, vegetation history tells something about cultural<br />
development during certain periods. If vegetation patterns in the past change suddenly, it<br />
might be possible that people influenced the landscape by cutting trees or cultivating plants.<br />
Human migration and vegetation are also related; people only settle on areas which are<br />
attractive for them.<br />
There is also a natural-historical interest; many scientists research vegetation changes<br />
through time. There are several hypotheses about the Dutch palaeo-landscape. One of the<br />
theories states that the Dutch vegetation history was characterized by an open landscape with<br />
a mosaic of grasslands, shrubs and small tree groups. An autochthon fauna of herbivores<br />
prevented that the open landscape turned into a closed forest (Vera, 1997).<br />
Others state that the natural vegetation was characterized by a closed deciduous forest with<br />
oak (Quercus sp.), elm (Ulmus), lime (Tilia), ash (Fraxinus excelsior), beech (Fagus<br />
sylvatica), hornbeam (Carpinus betulus) and hazel (Corylus avallana) being the dominant<br />
species. The autochthon fauna of herbivores did not have influence on the succession of the<br />
forests. Rejuvenation took place in forest gaps originated by the collapse of single trees,<br />
following the gap-phase model of Watt (1947).<br />
This research is focused on a restricted area called the Eem Valley. This area received<br />
attention of public and scientists when farmers found tree remains perfectly preserved, lying<br />
half-buried under the surface. Because nothing was known about the age of the trees, we<br />
decided to conduct a dendrochronological research and found out when they had been living.<br />
In addition, to gather information about the ecological context in which the trees were<br />
growing, soils samples were taken and pollen and plant-macrofossils were analyzed. This<br />
multi-proxy approach provides information and gives more certainty about the vegetation<br />
history of this area.<br />
This research has attracted the attention not only of the farmers but also of some volunteers<br />
who want to improve the cultural interest for this area. By inserting information signs along<br />
roads for walkers and cyclists they want to promote the historical and natural value of this<br />
area. Some results of this research have been included on their information signs.<br />
7
To reconstruct the palaeo-ecology of the research area some questions should be answered:<br />
1. How many growth phases are represented in the different tree trunks that are found?<br />
- When did the trees germinate and die?<br />
- Do the tree-ring patterns or the shape of the trees give information about the<br />
environment in which they lived?<br />
2. Is it possible to relate the tree growth of the different phases to soil development?<br />
- How did the soil develop during the different phases?<br />
- On what kind of soil did the vegetation grow?<br />
3. Is it possible to discern relevant changes in the tree-growth pattern? How did the<br />
climatic conditions of that period affect the tree growth?<br />
- Do the growth patterns of the trees growing together in the Eem Valley<br />
resemble each other?<br />
- Do these patterns match with trees growing at the same time at other locations?<br />
- Were there major significant climatic changes during the growth period of<br />
these trees, and if so, are they recognizable in the growth pattern?<br />
4. What was the distribution and diversity of the surrounding vegetation during the<br />
different phases?<br />
- How is the vegetation related to time?<br />
- What plant species were present according to the palynological analysis? Is it<br />
possible to describe local and regional vegetation?<br />
- Is the vegetation development influenced by anthropogenic factors?<br />
- Does the palynological analysis agree with the dendrochronological research?<br />
- Is the distribution and diversity of vegetation related to any of the hypothesis<br />
about the natural Dutch palaeo-landscape?<br />
This study will try to answer these questions. After this introduction I will describe the<br />
‘Material and methods’ in chapter 2. First I will present the geomorphological development of<br />
the Eem Valley (chapter 2.1), followed by a description of how we sampled the two sites<br />
(chapter 2.2). The chapter will end with the different methods of dendrochronology, pollen<br />
and macrofossil analysis and radiocarbon dating (chapter 2.3).<br />
Chapter 3 will present the results of the dendrochronological research (chapter 3.1), the pollen<br />
analysis (chapter 3.2), the macrofossil analysis (chapter 3.3) and the radiocarbon dating<br />
(chapter 3.4). The results are discussed in chapter 4 and the conclusions are given in chapter<br />
5. I will end this report with some suggestions of further research (chapter 6), the references<br />
used for this report (chapter 7) and finally the appendixes.<br />
8
2. Material and Methods<br />
2.1 Location and site description<br />
The Eem Polder is situated in the Eem Valley and it belongs to a nature reserve called<br />
‘Arkemheen- Eemland’. Figure 1 shows the current perimeter of the nature reserve and the<br />
two locations (1 and 2) where the vegetation remains were collected.<br />
1<br />
2<br />
Border National<br />
Landscape<br />
Open Area<br />
Location Groeneweg<br />
Location Lodijk<br />
Figure 1: The Eem Polder is situated in the central part of the Netherlands. The two studied sites are within the<br />
nature reserve.<br />
Geomorphological development of the Eem Valley<br />
The development of the Eem Valley can be traced back to about 370.000 years ago, during<br />
the Saalien period. At that time ice sheets covered large parts of Europe, including the<br />
northern part of the Netherlands, with its southernmost boundary between Haarlem and<br />
Nijmegen. The Saalien period is characterized by a succession of relatively warm and cold<br />
periods, also called the ‘Saalien Complex Stage’, after Litt and Turner (de Mulder, Geluk,<br />
Ritsema, Westerhoff and Wong, 2003), with the coldest phase at the end of the period (Figure<br />
2).<br />
9
Figure 2: The Saalien Complex Stage is characterized by a variation of warm and cold periods<br />
(after de Mulder et al., 2003).<br />
The start of the Saalien (Stadial I after Zagwijn, 1973) was characterized by a scarce<br />
vegetation cover with birch (Betula sp.) and pine (Pinus sp.) trees and a more abundant cover<br />
of grasses (Poaceae) (de Mulder et al., 2003). The uncovered parts of the landscape were<br />
influenced by wind which deposited aeolian sediments. During the next stage, the so-called<br />
Hoogeveen Interstadial (Zagwijn, 1973), trees became more abundant, but birches and pines<br />
still were the dominant species (de Mulder et al., 2003). After the Hoogeveen Interstadial<br />
period a short colder period occurred again, which is called Stadial II (Zagwijn, 1973) (de<br />
Mulder et al., 2003). During this short cold phase the vegetation became more open again,<br />
and local aeolian sediments were formed. This cold period was followed by a warmer period<br />
called the Bantega Interstadial after a pollen analysis on material collected in Bantega.<br />
During the last and coldest phase of the Saalien (185.000-130.000 yr B.P.), ice sheets reached<br />
their maximum movement to the south creating depressions and ice-pushed ridges in the<br />
landscape. Figure 3 shows the maximum expansion of the ice sheets during this period, with a<br />
large glacial depression representing the Eem Valley.<br />
10
Figure 3: During the coldest phase of the Saalien, a glacial<br />
depression formed the Eem Valley (after de Mulder et al., 2003).<br />
Ice-pushed ridges are hills formed by the pressure of glaciers on sediments deposited<br />
previously. In contradiction, moraines are built up by sediments transported by the glaciers<br />
from elsewhere (de Mulder et al., 2003). In the Netherlands the Utrechtse Heuvelrug is the<br />
best known example of an ice-pushed ridge.<br />
Depressions are caused when ice sheets produce large amounts of melt water that erodes the<br />
substrate underneath. The ice will sink into the depression, where it will grow thicker and<br />
cause even more erosion. When finally the ice sheets melt, the depression will be filled with<br />
melt water and as a result a lake is formed. The most important depressions surrounded by<br />
ice-pushed ridges at the Dutch landscape were formed during the latest phase of the Saalien<br />
period. Examples are the Eem Valley, the IJsseldal and Nordhorn in Germany (de Mulder et<br />
al., 2003).<br />
After the Saalien period the Late Pleistocene starts which is divided into the Eemien<br />
(130.000-115.000 yr BP) and the Weichselien (115.000-10.000 yr BP) period. The Eemien<br />
period is called after the river Eem, which has also given its name to the Eem Valley and<br />
which ends into the IJssel Lake. During this period global temperature increased, causing<br />
large ice sheets in Northern Europe and North America to melt. As a result sea level rose to a<br />
level even higher than it is at present (de Mulder et al., 2003).<br />
11
Figure 4: The Eem Valley was flooded during the transgression<br />
phase of the Eemien (after de Mulder et al., 2003).<br />
The map in figure 4 shows the transgression phase of the Eemien, with large parts of the<br />
Netherlands, including the Eem Valley, flooded because of the high sea level. The ice cover<br />
during the Saalien influenced the morphology of the landscape and determined the influence<br />
of this marine transgression. The transgression phase caused deposition of marine sediments<br />
and along these areas peat development occurred on a large scale, mainly in the Eem Valley<br />
and in the Hunzedal (province of Drenthe) (de Mulder et al., 2003).<br />
The Weichselien is the latest glacial period: glaciers covered parts of northern Europe without<br />
reaching the Netherlands. The early Weichselien (115.000-73.000 yr BP) is characterized by a<br />
cold climate with an open landscape where pine and birch trees dominated. This cold stage<br />
also consisted of warmer periods during which forests developed, still dominated by birch and<br />
pine but with some spruce (Picea abies) present. In the southern part of the Netherlands oak<br />
formed an important part of the vegetation development. The middle-Weichselien (73.000-<br />
13.000 yr BP) started with a decrease of the global temperature, causing expansion of land ice<br />
on higher latitudes and a lowering of the sea level. The forests disappeared due to the low<br />
temperatures and erosion influences the landscape.<br />
The late-Weichselien (13.000-10.000 yr BP) is characterized by a variation of short warm and<br />
cold stages. It started with a relative warm stage with an open landscape dominated by birch.<br />
After a cold phase the vegetation recovered and pine trees took over the pioneer vegetation of<br />
birch and juniper (Juniperus communis) (de Mulder et al., 2003). During the latest<br />
millennium of the Weichselien temperature decreased again.<br />
12
During the Holocene (10.000 yr BP - present), melting of ice sheets in North-America and<br />
Scandinavia caused sea level to rise again. The North Sea, which dried during the<br />
Weichselien, was flooded and inundated larger parts of the Netherlands than it did during the<br />
transgression phase of the Eemien (de Mulder et al., 2003). The fast rise of the sea level<br />
ended about 7000 yr BP (Figure 5). When sea level rises, subsequently groundwater levels<br />
rises too, triggering environmental changes and influencing vegetation development (Peeters,<br />
<strong>2007</strong>).<br />
Figure 5: The fast sea level rise caused by melting of ice sheets ended<br />
7000 years BP (after de Mulder et al., 2003).<br />
The Holocene period is the most relevant period to this study because it is the first period in<br />
which people influenced the landscape. Important changes in the vegetation development<br />
occurred due to changes in the precipitation regime (de Mulder et al., 2003), but also because<br />
of cultural development. At the start of the Holocene vegetation developed as a natural<br />
response to climate and environment and was reflected by the succession of tree species<br />
(Figure 6). During the Preboreal birch was the most common species, after which pines<br />
became dominant. During the Atlanticum a mixed oak forest developed (de Mulder et al.,<br />
2003). At the end of the Atlanticum culture developed towards permanent settlements and<br />
agriculture influenced the natural vegetation, which causes a relative decrease of trees and the<br />
first appearances of cultivated plants. During the Subboreal the decrease of trees and relative<br />
increase of herbs continued. During the Subatlatic the presence of natural vegetation came to<br />
an end and the natural landscape turned into a cultivated landscape.<br />
13
Figure 6: Vegetation development of the Late Weichselien and Holocene. This pollen diagram shows the general<br />
vegetation history of the Netherlands (after de Mulder et al., 2003).<br />
14
2.2 Material: Sampling sites (Location 1 and 2)<br />
The different sorts of vegetation remains used in this study were collected from various<br />
locations within the Eem Valley (Figure 1). Some of the remains were found in situ, buried in<br />
a grassland along the Groeneweg in Bunschoten (location 1). Other remains were gathered by<br />
a farmer along a highroad and on his grassland near the Lodijk in Bunschoten (location 2).<br />
The original location of the samples from the Lodijk is unknown.<br />
2.2.1 Location 1, Groeneweg<br />
Two large trees were identified in a grassland along the Groeneweg, because part of the<br />
trunks was visible in the adjacent ditch. These two trees are referred as GFE01 and GFE02<br />
(Groeneweg, Fossil oak, Eem Valley). The geographical coordinates (x;y) of GFE01 are<br />
155498.399; 472548.826 and those of GFE02 are 155473.385; 472559.495 (NedGraphics<br />
CAD/GIS). To determine the length and shape of the trees and the direction in which they<br />
were lying, metal sticks were inserted in the ground. With permission of the owner of the<br />
grassland a part of both trees was exposed by excavating (Figure 7a/b). In the exposed parts<br />
cross sections were cut with a chain saw. The thickest parts were cut into thinner samples at<br />
the laboratory (RACM) to make them fit under the microscope. From GFE01 two different<br />
sections were sawed, one closer to the roots and the other closer to the branches. From GFE02<br />
only one section was sawed. Because these trees were found in situ the surrounding sediments<br />
contain information about the environmental context of their growth. This information is<br />
present in the form of leaves, roots, seeds, pollen and spores. Right besides each stem, a<br />
sample of the profile was taken with a metal container (50x15x10 cm) (Figure 8a/b). These<br />
would cover the soil profile in which the trees were lying. The profiles are referred as GP01<br />
and GP02 (Groeneweg Profile) sampled near tree GFE01 and GFE02 respectively. Figure 9<br />
shows a schematic overview of the two trees and their location in the field. Their depth and<br />
the position of the profiles are shown in figure 10a and 10b. The two containers were stored at<br />
low temperatures to prevent decay of the material.<br />
Figure 7a: Tree GFE01.<br />
Figure 7b: Tree GFE02.<br />
Figure 8a: Profile GP01.<br />
Figure 8b: Profile GP02.<br />
15
Figure 9: A schematic overview of GFE01 and GFE02 in the grassland. Both trees were lying in the same<br />
direction. The black squares present the excavated parts with the sawed cross sections in white. The red lines<br />
present the location of the sampled soil profiles.<br />
Figure 10a: Schematic overview of GFE01 and its in situ location.<br />
Figure 10b: Schematic overview of GFE02 and its in situ location.<br />
16
Figure 11: LFE01 and LFE02 were collected by a farmer along<br />
a highway.<br />
2.2.2 Location 2, Lodijk<br />
At this location different trees were<br />
collected by a farmer. Two large trees<br />
were found along a highway where<br />
they became visible because of<br />
building activities (Figure 11). They<br />
are referred to as LFE01 and LFE02<br />
(Lodijk, Fossil oak, Eem Valley).<br />
LFE01 is 200 centimetre long and its<br />
diameter is 48 centimetre. LFE02 is<br />
255 centimetre long with a diameter<br />
of 37 centimetre. From LFE01 one<br />
cross section was taken, from LFE02<br />
two different ones. Other smaller tree<br />
fragments were gathered by the<br />
farmer on his land. Only the ones<br />
with enough rings for analysis were<br />
sampled. The original location of these samples is unknown because the trunks were moved<br />
prior to their discovery by the farmer and therefore no sediments were taken from this<br />
location.<br />
2.3 Methods<br />
To reconstruct the palaeo-ecology of the Eem Polder three different approaches were<br />
combined:<br />
1. The cross sections of trees were analyzed by means of dendrochronology<br />
2. In the soil sediments the content of pollen was analyzed<br />
3. A macrofossil analysis was performed by the determination of seeds, leaves, roots and<br />
small wood fragments<br />
2.3.1 Dendrochronology<br />
General<br />
The word dendrochronology is derived from the Greek words ‘dendron’ (tree), ‘chronos’<br />
(time) and ‘logos’ (the science of) and it refers to the science that studies annual tree-ring<br />
variables through time (Grissino-Mayer, <strong>2007</strong>). The main principles of dendrochronology are<br />
described in appendix 1 (Table 13). In temperate climate zones, tree species produce annual<br />
growth rings (a detailed explanation of tree growth and wood anatomy is described in<br />
appendix 2). Those tree rings will be wider or narrower depending on different factors such as<br />
precipitation, temperature, soil conditions and others. Tree-ring patterns of subfossil trees<br />
constitute a suitable record of proxidata for palaeoclimatological and ecological studies<br />
(Jansma, 2006).<br />
In Europe, oak (Quercus robur and Quercus petraea) is the most commonly used species for<br />
dendrochronological dating (Baillie, 1982). It is difficult to distinguish both species when<br />
remains of them are found. Only on living trees observation of the leaves and acorn cups may<br />
separate the two, but this can also be complicated because of the existence of hybrids (Baillie,<br />
1982). Despite of this, they can be used because of the common signal on environmental<br />
changes. When they are synchronized, the resulting calendars are useful for dating purposes.<br />
The growth pattern of oaks is influenced by climatic conditions that sometimes affected large<br />
areas. Because different oaks respond equally to annual climatic conditions, growth patterns<br />
17
of different oaks can be cross dated (Jansma, 1995). Cross dating is commonly used in<br />
dendrochronology. Tree-ring patterns of living trees can be matched with patterns of trees<br />
from earlier times, if a part of their growth pattern overlaps. Subsequently these growth<br />
patterns can be used to date older material (Figure 12). Doing so a dated chronology is<br />
developed and measurement series of an undated growth pattern can be dated with it.<br />
Figure 12: Wood can be dated by cross-dating<br />
if a part of their pattern overlaps (Grissino-<br />
Mayer, <strong>2007</strong>).<br />
The last ring of the measurement series (the closest<br />
one to the outer part of the tree) will provide the<br />
closest date to the felling or death of the tree<br />
(Jansma, 1995). Because absolutely dated oaks do<br />
not extend further back than 12,000 years the<br />
dendrochronological dating method is restricted<br />
when it is used for absolute dating. Also, it is not<br />
always possible to date wood with this method for<br />
several reasons such as short tree-ring patterns,<br />
growth disturbances, reduced reaction on climate<br />
changes and the absence of reference calendars.<br />
However, for approximately the last 8,000 years<br />
dendrochronology has an advantage over all other<br />
dating methods, such as stratigraphy or<br />
radiocarbon dating, because it provides absolute<br />
dates, exact to the year (Jansma, 1995).<br />
Calendars<br />
Dated chronologies (also called calendars) are useful to match undated growth patterns. The<br />
calendars are built up by detrending and averaging cross dated series. Individual tree growth<br />
is dependent on different environmental factors which influence the ring-width series. The<br />
principle of aggregate tree growth (Appendix 1, table 13) is used to standardize individual<br />
tree-ring patterns. The oldest dated annual ring in the Dutch dataset was formed in 6025 BC.<br />
Following the overview from the national research agenda for archaeology (Jansma, 2006)<br />
dendrochronological data replication is described for eight consecutive archaeological<br />
periods. The overview of Dutch calendars described below is based on data available at the<br />
RACM and the Ring Foundation. The detailed overview of these periods is shown in<br />
appendix 3.<br />
1. 6025 – 2850 BC (Mid and New Stone Age)<br />
2. 2850 – 2000 BC (New Stone Age)<br />
3. 2000 – 800 BC (Bronze age)<br />
4. 800 – 12 BC (Iron age)<br />
5. 12 BC – AD 400 (Roman time)<br />
6. AD 400 – 1000 (Early Middle Ages)<br />
7. AD 1000 – 1500 (Mid and late Middle Ages)<br />
8. AD 1500 – 1950 (Early modern and modern time)<br />
18
Figure 13: The dendrochronology calendar<br />
contains two hiatus (Jansma, 2006).<br />
Unfortunately the existing calendars contain a<br />
couple of hiatus, periods for which no data is<br />
available (Figure 13). Different factors cause these<br />
hiatus. First, in order to develop useful calendars at<br />
least ten dated samples covering the same period<br />
are needed. Especially from prehistorically periods<br />
often not enough samples are available. Second,<br />
poor conservation of wood, the local absence of<br />
forests and felling activities decrease the<br />
availability of wood.<br />
Sample preparation and the dating process<br />
Prior to the dating of wood, samples have to be prepared. First the surfaces have to be<br />
cleaned. A thin surface layer is removed with a razor blade, making the rings along the radii<br />
visible. It is important to follow the direction of the radii and measure the ring widths<br />
perpendicular to tree-ring boundary. After cleaning, white chalk is put on the sample to fill<br />
the vessels and to enhance the contrast between early spring vessels and the denser summer<br />
wood (Appendix 2). The width of each individual ring is measured, starting from the first ring<br />
near the pith towards the last ring at the outside, using a measuring device that registers the<br />
ring-width sequences (Past 4, <strong>2007</strong>). The tree-ring series are graphically displayed as a curve<br />
of ups and downs. The presence of sapwood rings or wane edge (Appendix 2) is registered for<br />
each sample. Whenever possible measurements are taken along different radii of the wood<br />
section and averaged into a mean curve per tree. The mean curve obtained from each tree has<br />
to be compared with existing chronologies for dating. Different statistical parameters describe<br />
the similarities between the curves; PV (Coefficient of Parallel Variation, also called<br />
Gleichläufichkeit (GL)), r (cross-correlation coefficient between two series) and t (the value<br />
resulting from a Students t-test on the highest value of r found when comparing the two<br />
series) (Jansma, 1995). PV calculates the fraction of ring widths that at a given position<br />
simultaneously show an increase or decrease relative to the previous ring width (Jansma,<br />
1995). The significance of this test depends on the overlap, a PV-value lower than 0.5 is not<br />
statistical significant. If the PV-value is 1 it means that the annual variation of ring widths is<br />
identical.<br />
‘The r is the correlation between the detrended series (x and y) and is estimated by:<br />
with n the number of values in each series, x i and y i the observed values in year i, x and y the<br />
mean of the series and S x and S y their standard deviation. To assess whether r xy deviates from<br />
zero, one uses the fact that the stochast has a Student distribution with (n-2) degrees of<br />
freedom’ (Jansma, 1995, pp 61).<br />
To produce a correct dendrochronological date the Students t-test must have a value of 5 or<br />
higher.<br />
19
Pointer years<br />
For studying the influence of climate on tree species, extreme ring widths in tree-ring series<br />
are a useful dendrochronological tool (Schweingruber, Eckstein, Serre-Bachet and Bräker,<br />
1990). Such extreme values are termed ‘event years’ and ‘pointer years’. An event year shows<br />
anomalous ring widths in a tree-ring series (Schweingruber et al., 1990) and a pointer year<br />
shows radial growth reduction or increase simultaneously replicated in several tree-ring series<br />
as a result of common forcing (Schweingruber et al., 1990). Several factors influence the<br />
appearance of event and pointer years, e.g. climate region, species, altitude (Desplanque,<br />
Rolland and Schweingruber, 1999).<br />
In this study pointer years of an existing dataset are compared to the ring patterns of the dated<br />
oaks (an overview of the datasets and their origin are shown in appendix 4). The dataset was<br />
selected because the chronologies of this dataset overlap in time with the dated trees of this<br />
research. The dataset contains several chronologies to make a statistical significant overview<br />
of pointer years. Their statistical comparison on pointer years was performed using the IPP<br />
Pointer Years Program, developed by Jansma (1991).<br />
The pointer years are compared on different significance levels;<br />
- 90% which requires a minimum of 5 ring sequences<br />
- 95% which requires a minimum of 6 ring sequences<br />
- 99% which requires a minimum of 9 ring sequences<br />
- 99.9% which requires a minimum of 13 ring sequences<br />
Using these four significance levels, the pointer years can then be visualized.<br />
2.3.2 Pollen<br />
General<br />
(Bennett and Willis, 2002) analyzed pollen content of sediments to investigate past climatic<br />
changes. At present, vegetation changes due to human impact, succession changes and biotic<br />
and abiotic factors also receive attention. Pollen analysis, more often referred as palynology,<br />
is concerned with the study of pollen, spores and other palynomorphs (plant and animal<br />
structures that are microscopic in size), both living as well as fossilized. Pollen grains are<br />
plants parts that contain a male nucleus for fertilization with the female nucleus in an ovule<br />
(Bennett and Willis, 2002). Spores are parts of ferns, mosses and fungi. Most pollen are<br />
dispersed either by wind, insects or water, plants with wind-dispersed pollen have the highest<br />
pollen production. Spores are wind-dispersed.<br />
That fraction of the pollen liberated into the atmosphere but not used for fertilization<br />
accumulates into sediments. Because the walls of pollen grains are very resistant against<br />
degradation, their preservation rate is high. The walls of pollen grains are built up of a<br />
mixture of cellulose and sporopollenin, both components are very resistant and chemically<br />
stable. The latter is resistant against most chemical and physical degradation, except oxidation<br />
(Bennett and Willis, 2002). This is an advantage because a large part of the pollen<br />
accumulated in anaerobic environments will be preserved. The other components of the<br />
sediments can be removed by chemical treatment in the laboratory, without dissolving the<br />
pollen. Birks and Birks (1980) described the main principles of palynology as follows: pollen<br />
and spores are produced in high abundances by plants and most of these produced<br />
palynomorphs accumulate in sediments. Pollen can be identified to various taxonomic levels.<br />
Because they are well mixed by atmospheric turbulence, the result is an uniform pollen rain in<br />
the sediments, where they are preserved in anaerobic environments. The pollen rain is a<br />
function of the composition of the vegetation. If pollen from sediments with a known age are<br />
20
identified, a vegetation reconstruction can be made in space and time. If many pollen samples<br />
are identified in a sequence of sediments, the vegetation changes throughout a period can be<br />
examined. And finally, using pollen spectra from different locations, it is possible to compare<br />
vegetation changes through time at different places.<br />
Sample preparation<br />
In preparation of the pollen analysis of the soil samples taken from the study sites at the<br />
Groeneweg, the surface of the profiles (GP01 and GP02) were cleaned to avoid contamination<br />
from surrounding sediments. Next, 12 samples at 4 centimetre intervals were taken from each<br />
container, starting from the bottom (Figure 14). Each subsample consists of 1 cm 3 material,<br />
taken from the middle part of the profile. All 24 samples were treated at the Laboratory of<br />
Palaeobotany and Palynology from Utrecht University to remove non-palynomorph organic<br />
matter and clastic sediments. Preceding this treatment a tablet with Lycopodium spore was<br />
added to each sample. The amount of Lycopodium in a sample can be used to estimate<br />
absolute pollen values. One Lycopodium tablet contains on average 18583 spores (Berglund<br />
and Persson, 2004) imbedded into the carbonate tablet. By adding hydrochloric acid (HCL)<br />
all carbonates (both of the tablet and of the sediments) were dissolved. Between all different<br />
steps of removing useless constituents with chemicals, the samples are homogenized,<br />
centrifuged and the added chemical is decanted. Some chemicals have more influence on the<br />
material when they are warmed up; therefore the samples were sometimes put in the dry-bath<br />
incubator to warm them. Homogenizing is important to ascertain that the chemicals have the<br />
same influence on the whole sample. It also avoids that heavier pollen grains sink down and<br />
escape the treatments. After removing the carbonates, potassium hydroxide (KOH) was added<br />
to remove organic debris. Unless pollen itself is also organic debris, this treatment will not<br />
solve the resistant wall of the grains. The remaining samples were sieved to remove the large<br />
material, and the sieved material, containing the pollen, was captured. From the remaining<br />
samples the silicates were removed using hydrogen fluoride (HF). Next to the silicates, the<br />
fluoride gels were removed with hydrochloric acid (HCL). Finally, the remaining organic<br />
matter was removed again using potassium hydroxide (KOH). The resulting samples mainly<br />
consist of palynomorphs. The samples were coloured with acetolysis. This makes it easier to<br />
recognize characteristics of certain pollen grains and removes finally the remaining organic<br />
matter except pollen grains. Before the samples were put on a microscope slide, glycerine was<br />
added, to prevent the samples from drying out on the slides. Due to the high density of<br />
palynomorphs one microscope slide for each sample was sufficient.<br />
Figure 14: Both profiles were divided into 12 samples from which the yellow labelled samples were selected for<br />
analysis.<br />
21
Counting and diagrams<br />
Due to time constraints not all 12 samples per container could be investigated. The selected<br />
samples for analysis are shown in figure 14. These samples were chosen to provide an<br />
overview of the whole profile and include the first and last sample. The palynomorphs were<br />
determined using a microscope at a magnification of 400x and for critical determinations a<br />
magnification of 1000x was used. During the analysis pollen and spores were classified and<br />
counted until a number of 300 tree pollen was reached. In addition, remarkable features such<br />
as fungi were registered. From a statistical point of view, counting until 300 tree pollen is<br />
enough to show the vegetation content of a sample (Bennett and Willis, 2002). After<br />
counting, the samples were scanned by O. Brinkkemper on pollen types not noticed before<br />
and these were included as rare types with ‘+’ (one specimen) or ‘++’ (more specimen) in the<br />
pollen diagrams. The amount of pollen and spores per samples were transformed to<br />
percentages and mapped in pollen diagrams against depth (Tiliagraph, 2001)<br />
Comparison with known data<br />
Pollen analysis can be useful to date sediments. In the Netherlands, a lot of research is done<br />
on vegetation development and a general vegetation history has been developed. This<br />
vegetation history is divided into several pollen zones, based on the first appearances or<br />
disappearances of certain species, e.g. the first appearance of beech in the Netherlands occurs<br />
around 2000 BC and the occurrence of beech pollen in a sample therefore indicates the date of<br />
the sediments after 2000 BC. The general overview is often used as a dating background.<br />
Three different samples from profile GP01 were radiocarbon dated and these dates were used<br />
to assign the vegetation to a time period. The radiocarbon dates were given as dates before<br />
present (BP) and were calibrated to dates before Christ (cal. BC) (Wincal25, Centre of Isotope<br />
Research at the university of Groningen).<br />
2.3.3 Macrofossils<br />
General<br />
Plant macrofossils are very helpful when reconstructing the vegetation history. In this study<br />
the term macrofossils refers to leaves, seeds, branches and roots. Compared to pollen, such<br />
remains are more sensitive to degradation, but when they accumulate into anaerobic<br />
conditions they will be well preserved. The advantage of macrofossil analysis over the<br />
analysis of pollen is that macrofossils represent the local vegetation because they have not<br />
become dispersed far from their source (Birks and Birks, 2006), whereas pollen can be blown<br />
in from elsewhere.<br />
Sample preparation<br />
To be able to relate the macrofossils to the pollen, the samples of macrofossils and pollen<br />
were taken at the same height in the profile. A slice of sediment was cut and the edges were<br />
removed to avoid contamination. To remove the material too small for macrofossil analysis,<br />
all samples were sieved using 250 and 125 μm mesh width. The volumes of the samples were<br />
measured before the samples were sieved. The resulting macrofossil samples were stored in<br />
demineralised water to avoid decomposition of the material.<br />
Analysis<br />
The macrofossil samples were analyzed on their content. The species in the fossil remains<br />
were determined and their number was estimated. Besides vegetation remains such as seeds,<br />
roots and leaves, also fragments of insects, fungi and charcoal were noted. Three samples<br />
from profile GP01 were analyzed by T. Vernimmen using a stereomicroscope at a<br />
22
magnification of 50x; (number 1 from the oldest part, number 6 from the middle part and<br />
number 12 from the youngest part).<br />
2.3.4 Radiocarbon dating<br />
General<br />
If wood can not be dated by means of dendrochronology, radiocarbon dating is a good<br />
12 13<br />
alternative. In nature, three different carbon isotopes exist, the two stable isotopes C and C<br />
14 14<br />
and the radioactive isotope C (Van der Plicht, 2006). C atoms are produced when cosmic<br />
particles from galaxies in outer space transform 14 N atoms in our atmosphere into radioactive<br />
14 C atoms, which then slowly decay (Ruddiman, 2002). Earth is protected against these<br />
cosmic rays by its magnetic field and therefore the atmospheric 14 C content is sensitive to the<br />
geomagnetic field strength and also to solar fluctuations (Kitagawa and Van der Plicht, 1998).<br />
14 C is added to the carbon cycle as CO 2 where it is used by living plants for their<br />
photosynthesis. Animals take up 14 C indirectly through their food. As soon as an organism<br />
dies, the exchange with CO 2 stops, and the remaining 14 C decays with a half lifetime of 5730<br />
years (Van der Plicht, 2006). By measuring the remaining 14 C content of a dead organism, the<br />
time that passes since its death can be estimated.<br />
Using the half lifetime of 14 C for dating is only possible if the half lifetime and the original<br />
concentration of 14 C are known, both being a problem within this method. When this method<br />
was first applied (1950) the half lifetime of 14 C was considered as 5568 year, but continuation<br />
of research showed that this half lifetime actually is 5730 years. Also, the atmospheric 14 C<br />
content is not constant because earth’s magnetic field and solar activity vary. Besides these<br />
two problems, isotope fractionation changes the 14 C content of organisms. To solve these<br />
problems a radiocarbon time scale was made, based on three principles (Van der Plicht,<br />
2006);<br />
- The original half-life time (from 1950) of 5568 years is used, instead of the actual half<br />
life time of 5730 years<br />
13<br />
- Fractionation is corrected using the stable isotope C<br />
14 14<br />
- C is measured relative to a standard, equal to the year 1950. C ages are given as BP<br />
(Before Present), where present is the year 1950.<br />
Because of these assumptions the created radiocarbon time scale is not equal to the calendar<br />
time scale, and the differences vary through time due to changing natural 14 C- content.<br />
Therefore the radiocarbon calendar has to be calibrated to an independent calendar.<br />
Radiocarbon calibration can be performed by 14 C dating of samples that also can be dated by<br />
an independent and preferably absolute dating method (Kitagawa and Van der Plicht, 1998).<br />
Tree rings are perfect samples for this purpose since they have an accurate annual time scale<br />
that can be dated by dendrochronology (Kitagawa and Van der Plicht, 1998).<br />
Dendrochronological calibration has been obtained back to 12.000 years. The calibration<br />
curve has been extended back to 26.000 BP, using dated corals (with 14 C and uranium<br />
isotopes) and laminated sediments from the Cariaco basin, in which the foraminifera (single<br />
celled organisms with calcareous tests) were dated (Reimer, Baillie, Bard, Bayliss, Beck,<br />
Bertrand, Blackwell, Buck, Burr, Cutler, Damon, Edwards, Fairbanks, Friedrich, Guilderson,<br />
Hogg, Hughen, Kromer, McCormac, Manning, Ramsey, Reimer, Remmele, Southon, Stuiver,<br />
Talamo, Taylor, Van der Plicht, and Weyhenmeyer, 2004). Recently different calibration<br />
curves have been obtained from varves, marine sediments and stalactites, extending the<br />
calibration curve back to 50.000 year BP. The only terrestrial 14 C samples used for this<br />
elongation are terrestrial remains such as roots, insects and leaves from laminated sediments<br />
of Lake Suigetsu in Japan (Kitagawa and Van der Plicht, 1998).<br />
23
Variations in the natural 14 C content cause fluctuations in the calibration curve. These<br />
fluctuations are called ‘wiggles’ (Van der Plicht, 2006) and although they complicate the<br />
calibration of 14 C dates, they can also be used to improve the accuracy of the radiocarbon<br />
dates. This is possible when several samples, with a known time period between them, are<br />
used. This method is called ‘wiggle-matching’ (Van der Plicht, 2006). This method is well<br />
suited when a piece of wood is not datable using dendrochronology. Because the time period<br />
between the rings is exactly known they from a chronology from which 14 C measurements can<br />
be obtained. Those measurements belong to a part of a calibration curve and are compared to<br />
the actual calibration curve to date the samples.<br />
Two different methods can be used; the conventional and the AMS-method (Van der Plicht,<br />
2006). The conventional method is based on radiometry and measures the radioactivity of the<br />
still present 14 C. For this method approximately 200 grams of carbon are needed and one<br />
measurement takes a few days. The AMS-method (Accelerator Mass Spectrometry) is based<br />
on mass spectrometry; elements of a sample are analyzed on their mass by giving them<br />
velocity, electrical charge and by separating them using a magnet. This method needs 1 mg of<br />
carbon, and measurements can be performed in less than one hour. Although the small sample<br />
fraction for the AMS-method is an advantage (it can date pollen, seeds and small<br />
macrofossils), the method also has some disadvantages; when using small samples the chance<br />
of contamination with other samples is larger and if the actual sample is not homogeneous, a<br />
small part of this sample might not be representative. Because of these disadvantages<br />
conventional measuring has the preference above the AMS-method when enough material is<br />
available (Van der Plicht, 2006).<br />
Sample preparation<br />
Besides the dendrochronological dating of wood, some wood and peat samples were dated<br />
using the radiocarbon method. From location Groeneweg, tree GFE02 was chosen, together<br />
with four samples from profile GP01. From location Lodijk four different pine samples were<br />
selected.<br />
Tree GFE02 was chosen because of<br />
its high number of tree rings, which<br />
enables dating the wood by<br />
‘wiggle-matching’. From this piece<br />
of wood every tenth ring was<br />
sampled, counting from the pith,<br />
ending with ring number 230<br />
(Figure 15) that represents<br />
sapwood. Exactly one ring per<br />
sample (with a minimum of 100<br />
milligrams) was cut out with a<br />
scalpel preventing the two<br />
surrounding rings being added to<br />
the sample.<br />
Figure 15: Subsamples of every tenth ring of tree GFE02 were radiocarbon dated.<br />
24
The pines did not contain enough tree rings for the ‘wiggle-matching’ method, but the AMSmethod<br />
could be performed. Because of the accuracy of the AMS-method, it was not<br />
necessary to measure exactly one ring. Counting rings from pine wood is more difficult than<br />
those of oaks because pines often form false rings. To select rings for radiocarbon dating one<br />
must be certain which rings were chosen. To make the dating more useful, some of the oldest<br />
rings were selected with a maximum of 5 rings. The rings of three different pine samples were<br />
subsampled for 14 C analysis:<br />
• LFE04: This sample contained at least 80 rings, but the oldest rings were difficult to<br />
measure and the exact amount was unclear. The subsample for 14 C dating contains<br />
ring 84 onwards with a maximum of 5 rings.<br />
• LFE07: This piece was characterized by a large part of reaction wood. Two different<br />
radii were measured, into two directions. Both of these radii correlate until ring 51,<br />
therefore rings 48 and 49 were chosen because these rings were obviously identifiable.<br />
• LFE08: Although this piece of wood contains 94 rings the outer rings were difficult to<br />
identify because they were very small and did not produce enough material for dating.<br />
Ring 70 till 74 were sampled because they were clearly visible.<br />
From profile GP01 four different samples were selected for radiocarbon dating. The first aim<br />
was to select macrofossils such as seeds, leaves or branches. But because most of the material<br />
was degraded, it was not possible to select complete macrofossils and three samples of bulk<br />
material had to be selected instead: sample 1 from the oldest part, sample 6 from the middle<br />
part and sample 12 from the youngest part. If roots were still recognizable, they were<br />
removed from the bulk material. Contamination of the bulk material with younger roots can<br />
cause 14 C ages to be too young (Yeloff, Bennett, Blaauw, Mauquoy, Sillasoo, Van der Plicht<br />
and Van Geel, 2006). Fortunately, in sample 12 also a branch was found, determined as<br />
Betula sp., and useful for radiocarbon dating.<br />
All samples were send to the Centre of Isotope Research at the university of Groningen<br />
(Figure 16) where they were dated using the AMS-method.<br />
Figure 16a: The Accelerator Mass Spectrometry<br />
measuring device.<br />
Figure 16b: The conventional measuring device.<br />
25
3. Results and Interpretation<br />
3.1 Dendrochronological results (Ring rapport number <strong>2007</strong>060)<br />
The two trees found in situ at location 1 (Figure 1) were both determined as oak (Quercus<br />
robur/petraea.). Unfortunately, synchronization of GFE01 with reference calendars did not<br />
give a significant result. Because both tree patterns have a significant comparison (α =<br />
0.0001), the dated tree GFE02 is used as a reference to date GFE01 (Table 1). The<br />
dendrochronological analysis showed that they were living in the same period: GFE01 from<br />
2719 until 2570 BC and GFE02 from 2882 until 2581 BC. The graphical comparison between<br />
GFE01 and GFE02 is shown in figure 17 and in more detail in figure 18, their common<br />
interval is 139 rings. The grey vertical areas represent the percentage of parallel variation<br />
(PV) (Chapter 2.3.1).<br />
Because GFE01 did not contain sapwood rings or wane edge, the exact date of death of the<br />
tree is unknown and can only be estimated. Date of death estimations (Jansma, <strong>2007</strong>) give 22<br />
± 7 sapwood rings for this tree. If the last heartwood ring of the sample would be measured,<br />
the date of death of the tree would be included in the interval –2570 + (22 ± 7). However, the<br />
number of heartwood rings missing towards the sapwood is unknown. Therefore it only can<br />
be said that the tree died after –2570 + (22 - 7) = - 2555<br />
Table 1: GFE01 is dated with GFE02 as a reference, both trees show a significant comparison.<br />
Tree Number Start End Date Sapwood THO GL Overlap Significance Reference<br />
of rings Year Year of estimation<br />
(α) chronology<br />
(BC) (BC) death<br />
(BC)<br />
after<br />
GFE01 150 2719 2570 2555 22 ± 7 6.461 69.780 139 0.0001 GFE02<br />
Figure 17: Graphical comparison of GFE01 (orange) and GFE02 (purple). The y-axis (logarithmic) shows the<br />
ring width (mm -2 ) and the x-axis gives the calendar year. The grey areas represent the PV.<br />
26
Figure 18: A detailed overview of the corresponding part of GFE01 (orange) and GFE02 (purple). The y-axis<br />
(logarithmic) shows the ring width (mm -2 ) and the x-axis gives the calendar year. The grey areas represent the<br />
PV.<br />
GFE02 is dated by synchronization with YPFKLSt (α = 0.0002) (Table 2). Table 3 shows the<br />
information of the used reference chronology. GFE02 contained 25 sapwood rings, but no<br />
wane edge. Calculations of the date of death (Jansma, <strong>2007</strong>) gave an sapwood estimation of 9<br />
± 8 rings. The last measured sapwood ring is from 2581 BC. The date of death of GFE02 can<br />
be calculated as: -2581 + (9 ± 8) = 2572 ± 8. This means that the date of death of GFE02 is<br />
after 2580 BC and because there are no heartwood rings missing, the tree died before 2564<br />
BC. The graphical synchronization between GFE02 and YPFKLSt is shown in figure 19, with<br />
the grey areas presenting the PV.<br />
Tree<br />
Table 2: GFE02 is dated with reference calendar YPFKLSt<br />
Number<br />
Number<br />
Start<br />
End<br />
Date of<br />
Sapwood<br />
of rings<br />
of<br />
Year<br />
Year<br />
death<br />
estimation<br />
sapwood (BC)<br />
(BC)<br />
(BC)<br />
rings<br />
GFE02 302 25 2882 2581 2580 /<br />
2564<br />
between<br />
THO<br />
GL<br />
Overlap<br />
Significance<br />
(α)<br />
Reference<br />
chronology<br />
9 ± 8 5.485 62.150 247 0.0002 YPFKLSt<br />
Table 3: Information about the used reference chronology<br />
Reference Country Province Number<br />
chronology<br />
YPFKLSt The Netherlands Zuid-<br />
Holland,<br />
Ypenburg<br />
Start Year End Year Author<br />
of rings (BC) (BC)<br />
247 2860 2614 E. Jansma, <strong>2007</strong>,<br />
based on data by<br />
Van Daalen<br />
(2001)<br />
27
Figure 19: Graphical comparison between GFE02 (green) and YPFKLSt (blue). The y-axis (logarithmic) shows<br />
the ring width (mm -2 ) and the x-axis gives the calendar year. The grey areas represent the PV.<br />
The tree-ring patterns of GFE01 and GFE02 were averaged into a site chronology (GFE1_2).<br />
The graphical synchronization of GFE1_2 (2882 – 2570 BC) with the reference calendar of<br />
YPFKLSt is shown in figure 20. The statistical results in table 4 shows that GFE1_2 and<br />
YPFKLSt synchronized together (α = 0.0002).<br />
Figure 20: Synchronization of the site chronology GFE1_2 (green) with the reference chronology YPFKLSt<br />
(blue). The y-axis (logarithmic) shows the ring width (mm -2 ) and the x-axis gives the calendar year. The grey<br />
areas represent the PV.<br />
Table 4: Results of the dendrochronological comparison of site chronology GFE1_2 with the reference<br />
chronology YPFKLSt.<br />
Site Number Start End Date of death THO GL Overlap Significance Reference<br />
Chronology of rings Year Year (BC) after<br />
(α) chronology<br />
(BC) (BC)<br />
GFE1_2 313 2882 2570 2555 5.992 62.550 247 0.0002 YPFKLSt<br />
28
The trees referred to as LFE01 and LFE02 from location 2<br />
(Figure 1) were also determined as oaks (Quercus<br />
robur/petraea). The trunks were characterized by the<br />
absence of branches; on both trees only one branch-growth<br />
was visible (Figure 21). Dendrochronological analysis<br />
indicates that these trees were living in the same period:<br />
LFE01 from 1678 until 1498 BC and LFE02 from 1688<br />
until 1476 BC. These dendrochronological dates are the<br />
result of the synchronization with reference calendar<br />
NLVEEN04 (α = 0.0001)<br />
Figure 21: The trunks showed<br />
one branch-growth.<br />
The ring pattern of LFE01 shows a normal growth pattern<br />
with juvenile wide rings at the beginning, followed by smaller rings due to the aging effect<br />
when the tree grows older. The ring pattern of LFE02 shows at the first beginning small ring<br />
widths where after a normal growth pattern appears, with the aging effect of small older rings.<br />
After 1582 BC the ring widths become remarkable smaller and subsequently it recovers from<br />
1566 BC onwards, showing wider rings.<br />
The ring widths of both trees were measured and compared to reference calendars. The best<br />
synchronization for both trees was found with the reference chronology NLVEEN04 (Table<br />
5). Table 6 shows the information of the reference chronology NLVEEN04. The samples did<br />
not contain sapwood or wane edge and the exact date of death of the trees is estimated.<br />
According to Jansma (<strong>2007</strong>) the sapwood estimation for sample LFE01 would be 25 ± 9<br />
rings. The number of heartwood rings missing towards the sapwood is unknown. Therefore it<br />
only can be said that the tree died after –1498 + (25-9) = -1482. Proceeding in the same<br />
manner, the sapwood estimations of LFE02 would be 28 ± 10. Again the number of<br />
heartwood rings till the first sapwood is unknown and it only can be stated that the tree died<br />
after –1476 + (28-10) = -1458.<br />
Table 5: Results of the dendrochronological comparison of LFE01 and LFE02 with the reference chronology<br />
NLVEEN04.<br />
Tree Number Start End Date of Sapwood THO GL Overlap Significance Reference<br />
of rings Year Year death (BC) estimations<br />
(α) chronology<br />
(BC) (BC) (after)<br />
LFE01 181 1678 1498 1482 25 ± 9 6.408 68.230 181 0.0001 NLVEEN04<br />
LFE02 213 1688 1476 1458 28 ± 10 6.887 65.730 213 0.0001 NLVEEN04<br />
Table 6: Information about the used reference chronology<br />
Reference Country Province Number Start Year End Year Author<br />
chronology<br />
of rings (BC) (BC)<br />
NLVEEN04 The Netherlands Zuid-<br />
Holland<br />
1601 2258 658 E. Jansma,<br />
1997<br />
29
Figure 22 shows the visual synchronization between LFE01 and LFE02. Their common<br />
interval is 181 rings, equal to the total amount of rings of LFE01. The graphical<br />
synchronization of LFE01 and LFE02 with the reference chronology NLVEEN04 is shown in<br />
figure 23 and 24 respectively. The site chronology LFE1_2 is the mean tree-ring pattern of<br />
LFE01 and LFE02. The results of the synchronization of LFE1_2 with NLVEEN04 are<br />
shown in table 7. The visual synchronization of this site chronology with the reference<br />
NLVEEN04 is shown in figure 25.<br />
Figure 22: The tree-ring pattern of LFE01 (orange) and LFE02 (purple) do synchronize together with a common<br />
interval of 181 rings. The y-axis (logarithmic) shows the ring width (mm -2 ) and the x-axis gives the calendar<br />
year. The grey areas represent the PV.<br />
Figure 23: Synchronization of the dated sample LFE01 (green) with the reference chronology NLVEEN04<br />
(blue). The y-axis (logarithmic) shows the ring width (mm -2 ) and the x-axis gives the calendar year. The grey<br />
areas represent the PV.<br />
30
Figure 24: Sample LFE02 (green) is also dated with the reference chronology NLVEEN04 (blue). The y-axis<br />
(logarithmic) shows the ring width (mm -2 ) and the x-axis gives the calendar year. The grey areas represent the<br />
PV.<br />
Table 7: Results of the dendrochronological comparison of site chronology LFE1_2 with the reference<br />
chronology NLVEEN04.<br />
Site<br />
Number<br />
Start<br />
End<br />
Date of THO<br />
GL<br />
Overlap<br />
Significance<br />
Reference<br />
Chronology<br />
of rings<br />
Year<br />
Year<br />
death<br />
(α)<br />
chronology<br />
(BC)<br />
(BC)<br />
(BC)<br />
(after)<br />
LFE 213 1688 1476 1458 7.671 66.430 213 0.0001 NLVEEN04<br />
1_2<br />
Figure 25: Synchronization of the site chronology LFE1_2 (green) with the reference chronology NLVEEN04<br />
(blue). The y-axis (logarithmic) shows the ring width (mm -2 ) and the x-axis gives the calendar year. The grey<br />
areas represent the PV.<br />
31
Besides LFE01 and LFE02 also some other small<br />
trunks were discovered but they were not useful for<br />
dendrochronology because the nearby root system<br />
deformed the tree rings. This root system was<br />
characterized by shallow roots (Figure 26) and the<br />
absence of deep taproots. Besides these trunks, many<br />
root systems were discovered in situ near the place<br />
where the oaks were found.<br />
Figure 26: A shallow root system.<br />
The small pieces of wood collected by the farmer on this land were determined as pine (Pinus<br />
sylvestris). Some of the collected pine trees showed remarkable growth of the stems. The<br />
bended forms (Figure 27) are also called ‘saber tooth growth’ (personal communication E.<br />
Jansma). Unfortunately none of these samples could be dated by dendrochronology. Table 8<br />
shows an overview of the pine trees which were measured.<br />
Table 8: Pine samples were<br />
measured on their tree rings. They<br />
could not be dated by means of<br />
dendrochronology.<br />
Tree Number Date<br />
of rings<br />
LFE04 97 Undated<br />
LFE07 96 Undated<br />
LFE08 94 Undated<br />
Figure 27: ‘Saber tooth growth’<br />
32
Pointer years<br />
The tree-ring patterns of the samples GFE01 and GFE02 were compared graphically with a<br />
pointer year graph obtained from 8 trees that lived during the same period but in a different<br />
area (Ypenburg, province of Zuid-Holland). These 8 trees were selected because they had a<br />
significant correlation with the oaks from location 1. The results are plotted in figure 28 and<br />
29 respectively. The pointer years require a minimum sample number (Appendix 4, Table 14).<br />
The results are given on four different significance levels (Appendix 4, Table15). The bars on<br />
the positive part of the graphic show the positive pointer years (parallel increased growth),<br />
whereas the negative part of the graphic shows the negative pointer years (parallel decreased<br />
growth). Table 16 (Appendix 4) shows which tree measurements from Ypenburg are<br />
included.<br />
99%<br />
95%<br />
90%<br />
99.9%<br />
99%<br />
95%<br />
90%<br />
90%<br />
95%<br />
99%<br />
99.9%<br />
Figure 28: The ring width pattern of GFE01 compared to pointer years of 8 trees. The bars show the pointer<br />
years on different significance levels (90, 95, 99 and 99.9 %). The shortest bars show the lowest significance<br />
level and the longest bars show the highest significance level.<br />
99.9%<br />
99%<br />
95%<br />
90%<br />
90%<br />
95%<br />
99%<br />
99.9%<br />
Figure 29: The ring width pattern of GFE02 compared to pointer years of 8 trees. The bars show the pointer<br />
years on different significance levels (90, 95, 99 and 99.9 %). The shortest bars show the lowest significance<br />
level and the longest bars show the highest significance level.<br />
The pointer years of 8 oaks from Ypenburg and two oaks from the Eem Polder show<br />
similarities but it is obvious that the statistical results of the comparisons are not highly<br />
significant. This is caused by the small amount of trees that is used for the analysis and not by<br />
an less significant event. During the growth of GFE01 only three pointer years were showed<br />
33
y the 8 trees from Ypenburg. This is not a statistical result, but the result of a small overlap<br />
between the tree-ring pattern of GFE01 and of the trees from Ypenburg.<br />
The tree-ring patterns of the samples LFE01 and LFE02 were compared graphically with a<br />
pointer year graph obtained from 32 trees from different locations in the province of Zuid-<br />
Holland. The tree-ring patterns of LFE01 and LFE02 are plotted with the pointer years of<br />
these 32 trees (Appendix 4, table 17) in figure 30 and 31 respectively. Table 18 (Appendix 4)<br />
shows the original location of these trees and table 19 (Appendix 4) shows which tree<br />
measurements are included.<br />
99.9%<br />
99%<br />
95%<br />
90%<br />
90%<br />
95%<br />
99%<br />
99.9%<br />
Figure 30: The ring width pattern of LFE01 (pink curve) compared to pointer years (bars) of 32 trees. The bars<br />
show the pointer years on different significance levels (90, 95, 99 and 99.9 %). The shortest bars show the lowest<br />
significance level and the longest bars show the highest significance level.<br />
99.9%<br />
99%<br />
95%<br />
90%<br />
90%<br />
95%<br />
99%<br />
99.9%<br />
Figure 31: The ring width pattern of LFE02 (pink curve) compared to pointer years (bars) of 32 trees. The bars<br />
show the pointer years on different significance levels (90, 95, 99 and 99.9 %). The shortest bars show the lowest<br />
significance level and the longest bars show the highest significance level.<br />
34
Table 9: Pointer years with a 99.9% significance level,<br />
three positive and four negative pointer years.<br />
Pointer years with a 99.9% significance level based on 32<br />
trees (province of Zuid-Holland)<br />
Positive Pointer Year (BC) Negative Pointer Year (BC)<br />
1660 1629<br />
1605 1628<br />
1582 1602<br />
1594<br />
Table 9 summarizes the positive and negative pointer years of the 32 trees from the province<br />
of Zuid-Holland on a 99.9% significance level. When only concentrating on these pointer<br />
years, the tree pattern of LFE01 shows some similarities. The first positive pointer year (1660<br />
BC) coincidence with an remarkable wider ring of LFE01. Also the positive pointer years<br />
from 1605 and 1582 BC show a sudden increase in ring width. The negative pointer years of<br />
1629 and 1628 are remarkable because they present two following years. LFE01 also shows<br />
smaller rings during both years. The negative pointer year from 1602 BC coincidence with a<br />
sudden decrease in ring width and although less obvious, also the ring width formed in 1594<br />
BC decreases. The tree-ring pattern of LFE01 also show some notable ring width changes that<br />
are not accomplished with pointer years of the 32 trees like a remarkable smaller ring at 1658<br />
BC and a remarkable wider ring at 1655 BC.<br />
The tree-ring pattern of LFE02 also shows many similarities with the pointer years of the 32<br />
trees from the province of Zuid-Holland. The positive pointer year from 1660 BC<br />
coincidences with a small ring width increase whereas the other pointer years are more similar<br />
with ring width changes of LFE02. Obviously, the tree-ring pattern of LFE02 shows a trend<br />
of smaller rings between 1580 and 1550 BC, only corresponding with one negative pointer<br />
year at 1569 BC.<br />
35
3.2 Pollen analysis<br />
The pollen analyses performed on the samples from Groeneweg (GP01 and GP02) resulted in<br />
two pollen diagrams (Figure 34a/b). The determined species are divided into several groups<br />
(upland trees, wetland trees, anthropogenic indicators, herbs, spores, scalariform perforation<br />
plate of wood, types and indeterminable pollen). Both figures show the lithology column of<br />
the profile with sand on the lower part and peat representing the major part of the profile.<br />
Some samples from profile GP01 were radiocarbon dated and these dates are visible in the<br />
figure. The radiocarbon dates were calibrated to calendar years, also plotted in the figure.<br />
The diagrams present percentages of the total amount of determined species, plotted against<br />
depth below NAP (Dutch Ordnance Datum). Two different graph styles are plotted in the<br />
diagram; the dark-grey colour shows the actual percentages of the graph, whereas the hatched<br />
graphs are 10 times exaggerated to make the presence of species visible. The first graph<br />
shows a cumulative diagram (percentages) of the upland trees, wetland trees and herbs. Pollen<br />
sums present the absolute total of upland trees, wetland trees and herbs. Rare types which<br />
were added to the diagram after scanning the samples are presented by the ‘+’ and ‘++’<br />
symbols. Appendix 5 shows the absolute concentration of upland trees, wetland trees and<br />
herbs together. This absolute curve is calculated by making use of the added and counted<br />
Lycopodium spores per sample relative to the counted pollen. The diagrams show no obvious<br />
vegetation changes, i.e. a sudden increase, decrease, appearance or disappearance of species,<br />
therefore the diagram is not divided into separate pollen zones.<br />
Description of pollen diagram GP01<br />
Tree species (both upland and wetland) represent the most dominant group in this diagram.<br />
The most dominant trees within the upland trees are hazel (Corylus avellana) and oak<br />
(Quercus sp.), with a minor dominance of birch (Betula). These three species have a rather<br />
constant appearance in the vegetation, although at the upper part hazel increases and oak<br />
decreases slightly. In the lowest part, lime (Tilia) has its highest appearance and after its<br />
amount decreases lime is present at a constant low amount. The same pattern applies for elm<br />
(Ulmus), although this species has even lower quantities. Pine (Pinus sylvestris) is constantly<br />
present with low numbers. Ash (Fraxinus excelsior), ivy (Hedera helix), holly (Ilex aquifolia)<br />
and mistletoe (Viscum album) all appear as rare types at the lower part and the first three<br />
appear along the profile with small quantities or as rare types. Pollen grains of ivy, holly and<br />
mistletoe are shown in figure 32.<br />
Figure 32a: Ivy<br />
(Hedera helix)<br />
Figure 32b: Holly (Ilex<br />
aquifolia). The left<br />
picture shows the polar<br />
view<br />
Figure 32c:<br />
Mistletoe<br />
(Viscum album)<br />
Woodbine (Lonicera periclymenum) is counted in the lowest sample and appears furthermore<br />
as a rare type several times. Finally, in the upper part a pollen grain of beech (Fagus<br />
sylvatica) was counted. Alder (Alnus sp.) is by far the most dominant wetland tree with a<br />
constant occurrence through the column. From alder and oak also some clumps of (immature)<br />
pollen were found (Figure 33). Willow (Salix) appears a couple of times.<br />
36
Herbs (including anthropogenic indicators) play a minor role in this vegetation. The most<br />
common herbs are heath (Ericales) and grasses. Spores were also determined in low<br />
quantities but present in all samples. Scalariform perforation plates appear in the middle and<br />
upper part. They are most probably from alder, but birch is also possible. One type of fungi<br />
(type 2, Gelasinospora) was found in this profile.<br />
Description of pollen diagram GP02<br />
Trees (upland and wetland) are also the most dominant group of<br />
species in profile GP02, with a higher abundance of upland trees in<br />
comparison with profile GP01. Hazel and especially oak are the most<br />
dominant upland species and they show an opposite appearance.<br />
Again some clumps of oak and alder pollen grains (Figure 33) were<br />
found. Birch is constantly present with a large increase in the upper<br />
part. Pine, lime and elm are constantly present, but in low quantities.<br />
Whereas beech occurred in GP01 only in the upper part, in GP02 it<br />
appears only in the lowest part. Ash, ivy and holly appear in low<br />
quantities along the profile, sometimes noticed as rare types.<br />
Figure 33: clumps of<br />
oak pollen.<br />
Woodbine and hornbeam (Carpinus) both were observed once as rare types. Alder represents<br />
again the most dominant wetland tree, with some occurrence of willow and bog myrtle<br />
(Myrica gale). The herbs (including anthropogenic indicators) are the minor group in this<br />
diagram, with again heath and grasses being the most frequent species. Some spores were<br />
determined, but they appear in low quantities along the profile. Scalariform perforation plates<br />
appear only in the upper part. Three different fungi types were determined, with a striking<br />
abundance of type 112 (Cercophora) in the highest part of the profile.<br />
37
Salix<br />
Ilex<br />
Asteraceae liguflorae<br />
Asteraceae tubuliflorae<br />
Alisma<br />
Apiaceae<br />
Caryophyllaceae<br />
Cyperaceae<br />
Ericales<br />
Vaccinium-type<br />
Galium-type<br />
Humulus lupulus<br />
Hydrocotyle<br />
Poaceae<br />
Potentilla-type<br />
Solanum dulcamara<br />
Sparganium erectum-type<br />
Typha latifolia<br />
Dryopteris-type<br />
Polypodium-type<br />
Pteridium-type<br />
Sphagnum<br />
Scalariform perforation plate<br />
v. Geel type 2<br />
Indet<br />
Indet. corroded<br />
Indet. folded<br />
Pollensum<br />
Artemisia<br />
Chenopodiaceae<br />
Plantago lanceolata<br />
Succisa<br />
Viscum album<br />
Alnus<br />
Tilia<br />
Tilia platyphyllos<br />
Ulmus<br />
Lonicera peridymenum<br />
Quercus<br />
Pinus<br />
Upland trees<br />
Wetland trees<br />
Herbs<br />
Betula<br />
Corylus<br />
Fagus<br />
Fraxinus<br />
Hedera<br />
Depth below NAP (mm)<br />
3875 ± 40 2350 BC<br />
2400 BC<br />
2600 BC<br />
4125 ± 35 2725 BC<br />
2800 BC<br />
3000 BC<br />
4475 ± 40 3185 BC<br />
1300<br />
1350<br />
1400<br />
1450<br />
1500<br />
1550<br />
1600<br />
1650<br />
1700<br />
1750<br />
1800<br />
Lithology<br />
Eemv alley, GFE01<br />
20 40 60 80 100<br />
20<br />
Humus Sand<br />
Figure 34a: pollen diagram of GP01<br />
20 40<br />
Upland Trees<br />
20 40<br />
Wetland Trees<br />
Anthropogenic indicators Herbs<br />
Spores Perf. Plate Types Indet.<br />
368<br />
352<br />
432<br />
350<br />
1360<br />
415<br />
Analysis: <strong>Marije</strong> <strong>Vlaar</strong>, <strong>2007</strong><br />
Radiocarbon Dates (BP)<br />
Ca lendar age
Myrica<br />
Salix<br />
Analysis: <strong>Marije</strong> <strong>Vlaar</strong>, <strong>2007</strong><br />
Plantago lanceolata<br />
Apiaceae<br />
Calystegia<br />
Cyperaceae<br />
Ericales<br />
Vaccinium-type<br />
Filipendula<br />
Galium-type<br />
Humulus lupulus<br />
Hydrocotyle<br />
Nuphar lutea<br />
Poaceae<br />
Potentilla-type<br />
Ranunculus<br />
c.f. Scr ophulariac eae<br />
Solanum dulcamara<br />
Sparganium erectum-type<br />
Typha latifolia<br />
Dryopteris-type<br />
Polypodium-type<br />
Pteridium-type<br />
Sphagnum<br />
Scalariform perforation pl ate<br />
Type 2 (Gelasinospora)<br />
Type 112 (Cercophora)<br />
Type 368 (Podospora)<br />
Indet<br />
Indet. folded<br />
Indet. corroded<br />
Pollensum<br />
1150<br />
1200<br />
1250<br />
1300<br />
1350<br />
1400<br />
1450<br />
1500<br />
1550<br />
1600<br />
20 40 60 80 100<br />
39<br />
Artemisia<br />
Asteraceae tubuliflorae<br />
Chenopodiaceae<br />
Eu-Rumex<br />
Alnus<br />
Tilia<br />
Ulmus<br />
Eemvalley, GFE02<br />
Upland trees<br />
Wetland trees<br />
Herbs<br />
Betula<br />
Carpinus<br />
Corylus<br />
Fagus<br />
Fraxinus<br />
Hedera<br />
Ilex<br />
Lonicera peridymenum<br />
Pinus<br />
20<br />
20<br />
Upland Trees<br />
20 40 60 80<br />
20<br />
Wetland Trees<br />
Anthropogenic indicators Herbs<br />
Spores Perf. Plate Types<br />
Indet.<br />
371<br />
368<br />
380<br />
619<br />
400<br />
393<br />
375<br />
20<br />
20<br />
20<br />
Quercus<br />
Humus Sand<br />
Dept h below NAP (mm)<br />
Lithology<br />
Figure 34b: Pollen diagram of GP02
3.3 Macrofossils<br />
The macrofossil analysis was carried out on three different subsamples from profile GP01.<br />
The selected subsamples were derived from the upper, middle and lower part of the profile<br />
(sample 12, 6 and 1 respectively). The estimation of the abundance of macrofossils per<br />
sample shows that the vegetation remains were not preserved well and that most of them were<br />
degraded. Some stem fragments were determined as common reed (Phragmites australis)<br />
although the determination was not described with certainty. Besides the macrofossil study on<br />
samples from the profile, also some pine wood was found in the field above the samples<br />
profile. The fossil content of the youngest sample (12) is the highest, containing also charcoal<br />
and insect remains. Table 10 shows the results of the macrofossil analysis of the three<br />
samples.<br />
Table 10: analysis of macrofossils<br />
Sample Volume Content Estimated number Determined<br />
number<br />
1<br />
(cm<br />
)<br />
39 Small root fragments 11-50<br />
Wood (root) 1 Betula sp.<br />
Small stem fragments 1-10<br />
Large stem fragments 1-10 Phragmites<br />
Cenococcum (fungi) 51-100<br />
6 41 Small root fragments 11-50<br />
Wood (root fragments) 1-10 Betula sp.<br />
Small stem fragments 1-10<br />
Small branch fragments 1-10<br />
Cenococcum (fungi) 1-10<br />
12 33 Small root fragments 1-10<br />
Small stem fragments 51-100<br />
Large stem fragments 1-10 Phragmites<br />
Wood (branch) 1-10<br />
Betula sp.<br />
Nodes 1-10<br />
Bark fragments 1-10<br />
Seeds 40 Betula sp.<br />
Beetle fragments 1-10<br />
Charcoal 1<br />
Cenococcum (fungi) 1-10
3.4 Radiocarbon dating<br />
In a first attempt both oaks GFE01 and GFE02 could not be dated by dendrochronology and<br />
therefore wood samples were sent for radiocarbon analysis. Both trees synchronized together<br />
(Figure 17) and by dating one of the trees, the age of the other tree can be assessed. From the<br />
23 samples of GFE02 which were sent to the Centre of Isotope Research, 6 were selected for<br />
AMS-dating. The dates were given as dates Before Present (BP) and by wiggle-matching<br />
these dates were transformed to dates Before Christ (BC). The dates (BP) were given with a<br />
40 year threshold for accuracy. Two different options were given as a result of the wigglematching<br />
(Figure 35a/b), depending on the goodness of fit (Chi 2 fitvalue).<br />
Figure 35a: Wiggle-match result, option 1.<br />
Chi 2 fitvalue: 0.895<br />
Figure 35b: Wiggle-match result, option 2.<br />
Chi 2 fitvalue: 0.980<br />
Table 11: The two wiggle-match options results in two different dates<br />
Ring number GrAnumber<br />
Age (BP) Calibrated calendar Age (BC) Calibrated calendar Age (BC)<br />
Wiggle-match 1, figure 34a Wiggle-match 2, figure 34b<br />
10 35608 4090 ± 40 2734 2815<br />
50 35609 4165 ± 40 2684 2765<br />
100 35610 4160 ± 40 2634 2715<br />
150 35613 4045 ± 40 2584 2665<br />
200 35616 4050 ± 40 2534 2615<br />
230 35615 4040 ± 40 2504 2585<br />
The two different options of the wiggle-match also resulted in two different dating for the tree<br />
as shown in table 11. According to the wiggle-match with the lowest Chi 2 fitvalue (option 1)<br />
the oldest measured ring (number 10) was formed in 2734 BC and the youngest measured<br />
ring (number 230) was formed in 2504 BC. The other wiggle-match result (option 2) showed<br />
that ring number 10 was formed in 2815 BC and ring number 230 in 2585 BC.<br />
The dendrochronological date of the last measured ring of GFE02 (ring number 302) showed<br />
that it was formed in 2581 BC. According to the two wiggle-match results ring 302 was<br />
formed in 2432 BC (option 1) or in 2513 BC (option 2).<br />
41
From profile GP01 four samples were dated by the AMS- radiocarbon dating method. The<br />
results show that the profile was formed between 4475 ± 40 and 3875 ± 40 year BP. The<br />
radiocarbon dates (BP) were calibrated to calendar years (BC) (Wincal25, Centre of Isotope<br />
Research at the University of Groningen). First a 2σ range was calculated from which the<br />
average age (cal BC) of the samples was calculated (Table 12). The depth-age relation of this<br />
profile is given in figure 36. This rather linear relation (0.6 mm/ year) is used to assume a<br />
constant accumulation rate of sediments through time and to deduce calendar ages for the<br />
whole profile.<br />
Table 12: 14 C dates of profile GP01<br />
Sample GrAnumber<br />
Age (BP) 2σ range (cal Average age (cal Material<br />
BC)<br />
BC)<br />
GP01-1 35440 4475 ± 40 3026-3342 3185 Peat<br />
GP01-6 35441 4125 ± 35 2579-2871 2725 Peat<br />
GP01-12 35444 3875 ± 40 2208-2476 2350 Peat<br />
GP01-12 35445 3855 ± 40<br />
Branch (Betula sp.)<br />
Age (cal. BC)<br />
-2200<br />
-2400<br />
-2350<br />
-2600<br />
-2725<br />
-2800<br />
-3000<br />
-3185<br />
-3200<br />
1250 1305 1360 1415 1470 1525 1580 1635 1690 1745 1800<br />
Depth below NAP (mm)<br />
Figure 36: Depth-age model of the analyzed profile.<br />
42
4. Discussion<br />
Dendrochronology<br />
The oaks found at the Lodijk (LFE01, 1678-1498 BC, date of death after 1482 BC and<br />
LFE02, 1688-1476 BC, date of death after 1458 BC) lived during a more recent period than<br />
the oaks found near the Groeneweg (GFE01, 2719 – 2570 BC, date of death after 2555 BC<br />
and GFE02, 2882 – 2581 BC, date of death between 2580 and 2564 BC). The shape of the<br />
trunks and the tree-ring patterns of the trees contain information about their growth<br />
conditions. The branchless character of LFE01 and LFE02 gives information about their<br />
direct surrounding: diminished branching and enhanced stem elongation is typical of crowded<br />
trees to avoid shading (Weinig, 2000), which indicates that both trees were growing in a<br />
closed forest with competition for light and space. The tree-ring pattern of LFE01 did not<br />
show any remarkable growth conditions, but its pattern displays the aging effect<br />
(Schweingruber, 1996). The tree-ring widths of LFE02 show some remarkable patterns which<br />
can be caused by stressed growth conditions. The stress periods can be due to different factors<br />
such as temperature changes, soil characteristics and fungi or viruses (Pedersen, 1998).<br />
The discovery of several root systems near the location where LFE01 and LFE02 were found<br />
might indicate that more trees were growing during this period (ca. 1700-1450 BC). The<br />
absence of taproots on some pine trunks indicates a high groundwater level. The development<br />
of deep taproots versus wide spreading roots depends on the depth of the groundwater table<br />
(Callaway, 1990). This implies that these trees were growing on a saturated soil where<br />
taproots cannot develop. The branchless trunks and the abundant root-systems in the field<br />
indicate that the trees possibly were once part of a closed forest. The preservation conditions<br />
must be perfect to conserve these trees such a long time. Aerobic conditions speed up the<br />
decomposition process (Kristensen, Ahmed and Devol, 1995). These trees were preserved in<br />
oxygen-free conditions and therefore hardly degraded.<br />
The dendrochronological analysis of GFE01 and GFE02 was complicated because the trees<br />
did not date with any of the existing reference calendars. The reference calendars used to date<br />
Dutch oaks are compiled from measurements on oaks that were mostly growing on humid<br />
soils such as bogs. These oaks are also known as bog oaks and their ring patterns often<br />
contain obvious growth depressions lasting for a decade or more (personal communication E.<br />
Jansma). The ring-width patterns of GFE01 and GFE02 show a rather constant growth,<br />
without remarkable growth reductions. The fact that the oaks from the Groeneweg could not<br />
be dated with the existing references calendars based on bog oaks could indicate that the oaks<br />
from location 1 were growing on less saturated substrates. The local growth conditions of the<br />
oaks from the Groeneweg were different from the environmental factors that steered the<br />
growth of the oaks that served as building blocks of the existing ‘wet’ reference calendars.<br />
In a second attempt to date the undated trees from the Groeneweg, a new chronology was<br />
compiled, using data from previous research on an ancient forest in Ypenburg (62 oaks<br />
growing from 3000 until 2500 BC; Van Daalen, 2001). These eight trees were chosen because<br />
they gave a significant date with the oaks from location 1. If more measurement series of oaks<br />
from Ypenburg were inserted into the calendar, the Ypenburg material became dominant in<br />
the calendar and the match with the series from location 1 decreased. The obvious reason is<br />
that the local growth at Ypenburg was not identical to those at location 1. The dated oaks<br />
from Ypenburg presumably grew on dunes and although peat development occurred on a<br />
large scale during that time, it did not reach these topographic highs (Van Daalen, 2001). The<br />
fact that the two oaks from location 1 could not be dated with bog oak chronologies, but<br />
43
synchronize well with the Ypenburg oaks that grew on less saturated soils, indicates that the<br />
oaks from location 1 also grew on unsaturated soils. Although the growth environment of the<br />
oaks clearly cannot be characterized as a bog, the trunks were preserved in peat. Both buried<br />
trees were very large and GFE02 also contained easily degradable sapwood rings, which<br />
means that the burying of the trunks and the development of the peat must have been taken<br />
place before the trees died or soon after their death. The absence of plant species that<br />
characterize peat development is in contrast with the peat layers in which the trunks were<br />
found. The trees were probably growing on a very humus-rich soil and somehow the<br />
environment changed and increased the conservation of the trunks. The tree-ring patterns do<br />
not show any sign of coming death (illness or death by old age, expressed in narrow ring<br />
widths), so the trees probably died because of a sudden external event. The results of this<br />
study do not shed any light on the chronological relationship between the start of the peat<br />
development on this location on the one hand, and the dying-off of the trees and their<br />
subsequent sedimentation in the soil on the other hand.<br />
Compared to dendrochronology the radiocarbon dates of the wood, which were performed<br />
before the oaks could be dated with the reference calendar from Ypenburg, give a different<br />
date for the oaks. Two different wiggle-match options show a difference of 149 (± 40) years<br />
with the dendrochronological date and the other option gives a difference of 68 (± 40) years.<br />
It is difficult to explain this large difference between the radiocarbon and dendrochronological<br />
dates. The statistical results of the dendrochronological dates are very significant and<br />
therefore there are no doubts about their accuracy. The software used for wiggle-matching<br />
lacks proper statistical foundation. It does not make clear if the match fits within the<br />
dendrochronological dating. Probably radiocarbon measuring with a higher resolution, e.g.<br />
every tenth ring instead of every fiftieth ring, will increase the statistical certainty about the<br />
measurement and may even change the wiggle-match results.<br />
Regarding the dating of pinewood, the main problem is the absence of good reference<br />
calendars in the Netherlands. A problem with the pines is their relative young age, which<br />
results in measurement series containing too few rings for statistical comparisons and cross<br />
dating. The remarkable stem form (‘saber-tooth growth’) of some pines is the result of<br />
stabilization after subsiding of the unstable substrate. Because the pines could not be dated by<br />
dendrochronology, samples were sent to the Centre of Isotope Research at the University of<br />
Groningen for radiocarbon dating. Unfortunately the results are not yet available due to the<br />
busy schedule of this research centre. The ecological past of the pines and their relation to the<br />
oaks from location 2 are still unknown.<br />
Pointer Years<br />
If tree-ring patterns from a certain location follow the pattern of pointer years from trees<br />
growing elsewhere, it may imply that the same environmental factors have influenced the<br />
vegetation on a large scale. It is difficult to explain why the pointer years of oaks growing in<br />
the Eem Valley correspond with the pointer years of oaks growing in the province of Zuid-<br />
Holland. There is little known about environmental or climatic changes in the period during<br />
which these trees were living. However, some discussion exists about the pointer years from<br />
1629 and 1628 BC. There is some evidence that a volcanic event influenced at least the<br />
northern hemisphere during this time. Ice cores from Greenland, bristlecone frost rings from<br />
North-America and Irish oak minimum-growth events show that there was a large volcanic<br />
dust-veil event in 1628 BC (Baillie, 1989). Some scientists suggest that these phenomena are<br />
related to the dating of the Thera (Santorini) eruption (Baillie, 1991). Thera is a volcanic<br />
island in the Mediterranean and belongs to Greece. However, none of the evidences specifies<br />
44
which volcanic eruptions were actually involved and therefore it cannot be proven that the<br />
1628 BC event was connected to the (undated) Thera eruption (Baillie, 1989). One of the<br />
problems to be solved is the link between a volcanic dust event and reduced oak growth in<br />
northern Europe. These oaks are hardly affected by cooling or low pressure anomalies, both a<br />
result of major eruptions, but it is possible that the oaks were responding to raised water<br />
tables (Baillie, 1989). Because of a time lag between volcanic eruptions and the resulting<br />
environmental changes elsewhere, it may be possible that the eruptions took place a couple of<br />
years before the occurrence of growth responses in European oaks (Baillie, 1989). Volcanic<br />
eruptions in the past may have caused large-scale environmental changes on which trees<br />
throughout the world reacted. The main problem is that it is unknown what kind of eruption<br />
influenced the 1628 BC event and if this was really the Thera-eruption.<br />
At 1629 and 1628 BC the 32 trees from Zuid-Holland and the two oaks from location 2 (Eem<br />
Polder) show negative pointer years. If the volcanic eruption from 1628 BC has influenced<br />
the growth of these trees, it is surprising that the preceding year (1629 BC) already shows a<br />
negative pointer year. This pointer year can not be related to an event that happened a year<br />
later. In addition, one would expect the reduced growth to last after the volcanic event of 1628<br />
BC, because of lagged environmental response and the fact that tree growth would need some<br />
time to recover. This lagged response does not show up in the material. There must have been<br />
a general environmental impact (or combination of factors) that for two years negatively<br />
influenced the growth of the 32 trees from Zuid-Holland and the two oaks from the Eem<br />
Polder, but the reason for this reduced growth can not be rashly ascribed to the Thera eruption<br />
of 1628 BC.<br />
Pollen and Macrofossil analysis<br />
To reconstruct the appearance of the surrounding vegetation near the oaks which were<br />
discovered in situ, pollen and macrofossil analysis were performed. General pollen zones of<br />
the Netherlands are used to relate the vegetation history to time. The exclusion of local<br />
circumstances within these overviews should be considered when these overviews are used as<br />
a dating method. During the pollen analysis of this study several species were noticed as a<br />
time indicator. According to chronostratigraphic pollen zones of (Bastiaens, Brinkkemper,<br />
Deforce, Maes, Rovekamp, Van den Bremt and Zwaenepoel, 2006) the pollen samples from<br />
the Eem Polder presented a vegetation development from the Early-Atlantic period (7000-<br />
5500 yr BC). The division of pollen zones from (Bastiaens et al, 2006) are based on a<br />
previous developed overview from (Berendsen and Zagwijn, 1984). The addition of new<br />
species as time indicators to the overview of (Berendsen and Zagwijn, 1984) are based on<br />
more recent research on pollen, but also on the analysis of seeds and wood remains (Bastiaens<br />
et al, 2006).<br />
According to the pollen zones of (Bastiaens et al, 2006) oak and alder are the most dominant<br />
species during the Early-Atlantic. Lime and elm are also abundant and ivy, mistletoe and<br />
holly are present. Furthermore, pine trees are present but in low quantities (Bastiaens et al,<br />
2006). The transition to the Late-Atlanticum is characterized by the first appearance of<br />
cultivated species. During the next stage, the Subboreal, European yew (Taxus baccata)<br />
appears, lime and elm are no longer abundant, hazel increases in appearance and cultivated<br />
plants are common. The Subboreal ends with the appearance of beech (Bastiaens et al, 2006).<br />
Another general pollen zone overview is presented by de Jong (1983). This overview<br />
describes the pollen zones with less detail. The Atlanticum is characterized by the occurrence<br />
of more than 5% elm, alder and oak are dominant and pine appears in low quantities. The<br />
45
transition to the Subboreal is defined by a decrease of elm and the appearance of cultivated<br />
species. The first occurrence of beech represents the end of this period.<br />
Because of some important indicators such as ivy, mistletoe and holly and because of the<br />
absence of cultivated plants and European yew, the results of the pollen analysis of this study<br />
were interpreted as pointing to the Early-Atlanticum (7000-5500 yr BC). However, the<br />
radiocarbon analysis of three soil samples from the profile points out that the profile dates<br />
from 3185 to 2350 yr BC (Subboreal). Although it might be possible that the radiocarbon<br />
dates are too young due to the presence of roots in the bulk material, after these results the<br />
date based on the general pollen zones was regarded with more caution. According to<br />
(Bastiaens et al, 2006) the presence of European yew and cultivated plants is an important<br />
characteristic of the Subboreal. Both indicators were not noticed in the studied pollen samples<br />
but their absence can be the result of the local influences. European yew pollen is a difficult<br />
species for analysis because it is sensitive to folding and fracturing. Because this species does<br />
not appear in high quantities, local factors may influence its appearance. A detailed palaeoecological<br />
study on micro- and macrofossils from the Borchert (North-Eastern part of The<br />
Netherlands) resulted in a pollen diagram dated from 10.000 to 1300 yr BP (Van Geel,<br />
Bohncke and Dee, 1980). In this diagram European yew is not present during the whole<br />
period. Another palaeo-ecological study on sediments from Ilperveld (Western Netherlands)<br />
resulted in a pollen diagram dated from about 3000 to 1000 yr BC (Bakker and Van<br />
Smeerdijk, 1981). Again European yew is not found during the whole period. This may<br />
indicate that this species should not be put out as an important time indicator. There are also<br />
examples in the Netherlands of pollen analysis that contain pollen of European yew. A pollen<br />
diagram from the ‘Dommelvalley’ (North Brabant) contains European yew pollen during the<br />
Subboreal period (Janssen, 1972). They appeared during clearance phases, at a time when the<br />
forests were less dense (Janssen, 1972). Other time indicators according to (Bastiaens et al,<br />
2006) for the Early-Atlanticum are ivy, mistletoe and holly which are not classified at the<br />
Late-Atlanticum and the Subboreal. However, the palaeo-ecological study from the Borchert<br />
(Van Geel et al, 1980) showed that these species continued their appearance through these<br />
periods and therefore could also be present in the studied samples from the Eem Polder. The<br />
appearance of elm in the profile from the Eem Polder agrees with the Subboreal period of de<br />
Jong (1983). According to his pollen zones elm should be present with more than 5% during<br />
the Atlanticum and with low quantities during the Subboreal. In both profiles from the Eem<br />
Polder elm appears with low quantities, indicating the Subboreal period. Finally, according to<br />
both pollen zones of (Bastiaens et al, 2006) and de Jong (1983) the Subboreal is characterized<br />
by the presence of cultivated plants. The profiles from the Eem Polder show some<br />
anthropogenic indicators, but in very low quantities. This could also be caused by local<br />
factors such as the absence of people in this area during this period. The discussed time<br />
indicators and the comments on the general pollen zones assume that the three radiocarbon<br />
dates of the profile are correct.<br />
Both pollen diagrams show a very high amount of trees (upland and wetland) relative to the<br />
herbs. This implies that the abundant trees formed a closed forest in which herbs could not<br />
develop because of shading. The most dominant tree species in this forest were hazel, oak and<br />
alder. Pollen dispersal in trees is a combination of local dispersion and long-distance transport<br />
(Streiff, Ducousso, Lexer, Steinkellner, Gloessl and Kremer, 1999). The clumps of<br />
(immature) pollen constitute the gravity component and indicate local occurrence of the<br />
corresponding taxa. Clumps of pollen seem to be associated with local or extra-local pollen<br />
deposition (Janssen, 1986) implying local appearance of oak and alder in this site. The pollen<br />
analysis from the ‘Dommelvalley’ (Janssen, 1972) shows that European yew only appears<br />
during clearance phases when the forest was less dense. The density of the forest from the<br />
46
Eem Polder might explain why European yew did not appear in this vegetation. Besides very<br />
abundant trees as hazel, oak and alder, also pine appears constantly in this pollen diagram,<br />
although with very low quantities. Pollen of pine trees is susceptible to long-distance dispersal<br />
and due to the high pollen production of pines it is difficult to determine the presence of pine<br />
on a local scale (Eide, Birks, Bigelow, Peglar, and Birks, 2006). Hence, plant macrofossils of<br />
pines are useful to refine and correct the pollen analysis of pine trees (Birks and Birks, 2000).<br />
Because the pollen percentages of pines are very low in both diagrams, pine trees were not<br />
abundant in the forest. However, remarkable pine macrofossils were found above the sampled<br />
profile. This means that pine trees were locally present and probably their dominance<br />
increased in time. Only three samples from profile GP01 were analyzed. The macrofossil<br />
remains determined as birch indicate that this tree was locally present. The older samples (1<br />
and 6) only contain birch roots which could be grown in after the peat formation and are no<br />
proof the local appearance of this species. If macrofossil remains of common reed occur in the<br />
profile it is difficult to explain why this species is not noticed during the pollen analysis. This<br />
species was also found above the sampled profile and could indicate the formation of peat and<br />
the presence of a high groundwater table. Charcoal could indicate ancient anthropogenic fires<br />
although clearance of vegetation is not visible in the pollen diagrams. Analysis of more<br />
samples may not increase the knowledge of the macrofossils because these three samples<br />
show that the material was largely degraded.<br />
Some of the herbs are grouped as anthropogenic indicators according to (Behra, 1981). Both<br />
pollen zones overviews of (Bastiaens et al, 2006) and (Jong, 1983) show that the Subboreal is<br />
characterized by the presence of cultivated species. The anthropogenic indicators from the<br />
Eem Polder appear in very low quantities and important indicators for agriculture such as<br />
cereals are absent. If the pollen diagrams of the Eem Polder present a closed forest it might<br />
explain why anthropogenic indicators such as cereals are missing. This would not necessarily<br />
mean the cereals were not growing in the surrounding environment; they release little pollen<br />
and the dispersal zone is very limited (Edwards, Whittington, Robinson and Richter, 2005). It<br />
is possible that people influenced the area around the forest by cultivating plants and that<br />
pollen of this species could not reach the forest area because of its density. The pollen<br />
diagrams most likely show an extra-local vegetation history of a closed forest where herbs<br />
could not develop and in which surrounding vegetation could not disperse their pollen.<br />
The peat layer from which the vegetation remains were sampled prevented decaying of the<br />
organic material and stored it for many years. As the tree-ring patterns of GFE01 and GFE02<br />
and their synchronization with the oaks from Ypenburg shows, the substrate on which the<br />
oaks were growing was not very saturated. The absence of water plants and the low<br />
abundances of typical peat bog species such as peat moss (Sphagnum), confirm that the forest<br />
grew on a relatively dry substrate. Every year the abundant trees dropped their leaves and the<br />
substrate probably turned into a very humus-rich soil. The high absolute number of pollen in<br />
this profile indicates either a slow accumulation rate of soil layers or an extraordinary pollen<br />
release from the vegetation. The calibrated radiocarbon dates can be used to estimate the<br />
accumulation rate of the sediments in the profile. From these dates it is assumed that the<br />
accumulation rate was rather constant (± 0.6 mm/year). At a certain time the humus rich soil<br />
turned into a peat bog, but the results of the pollen analysis do not show when exactly this<br />
happened.<br />
The growth of the oaks GFE01 and GFE02 occurred during the same period as the<br />
development of the profile in which they were found. The oak pollen curve of GP01 does not<br />
show a sudden decrease and this might indicate that there was no event that influenced the<br />
47
death of many trees. The profile of GP02 does show a decrease of oak pollen in the upper<br />
part. At the same part the amount of Cercophora increases. This fungus does not only appear<br />
on dung, but also on decaying wood (Van Geel et al, 1980). Profile GP02 is not dated but<br />
according to the NAP-depths relative to GP01 the decrease of oak pollen and the increase of<br />
the fungi probably happened after the death of the studied oaks.<br />
The results of this study are relevant to our understanding of the natural vegetation in the<br />
Netherlands. The evidence of branchless tree trunks, abundant root systems and a very high<br />
tree ratio relative to herbs or anthropogenic indicators points at the previous existence in the<br />
Eem Polder of a closed deciduous forest. At least locally this finding disagrees with the<br />
hypothesis stated by Vera (1997) that the environment in the Netherlands was characterized<br />
by an open landscape. Although both phases of tree growth identified in the Eem Polder<br />
indicate that a closed forest existed, it is possible that the results only have a restricted<br />
geographical bearing. Therefore the general assumption that the Dutch natural vegetation was<br />
characterized by an open landscape with widely distributed single trees cannot be refuted<br />
based on the results presented here.<br />
48
5. Conclusions<br />
The different trees that were studied during this research have been living during two different<br />
time periods. The oaks from location Groeneweg (GFE01 and GFE02) died after 2555 and<br />
between 2580 and 2564 BC respectively, while the oaks from the Lodijk (LFE01 and LFE02)<br />
died after 1482 and 1458 respectively.<br />
Groeneweg<br />
The studied oaks at Groeneweg grew on a relatively dry substrate. The oaks could not be<br />
dated by synchronization with oaks that grew in wet circumstances and have been preserved<br />
in Dutch bogs, but they gave a significant correlation with oaks from Ypenburg which also<br />
grew on an unsaturated soil. The absence of waterplants or plant species that characterizes<br />
peat formation confirmed that the trees at Groeneweg were growing on an unsaturated soil.<br />
Their synchronization with oaks from Ypenburg showed that the growth patterns of the trees<br />
reflect regional growth conditions, while the pollen analysis showed local conditions.<br />
According to the pollen, the oaks were growing in a closed forest with a dominance of oak,<br />
hazel and alder. Herbs could not develop because of shading by the closed canopy. In<br />
addition, anthropogenic indicators were almost absent, which indicates a minor influence of<br />
people on this area. The pollen diagrams did not show a remarkable decrease of the<br />
appearance of species, implying that the forest did not suffer any dramatic changes or<br />
clearances. This vegetation reconstruction showed the local appearance of plant species and<br />
from this study the regional vegetation development or the geographical distribution of the<br />
forest is unknown. The tree patterns of the oaks showed that the trees died because of a<br />
sudden event, whereas the constant amount of tree pollen indicated that this event did not<br />
influenced the whole forest or even a substantial part of it. Because the trees were preserved<br />
well, one of them even still retaining its sapwood rings, it is possible that peat development<br />
started before the death of the trees. This study does not clarify how the change from an<br />
unsaturated soil into peat development took place, when this occurred and on what scale it<br />
affected the forest near the Groeneweg.<br />
Lodijk<br />
The studied oaks sampled at the Lodijk also were growing in a closed forest, which is<br />
indicated by the fact that the trunks were relatively branchless. Also, near the original location<br />
of the studied oaks more root systems were found, implying that at some point in time more<br />
trees were growing at this location. Because the age of these root systems is unknown, it is<br />
unclear if these trees were growing during the same period as the studied oaks. The<br />
morphology of some root systems indicates a high groundwater level and a saturated soil. The<br />
tree-ring patterns of the studied oaks show that they grew unhindered by any remarkable<br />
stress. Their tree-ring patterns show wider and narrower rings that occur synchronously with<br />
pointer years in the patterns of trees which grew in the province of Zuid-Holland. Remarkable<br />
reduced growth occurred during 1629 and 1628 BC. From this study it is not clear if these<br />
growth reductions are caused by the eruption in 1628 BC of the volcano Thera.<br />
The results of the radiocarbon dating of the pines are not accomplished yet. The relation<br />
between the pines and the oaks from the Lodijk is still unknown. Some pines showed ‘saber<br />
tooth growth’, indicating that they were growing on an unstable substrate.<br />
49
6. Further Research<br />
After finishing this research on the vegetation history of the Eem Polder some questions still<br />
remain. Therefore it is recommendable to perform further research in this area to extend the<br />
knowledge about vegetation history. In both locations near the Groeneweg and near the<br />
Lodijk further research questions can be proposed, but also more knowledge is required about<br />
a possibly relationship between the two locations.<br />
The reconstruction of the vegetation history indicated that a closed forest dominated the<br />
landscape at both study sites, but little is known about the geographical distribution of these<br />
forests. It is likely that the vegetation reconstruction presented the local occurrence of the<br />
vegetation development and therefore it would be interesting to find out on what scale the<br />
forests developed. It might be possible that the forests were distributed over a large area and<br />
maybe they migrated through the landscape. By searching for more tree remains at both<br />
locations the distribution of the forest can possibly be studied. The forest near the Groeneweg<br />
grew about 1000 years earlier than the forest near the Lodijk but possibly there is a relation<br />
between the two. Because the tree remains from the Lodijk used during this study were not<br />
found in situ, nothing is known about the surrounding vegetation development in this area. If<br />
in situ tree remains can be localized, a pollen analysis should be performed at this location<br />
too. Furthermore, the relation between the collected pines and oaks from the Lodijk is still<br />
unknown. By dating the pines this relationship can be revealed and the remarkable growth<br />
forms of the pines can also explain the growth conditions of the oaks.<br />
The pointer year analysis indicated that the tree-ring patterns of the studied trees show<br />
similarities with trees growing in the province of Zuid-Holland. Maybe it is possible to reveal<br />
what kind of growth conditions were responsible for these similarities and if remarkable<br />
environmental changes occurred during this period.<br />
There are also some questions left about the downfall of the forest discovered near the<br />
Groeneweg. According to the pollen analysis no decrease in tree pollen occurred during the<br />
studied period and no remarkable event influenced the whole forest. But because both trees<br />
were preserved very well the burying of the trees must have happened soon after their death.<br />
This implies that peat development already started and it would be logical that more trees<br />
suffered from the increased saturation of the soil. Maybe it is possible that on a local scale no<br />
peat developed and the humus-rich soil characterized the substrate. Information about this can<br />
be obtained by biochemical research on the soil profile. More information about the soil<br />
development and the relation with the downfall of single trees or the whole forest is needed to<br />
discover the sequence of events. To investigate the influence of changes in soil development<br />
on the forest a pollen analysis should be performed on younger soil samples to discover a<br />
decrease of trees. Also dendrochronological research on other tree remains might answer<br />
questions about environmental changes and other tree-ring patterns may indicate stress<br />
conditions of growth.<br />
The difference between the radiocarbon and the dendrochronological dating of the oaks from<br />
the Groeneweg should be explained in more detail. Maybe it is possible to revise the<br />
radiocarbon dates by studying the significance of the statistics used with this dating method.<br />
To two locations in the Eem Polder are very interesting for further research and hopefully<br />
more details about the vegetation history of this area can be discovered. Also, besides these<br />
two locations it is likely that more tree remains can be collected from many other locations<br />
50
within the Eem Polder. Further research can hopefully complete the reconstruction of the<br />
palaeo-ecology of the Eem Polder.<br />
7. References<br />
Baillie, M.G.L. (1982). ‘Tree-ring dating and Archaeology’. In Varley, G.C., and Gradwell,<br />
G.R. 1962. The effect of partial defoliation by the Caterpillars on the Timber<br />
Production of oak trees in England. XI Internationaler Kongress fur Entomologie<br />
Wien 1960, Sonderdruck aus den Verhandlungen, Bd II, pp. 211-14.<br />
Baillie, M.G.L. (1989). ‘Proceedings of the third International Congres, Santorini’.<br />
Baillie, M.G.L. (1991). ‘Marking in Marker Dates: Towards an Archaeology with historical<br />
Precision’. World Archaeology, volume 23, issue 2, pp 233-243.<br />
Bakker, M. and Smeerdijk, D.G. van. (1981). ‘A palaeo-ecological study of a late Holocene<br />
section from ‘Het Ilperveld’, Western Netherlands’. Review of Palaeobotany and<br />
Palynology, 36, pp 95-103.<br />
Bastiaens, J., Brinkkemper, O., Deforce, K., Maes, B., Rovekamp, C., Bremt, P van den,<br />
Zwaenepoel, A. (2006). Inheemse bomen en struiken in Nederland en Vlaanderen,<br />
herkenning, verspreiding, geschiedenis en gebruik. Amsterdam, Boom.<br />
Behra, K. (1981).’The interpretation of anthropogenic indicators in pollen diagrams’. The<br />
Journal of Ecology, volume 76, issue 3, pp 903.<br />
Bennett, K.D. and Willis, K.J. (2002). Pollen. In, J.P. Smol, H.J. B. Birks and W.M. Last<br />
(eds.) Tracking Environmental Change Using Lake Sediments, Volume 3: Terrestrial,<br />
Algal, and Siliceous Indicators. Kluwer Academic Publishers, Dordrecht.<br />
Berendsen, J.H.A, & W.H. Zagwijn. (1984). ‘Some conclusions reached at the symposium on<br />
geological changes in the western Netherlands during the period 1000-1300 AD.’<br />
Geologie en Mijnbouw 63, pp 225-229.<br />
Berglund, B. And Persson, T. (2004). Information on spore tablets Retrieved from<br />
http://www.geol.lu.se/personal/tsp/Ly483216Eur.pdf.<br />
Birks, H.H and Birks, H. J. B. (2000). Future Uses of Pollen Analysis Must Include Plant<br />
Macrofossils. Journal of biogeography, volume 37, issue 1, pp 31-35.<br />
Birks, H.J.B. and Birks, H.H. (1980). Quaternary palaeoecology. New Jersey,<br />
The Blackburn Press.<br />
Birks, H.J.B. and Birks, H.H. (2006). ‘Multi-proxy studies in palaeolimnology’. Vegetation<br />
history and Archaeobotany, volume 15, issue 4, pp 235-251.<br />
Callaway, R.M.(1990). ‘Effects of soil water distribution on the lateral root development of<br />
three species of Californian oaks’. American Journal of Botany, volume 77, issue 11,<br />
pp 1469-1474.<br />
Cook, E.R. and Kairiukstis, L.A. (1990). Methods of dendrochronology. Applications in the<br />
51
environmental sciences. Dordrecht, Kluwer Academic Publishers.<br />
Daalen, S van. (2001). Een reconstructie van het moerasbos te Ypenburg. MSc thesis,<br />
University of Wageningen<br />
Desplanque, C., Rolland, C. and Schweingruber, F.H. (1999). ‘Influence of species and<br />
abiotic factors on extreme tree ring modulation: Picea abies and Abies alba in<br />
Tarentaise and Maurienne (French Alps)’. Trees, volume 13, issue 4, pp 218-227 .<br />
Edwards, K.J., Whittington, G., Robinson, M. and Richter, D. (2005). ‘Palaeoenvironments,<br />
the archaeological record and cereal pollen detection at Clickimin, Shetland,<br />
Scotland.’ Journal of archaeological science, volume 32, issue 12, pp 1741-1756.<br />
Eide, W., Birks, H.H., Bigelow, N.H., Peglar, S.M. and Birks, H.J.B. (2006). ‘Holocene forest<br />
development along the Setesdal valley, southern Norway, reconstructed from<br />
macrofossils and pollen evidence.’ Vegetation history and archaeobotany, volume 15,<br />
issue 2, pp 65-85.<br />
Fritts, H.C. (1976). Tree rings and climate. London, Academics.<br />
Geel, B van, Bohncke, S.J.P. and Dee, H. (1980). ‘A palaeo-ecological study of an upper late<br />
glacial and Holocene sequence from ‘De Borchert’, The Netherlands.’ Review of<br />
Palaeobotany and Palynology, 31, pp 367-448.<br />
Henri D. Grissino-Mayer, Department of Geography, The University of Tennessee, from:<br />
http://web.utk.edu/~grissino (<strong>2007</strong>).<br />
Jansma, E. (1995). Rememberings. The development and application of local and regional<br />
tree-ring chronologies of oak for the purpose of archaeological and historical<br />
research in the Netherlands. PhD dissertation, University of Amsterdam.<br />
Jansma, E. (2006). Dendrochronologie, NoaA Chapter 3 (versie 1.0), (www.noaa.nl), 1-40.<br />
Jansma, E., (<strong>2007</strong>). Datering, herkomst en bouwvolgorde van De Meern 4. In: T. de Groot &<br />
J.-M.A.W. Morel (red.), <strong>2007</strong>: Het schip uit de Romeinse tijd De Meern 4 nabij<br />
boerderij de Balije, Leidsche Rijn, gemeente Utrecht. Waardestellend onderzoek naar<br />
de kwaliteit van het schip en het conserverend vermogen van het bodemmilieu.<br />
Amersfoort, RACM.<br />
Janssen, C.R. (1972). ‘The palaeo-ecology of plant communities in the Dommel Valley,<br />
North Brabant, Netherlands.’ The Journal of Ecology, volume 60, issue 2, pp 411-437.<br />
Janssen, C.R. (1986). ‘The use of local pollen indicators and of the contrast between regional<br />
and local pollen values in the assessment of the human impact on vegetation.’<br />
Quarternary research,16, pp 80-96.<br />
Jong, J de. 1983. Chronostratigraphic subdivision of the Holocene in the Netherlands.<br />
52
Kitagawa, H. and van der Plicht, J. (1998). ‘Atmospheric Radiocarbon Calibration to 45.000<br />
yr B.P.: Late Glacial Fluctuations and Cosmogenic Isotope Production.’ Science,<br />
volume 279, issue 5354, pp 1187-1190.<br />
Kristensen, E., Ahmed, S.I. and Devol, A.H. (1995). ‘Aerobic and Anaerobic decomposition<br />
of organic matter in marine sediment: Which is fastest?’ Limnology and<br />
oceanography, volume 40, issue 8, pp 1430-1437.<br />
Mauseth, J.D. (1998). Botany, an introduction to plant biology. University of Texas, Austin.<br />
Sudbury, Jones and Bartlett Publishers.<br />
Mulder, E.F.J. de, Geluk, M.C., Ritsema, I.L., Westerhoff, W.E. and Wong, T.E. (2003). De<br />
ondergrond van Nederland. Groningen, Wolters Noordhoff.<br />
Pedersen, B.S. (1998). ‘The Role of Stress in the Mortality of Midwestern Oaks as Indicated<br />
by Growth prior to Death;. Ecology, volume 79, issue 1, pp 79-93.<br />
Peeters, J.H.M. (<strong>2007</strong>). Hoge vaart-A27 in context: Towards a model of Mesolithic –<br />
Neolithic land use dynamics as a framework for archaeological heritage management.<br />
PhD dissertation, University of Amsterdam.<br />
14<br />
Plicht, J. van der. (2006). De C methode, NoaA Chapter 4 (versie 1.0), (www.noaa.nl), 1-25.<br />
Reimer, P.J., Baillie, M.G.L., Bard, E., Bayliss, A., Beck, J.W., Bertrand, C.J.H., Blackwell,<br />
P.G., Buck, C.E., Burr, G.S., Cutler, K.B., Damon, P.E., Edwards, R.L., Fairbanks,<br />
R.G., Friedrich, M., Guilderson, T.P., Hogg, A.G., Hughen, K.A., Kromer, B.,<br />
McCormac, G., Manning, S., Ramsey, C.B., Reimer, R.W., Remmele, S., Southon,<br />
J.R., Stuiver, M., Talamo, S., Taylor, F.W., Van der Plicht, J. and Weyhenmeyer, C.E.<br />
(2004). ‘IntCal04 terrestrial radiocarbon age calibration, 0-26 cal. KYR BP.’<br />
Radiocarbon 46, pp 1029-1058.<br />
Ruddiman, F.W.(2002). Earth’s climate: past and future. New York, W. H. Freeman and<br />
company.<br />
Ruelle, J., Yamamoto, H. and Thibaut, B. (<strong>2007</strong>). ‘Growth stresses and cellulose structural<br />
parameters in tension and normal wood from three tropical rainforest angiosperm<br />
species.’ Bioresources, volume 2, issue 2, pp 235.<br />
Schweingruber, F.H. (1996). Tree rings and environment dendroecology. Berne, Paul Haupt<br />
Publishers.<br />
Schweingruber, F.H., Eckstein, D., Serre-Bachet, F. and Bräker, O.U. (1990). ‘Identification,<br />
presentation and interpretation of event years and pointer years in dendrochronology.’<br />
Dendrochronologia 8, 9-38.<br />
Streiff, R., Ducousso, A., Lexer, C., Steinkellner, H., Gloessl, J. And Kremer, A. (1999).<br />
‘Pollen dispersal inferred from paternity analysis in a mixed oak stand of Quercus<br />
robur L. and Q. petraea (Matt.) Liebl.’ Molecular Ecology, volume 8, issue 5, pp 831,<br />
11p.<br />
53
Vera. F.W.M. (1997). Metaforen voor de Wildernis. Eik, Hazelaar, rund en paard. ’s-<br />
Gravenhage, Ministerie van Landbouw, Natuurbeheer en Visserij.<br />
Watt, A.S. (1947). ‘Pattern and process in the plant community’. Journal of Ecology,<br />
35, 1-22.<br />
Weinig, C. (2000). ‘Differing selection in alternative competitive environments: shadeavoidance<br />
responses and germination timing.’ Evolution, volume 54, issue 1, pp 124-<br />
136.<br />
Yeloff, D., Bennett, K.D., Blaauw, M., Mauquoy, D., Sillasoo, U., van der Plicht, J. and van<br />
Geel, B. (2006). ‘High precision 14C dating of Holocene peat deposits: A comparison<br />
of Bayesian calibration and wiggle-matching approaches.’ Quaternary geochronology,<br />
volume 1, issue 3, pp 222-235.<br />
Zagwijn, W.H. (1973). ‘Pollenanalytic studies of Holsteinian and Saalien beds in the northern<br />
Netherlands’. Mededelingen Rijks Geologische Dienst, Nieuwe serie 24, pp 139-156.<br />
54
Appendix<br />
Appendix 1<br />
Table 13 (Grissino-Mayer, <strong>2007</strong>)<br />
Principle<br />
The Uniformitarian Principle<br />
The Principle of Limiting Factors<br />
The Principle of Aggregate Tree Growth<br />
Explanation<br />
‘The present is the key to the past’, stated by<br />
James Hutton in 1785; by knowing the<br />
relationship between tree-rings and climate at<br />
present, it is possible to make climate<br />
reconstructions of the past using tree-rings of<br />
old, long-lived trees.<br />
The environmental variable that is most<br />
limiting for plant growth, determines the<br />
conditions of the plant; if tree growth is<br />
limited by precipitation, trees will not grow<br />
faster than allowed by the amount of rainfall,<br />
and ring-width is a function of precipitation<br />
(Fritts, 1976).<br />
Individual tree growth can be described as a<br />
function of different environmental factors,<br />
both human and natural, that effect the tree<br />
growth. This is given as the function:<br />
R t= A t + C t + δD1 t + δD2 t + Et<br />
R t = the observed ring-width series<br />
A t = the age-size related trend in ring width<br />
C t = the climatically related environmental<br />
signal<br />
D1 t = the disturbance pulse caused by a local<br />
endogenous disturbance<br />
D2 t = the disturbance pulse caused by a<br />
standwide exogene disturbance<br />
E t = the largely unexplained year-to-year<br />
variability not related to other signals<br />
(Cook and Kairiukstis, 1990).<br />
The Principle of Ecological Amplitude<br />
Using this formula we assume that A t, Ct and<br />
E t are always present, while D1 t and D2t are<br />
only present (δ=1) if a disturbance had<br />
occurred at a certain year, and not present<br />
(δ=0) if no disturbances occurred (Cook and<br />
Kairiukstis, 1990).<br />
The ecological amplitude refers to the area in<br />
which tree species can grow and reproduce.<br />
This principle is used within<br />
dendrochronology because the most useful<br />
trees are found near the margins of their<br />
natural range, determined by latitude,<br />
longitude and elevation.<br />
55
Appendix 2: Tree growth and wood anatomy<br />
There are two fundamental types of plant bodies: the primary plant body (herbaceous plants<br />
that only undertake primary growth) and the secondary plant body (trees or shrubs with<br />
secondary growth that produce wood and bark) (Mauseth, 1998).<br />
Plants contain two types of vascular tissue, xylem and phloem, which both are types of<br />
transport tissue. Xylem conducts water and minerals and phloem distributes sugars and<br />
minerals (Mauseth, 1998). Water and minerals enter the xylem in the roots and are conducted<br />
upwards to the leaves and stems, travelling through dead cells. Once it enters in the right<br />
location, the water evaporates and the minerals are absorbed (Mauseth, 1998). The phloem<br />
consists of living cells that pick up sugars from areas where they are abundant and transport<br />
them to areas where they are needed. Both types of vascular tissues exist in primary and<br />
secondary plant bodies. The vascular cambium is a thin layer of cells located between the<br />
wood and the bark (Figure 37). When active (during the growing season) it produces xylem<br />
towards to inside and phloem or bark towards the outside (Baillie, 1982). The growing season<br />
on temperate climates ranges from March to late September.<br />
Figure 37: Anatomy of a tree stem.<br />
56
At the beginning of the growing season the cambium of deciduous trees produces large<br />
vessels to provide the newly grown leaves with transport canals; these vessels are also called<br />
spring vessels (Figure 38). As the growing season proceed and the tree is full of adult leaves,<br />
stability becomes the major priority and the cambium produces small and compact cells that<br />
will provide strength to the whole structure. This type of wood is called summerwood (Figure<br />
38) (Mauseth, 1998). At the end of the growing season cell division stops; the last cells are<br />
formed as heavy fibers with thick secondary walls. The spring vessels and the summerwood<br />
together represent the growth of one year, also called an annual ring (Mauseth, 1998).<br />
Spring vessels<br />
Summerwood<br />
Spring vessels<br />
Figure 38a: Spring vessels<br />
(white cells) and<br />
summerwood (brown bands)<br />
together represent one annual<br />
ring<br />
Figure 38b: a detailed picture<br />
of spring vessels and<br />
summerwood<br />
The width of individual growth rings may vary greatly from year to year as a function of<br />
environmental factors such as light, temperature, rainfall and available soil. The forming of<br />
spring vessels is dependent on the reserves of the tree built up in previous years and the<br />
summerwood is dependent on the food supplies of the particular year during which the<br />
summerwood formation takes place. The width of the spring vessels varies little from year to<br />
year, which means that the variation of ring widths in deciduous trees in fact represents the<br />
variation in the summerwood width.<br />
The obvious ring pattern, the long lifetime and the annual character of oaks are advantages for<br />
dendrochronology. In the cross sections of oaks, two zones can be distinguished: the centre,<br />
almost always darker in colour than the outer wood, is called heartwood. The outer section,<br />
usually moister and lighter in colour is called the sapwood. These two regions exist because<br />
vessels can break as the result of freezing, wind vibration, insects and other factors. If a<br />
column is broken, it makes no sense to pull water in for transport and the vessel usually never<br />
conducts water again. Therefore surrounding parenchyma cells push protoplasm through the<br />
pit into the vessels, forming a plug called tylosis that will seal them off (Mauseth, 1998).<br />
These changes make the heartwood more durable, less penetrable and more protected against<br />
decay by organisms. The sapwood contains living parenchyma and waterfilled conducting<br />
vessels, full of xylem sap. A new layer of sapwood is formed every year, and one annual ring<br />
is converted into heartwood on average each year (Mauseth, 1998). The thickness of the<br />
sapwood of oaks is more or less constant, making it useful for dendrochronological research.<br />
The wane edge is the last ring under the bark, presenting the last year of growth. It is possible<br />
to date the death of a tree exactly if the wane edge is still present. Despite the fact that it is<br />
very unusual to find complete oak trunks from the past, including all sapwood rings and bark,<br />
it is possible to estimate the original amount of sapwood rings, and thus to estimate the death<br />
of the tree.<br />
57
Coniferous trees may not form a ring along the total circumference of the stem if they are<br />
influenced by extreme weather (Fritts, 1976). When this happens, the ring is missing along<br />
the radii of a tree-ring series. It is also possible that growth conditions temporarily become<br />
poor before the end of the growing season. During this period vessels grow smaller with<br />
thicker walls. When conditions improve again during the same season, cells are formed larger<br />
with thinner walls. These intra-annual growth bands are also referred as false rings (Fritts,<br />
1976). Trees are not only influenced by climatological factors, but also on abiotic ones. Trees<br />
react on mechanical stress such as wind exposure, soil movement or shading. If the substrate<br />
on which a tree grows moves, the tree will respond to restore verticality. Trees are able to<br />
bend their trunk by the formation of stressed wood, called reaction wood, to stabilise or reach<br />
for light (Ruelle, Yamamoto and Thibaut, <strong>2007</strong>). If ring widths of coniferous trees are<br />
measured for dating or climatological research one should be aware of missing rings, false<br />
rings or reaction wood.<br />
58
Appendix 3: Overview of archaeological periods describing dendrochronological calendars,<br />
based on dendrochronological research until 2004 (Jansma, 2006).<br />
6025 - 2850 BC<br />
The dendrochronological dataset of this period contains 95 measurements of oaks from<br />
natural habitats in ancient peat bogs and river areas. It also contains 7 measurements on oaks<br />
from a road and a bridge made of wood through a peat bog. From cover-sand areas at the<br />
eastern and southern Netherlands no wood data are available. Chronologically, the data set of<br />
this period contains two hiatus: 5467 - 5121 BC and 4557 - 3643 BC.<br />
2850 - 2000 BC<br />
The dendrochronological dataset of this time interval contains 115 measurements of oaks<br />
from natural habitats in ancient peat bogs and 10 measurements from archaeological sites.<br />
The dataset covers the whole period, without any hiatus. The oaks of this dataset grew on<br />
different soils (also sand-cover areas) but their geographical distribution is restricted; 66<br />
series are derived from Ypenburg and a cluster of oaks from the southern area of the province<br />
of Zuid-Holland.<br />
2000 - 800 BC<br />
This dendrochronological dataset contains 83 measurements of oaks from natural habitats<br />
(peat bogs) mostly from the western past of the Netherlands. It contains 5 measurements of<br />
oaks from archaeological sites. Except for the western part of the Netherlands, this period<br />
must yet be described dendrochronologically.<br />
800 - 12 BC<br />
Measurements of 103 oaks from natural habitats and 22 oaks from archaeological sites cover<br />
this whole period. Within this dataset especially the northern and eastern part of the<br />
Netherlands and the province of Zeeland are not presented adequately.<br />
12 BC - AD 400<br />
For this period an extraordinary amount of dendrochronological measurements are available;<br />
296 measurements from native Roman context (5 silver fir, 14 ash, 277 oak), 290 from<br />
military context (26 elm and ash, 264 oak), 16 from civil context (oak), 33 from religious<br />
context (oak) and 62 from natural context (26 beech, 36 oak). The northern part of the<br />
Netherlands is not presented well in this dataset.<br />
AD 400 - 1000<br />
This period contains 179 measurements from archaeological context (oaks) and 62<br />
measurements from natural context (12 beech, 50 oak). The dataset covers the whole period,<br />
but the native material found in situ dates no later than AD 570. The Dutch calendar of native<br />
oaks found in situ here reaches its chronological end. In this dataset the northern part of the<br />
Netherlands and the province of Zeeland are not presented well.<br />
AD 1000 - 1500<br />
From this period 464 archaeological measurements are available (2 pine, 3 ash, 20 silver fir,<br />
439 oak). The eastern part of the Netherlands is not presented well in this period.<br />
During this period large amounts of wood were imported, and therefore the material is not<br />
useful to develop native calendars. The data are useful to determine the origin of transported<br />
wood and to develop regional import calendars.<br />
59
AD 1500 - 1950<br />
The archaeological dataset of this period contains 245 measurements (2 beech, 2 pine, 8 silver<br />
fir, 233 oak). Besides this a lot of data are available from monumental buildings still present.<br />
Some of these data overlap with the time span represented by the growth rings in living oaks.<br />
But the data are in many cases not native and therefore not very useful for indigenous<br />
calendars and climatological research directed at the Netherlands.<br />
60
Appendix 4<br />
Table 14: Minimum sample amount vs. significance<br />
n 99.9 99 95 90<br />
1 - - - -<br />
2 - - - -<br />
3 - - - -<br />
4 - - - -<br />
5 - - - 97<br />
6 - - 98 92<br />
7 - - 94 88<br />
8 - - 91 85<br />
9 - 98 88 83<br />
10 - 96 86 81<br />
11 - 93 84 79<br />
12 - 91 82 78<br />
13 100 89 81 77<br />
14 98 88 80 75<br />
15 96 87 79 75<br />
16 94 85 78 74<br />
17 93 84 77 73<br />
18 92 83 76 72<br />
19 90 82 75 71<br />
20 89 81 74 71<br />
21 88 80 74 70<br />
22 87 80 73 70<br />
23 87 79 73 69<br />
24 86 78 72 69<br />
25 85 78 72 68<br />
26 84 77 71 68<br />
27 84 77 71 68<br />
28 83 76 70 67<br />
29 82 76 70 67<br />
30 82 75 70 67<br />
31 81 75 69 66<br />
32 81 74 69 66<br />
33 80 74 69 66<br />
34 80 74 68 66<br />
61
Table 15: Pointer years of 8 trees from Ypenburg<br />
JAAR No % Min/Max a<br />
-2794 5 100% -<br />
-2791 6 100% + **<br />
-2790 6 100% + **<br />
-2776 6 100% - **<br />
-2774 6 100% + **<br />
-2772 7 100% - **<br />
-2770 7 100% + **<br />
-2764 7 100% + **<br />
-2762 7 100% - **<br />
-2759 7 100% - **<br />
-2758 7 100% - **<br />
-2754 7 100% + **<br />
-2746 7 100% + **<br />
-2744 7 100% + **<br />
-2738 8 88% +<br />
-2737 8 88% -<br />
-2735 8 100% + **<br />
-2734 8 100% + **<br />
-2721 7 100% - **<br />
-2715 7 100% + **<br />
-2697 7 100% + **<br />
-2680 5 100% +<br />
****** : z = 3.30 and certainty 99.9%<br />
**** : z = 2.57 and certainty 99%<br />
** : z = 1.96 and certainty 95%<br />
- : z = 1.64 and certainty 90%<br />
90.0% certainty: 14 maxima and 8 minima<br />
95.0% certainty: 12 maxima and 6 minima<br />
99.0% certainty: 0 maxima and 0 minima<br />
99.9% certainty: -2 maxima and -2 minima<br />
Time span =-2889 - -2658<br />
62
Table 16: A total of 10 trees were compared to determine pointer years.<br />
Dated studied<br />
chronologies<br />
GFE010<br />
GFE020<br />
Ypenburg<br />
ypf201<br />
ypf602<br />
ypf061<br />
ypf131<br />
ypf221<br />
ypf231<br />
ypf241<br />
ypf733<br />
63
Table 17: Pointer years of 32 trees from the province of Zuid-Holland<br />
JAAR No % Min/Max a JAAR No % Min/Max a<br />
-1699 7 100% + ** -1585 16 88% - ****<br />
-1689 10 90% + ** -1584 16 75% -<br />
-1687 11 91% - ** -1582 16 94% + ******<br />
-1686 11 82% - -1569 13 92% - ****<br />
-1683 11 82% - -1557 9 89% + **<br />
-1679 11 91% + ** -1549 8 88% -<br />
-1675 11 91% + ** -1544 6 100% - **<br />
-1673 14 86% + ** -1543 6 100% - **<br />
-1670 14 93% - **** -1538 6 100% + **<br />
-1669 14 93% + **** -1534 6 100% + **<br />
-1668 14 79% - -1522 7 100% + **<br />
-1667 14 93% - **** -1511 7 100% - **<br />
-1664 14 86% + ** -1509 8 88% +<br />
-1660 14 100% + ****** -1503 8 88% +<br />
-1653 14 86% - ** -1495 8 88% -<br />
-1652 14 86% + ** -1485 7 100% - **<br />
-1651 14 86% + ** -1470 6 100% - **<br />
-1648 14 86% - ** -1458 9 89% + **<br />
-1646 14 93% - **** -1455 9 100% - ****<br />
-1643 14 86% + ** -1453 9 100% - ****<br />
-1640 15 80% - ** -1446 9 89% - **<br />
-1638 15 80% + ** -1445 9 89% + **<br />
-1637 15 80% - ** -1435 10 100% - ****<br />
-1633 16 81% - ** -1434 10 100% - ****<br />
-1632 16 75% + -1429 10 90% + **<br />
-1629 16 100% - ****** -1425 8 88% +<br />
-1628 16 94% - ****** -1423 8 88% -<br />
-1625 15 87% + **** -1420 8 100% - **<br />
-1624 15 80% + ** -1417 8 88% -<br />
-1622 15 93% + **** -1414 8 88% -<br />
-1621 15 93% - **** -1413 8 88% -<br />
-1618 16 75% + -1407 8 88% -<br />
-1617 16 75% - -1404 10 100% + ****<br />
-1615 16 75% - -1392 9 89% + **<br />
-1613 16 88% + **** -1390 8 100% + **<br />
-1611 16 75% - -1385 9 100% - ****<br />
-1608 16 81% - ** -1383 10 90% - **<br />
-1605 16 94% + ****** -1381 11 91% + **<br />
-1604 16 81% + ** -1373 13 92% - ****<br />
-1602 17 94% - ****** -1372 13 77% -<br />
-1601 17 76% - -1368 13 85% + **<br />
-1599 17 82% + ** -1367 13 77% -<br />
-1597 17 88% + **** -1361 13 77% +<br />
-1594 17 100% - ****** -1360 13 77% -<br />
-1593 17 76% - -1358 13 85% - **<br />
-1591 17 82% + ** -1357 13 92% + ****<br />
-1588 16 81% + ** -1356 13 85% + **<br />
-1587 16 88% - **** -1355 13 77% +<br />
-1353 13 77% -<br />
64
****** : z = 3.30 and certainty 99.9%<br />
**** : z = 2.57 and certainty 99%<br />
** : z = 1.96 and certainty 95%<br />
- : z = 1.64 and certainty 90%<br />
90.0% certainty: 44 maxima and 53 minima<br />
95.0% certainty: 37 maxima and 33 minima<br />
99.0% certainty: 10 maxima and 17 minima<br />
99.9% certainty: 1 maxima and 2 minima<br />
Time span = -1701 - -1350<br />
Table 18: The original location of the 32 trees from the province of Zuid-Holland<br />
Site name Start Year End Year Locations<br />
(BC) (BC)<br />
Site1c_b -1748 -767 Ouderkerk a/d IJssel<br />
Papendrecht<br />
Site 1b -2256 -1141<br />
Papendrecht<br />
Papendrecht<br />
Sliedrecht<br />
Alphen a/d Rijn<br />
Bodegraven<br />
Gouderak<br />
Hazerswoude-Rijndijk<br />
Hazerswoude<br />
Hazerswoude<br />
Oudewater<br />
Zoeterwoude dorp<br />
65
Table 19: A total of 34 trees were compared to determine pointer years,<br />
32 tree measurements from the province of Zuid-Holland were included.<br />
Dated studied Site1c_b Site 1b<br />
chronologies<br />
LFE010<br />
LFE020<br />
OKf082<br />
PAf041<br />
PAf071<br />
PAFB20<br />
PAFo31<br />
PAFo81<br />
PAF141<br />
PAF152<br />
SCf161<br />
afe142<br />
bof021<br />
bof041<br />
gof011<br />
gof020<br />
gof031<br />
gof051<br />
gof061<br />
gof071<br />
gof081<br />
hrf030<br />
hwf030<br />
hwf041<br />
hwf111<br />
hwf181<br />
hzfn11<br />
ouf022<br />
ouf032<br />
ouf040<br />
ouf061<br />
zdf040<br />
zdf070<br />
zwf042<br />
66
Appendix 5: absolute pollen diagrams<br />
Eem Eemvalley, Valley, GFE01 GP01<br />
1300<br />
Total concentration<br />
Eemvalley, GFE02<br />
Eem Valley, GP02<br />
1150<br />
Total concentration<br />
1350<br />
1200<br />
1400<br />
1250<br />
1450<br />
1300<br />
1500<br />
1350<br />
Depth<br />
1550<br />
Depth<br />
1400<br />
1600<br />
1450<br />
1650<br />
1500<br />
1700<br />
1550<br />
1750<br />
1800<br />
2000 4000<br />
(* 1000)<br />
1600<br />
500 1000 1500<br />
(* 1000)<br />
67