A Mega-Tsunami in the Baltic Sea 1171 BC: Geological Records with Special Reference to the Lake Mälaren Area in Sweden

At about 3000 C14-year BP or 1200 cal. yrs BC, the Baltic Sea experienced a mega-tsunami with a wave-height of 10 m or more, and a run-up height of up to 16.5 m. This event had significant geological and archaeological effects. We explore the records from the Lake Mälaren area in Sweden. The tsunami event is linked to seismic ground shaking and methane venting tectonics at several sites. The triggering factor is proposed to be the Kaali meteor impact in Estonia of the same age. The documentation of a mega-tsunami in the middle of the Bronze Age has wide implications both in geology and in archaeology. The archaeological key sites at Annelund and Apalle are reinter-preted in terms of tsunami wave actions remodelling stratigraphy. By extensive coring, we are able to trace the tsunami effects in both off-shore and on-shore environment. At the time of the event, sea level was at +15 m (due to isostatic uplift). The tsunami wave erosion is traced 13.5 m below sea level. The tsunami run-up over land is traced to +29.5 m to +31.5 m (occasionally even higher), implying a run-up of 14.5 - 16.5 m. In Ångermanland, the tsunami event was absolutely dated at 1171 varve years BC. Archaeologically, the tsunami event coincides well with the transition between Periods II and III of the South Scandinavian Bronze Age. Period III has traditionally been difficult to identify in the cultural materials of the Lake Mälaren region.


