Abstract
Keywords
Introduction
Lina myr, a drained and cultivated area of fenland in central-eastern Gotland, is an area of archeological significance with a complex history of shoreline displacement. Archeological excavations have indicated that people have been active in the local area since the Mesolithic period (Lund, 1996; Österholm, 1989; Seving, 1986; Wallin, 2010; Welinder, 1975) and during the Neolithic period (Hallström, 1971; Sundberg, 2008). Nowadays these archeological sites are located inland but were perhaps once situated next to the Littorina Sea shorelines. Archeological excavations at Gothemshammar, 7 km northeast of Lina myr (Figure 1), unearthed the foundations of a rampart, consisting of a 500-m-long cavity wall. Radiocarbon dating of the domestic materials found intermixed with the wall foundations indicate that the rampart was built 2850–2540 cal. yr BP (Wallin, 2010; Wallin and Martinsson-Wallin, 2018). A large Bronze Age stone cairn and stone ship setting (a type of burial custom common on Gotland) 5 km northeast of Lina myr (Martinsson-Wallin and Wehlin, 2017) also confirms the historical importance of the area for Bronze Age inhabitants of Gotland. Around 400 stone ship settings are found on Gotland and are usually located along the Ancylus Lake and Littorina Sea shorelines of Gotland (Hansson, 1927). According to Wehlin (2010), around 15% of the stone ships on Gotland can be found around Lina myr. Martinsson-Wallin and Wehlin (2017) stated that several of these stone ships are located along the northern side of the Gothem River and Lina myr and are positioned in an eastwest direction. This differs from the majority of stone ships settings on Gotland, which are found in coastal locations around the island placed with the bow toward the south. According to Ohlsson (1984) at least up to medieval times, a historical waterway, comprising the Gothem River and some areas of open water, could have been used for transport through Lina myr, into central Gotland. This important waterway could explain the difference in the orientation in the ship settings. The outlet of this water system is located at Vitviken bay, northeast of Lina myr. This sea connection would have made Lina myr an important center for transport, trade, and settlement.