Introduction
Today, Sweden has a low to moderately low seismic activity. For a long time, this was also considered to be the case for the entire postglacial period of Sweden and Fennoscandia. This idea totally changed in the late 20 th century, when it became clear that Sweden was a high-seismic region in deglacial time when the rate of glacial isostasy peaked [1]. Thanks to intensive field work this was later elevated to a new paradigm [2], documenting 52 separate paleoseismic events in the last 12,000 years. Each event was documented by multiple criteria (lending internal support and control to each of the events; [3]. Tsunami events were recorded at an increasing number of sites [2] [3] [4] [5] [6] up to 17 events [7] [8] [9]. The tsunami wave-heights observed provided an independent method of estimating the seismic magnitude of ground shaking [9]. This implies that tsunami events have become a natural process in the postglacial evolution of the Baltic and the Kattegat Seas. Originally, all the events at about 3000 -4000 BP were believed to represent separate paleoseismic events generating local tsunamis [2]. The first thing pointing to something different was that the events recorded at Hudiksvall and Forsmark could be combined into one single event occurring 2900 C14-yrs BP [8]. Not until recently, did it become obvious that there was one big tsunami event at 3000 BP or 1200 cal. yrs BC (Figure 1), which was Baltic-wide and was initiated by the Kaali impact event [10] [11]. This is where the present paper starts.
The observed tsunami wave-height (or run-up) is in the order of 10 -16 m (>6 to 20.5 m), which corresponds to a seismic ground shaking of magnitude ranging from M 7.8 to 8.6. The Swedish data set is given in Figure 2 on the basis of observed relation between tsunami wave height and earthquake magnitude [9] (here redrawn with one Baltic-wide event at around 3000 BP).
The methods applied are sediment coring, sediment characteristics, radiocarbon dating (the AMS technique) and analysis of relevant literature.
Sediment coring was performed at 24 sites (totalling 52 cores): with the majority of 24 sites (A-Y) in the Enköping-Bålsta area. The sediment characteristic is a prime factor for the identification of tsunami events (e.g. [27]). In this case, it was important to observe, document and interpret both the on-land signals and the off-shore signals. The core taken above the +15 m shore level at 3000 BP usually contained tsunamites starting with a distinct erosional surface to the clay below (sharp and irregular). Often erosion could have gone down into the varved clay from the deglaciation time, implying that 7000 -8000 years of sediment was removed. The tsunamite itself consists of a sandy layer of graded bedding (fining upwards) or unsorted beds with grains lacking rounding. The highest altitude of their occurrence provides a measure of the tsunami wave-height and run-up capacity ( Figure 1). The tsunamite found in the off-shore environment below +15 m had a quite different appearance. Usually, it was a bed of pure fine clay from the settling of material set in suspension by the base of the tsuna-mi wave hitting the seabed and generating erosion. This was understood by the occurrence of lagoonal gyttja-clay below the fine clay. This surprising stratigraphy was observed by [17] [18] [19] [20], but not understood. The tsunamite characteristics are illustrated in Figure 3 and were used as a key to the understanding of the stratigraphy observed in the sediments cores collected (cf. Figure 5). Figure 3. Characteristics of tsunami beds (tsunamites) in the on-shore area (left) and off-shore area (right) with respect to the +15 m sea level at the time of the tsunami event.

The Lake Mälaren Valley
Lake Mälaren Valley (Swed: Mälardalen) is a broad lowland valley stretching 200 km inland from Stockolm and including both Lake Mälaren and Lake Hjälmaren. The location of the sites (3.1 to 3.10) here discussed are given in Figure 4.
Our collaboration on the wider documentation of the 3000 BP tsunami in geological and archaeological data started with the Tolan-project in Upplands-Väsby [24]. The first target was the archaeological investigation at Annelund at Enköping (15), because at this site a "silt-layer" of unknown origin had caused much confusion. It seemed natural to suspect that the "silt-layer" at Annelund might be a tsunami layer. This idea initiated an extensive coring campaign ( Figure 5).      Sävsta Kärr (coring site B) refers to a swamp and bog at an elevation of +18 m ( Figure 6). It occupies a depression ( Figure 5) with an outlet threshold at about +20 m. Coring was performed at two nearby sites. A quite typical tsunamite was retrieved between lacustrine gyttja and Baltic clay. The tsunamite is 2.8 cm thick and consists of a clayey silty sand with a strong erosional contact to the clay below, which is a varved clay from the deglaciation period, implying that about 7000 years' Baltic sediments were removed by the erosion of the tsunami event ( Figure 7).
From these two sites we got stratigraphic evidence of a tsunami ( Figure 6) and [24] and time coincidence with the Kaali impact in Estonia [10] [11]. It also worked as the incitement of the Enköping-Bålsta investigation here presented.

The Annelund Archaeological Site at Enköping
In the eastern part of the city of Enköping, an archaeological site ( Figure 3, site 1 in Figure 4) was discovered in connection with the construction of a new highway at Annelund, and documented [15]. Above archaeological material from the Early Bronze Age and partly also above a gallery grave from the Late Neolithic, but below archaeological material from the Late Bronze Age and the Iron Age, there is a layer termed the "silt layer" [15]. It has a wide lateral extension ( Figure   8). It is 5 -30 cm thick and reaches up to about 28.5 m. Occasionally, it contains lenses of gravel and stones of 5 to 20 cm size. It exhibits flow structures and includes Baltic diatoms. Obviously, it is a Baltic Sea deposit accumulated under quite strong water flow forces. Its origin remained unsolved, however, which posed major interpretational problems for the stratigraphy and chronology [15].
We now propose that the so-called "silt layer" (red in Figure 8) constitutes a sedimentary bed laid down during a tsunami event. If so, it is a "tsunamite" [27]. The fact that the "silt layer" contains lenses of gravel and stones ranging from 5 to 20 cm is consistent with a rapid and violent mode of deposition. Figure 8. Two profiles from the archaeological site Annelund: "silt-layer" red, cultural layer dark, sand yellow and clay blue. The silt-layer reaches up to +28.5 m in (a), is quite extensive in (b), and rests upon cultural layers from the late Neolithic and early Bronze Age (redrawn from [15]. At about 3000 BP sea level was at about +15.0 m. As the "silt bed" reaches up to about +28.5 m in Annelund, the tsunami wave-height plus tsunami run-up must have been at least 13.5 m as illustrated in the local sea level curve ( Figure 9). Figure 9. Sea level changes observed. At about 3000 BP, sea level was at +15 m, and the area was struck by a mega-tsunami with w run-up of 13.5 to 16.5 m in the Enköping-Bålsta area.
By an intensive coring program (52 sediment cores at 24 sites, Figure 5), the tsunami event was well documented by the occurrence of "tsunamites" ranging environmental conditions. This level is dated at 3060 ± 100 C14-years BP, which makes a correlation with the tsunami event possible.
The Annelund archaeological site [15] deserves a deeper analysis in view of our suggestion that the "silt layer", in fact, is a tsunami layer deposited at a Baltic-wide tsunami event occurring at about 3000 BP.
The silt layer overlies the Neolithic cultural layer including a gallery grave of Late Neolithic age and a "bone-pit" (A6) at +27.0 -29.1 m. The silt layer extends at about +26 -27 m, and occasionally goes up to 28.0 -28.5 m. The silt layer is said to become sandier to the east.
The gallery grave is at its edge covered by a relatively thick layer of littoral sand and gravel (according to [15]). Its level is about +28.5 m. Because the underlying gallery grave represents land conditions in Late Neolithic time, the littoral sand must represent an extreme storm episode or, in our view, a tsunami event. In the gallery grave, there is a gravel filling (including some bones) above the original land-surface at +28.5 m.
The littoral sand and gravel separates the gallery grave from a covering mound of brittle burned stones (A1). On the whole there are four, maybe five, mounds of brittle-burned stones. There are also three house-foundations marked by postholes.
Besides the identification of the silt-layer as a tsunamite, the Annelund stratigraphy calls for an additional archaeological revision beyond the scope of this paper.
There are a large number of irregular pits scattered over the investigation area.
They are hard to explain and do not seem to fit any known function. Therefore, it has been suggested that they may represent the remains from root systems of removed trees. This is very interesting in view of the tsunami event proposed.
An instantaneous flooding of the ground during a tsunami, may have left standing dead trees on the ground, which the people removed, hence creating the irregular pits scattered all over the area (in Fanna Mosse, there are palynological evidence [15] [28] that light demanding trees and spruce decreased just at this time, and grass increased; cf. above).
Unfortunately, we have no C14-dates from the littoral sand and gravel layer covering the gallery grave or from the gravel filling in the gallery grave, which contained skeleton remains from a few persons. Stratigraphy and sediment characteristics clearly show that the area was over-washed by a tsunami in the mid Bronze Age, however.   Figure 12).

The Boglösa Area
The Boglösa area lies 5 km to the south of Annelund (Figure 4). In a small bog situated at +25.5 m our third coring was performed (Site C) and a very clear tsunami layer was retrieved: a 3.2 cm layer of sand and gravel with fine sand in the uppermost 0.5 cm (i.e. a graded bedding as shown in Figure 3 left) between Baltic clay and lacustrine gyttja (only 6.5 cm thick) and peat. The contact to the clay was strongly erosive and irregular, and cut down into the glacially varved clay (implying the removal of 7000 -8000 years of sediment). The grains were angular and badly sorted precluding a littoral shore formation. This very clear tsunami record initiated extended work: with a total of six coring points in Site C, and additional Sites D, E and F ( Figure 10).
The tsunami is recorded in all six cores in Site C. It starts with an irregular erosion surface reaching down into the glacially varved clay (implying an extensive hiatus). The tsunami layer is of different character and thickness in the six cores: C1: 3.2 cm sand and gravel to fine sand in graded bedding. C2: erosional hiatus. C3: 2.3 cm gravel followed by 18.2 cm of alternating organic and silty layers. C4: 5.5 cm gravelly sand. C5: 5.7 cm gravelly sand. C6: and 1.7 -3.2 cm of silty sandy gravel as shown in Figure 11. Core D at +24.0 m showed a strongly weathered clay surface (i.e. land surface) covered by a thin layer of brushwood peat. Core E at +29.0 m had a strong erosional contact at the top of a hard clay. Core F at +29.5 m (two core sites) shows a basal clay with a strong erosional surface at its top, 2.8 cm of sandy gravel with angular stones (the tsunami bed), 27.0 cm of lake gyttja, and 34.0 cm of peat. The Core F stratigraphy indicates that the tsunami wave reached up to, at least, 29.5 m (i.e. 1.0 m higher than the silt-layer in Annelund). Two C14-dates were performed: black wood just in top of the tsunami layer gave 335 ± 29 BP (1560 ± 90 AD) and 1.2 cm clay gyttja directly above gave 1343 ± 30 BP (667 ± 17 AD), indicating far too young ages (and in reversed stratigraphic order) due to roots and humic acids.
In the Boglösa area, there is a famous and quite unique rock carving from the Bronze Age [29] [30], known as "Brandskogsskeppet". The ship is 4.2 m long.
Judging from the size of the men (0.7 -0.8 m above their waistlines) and their paddles (about 1 m long), the ship is likely to be drawn more or less in natural size ( Figure 12). The elevation of the rock carving was leveled at "between 28 and 29 m above sea level" [29] (or about +28.5 -29.5 m in the present reference system; p. 31). Consequently, it matched the run-up level of the tsunami wave (+28.5 m in Annelund and +29.5 m in Boglösa) perfectly well. This may be a coincidence (BA & PAA) or a causal linkage to the tsunami event (view here proposed by NAM). In the second case, the place may mark the place of wreckage and death, or the place of a remarkable salvage of ship and men. From a tsunami dynamic point of view, wreckage is likely and salvage possible.
If the men were out fishing in the open water to the south when the tsunami hit the area, the ship may just have been lifted up (as was the case in the Mal-dives at the 2004 tsunami; [7]) and moved ashore at the limit of the run-up of the tsunami wave at the hill slope just at or close to where the carving was cut ( Figure 13). Consequently, the carving may illustrate a miraculous rescue from the tsunami disaster. Figure 12. Brandskogsskeppet, a 4.2 m long rock carving from the Bronze Age (based on [30]). The carving was leveled by [29] at between +28.0 -29.0 m, which implies that the elevation corresponds to the tsunam run-up level. This may be a coincidence or rather a reason to seek a causal connection. Figure 13. Synthesis of the Enköping tsunami data. Due to the glacial isostatic uplift the shore is continually displaced downwards. When sea is at about +15 m at about 3000 BP, the mega-tsunami strikes the area, and the tsunami wave instantaneously floods the area up to about +28.5 -29.0 m. The Brandskogsskeppet rock carving lies right at the tsunami run-up level and may commemorate a wreckage or rather the miraculous salvage of ship and men (provided they were out on the open water to the south, from where they may have been "lifted" ashore).  Figure 13 gives a summary of the tsunami records obtained from Sites A to F, Annelund and the "Brandskogsskeppet" rock carving. The shore was successively displaced downwards due to the glacial isostatic uplift of the area. At about 3000 BP, sea was at +15 m. At that moment, the region was struck by a mega-tsunami with a wave-height plus run-up height of about 13.5 m (Figure 9).
In the off-shore environment south and southeast of Boglösa eight sediment cores were taken ( Figure 5; G-M, S). Core J at +14.5 m shows a 15.0 cm tsunamite of silt with angular stones between 48 cm of weathered fine clay above and >30 cm of weathered glacial clay below. This means that it has a characteristic of an on-shore to shore tsunamite ( Figure 3). The other cores have typical off-shore tsunami characteristics (Figure 14), where the tsunami is represented by a layer of fine clay from the settling of particles set in suspension by the erosion of the clayey seabed surface ( Figure 3, right side). It implies that the erosion by the base of the passing tsunami wave has gone down, at least to +0.5 m (H) to +1.5 m (G), i.e. 14.5 m to 13.5 m below the sea level at +15 m. This is as much below sea level ( Figure 14) as was recorded above sea level ( Figure 12), which fits well with the theoretical symmetry of a tsunami wave (cf. 5.1, below). Figure 14. Tsunami records from the off-shore environment, indicating basal wave trimming of the seabed at least down to +0.5 to +1.5 m, i.e. 14.5 to 13.5 m below sea level, and generating a submarine "cloud" of suspended matter, deposited as a fine clay (violet) on top of shallow-water gyttja-clay (greenish grey).
The off-shore characteristics of a tsunami event in the form of erosion of the seabed ( Figure 3, right hand image) fit the model of [7] quite well.
In the local map descriptions [17] [18] [19] [20], it had been observed that there occasionally occurred fine clay on top of gyttja clay. They noted it as "unusual" and were not able to present any tenable interpretation. They undertook a number of grain-size analyses, however, which are quite informative. The covering fine clay consists of >50% clay, <50% silt and organic matter <2%, i.e. a composition fitting the deposition of suspended matter. The gyttja clay below consists of 4% -5% sand (even course sand particles), <50% silt, <50% clay and with organic matter >2% (up to 4% -5%), implying a combination of bed-load transport and suspended matter flocculation (in agreement with our coring data). The C14-dates obtained from Sites Q and R ( Figure 14) give far too old ages, which implies that they both represent reworked and redeposited organic matter (the 13 C values indicate brackish environment). The C14-date from Site S was 2396 ± 29 C14-years BP or 455 ± 50 cal. yrs BC, an age which corresponds to the climatic changes at the Sub-Boreal to Sub-Atlantic boundary; viz. a short period of intensive drought [31] followed by a general change to cooler and wetter climate conditions [32]. The site was isolated from the Baltic at about 1050 AD. A 2-point sedimentation rate graph is presented in Figure 15. It proposes that the tsunami ended at 1175 BC at the fine clay/gyttja clay boundary, an age which corresponds well with other dates of the mega-tsunami event (Appendix).

The Håbo Area
In this paper, the Håbo area refers to the hilly area between the Ekolsundsviken Bay and the Ullfjärden lakes hosting the villages of Övergran, Yttergran and Bålsta ( Figure 16). The Uppsala Esker passes on its eastern border. Two cores were taken in the southern lowland (Q and R; Figure 14), four sites were cored in the central valley (T-X) and one (Y) on the northern outlet side at the tsunami event ( Figure 17). The area seems to have been quite densely populated in the Bronze Age with an archaeological key-site at Apalle (point 3 in Figure 16). International Journal of Geosciences   The base of the central valley is at about +13.5 m, implying that a shallow lake was isolated in the area around Site T. This former lake stage is here called "Ancient Lake Gran" (Figure 18). It was more or less level with the sea at +15 m at 3000 BP. The tsunami invaded the valley basin from the west (with threshold at about +23 m), the south (with threshold at about +26 m), and the east (with threshold in the order of +15 -20 m). The tsunami seems to have over-washed the threshold in the north (at about +27 -28 m). The main outlet of the collapsing water masses trapped within the basin must have escaped via the eastern outlet ( Figure 16). Thick layers of fine clay in Sites T, U and X represent the settling of fine particles set in suspension by the tsunami wave. In Site U, the erosion has gone down into the glacial clay, leaving a hiatus of about 7000 years, indicating that the valley floor was swept clean of vegetation, soil and a substantial package of sediments.
Core T was taken in the middle of Ancient Lake Gran ( Figure 18) and exhibits a key record in stratigraphy and chronology ( Figure 19). It begins with a banded Baltic clay (probably a Littorina Sea clay, and not a glacial varved clay), which at its upper limit has an erosional surface cutting down into the clay. It is covered by 2.5 cm of gyttja silt representing the base of a tsunamite. It is covered by 19.5 cm gyttja and clay gyttja, and finally 127.5 cm of fine clay (from the settling of suspended particles in later phase of the tsunami event). The question is what the gyttja bed represents; 1) a normal lake sequence or 2) the final part of an initial tsunami phase of graded bedding starting with an erosion surface and followed by gyttja silt, gyttja and clay gyttja. The second alternative is here proposed. Figure 18. A local lake, termed "Ancient Lake Gran", was isolated between the valley and the esker when sea level was at +15 m at about 3000 BP, when the tsunami event took place.  The tsunami at Site T ( Figure 19) starts with erosion of the seabed surface followed by the deposition of 2.5 cm silt (bedload), 19.5 cm of gyttja and clay (a mixture of re-deposition of organic matter and setting of suspended particles) and 127.5 cm of fine clay (representing the settling of suspended particles). Gyttja 2 -5 cm above the silt layer was C14-dated at 5113 ± 31 BP or 3882 ± 82 cal. yrs BC (indicating the re-deposition of older reworked material). This means that the upper 149.5 cm, in fact, represents one huge tsunamite divided into 3 parts of bedload transport and settling of suspended particles in a mega-graded-bedding unit (black arrow in Figure 19).
Site X refers to an overgrown dead-ice basin, where 3 cores were taken (     In the north, the overflow threshold is at about +27 -28 m (Figure 22), and on the northern side the geological map [17] shows a wedge like layer of fine sand (violet field in Figure 22), which in view of the tsunami episode, might be the deposition of material from a violent overflow episode from the south (blue arrow in Figure 22). Core Y was therefore taken in the wedge (Figure 23), and, indeed, it shows 1.25 m of alternating fine sand (from tsunami bedload) and fine clay (from the tsunami suspended matter). Figure 23. Core Y from the deltaic fine sand wedge associated with the +15 m shore level just north of the tsunami overflow area to the north ( Figure 22). The deltaic deposits are 125 cm thick and consist of fine clay (from the particles in suspension) and fine sand (from the bedload at the overflowing of the outlet threshold). The main +15 m outlet of the valley in the east is a deep erosional gauge ( Figure 24). Prominent depositional terraces were leveled at +14.5 m, +11.2 and +3.5 m. The +14.5 m terrace fits with a sea level at +15 m at the time of the tsunami, which must have had a violent outflow here.
The western inlet area (at the village of Krägga), the strong forces of the tsunami wave are seen in the eastward fingering of the +25 m contour (Figure 16). One valley points directly to the Apalle area, where traces of the tsunami wave are seen up to +31.5 m.

The Apalle Archaeological Site
Apalle is an archaeological site from the Bronze Age, about 350 × 70 m (NW-SE), with an unusually thick, complex cultural layer [16] in the cultivated fields surrounding a crest with bedrock outcrops that reach above +31 m. This site is interesting to reconsider, both geologically and archaeologically.
Excavations in the later 1980s revealed a cultural layer over, under and between 79 house foundations, with post holes, hearths, clay flooring, two cultivated plots and eleven artificial wells. The settlement existed between the late Neolithic and the early Iron Age, with two main phases: the early Bronze Age and the latest Bronze Age-earliest Iron Age [16]. We note that the tsunami even occurred right in between those two main phases.
The Apalle monography [16] did not use stratigraphic units (i.e. chrono-stratigraphy) but "types of layers" (i.e. layers without stratigraphic order) as illustrated in Figure 25. This leads to dialectics between the original description [16] and the character of some layers in Figure 25. We have not been able to resolve this even among ourselves, so we leave it for a future revision. In Figure  25, the geological suggestions (NAM) are expressed, however, for the sake of possible tsunami identifications.
The stratigraphy begins with a normal regression sequence from Baltic clay (layer 9) to littoral shore sand (8) and a land-surface with remains of vegetation (6) dated at 3140 ± 110 BP or 1390 ± 130 cal. yrs BC. This date surely pre-dates the tsunami event. Figure 25. Stratigraphy at Apalle [16] with geology-based reinterpretation in terms of the re-deposition at two successive tsunami waves (marked yellow and tan). The Layer 6 land surface (1390 ± 130 cal. yrs BC) pre-dates the tsunami. Layer 11 is a "till" containing "gravel and cultural remains". This indicated re-deposition under strong currents: we propose the tsunami. Layer 7 at the top of Section B is our problem. Clay mixed with pebbles and stones sounds like a typical tsunami deposit (the later phase with fine clay deposition). Chronologically it is held to be younger, however [16]. The material covering the old land-surface is very complex (Figure 25). The tsunami event seems to open new and logical interpretations, however. But there remains a lot of work to be done with respect to chronology [16] and strict correlation among radiocarbon dates and layer types. This will become an interesting challenge for future work.
At this site, the tsunami run-up must have reached +31.5 m (Figure 25), which is higher than at Annelund (+28.5 m) and Boglösa (+29.5 m). The run-up level of a tsunami wave is not something fixed in elevation, rather it varies greatly with respect to exposure, shore configuration and possible interference of separate waves (in this case waves from the west and from the south). Therefore, it is concluded (NAM) that the tsunami wave must have reached up to +31.5 m here at Apalle.
Several houses and wells were excavated [16] [33] (Figure 26). The wells were dug in clay. Tree branches were pushed down around the edges to keep the clay in position. Three of those were dated: 3110 ± 30 BP or 1379 ± 42 cal. yrs BC in well A1992, 3090 ± 65 BP or 1354 ± 73 cal. yrs BC in well A1640, and 3045 ± 55 BP or 1303 ± 77 cal. yrs BC in well A1514. All the 3 dates should pre-date the tsunami. The 3 dates overlap in the zone of 1338 to 1380 cal. yrs BC (1351 ± 21 BC), which provides an important date of the period preceding the tsunami event. Figure 26. Location of houses (violet) and wells (blue) at Apalle (redrawn from [32]).  The wells are filled by "clay and anthropogenic material" according to [33]. Figure 27 gives the stratigraphy of well A1640. Layer 14 (the first layer to fill the well) is "clay rich in wood remains and with occasional stones". Layer 13 is "clay with brittle-burned stones, charcoal and wood remains". Layer 12B is "a clay layer with occasional stones and cultural remains". Layer 12 is "clay with cultural remains, stones and bones" (including a human bone). This looks far from a natural filling. The alternations between more clayey layer (14 and 12B) and layers of more cultural remains and stones are proposed to represent the overflowing of the well by two successive tsunami waves both including a running-up phase (14 and 12B) and a swashing-back phase (13 and 12) drawing down material from the surroundings (including one human bone together with bones of cow and pig). Layer 11 post-dates the tsunami.
From the elevation (Figure 27), it is clear that the tsunami must have reached, at least, up to +28.5 m (in order to flood well A1640), which is in full agreement the records from Annelund, Boglösa and the coring in the Håbo area.
It seems significant that pollen analysis of the well A1640 sequence lack signs of progressive development, unlike the continual development recorded in a progressively infilled well from the Iron Age [34] [35].
Ullén [16] recognized 54 houses (43 separate houses according to [36]), which she arranged in 5 phases of housing with related activities (wells, burnt mounds, farming, etc.). With a contemporaneous tsunami event, things become complicated, however. Below follows a geological based (NAM) interpretation of the material presented [16].
Phase 1 includes 2 houses. There are two dates from house K10: soot dated 3055 ± 60 BP (1310 ± 79 cal. yrs BC) and a food crust on a jar fragment dated 2835 ± 50 (1010 ± 72 cal. yrs BC). This means a difference of about 300 years, which seems questionable. The jar fragment with a food crust was obtained below "the clay floor", considered to be human made. If the clay comes from the tsunami, however, things seem more logic. There was an early house with an age of about 3055 ± 60 BP or 1310 ± 79 cal. yrs BC. Then came the tsunami and flooded the area and deposited the jar fragment in a layer of clay, which would give an age of the tsunami of about 2835 ± 50 or 1010 ± 72 cal. yrs BC (Appendix).
Houses K9 and K27 overlap each other. Between the house remains there is a clay floor of 4 m (double the size of the houses), which may represent a tsunami deposit (Figure 31 in [16]). Unfortunately, the houses are undated.
Phase 2 includes 10 houses and 5 dates ranging from 3110 to 2970 BP (mean: 3030 ± 40 BP) or from 1360 to 1195 BC (mean: 1268 ± 80 cal. yrs BP). This house phase is likely to pre-date the tsunami.
In conclusion, this means that the Apalle Sections A and B establish a tsunami run-up of about 16 meters, based on a contemporary shoreline level +15 m. The evidence seen in the wells, establishes a necessary run-up of a minimum of 13.5 meters. This coincides well with the data from the Annelund site, the Boglösa area and the coring in the Håbo area. A comparison of available C14-dates is given in (Appendix).

The Western Part of the Lake Mälaren Area
Theoretically, the tsunami might have reached all the way in to Örebro ( Figure  4). With the new tool of a mega-tsunami event at about 3000 BP, it seems urgent also to revisit geological and archaeological sites in the entire Lake Mälaren region ( Figure 4). Lake Hjälmaren has a present level of +21.8 m (drained by about 1.3 m in 1878-1887), so even this region might have been affected by the tsunami. We have investigated the map description and relevant literature of the area without finding any firm records of the Lake Mälaren 3000 BP mega-tsunami here described.
The tsunami wave-height is likely to decrease with distance inland. In the Enköping-Bålsta area a run up of 13.5 to 16.5 m (the wave reaching 28.5 to 31.5 m) is well established. Further westwards, the situation is still unexplored, however.
Southeast of Västerås [19] and between Strängnäs and Eskilstuna [20] there are records of fine clay resting upon gyttja clay including sand from bedload deposition as mentioned above (and marked with red dots in Figure 4).
Robertson [37] made a study of the shore-level displacement based on the isolation of lake basins at different elevations in the Eskilstuna area. She met problems in the dating of isolation basins in the +24 to +35.7 m level. The isolation of Fågelmossen at +35.7 m was dated (C14-ages) at 3850 ± 130 BP, Fräksmyren at +31.1 m at 4135 ± 70 BP, Janslunda at +29.3 m at 4225 ± 110 and right below 4045 ± 70 BP (i.e. in reversed order), Stormossen at +28 m at about 4400 BP, and Borsökna sjön at +24 m at 4290 ± 195. These are, of course, impossible dates as it would imply that basins at 12 m difference in elevation would all have been isolated at about 4100 BP and in a nearly opposed sequential order. Besides, all dates were old "whole-sediment" dates before the AMS-dating technique.
The Janslunda bog (at +29.3 m) was cored in May, 2020 (our Site Z). The original paper [37] included some points that seemed open for reconsideration, viz.: 1) a 3 cm "sandy gravel" at the boundary between gyttja clay and Baltic clay (which sounded like a tsunamite), 2) the occurrence of a 51% peak of planktonic brackish-marine diatoms in the gyttja clay above (which usually is a criterion to identify a tsunami bed [7]), and 3) the inconsistencies of the C14-dates which opened for a revision.
The new core taken ( Figure 28) records a deep erosion (irregular and sand filled) going down into grey clay (maybe of the Ancylus lake stage), implying a hiatus of about 6000 years.  The irregular erosions surface ( Figure 29) is covered by 1.0 -2.5 cm of sand including grits of fine gravel and coarse gravel that lack signs of rounding excluding a littoral deposition; rather an instantaneous and violent tsunami deposition. This bed has all the characteristics of a tsunami bed. The covering 3.3 cm consists of gyttja-clay (fine clay) with silt and angular grits of sand and gravel (even coarse gravel). This unit can only have been formed due to a mixture of violent bedload deposition and settling of suspended material. The covering 0.7 -1.0 cm of fine-clay-gyttja seems to represent the final part of the unit with only the settling of suspended gyttja and clay matter.
The diatom flora (as measured at 6 levels by [37]) can now be understood in other terms: 1) 60% Fragilaria sp. + 28% brackish to brackish-marine species represents a mixed flora dominated by reworked material, 2) 51% of Melosira westii, a planktonic brackish-marine species, fits the characteristics of tsunami deposits [7], 3) a mixture of 50% of both brackish and freshwater species seem to mark the later phase of the tsunami, 4) 85% Fragilaria sp. implies final setting of suspended material close to the isolation contact, and (5 -6) represent full lacustrine conditions. The recording of the mega-tsunami tsunami event seems now firmly fixed in Figure 29. It represents a tsunami wave height of at least 14.5 m, which implies that Lake Hjälmaren (+24 m) must also have been invaded by the tsunami. A record of this seems to occur in Ullavisjön (Magnusson, 1970), but we were not able to document this by coring at the site. Figure 30. The Långaedet area and the Lower and Middle Lake Marviken. Three sediment cores were taken in a small bog at +29.5 m (red dots) and the core obtained (right) with a clear tsunami layer. In the valley side to the west, there is a big earth slide and an erosional gully ending in an alluvial cone (yellow) with its inner apex at +34 m, its break-in-slope at +26 m and its foot at +17.7 m, just above the shore at +15 m at the time of the tsunami.