(a) The location of Gotland, (b) Lina myr, and (c) surface geology (SGU, 2017).
Previous paleoecological studies have used fossil pollen to research the vegetation history of Gotland (Eriksson, 1992; Munthe et al., 1927; Österholm, 1989; Påhlsson, 1977; Pettersson, 1958; Sernander, 1894). A few studies have focused on the paleoecology of Lina myr. Thomasson (1927) described the diatom stratigraphy of Lina myr; Lundqvist (1928) described diatoms and mollusks in a layer of gyttja, which he attributed to Littorina Sea deposits. Svensson (1989) investigated the shore displacement of Gotland using pollen stratigraphy as a method of correlating sites all over the island. However, there have been no studies directly aimed at investigating the prehistoric human impacts on vegetation at Lina myr.
In order to understand when Lina myr may have been connected to the sea, and therefore possibly used as an inland waterway for transport, it is important to explore the history of shore displacement in this area. However, the history of sea-level change and isostatic uplift around the Baltic basin is complex. Discrepancies between shore displacement curves from different parts of the Baltic basin have previously been discussed by Björck and Svensson (1992) and Uścinowicz (2003). There are also apparent inconsistencies around the timing of the Littorina transgression and number of sub-phases.
Based on studies in Blekinge, southeastern Sweden, Berglund (1971), Berglund and Sandgren (1996), Berglund et al. (2005), and Risberg et al. (2005) suggested that the Littorina transgression onset occurred at about 8500 cal. yr BP and distinguished five separate phases. Several transgressive sub-phases have also been seen in Denmark (Christensen and Nielsen, 2008; Wohlfarth et al., 2008). However, in Finland only one uniform transgression has been identified (Eronen, 1974; Miettinen et al., 2007). Similarly, only one main transgression has been detected in Estonia (Raukas et al., 1995). At Vääna lagoon, Estonia, the Littorina transgression began about 8300 cal. yr BP and continued until about 7000 cal. yr BP. A similar timing was observed on Saaremaa Island, Estonia, where the transgression began about 8300–8200 cal. yr BP and lasted until 7300 cal. yr BP (Saarse et al., 2009a, 2009b). A twofold transgression started at 8400 cal. yr BP on the Karelian Isthmus is suggested by Miettinen et al. (2007). Rosentau et al. (2013) proposed that the Littorina transgression occurred as a single event and started 8300 cal. yr BP based on coastal sediments in the eastern Gulf of Finland.
The detailed shoreline displacement and vegetation history in relation to archeological sites, climate, and uplift is poorly understood on Gotland. Our overall aim was to investigate the interplay between sea-level change, the development of waterways, human settlement, and vegetation changes on Gotland. Sea-level fluctuations were determined using diatom and geochemical proxies from new sediment cores from Lina myr. The timing and impact on vegetation by humans was investigated using fossil pollen.
Materials and analysis
Study site
Lina myr (57° 34′ 31.764″ N, 18° 37′ 59.124″ E), situated 6.3 km southwest of Åminne and 21 km east of Visby, is a 9-km2 area of cultivated, drained fenland. The area stretches for 6 km from southwest to northeast and 2 km in an east to west direction. Lina myr drains in a northeast direction toward the east coast of Gotland via the Gothem River. The Gothem River was straightened during the draining of the fen in the 1940s. Today, the elevation of the Lina myr fen is about 9 m a.s.l. The original basin threshold is estimated to be 9.5–9.9 m a.s.l., based on a cross-section established from modern digital evaluation models (Barliaev, 2017). To the west of the fen the surface geology consists of marl and further to the west of glacial till (Figure 1c). To the east, limestone bedrock is exposed. The island of Gotland consists of Silurian limestone (Hede, 1925). The bedrock of Lina myr is classified as the Halla Formation (Laufeld, 1974).
Coring, sampling, and stratigraphy
Fieldwork was conducted in August 2016 and coring was carried out along three transects (Figure 2) to investigate the stratigraphy of the fen. A representative master core was collected at site 7 in section B. A Russian peat corer, 1 m in length and 4.5 cm in diameter, was used for sampling. Supplementary cores were also required in order to provide material for dating. Therefore, 15 parallel cores were retrieved as overlapping sequences in a 2-m2 grid surrounding the master core. A Russian peat corer, 70 cm long and 7 cm in diameter, was used to retrieve this material. The cores were stored in a cold room with a constant temperature of 4–5°C. The stratigraphy of the master core was investigated in the field, combined with a more detailed description in the laboratory. The description was based on analysis of the material under a stereo microscope with ×25 magnification using Munsell color charts, combined with loss on ignition.