The Åker-Långaedet Area
The area is traversed by an old fracture zone, now occupied by a system of lakes: the Lake Marviken lake system. The word "marviken" goes back to "fractured" [25], and, indeed, the valley is characterized by fractured bedrock, fault movement, earth slides and methane venting tectonics [2] [22] [23]. Sediments cores in the lake record a tsunami bed over a distance of at least 6 km. It is dated at 3055 ± 70 BP or 1303 ± 91 cal. yrs BC. Previously, it was interpreted as a local lake tsunami generated by a slide [23]. Now, it seems obvious that it was contemporaneous with the mega-tsunami event recorded in the whole of the Lake Mälaren valley at about 3000 BP. Therefore, the site was revisited in 2020. In a small bog at +29.5 m ( Figure  30 There are also morphological criteria in the valley side to consider ( Figure  31). In the valley slope, there is an erosional gully with back-erosion up the slope and with a sub-areal alluvial cone in front ending at +17.7 m just above the shore of the +15 m sea level at 3000 BP ( Figure 30). The apex of the alluvial cone is at +34 m and the outer break-in-slope is at +26 m. The gully is assumed to have been formed by back erosion initiated by the tsunami run-up at about +34 m due to narrowing of the valley and maybe also the interference of waves coming both from the southeast and the northwest. There is also a land slide ending at +29 m.  The new data indicate that the tsunami was not local but a part of the Baltic-wide event occurring at about 3000 BP, and that the run-up was at least 13 m (at the core site) and maybe as much as 19 m at the gulley ( Figure 30).

Upplands-Väsby (Tolan)
As a part of an archaeological investigation, we also undertook a geo-tectonic analysis of the Tolan site and Sandaskogen area ( [24] and supplementary material to [11]). Obviously, the area had suffered severe ground-shaking. In lake Sandakärret at +35.5 m, a tsunami bed was recorded (in 3 cores) consisting of unsorted and angular grains. It was dated at 3073 ± 33 C14-yrs BP or 1315 ± 85 cal. yrs BP (in full agreement with other dates of the mega-tsunami). In order to reach into the lake basin, the tsunami-wave need to have had a run-up of 20.5 m (which is the highest value recorded).
Some 5 km to the south, Henkel [38] has described an earth slide over a clay field and deformed bedrock, which seems to correlate to the Tolan deformations.

Saltsjöbaden
A tsunami bed in Lake Dammsjön at +24 m indicates a tsunami height of at least 9 m (rather 10 to reach over the outlet threshold). The site is described in the supplementary file to [11]. Cores taken at +28 -29 m did not have any visible tsunami layers. The bedrock steepness outside the threshold suggests that the run-up above tsunami wave-height was small to absent.

Södertörn (Fituna-Kvarnberget)
Kvarnberget is a type-site for methane venting tectonics (MVT) [11] [22] [26]. It is linked to the Fituna delta at +15 m [11] [22] deposited at about 3000 BP as a consequence of the tsunami event, the wave-height of which was estimated at a minimum of 10 m.

Tystberga
Tystberga refers to some old sites investigated by Florin [39] with respect to sea level changes and archaeology. It is clear that his high-amplitude sea level changes were in need of revision. The tsunami event at Fituna on the Södertörn peninsula fits both in time and elevation [11] [22]. At Rogsta Hage in Tystberga, there is a section which seems to record a tsunami layer extending up to +31 m (supplementary material to [11]). Therefore, we seem to be dealing with a 16 m tsunami run-up.

Varve Dating and Tsunami Records in Ångermanland
At Umeå, the tsunami height was established at 11 m and at Hudiksvall at 12 m ( Figure 1). As the River Ångermanälven area lies in between these two points, one might expect to also record the tsunami at River Ångermanälven. This area is interesting because here we have a chronology based on annual varve counting [40]. The original varve counting by Lidén [40] was not connected to the present, however. This connection was later successfully established by Cato [41] [42] [43]. The dataset consists of 16 varve dated deltas ranging in elevation from +134 m to +10 m and in time from 6072 BC to 555 AD. After a careful investigation of available data, it now seems possible to identify a significant tsunami event and provide an absolute date in terms of the Swedish varve chronology [42]. The effects of a mega-tsunami are revealed in the following three facts (as illustrated in Figure 32): 1) Site 16 is a terrace with its last varve dated 1165 BC, though it would have been at sea level at 600 BC (i.e. ca 665 varves are missing). The terrace lies 7 m below the sea level at that time. It is covered by 7 sandy varves, which seems to represent a stratified tsunamite of 118 cm thickness (composed of two waves).
2) Site 9 is a terrace at +46 m and has its last varve dated at 3086 BC, though it would have been at sea level at 2100 BC (i.e. ca 986 varves missing). The terrace lies 14.5 m below the sea level at the time of the last varve. It is covered by a potential tsunamite.
3) The varve dates of Sites 10 to 13 are off-set by 120 -200 years from the straight line of shore displacement (blue line in Figure 32). Figure 32. Ångermanland elevation/time plot based on data from [42]: blue line gives shore-level displacement, red plus-signs refer to varve dated deltas, blue plus-signs refer to terraces, red dot in a U refer to varve dated lake records [44] [45], purple line gives sea level at tsunami event (+32.5 m) and terrestrial run-up (Site 9), submarine bedload (Site 16), and run-up erosion at Sites 10 -13. The tsunami seems to have had a total run-up of about 13.5 -18 m and submarine erosion/deposition effects down to −7 m. I think we can now provide quite a firm interpretation. In the varve year 1171 BC, the area was struck by a mega-tsunami. It deposited a 118 cm thick stratified tsunamite on top of the Site 16 terrace implying a submarine bedload deposition at a depth of −7 m. The same event had a run-up of at least 13.5 m, eroding the top surface of the Site 9 terrace and depositing a stratified tsunamite of 688 cm (previously counted as 8 varves [42]) or (as shown in Lidén's original stratigraphic drawing, Figure 9 in [42]), a unit of 4.7 m consisting of 4 beds of gravel, which might represent a sequence of 4 tsunami waves (in analogy with the 10,430 BP tsunami in southern Sweden consisting of 6 cycles of graded bedding in Figure 4   The +33 m level has a very special shore-morphology with ridges, shore-spurs and erosion channels ( Figure 34). We cored the Lintjärnen bog and retrieved the tsunami bed covered by 13 varves and a peat cover. Gyttja in varve 13 was C14-dated at 1306 ± 87 cal. yrs BC. The date represents redeposited organic matter from the tsunami run-up. International Journal of Geosciences

Discussion: Geology
There are two points of general implication that need to be mentioned before our general discussion; that is the C14-datings and the reference level of elevation.
All our C14 dates refer to AMS-dates with absolute age calibration according to Reimer et al. [46]. Older dates referred to in connection to archaeological material may be calibrated with older standards. The differences are small, however.
In a few cases [37] the dates refer to the old pre-AMS technique. The ages given in the text are either C14-ages in years BP (i.e. before 1950) or in calibrated yeas BC, "cal. yrs BC", (in a few cases in cal. yrs AD). has increased, but for the last 3000 years it seems to have remained stable [47].
We may therefore, conclude that sea level 3121 years before HY200 is likely to have been at +15 m to +16 m. In this paper we used +15 m.

General Discussion
At about 3000 BP or 1200 cal. yrs BP something both unique and revolutionary happened: the Kaali impact occurred in Estonia (e.g. [49]) and initiated violent ground shaking, methane venting tectonics and a mega-tsunami affecting the whole of the Baltic [2] [10] [11] [22]. In this paper, we have explored the tsunami effects in general and specifically in the Enköping area. We have combined geological and archaeological data, and, indeed, where ever we work, we find evidence of a major tsunami event at about 3000 BP when sea level was at about +15 m in the Stockholm area.
We found the following tsunami run-up values ( Figure 1 The documentation is extensive. The run-up values differ from about 10 m to about 15 ± 1.5 m, maybe even 19 -20.5 m. This is not at all surprising because the run-up level is not fixed but varies with several parameters like exposure, off-shore slope, land topography, interference of waves [12]. In our case, we also have records of the off-shore effects on the seabed ( Figure  3), with erosion going down 13.5 -14.5 m below the sea level 3000 years ago: at +15 m ( Figure 14). With such a deep off-shore effect, it is quite natural to have the run-up wave over land, at least, as high as this and even somewhat higher. Figure 35. Core records of the tsunami wave submarine seabed surface trimming below sea level at the tsunami event (+15 m at 3000 BP) and the on-land run-up and tsunamite deposition in the Enköping-Bålsta area (black dots refer to elevation of cores retrieved). The tsunami amplitude is +15 m ± 1.5 m (pink arrows). In the land area, the run-up may occasionally go up to 19 -20.5 m, however.
Because the Enköping-Uppsala-Stockholm area is a lowland characterized by clay deposition, the tsunami wave sent up clouds of fine clay particles in suspension. This led to the deposition of fine clay beds in the off-shore area. The tsunami run-up over land, usually led to severe erosion; sometimes the erosion went down all the way into the glacially varved clay, leaving a hiatus of about 7000 -8000 years below the tsunamite. Figure 35 graph is quite unique in recording both the on-shore run-up values, and the submarine erosional/depositional effects. The tsunami deformations are symmetrically ordered in ±13.5 m around the +15 m sea level at the time of the mega-tsunami event. This implies quite symmetrical waves ( Figure 36) before the tsunami hits the coast and runs-up over the land surface causing erosion, deformation and re-deposition (as observed by coring, e.g. Figure 13, and reinterpretation of archaeological data; e.g. Figure 25, Figure 27). International Journal of Geosciences

Additional Data
It seems where ever we go, we find new examples of the Mid Bronze Age tsunami event. At Björklinge (20 km north of Uppsala; [50]), hearths with artefacts and bones were found at +36 m and +34 m (site Tibble II) covered by 10 -15 cm of gravelly sand with occasional stones; earlier considered "remarkable", now perfectly well understood in terms of the 3000 BP tsunami event. At a close-by site (Tibble III), there was a hearth at +30.6 m including a flint tool from the Bronze Age. The hearth was covered by 15 cm of gravelly sand; obviously from our 3000 BP tsunami event.
The thin layers of gravelly sand at Björklinge are typical tsunami deposits, like the "silt layer" at Annelund (above). Similar layers are recorded at former tsunami events all over to globe. They were also observed at the huge Indonesian tsunami in December, 2004 [51] [52] [53].
Even a thin layer of sand and silt may alter the permeability of a tsunami-flooded bed. This may change the drainage conditions and, for example, initiate peat formation. Deposition of gyttja may act as a fertilizer.
A tsunami may also lead to severe costal destruction [14] [27]. Ridges and mounds may become remodelled to flat surfaces, and heavy objects like boulders may become relocated [51].

Gillberga Gryt in NE Uppland Revisited
Gillberga Gryt refers to a violent bedrock deformation, where a previous "roche moutonnée" has been deformed into a heap of huge angular blocks, often with one surface representing the old glacially polished surface, indicating that the deformation took place in postglacial time [2] [17] [54]. It was interpreted in terms of a high-magnitude earthquake occurring about 10,000 BP [2] [54]. The bedrock deformation is bounded to the north by an extensive peat bog at an elevation of +20 m. The level implied that the tsunami event was likely to be recorded in the bog stratigraphy.
We undertook an extensive coring program in the bog ( Figure 37). The profile is 130 m long and contains 7 coring sites. In sites 1 and 2 additional cores were taken at 1 m interval to record the topography of the rock surface that stopped deeper coring. It reveals a heavily irregular surface indicating the presence of blocks 14 m (1) and 30 m (2) from the present foot of Gillberga Gryt (angular blocks). From sites 5 to 8, a typical tsunami layer extends over at least 86 m (just as the layers we recorded at so many places in the Lake Mälaren region). The tsunami bed consists of fine sand and silt and reworked gyttja with occasional pebbles. It exhibits graded bedding in cores 5 and 6.
Gillberga Gryt lies just to the south of the Forsmark area from where a tsunami was recorded at several sits ( [7] and Supporting Material to [11]).

Archaeological Perspectives
This paper presents geological evidence suggesting the occurrence of a mega-tsunami in eastern Sweden around 3000 BP, more exactly 1171 BC. A primary question is whether this catastrophe is also traceable in the archaeological material of the Lake Mälaren region. For archaeologists, there are several key and methodological problems, not least: what is there to look for, where is such evidence found and what is important to document and analyze?
The above discussions have looked at three examples of anomalies documented in earlier archaeological excavations that have been described as difficult to understand or unusually explained archaeologically: at Annelund, Apalle and Björklinge. These have all followed contemporary archaeological practices and praxes. The first two took place in the mid-1980s during a dynamic period of archaeological development and administrative change and an intensive period of field work involving excavations and ancient monument surveys. This led to new archaeological understandings of settlement variations and the cultural materials of Bronze Age and their results became the focus of many studies (J. Anund in [16]). Previously, very few excavations had been carried out in the Trögden and Håbo areas, and discovery of evidence suggesting a tsunami catastrophe was still a decade farther into the future.
Here, attempts have been made to include perspectives of two complex and independent disciplines, each of which works with complicated and independent processes: natural geological and human activity. Geological raw materials are undeniably mutual to both disciplines-soils and rocks-where the distinction between the two is not wholly distinct. When do geological materials and evidence of natural processes become part of the scope of archaeological processes as evidence reflecting human activity: artifacts, ancient monuments, humanly created cultural layers and stratigraphy? In many instances, these are inextricably intertwined. From an archaeological perspective, this is not an obvious question and there is no standardized answer (Anderson Ambrosiani, forthcoming).
Since the early 1900s, principles of the natural processes of isostasy and eustasy have been used in archaeological fieldwork. Isometric information was a key factor combined with excavation results in early Swedish archaeological settlement studies [55]. Archaeologists have continued to use such landscape information as a tool in fieldwork: in surveys and excavations, in interpreting ancient monument distributions and in developing and discussing models of settlement development.
However, it is as obvious to archaeologists that archaeological phenomena cannot be fully explained by geology as the opposite is to geologists. At the same time, experiences from this exploratory project and paper suggest that knowledge of such an event can be potentially interesting to archaeology with new possibilities for archaeological research. Adding such a perspective to theoretical and archaeological practices and praxes is a challenge to archaeological fieldwork and explanation alike. For archaeology, exploration and interpretation of sites and materials must be carried out source critically with a foothold in archaeology and respectfully, with the goal of maintaining the authenticity and integrity of the archaeological materials.
Future ancient monument studies and archaeological excavations, particularly in the Lake Mälaren region, that deem these possibilities interesting and worthwhile to explore further, may formulate research questions and create strategies for identification and documentation already in the planning stages. This event may create a postand ante quem for a long prehistoric period that until now has lacked a fixed chronological point in time. Cultural materials from the Lake Mälaren region, roughly from Period III, the time of the tsunami and the transition between the early and late Bronze Age, have long been notably scarce.

Conclusions
The tsunami records from Umeå, Hudiksvall, Forsmark, Upplands-Väsby, Saltsjöbaden, Åker, Södertörn, and Tystberga, all occurring at about 3000 BP, were in this paper put to a detailed test and analysis with respect to geology and archaeology in the Enköping-Bålsta area. The new data provide evidence that, indeed, the area was struck by a violent tsunami at about 3000 BP or about 1200 cal. yrs BC.
The on-shore tsunami run-up was determined at +13.5 m at Annelund, +14.5 m at Boglösa, +13.5 m in the Håbo area and +16.5 m at the Apalle archaeological site.
In the off-shore area traces of erosion and deposition of fine particles in tsunami suspension were traced 13.5 -14.5 m below sea level ( Figure 35).
The tracing of a tsunami wave, both in the on-land run-up phase and the off-shore erosional phase, is unique and acts as a confirmation of the tsunami episode.
The exact timing of the tsunami differs somewhat. The dates referring to the tsunami event itself fall within 2900 -3070 BP or 1100 -1300 cal. yrs BC (i.e. 1200 ± 100 BC). There are a number of dates pre-dating the tsunami; their age falls within 3400 -3000 BP or 1650 -1300 cal. yrs BC. There is also a group of dates post-dating the tsunami; their ages are <2900 BP or <1100 cal. yrs BC. The dates are compiled in Appendix.
At a late stage (April-June, 2020), the Ångermanälven analysis (Section 5) was performed. Evidence of the occurrence of a mega-tsunami was established, and it was found that the tsunami must have occurred in varve-year 1171 BC or 3121 BP, an age which fits perfectly well with previous C14-dates of the tsunami event in Sweden (Appendix, [10] [11] and the age of the Kaali impact [10] [11] [49]. The occurrence of a tsunami wave (or sequence of two or more waves) with a run-up amplitude of 10 m to about 15 m ± 1.5 m, maybe even 19 -20.5 m, in the middle of the Bronze Age will have wide, maybe even revolutionary, effects on our understanding of the archaeological evolution during the Bronze Age. We have analysed the records from Annelund and Apalle, and been able to propose significant reinterpretations. We feel that this is just the beginning, however.
The Brandskogsskeppet rock-carving is also considered in relation to the tsunami. The level of the rock surface at +28.5 m coincides with the tsunami run-up in the Boglösa area.
The western part of the Lake Mälaren region and the Lake Hjälmaren region are in need of extended analyses, like the whole of the South Baltic area. It is our sincere hope that our paper will inspire such studies.