The stratigraphy of Lina myr fen and coring locations.
Radiocarbon dating
In the laboratory, the 15 parallel cores were visually correlated by aligning stratigraphic boundaries and color variations. Subsamples were soaked overnight in 10% KOH and were then sieved using a 0.25-mm sieve. Macrofossils were identified using a stereo microscope under ×25 magnification and extracted with soft tweezers and radiocarbon dated. Eight dates were based on terrestrial macrofossils, which were combined with five bulk sediment samples (Supplemental Table 1, available online). Bulk dates were used where no macrofossils were available. To estimate the reservoir ages for the basin, a bulk sample was dated at the same depth as a terrestrial macrofossil sample (201 cm depth). Calibration of 14C ages and the construction of the age–depth model, which was a basic linear regression between neighboring levels, were carried out in Clam 2.2 and R Workspace using IntCal13 (Reimer et al., 2013). The reservoir age was subtracted from the bulk sediment 14C ages prior to calibration. Ages are reported as cal. yr BP (before 1950 AD).
Organic matter content
%LOI550 (loss on ignition 550) and % carbonate content (%LOI950) determinations (Dean, 1974) were made to estimate TIC and TOC (total inorganic and organic carbon, respectively). For sediments with low carbon values, a multiphase carbon analyzer Leco RC-512 was used.
Geochemistry
The cores were analyzed using an Itrax x-ray fluorescence (XRF) scanner (Croudace et al., 2006). The master core was scanned, using a molybdenum tube set at 30 kV voltage and current of 50 mA. Measurements were taken at 0.2 mm resolution for 15 s each. The obtained dataset includes a total of 15 elements (Br, Ca, Cu, Fe, K, Mn, Ni, Pb, Rb, S, Si, Sr, Ti, Zn, Zr) along with Mo incoherent and coherent scattering signals. Itrax data are uncalibrated (counts or count-rate) and therefore all the values were normalized by Ti which is considered to be a conservative element (Weltje and Tjallingii, 2008). Furthermore, normalization helps to compensate for the influence of the dilution effect of organic matter and physical defects (Löwemark et al., 2011). Average values were calculated for 1 cm intervals for smoothing the graphs.
Diatom analysis
Samples were treated according to the standard method compiled by Battarbee et al. (2002). Counting of diatoms was performed under light microscope using ×1000 magnification and immersion oil. Identifications of diatom valves were mainly based on Cleve-Euler (1951, 1952, 1953, 1955), Mölder and Tynni (1967–1973), Tynni (1975, 1976, 1978), Krammer and Lange-Bertalot (1986, 1988, 1991a, 1991b), Snoeijs (1993), Snoeijs and Balashova (1998), Snoeijs and Kasperoviciene (1996), Snoeijs and Popatova (1995), Krammer (2000, 2002, 2003), Lange-Bertalot (2001), Levkov (2009), Lange-Bertalot et al. (2011), and Levkov et al. (2013, 2016). The aim was to count 300 valves; however, lower basic sums were accepted where diatoms were sparse (lowest 85.5). Identified diatom species were grouped according to their salinity preferences (coastal brackish-marine taxa, planktonic brackish-marine taxa, halophilous taxa, lagoonal taxa, indifferent taxa, freshwater taxa), as well as water depth (coastal, planktonic, and benthic). The diatoms were grouped according to the following literature: Miller (1973), Lie et al. (1983), McQuoid and Hobson (1988), Zong (1997), Laffaille et al. (2002), and Hedenström and Risberg (2003). Since the Baltic Sea has not experienced fully marine conditions, that is,
Tilia 2.0.41 software was used for plotting diagrams. CONISS cluster analysis was employed as an additional tool to eye-matching for defining diatom zones (Grimm, 1987).
Carbon and nitrogen determination
The C/N ratio is useful for indicating the type and source of organic matter (Rice and Hanson, 1984). Variations in C/N ratio can occur because algae are enriched in nitrogen and depleted in carbon compared with vascular plants (Tyson, 1995). C/N ratios of <8 typically indicate marine sediments (Bordovskiy, 1965) and according to García-Alix et al. (2012), a C/N ratio of <10 indicates aquatic organic matter. A C/N ratio of 10–20 indicates a mix of terrestrial and aquatic organic matter inputs and a C/N ratio of >20 signifies that terrestrial organic matter was predominant (Jones et al., 2013). Due to these variations in the C/N ratio of different sediments, C/N ratios can be used as proxy for relative sea level in isolation basins (Mackie et al., 2005). Forty subsamples spaced 7 cm apart at around 1 cm3 volume were subsampled from the gyttja and calcareous gyttja. No subsamples were taken from the fen peat as the source of the organic material was obviously terrestrial. No subsamples were analyzed from the clay as it can be assumed this formed in deep water and to have low organic matter content. The samples were then treated with 10% HCl and stirred in order to break up calcium carbonates, then dried overnight at 50°C in an oven and ground in a marble pestle and mortar. A microbalance was used to weigh 100 μgC–1 mgC into LUDI SWISS ϕ9.0-mm height 10-mm tin foil capsules and sealed with flat-edged tweezers. The C/N ratio measurement samples were made using a CHNS analyzer. The overall uncertainty in determining C/N (
Pollen analysis
A total of 28 samples were prepared and counted with higher resolution in the upper part of the core. A standard procedure for pollen preparation was used following Berglund et al. (1986). Hydrofluoric acid was used on all samples in order to remove silicates. To calculate pollen concentration rates, tablets with a known number of
Results and interpretations
Stratigraphy
In general, the Lina myr basin is filled with glacial and postglacial clay. A thin sandy layer lies on top of the clay, covered by gyttja and fen peat (Figure 2). The gyttja was divided into three layers: two layers of calcareous gyttja separated by a layer of gyttja. In the master core (core 7), the following eight sedimentary facies were identified: clay and sand (0.15 %LOI550 and 1.42 %LOI950 400–339 cm), calcareous gyttja (6.8 %LOI550 and 28.28 %LOI950 339–307 cm), gyttja (26.58 %LOI550 and 3.46 %LOI950 307–79 cm), calcareous gyttja (21.38 %LOI550 and 30.95 %LOI950 79–57 cm) and fen peat (84.32 %LOI550 and 2.33 %LOI950 57–40 cm). The upper 40 cm of the sequence was not used for analysis since the top surface of the fen was reworked due to plowing.
Chronology
Thirteen radiocarbon dates, ranging between 9102 and 2357 cal. yr BP (Supplemental Table 1, available online), were used to construct the age-depth model (Figure 3). At 201 cm depth, the date of a terrestrial macrofossil (5567 ± 51 14C yr BP) was compared with the date of a bulk gyttja sample (6018 ± 32 14C yr BP) and the 451-year age difference between the samples was used to estimate the reservoir effect. This reservoir age was subtracted from all bulk sediment dates prior to calibration. Two dates (Ua-54786 and Ua-54320) were rejected from the age-depth model since they were regarded as outliers (see Supplemental material, available online, for ages).

Age-depth model Clam 2.2 and R Workspace IntCal 13.14C atmospheric curve for the Northern hemisphere (Reimer et al., 2013). The red points indicate outliers and the blue boxes indicate bulk dates which have had the reservoir age subtracted to them prior to calibration.
XRF geochemistry
When Br is compared with Ti, which is considered a relatively depleted element in the coarser-grained fractions of sediments, Br/Ti can be used as a proxy for marine organic matter (de Boer et al., 2014; Ziegler et al., 2008). Br/Ti has been used to indicate the ingress of marine water into coastal lakes caused by cyclone activity (Oliva et al., 2018). It could be argued that longer term changes, such as marine transgressions, could also be detected by fluctuations in the Br/Ti ratio. There were three apparent peaks in the Br/Ti for Lina myr – between 275 and 245 cm depth, culminating around 265 cm depth, between 165 and 155 cm and a sustained increase in the calcareous gyttja between 57 and 79 cm (Figure 4). The Ca/Ti ratio can be used to show relative changes in calcium carbonate content (Piva et al., 2008). The Ca/Ti ratio peaks in the upper and lower calcareous gyttja and mirrors LOI950.

The master core stratigraphy, loss on ignition, Br/Ti, Ca/Ti, C/N, brackish-marine planktonic diatom species share (including
Carbon-to-nitrogen ratios
The C/N was relatively high for the lower calcareous gyttja (C/N ratio between 11 and 17) indicating mixed sources of organic matter. The ratio started to shift toward lower values (C/N ratio of 9) after the transition to gyttja; however, the ratio did not fall below 7, which would indicate predominance of marine organic matter (Emerson and Hedges, 1988). After the initial shift to lower C/N ratios in the gyttja, the ratio shifted to a higher value of 12 and then again to a lower value of 9 a total of three times. These three shifts to lower values at 208, 173, and 117 cm depth in the stratigraphy may indicate increases in marine influence. The C/N ratio increased again to 25 at 89 cm depth with the boundary of the gyttja and upper calcareous gyttja deposits but then decreased again to 11 indicating a shift from terrigenous organic matter in the lake to a mixed source of organic matter.
Pollen record
The pollen diagram (Figure 5) was dominated by arboreal taxa such as
PAZ 1 (330–260 cm, 8850–7440 cal. yr BP, 6900−5490 BC)
Arboreal taxa of
PAZ 2 (260–190 cm, 7440–6210 cal. yr BP, 5490–4260 BC)
PAZ 3 (190–97 cm, 6210–4680 cal. yr BP, 4260–2730 BC)
The most dominant taxa in this zone were
PAZ 4 (97–40 cm, 4680–2350 cal. yr BP, 2730–400 BC)
The onset of PAZ 4 was defined by the introduction of the human impact indicator taxa

Pollen percentage diagram for Lina myr.
Diatom record
A total of 163 diatom taxa were identified and seven diatom assemblage zones (DAZ 1–7) were established (Figure 6). No diatoms were observed between 400 and 338 cm where the sediment consisted of clay and sand.
DAZ 1 (339–307 cm, 9300–8500 cal. yr BP, 7350–6550 BC)
The diatom composition was characterized by freshwater benthic taxa. Dominant genera are
DAZ 2 (307–290 cm, 8500–8050 cal. yr BP, 6550–6100 BC)
This was based on one diatom assemblage and was dominated by taxa with a wide range of salinity preferences, for example,
DAZ 3 (290–275 cm, 8050–7700 cal. yr BP, 6100–5750 BC)
This was characterized by a predominance of the coastal brackish-marine diatom
DAZ 4 (275–202 cm, 7700–6400 cal. yr BP, 5750–4450 BC)
This displayed even shares of brackish-marine and halophilous taxa. However,
DAZ 5 (202–150 cm, 6400–5650 cal. yr BP, 4450–3700 BC)
This was dominated by coastal brackish-marine taxa reaching
DAZ 6 (150–79 cm, 5650–3800 cal. yr BP, 3700–1850 BC)
This was dominated by coastal brackish-marine taxa reaching
DAZ 7 (79–53 cm, 3800–2600 cal. yr BP, 1850–650 BC)
This was characterized by low diatom abundance. Freshwater benthic species, for example,

The diatom summary diagram. Dominating taxa for the various groups are:
Discussion
Shore displacement reconstruction
The shore displacement curve shown in Figure 7 illustrates nine phases based on stratigraphy and diatoms combined with geochemical parameters. The aim of this curve is to tentatively reconstruct past relative sea-level fluctuations. The shore displacement curve is based on the elevation of the Lina myr threshold and the highest shoreline for the Littorina Sea (22 m). This highest shoreline is based on the occurrence of beach ridges and wave washed material found 21–22 m a.s.l. close to Lina myr.

The tentative shore displacement curve constructed from investigations of the master core from Lina myr and comparisons for other studies around the Baltic Sea (modified from Barliaev, 2017). L1, L2, and L3 represent transgressional sub-phases of the Littorina Sea.
Stratigraphic data from phase I (prior to 8980 cal. yr BP) indicated that Lina myr was submerged beneath the Ancylus Lake. Thin layers of sand above the clay indicated reworking of glacial till through wave action and possible erosion. According to Svensson (1989), the Ancylus regression was rapid in the beginning (5–10 m) followed by a slowdown. Svensson (1989) also stated that the first isolation of the Lina myr basin, after the Ancylus Lake stage, occurred around 8600 radiocarbon yr BP, that is,
Phase II (8980–8500 cal. yr BP) represented an isolation of the basin from the Ancylus Lake allowing the formation of calcareous gyttja, which was accumulated in a small freshwater lake, as indicated by
Diatom analysis from phase III (8500–8050 cal. yr BP) suggested that the lower portion of the gyttja represents a lagoonal stage with brackish conditions (Miller, 1986). This was supported by a decrease in the C/N ratio, which indicated a shift toward marine influence. The water depth and salinity were still low, as non-brackish and benthic diatoms dominated; this was also reflected in the relatively low Br/Ti ratio. This interpretation is supported by a peak in lagoonal taxa, for example,
The timing of the onset of the Littorina transgression is in agreement with L1 in Blekinge (Berglund et al., 2005). It also coincides with the initiation of the transgression at Narva-Luga lowland and Pärnu area (Rosentau et al., 2013; Sandgren et al., 2004; Veski et al., 2005).
The initial period of phase IV (8050–7700 cal. yr BP) was characterized by the predominance of
Throughout phase V (7700–6400 cal. yr BP), C/N ratios were relatively low, indicating a stronger marine influence.
The L2 sub-phase was followed by a decrease in planktonic diatoms and in the Br/Ti ratio. This corresponded to a possible regression as the sea-level rise was overtaken by postglacial isostatic uplift. We find that a regression took place between 7400 and 7000 cal. yr BP, which coincides with still stand conditions at the Karelian Isthmus (Miettinen et al., 2007) and a regression in Blekinge (Berglund et al., 2005) and the eastern Gulf of Finland (Rosentau et al., 2013).
After L2, low proportions of planktonic brackish-marine taxa indicate a low sea level. The relatively low Br/Ti ratios support a decrease in marine influences. It is, however, impossible to estimate the absolute water depth lowering.
The dominance of brackish-marine planktonic diatoms and a peak in the Br/Ti ratio in phase VI (6700–6200 cal. yr BP) may indicate an increase in water salinity and an increase in relative water depth around 6100 cal. yr BP. C/N ratios were relatively low at this time supporting the premise that aquatic organic matter input was relatively higher. This increase in salinity is in agreement with Westman et al. (1999), who suggested the highest salinity conditions (about 20‰) of the Littorina Sea between 6500 and 5000 cal. yr BP. This peak was caused by an increased volume of saline water inflow from the Atlantic as a result of the enlarged cross-section of the Danish straits (Gustafsson and Westman, 2002; Westman et al., 1999). As the study lacks any proxy for direct estimation of past water depths, it is difficult to assess whether the Littorina limit (i.e. 22 m a.s.l. at Väänä Lagoon NW Estonia, (Saarse et al., 2009a) was reached during the L2 or L3 sub-phase.
The L3 sub-phase, recorded in the older part of phase VI, was followed by a regression, which led to the eventual isolation of the basin. No evidence of minor transgressions after 6100 cal. yr BP was found, despite that indications of minor transgressions have been detected in eastern Svealand and Blekinge, Sweden (Berglund et al., 2005; Risberg et al., 2005). Planktonic brackish-marine diatoms were gradually replaced by coastal brackish-marine taxa indicating salinity remained relatively high, whereas water depth was falling. This supports the premise that the global eustatic sea level was stable, while isostatic land uplift became the main factor in shoreline displacement on Gotland.
Phase VII (5700–3800 cal. yr BP) was characterized by the dominance of coastal brackish-marine diatom taxa as a result of isostatic uplift and basin infilling. The final isolation from the Littorina Sea occurred 3800 cal. yr BP when the sea level dropped below the Lina myr basin threshold, which had probably been eroded and was lower than during the previous isolation event. The younger freshwater stage, phase VIII (3800–2700 cal. yr BP), was characterized by a lower water level due to erosion of the threshold and further sediment infilling of the basin. The calcareous gyttja deposits which accumulated during phase II were thicker and geographically more widespread than the more recent calcareous gyttja deposits. This indicated that the older isolated freshwater lake was larger than the younger isolated phase. The scarcity of diatoms frustules circumstantially indicates a high pH environment since diatoms are often poorly preserved in high pH environments (Flower and Ryves, 2009). An increase in the C/N ratio indicates a shift to higher inputs of terrestrial organic matter. The relatively high Br/Ti ratio may be caused by an increase in organic matter (Kalugin et al., 2013) rather than marine water input. As the basin shallowed over time, due to continuing input of sediments, the lake became overgrown. Accumulation of fen peat began 2700 cal. yr BP at the sampling site, which is located in the middle of the fen.
Phase IX (2700 cal. yr BP until present) is represented by the formation of a fen dominated by

Summary of findings. L1, L2, and L3 represent transgressive sub-phases of the Littorina Sea.
Vegetation reconstruction and human impacts
Early humans and the environment
The first human settlements on Gotland are from about 9200 cal. yr BP (Apel et al., 2018). Population numbers are interpreted to be generally low during the Mesolithic, with periods of variation in settlement intensity and probably a marked decline during the Late Mesolithic. During part of this period, there is evidence of hunting for marine mammals (Apel and Storå, 2017). One of few Mesolithic sites found on Gotland is at Svalings, close to Lina myr. The site was dated to about 7000 cal. yr BP by radiocarbon dating of seal bones (Apel et al., 2018). Seal bones are prone to having an inbuilt marine reservoir age; the study by Apel et al. (2018) attempted to correct for this by subtracting 90 years from the uncalibrated ages. As seen from the C/N ratio and diatom data in the present study, part of the Mesolithic period is characterized by increased marine influence (diatom phase V: 7700−6400 cal. yr BP), and a sub-phase of the Littorina transgression. Hunter-gatherer populations are difficult to detect in pollen records since they generally have low impact on their surrounding vegetation. The first Early Neolithic finds on Gotland are dated to about 5800 cal. yr BP; however, the Early Neolithic culture could have been introduced about 200 years earlier, around 6000 cal. yr BP (Apel et al., 2018). The Early Neolithic finds belong to the Funnel Beaker Culture (FBC), which is the first farming society in Northern Europe and Scandinavia. The FBC complex, typical for the Early Neolithic period, lasts until about 5000 cal. yr. BP on Gotland (Apel et al., 2018). Studies from Scania, southern Sweden (Gron et al., 2015), and Mälardalen, eastern central Sweden (Isaksson and Hallgren, 2012), show that cattle were used for dairy farming during this period. The earliest evidence for cattle farming in Sweden comes from Almhov, Scania, and is dated to 5900−5450 cal. yr BP (3950−3500 cal. BC). It is believed that the Early Neolithic culture on Gotland had connections with Southern Scandinavia based on finds such as imported flint axes, amber, and possibly even pottery (Österholm, 1989). During the Early Neolithic, pigs were brought to Gotland; however, bones from sheep/goat and cattle are first found during the Middle Neolithic (Apel et al., 2018).
Gotland’s first collective burial site, from the late Early Neolithic is located at Ansarve in eastern Gotland. Analyses of bone material dated to 5300−4900 cal. yr BP indicates mainly terrestrial food intake for the studied individuals (Apel et al., 2018). Terrestrial diets may have been more important than marine at this time; however, impacts on vegetation are not evident in the pollen record.
Human impacts on vegetation were first seen clearly in the Lina myr pollen diagram at 4580 cal. yr BP (lower part of PAZ 4). This is after pigs are thought to have been introduced to Gotland and after the first identification of dairy farming in Scandinavia (Gron et al., 2015). Since the pollen record from Lina myr is interpreted to represent general regional vegetation change for large parts of Gotland, this indicates that humans began to change the general vegetation patterns on the island from the Middle Neolithic period (5000 cal. BP to 4200 cal. yr BP), associated with the Pitted Ware Culture (PWC). According to frequency distribution of dates for the PWC culture, this culture is linked to a marked population intensification on Gotland (Apel et al., 2018). The most important pollen taxa indicating human impact in PAZ 4 is
Within PAZ 4,
Lina myr was connected to the Littorina Sea until the Late Neolithic when it became a lake at about 3820 cal. yr BP. It is possible that Lina myr was once part of an inland water system; however, this must then have been before 3820 cal. yr BP and hence during the Mesolithic and Neolithic BP and hence during the Mesolithic and Neolithic periods. Gothemshammar (dated to around 2850−2450 cal. yr BP) and other Bronze Age finds in the area such as Tjelvar’s grave (dated to 3056−2941 cal. yr BP) and a location named ‘Killingborgen’ close by were likely used as ritual sites. The large cairn Majsterrojr just to the south of Tjelvar’s grave is a type of structure built in the Early Bronze Age (around 3500 cal. yr BP) but probably served as a ritual site for several thousand years (Martinsson-Wallin and Wehlin, 2017). These structures were built after the isolation of the lake from the Littorina Sea. After the isolation, the basin was a lake for about 1000 years, until the Late Bronze Age, when it became a fen. It is likely that the relatively large lake that filled the Lina myr basin was a good resource for fish and freshwater for the Bronze Age population around Lina myr and the Gothem River. The Gothem River connected the lake with the Sea to the east and the Hörsne River. During the early part of the Middle Bronze Age, at 2920 cal. yr BP, pollen grains of
Conclusion
In the earliest part of the Lina myr record (9000−8980 cal. yr BP), the basin was submerged by the Ancylus Lake. At about 8980 cal. yr BP, the basin threshold located at the minimum altitude of 12 m a.s.l. emerged and the basin was isolated for the first time. Between 9500 and 8500 cal. yr BP, a freshwater lake with a maximum depth of 5 m and highly reducing conditions existed in the Lina myr basin. The basin became a bay of the Littorina Sea during a transgression, which initiated at about 8500 cal. yr BP. Two major sub-phases of the Littorina transgression could be identified. The first lasted from 8100 to 7500 cal. yr BP (L2), and the second from 6500 to 6000 cal. yr BP (L3). The highest Littorina shoreline was formed during one of these sub-phases. The highest salinity was reached during the second sub-phase (L3). This study agrees about the timing of the Littorina transgression with other studies with similar patterns of postglacial isostatic uplift. These findings suggest that these events were trans-Baltic. From 6000 cal. yr BP, the area experienced a regression that eventually led to the second isolation of the basin at about 3820 cal. yr BP. The following lake phase lasted for about 1000 years. During this time the lake slowly became shallower and at 2700 cal. yr BP also the central parts of the basin had transformed into a fen. Since the Lina myr basin is relatively large, the pollen record from the Lina myr core mainly represents the vegetation of most of Gotland. Prior to the Neolithic, most vegetation changes were due to ecological factors such as adaptation of plants to changing environmental conditions due to, for example, changing sea levels in connection with the Littorina transgressions. At 4580 cal. yr BP, in the Middle Neolithic, barley, a clear indicator of cultivation, is observed for the first time in the pollen record. Barley continues to be present in the pollen record together with
It is clear from the archeological record that sites surrounding Lina myr were of great importance to Gotland’s prehistoric inhabitants. Perhaps this importance was due to its coastal location at this time which may explain why stone ship settings are orientated in a different way to stone ship settings on the rest of Gotland; the bow of many of these ships would have been oriented toward the modern day location of Lina myr. Although people have been in this area for around 7000 years, no obvious human impacts were seen in the vegetation record in the Mesolithic and Early Neolithic. Compared with the known ages of archeological sites, there was a lag time between the arrival of pigs to Gotland (Early Neolithic) and farming or animal husbandry to the opening up of the landscape at Lina myr. The intensification of land-use in the Middle and Late Neolithic may be due to increased populations and the pressure to produce food.
Supplemental Material
Diatom_diagram_supplementary_data-01-01 – Supplemental material for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP
Supplemental material, Diatom_diagram_supplementary_data-01-01 for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP by Nichola Ann Strandberg, Aleftin Barliaev, Helene Martinsson-Wallin, Jan Risberg, Martina Hättestrand, Ian Croudace, Malin Kylander and Yusuke Yokoyama in The Holocene
Supplemental Material
Diatom_groupings_with_reference_list – Supplemental material for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP
Supplemental material, Diatom_groupings_with_reference_list for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP by Nichola Ann Strandberg, Aleftin Barliaev, Helene Martinsson-Wallin, Jan Risberg, Martina Hättestrand, Ian Croudace, Malin Kylander and Yusuke Yokoyama in The Holocene
Supplemental Material
Supplemenary_documents_most_common_diatom_taxa-01 – Supplemental material for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP
Supplemental material, Supplemenary_documents_most_common_diatom_taxa-01 for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP by Nichola Ann Strandberg, Aleftin Barliaev, Helene Martinsson-Wallin, Jan Risberg, Martina Hättestrand, Ian Croudace, Malin Kylander and Yusuke Yokoyama in The Holocene
Supplemental Material
Table_1 – Supplemental material for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP
Supplemental material, Table_1 for Landscape development at Lina myr fen, Eastern Gotland, 9000−2500 cal. yr BP by Nichola Ann Strandberg, Aleftin Barliaev, Helene Martinsson-Wallin, Jan Risberg, Martina Hättestrand, Ian Croudace, Malin Kylander and Yusuke Yokoyama in The Holocene
Footnotes
Funding
Supplemental material
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
