Introduction

Disturbances are short-term, temporal variations of environmental conditions and force impulses in changes of the nature of forest ecosystems (Seidl et al. 2017; Dale et al. 2001). Factors such as classified in natural (e.g. orbital forcing, volcanic eruptions, drought, fire, frost, windstorms, landslides, animal attacks; Esper et al. 2012; Esper et al. 2015; Hartl-Meier et. al. 2014a; b; Attiwill 1994) and anthropogenic (e.g. settlement, harvesting, land use, pollution; Neumann 1991; Popradit et al. 2015; Hartl et al. 2019) can affect plant physiology (Attiwill 1994; Büntgen et al. 2009) which can have negative effects on tree growth and vitality or even lead to tree mortality due to their intensity (Seidl et al. 2014). However, the impacts of disturbance factors on forest dynamics do not only depend on the intensities of the factors but also on the characteristics of the tree’s site and, particularly, on the tree species (Kollas et al. 2013).

In our study, we examined trees of the species Betula pubescens Ehrh. (downy birch) at the latitudinal tree-line at the Kåfjord next to Alta in northern Norway. Downy birches are the predominant tree species in this region (Tømmervik et al. 2009). B. pubescens is a mid-sized, broadleaved tree species that occurs naturally all over Europe. Its ability of growing rapidly on well-drained soils as well as on poorly developed sites predestinate downy birch trees to be a pioneer tree species, even in boreal regions (Seppä et al. 2002; Beck et al. 2016). Thus, downy birches are well adapted to the environmental conditions at the latitudinal tree-line, but, nevertheless, there are also growth-limiting disturbances, especially for rapid changes in climate conditions (e.g. Opala et al. 2015; Levanic and Eggertsson 2008), moth outbreaks (e.g. Karlsen et al. 2013; Jepsen et al. 2011) and, particularly around the Kåfjord, short-term artificial pollution (Hartl et al. 2019).

At extreme sites like altitudinal or latitudinal tree-lines, tree species respond generally sensitively on changes of climate conditions (Körner 2012; Kullman and Öberg 2009). Whereas precipitation shortage can lead to a negative growth response at relatively dry sites (Beck et al. 2016; Lyu et al. 2019), the northern tree-line is typically controlled by temperature (Kullman 2003; Körner 2012). In a strict sense, long-term changes in temperature belong to climate variability and are no disturbances. However, especially exceptional, short-term changes in temperatures can lead to a disturbance of tree’s growth period. In fact, unusual cold temperatures at the beginning of the bud-burst period of B. pubescens in May can interrupt the tree’s unflowering and shift the onset of the growing season to a later time (Karlsson et al. 2003; Delpierre et al. 2018; Eckstein et al. 1991). Low temperatures in late summer can terminate the tree growth earlier as expected (Klaveness and Wielgolaski 1996). Both phenomena can shorten the time of growth severely and disturb the current tree’s growth and even the subsequent years that can lead, in serious cases, to a tree’s dieback (Karlsson et al. 2004; Eckstein et al. 1991). Mean July temperature with additional effects of June and/or August as expressed by JJA mean temperatures is a good indicator for the general growth conditions in northern Fennoscandia (Karlsson et al. 2004; Kirchhefer 2001; Esper et al. 2014).

The most common biological disturbance factors of birch trees in Fennoscandia that have severe negative impact on tree growth are insect outbreaks of autumnal moths (Epirrita autumnata Bkh.) and winter moths (Operophtera brumata L.) (Tenow 1972; Karlsen et al. 2013; Jepsen et al. 2008, 2011). The life cycle of the moths begins with the egg hatch between Mid-May and late June. Their larvae feed on expanding birch leaves (Fig. 1e) and can cause massive defoliation and diebacks of a whole site during vegetation period until the moth population collapses (Babst et al. 2010b; Beck et al. 2016). Moth outbreaks can stay over 2 years at the same site, but they occur most of the time synchronously in a 10-year cycle barring that the winter moth pests delay sometimes by 1–2 years. Thus, there is a negative, periodically occurring effect on growth of B. pubescens that is not only limited on the first year of disturbance event (Bylund 1999; Neuvonen et al. 1999; Hagen et al. 2007).

Fig. 1
figure 1

© Tirpitz Museum), and e a foliation feeding moth (here E. autumnata, © Photo: Moritz Klinghardt, Norwegian Institute for Nature Research)

a Map of Fennoscandia showing the location of the study site (orange dot) in northern Norway. b The sample site is a pure birch site at the southeast coast of the Kåfjord. c Climate diagram from Alta Lufthavn considering data from 1964 to 2017. d Pollution through smoke screens hiding the battleship Tirpitz in 1944 (

Growth response of trees on long-term artificial air pollution was a topic for many studies of the last decades (e.g. Elling et al. 2009; Malik et al. 2012; Däßler and Lux 1984; Haselhoff and Lindau 1903). However, the kind of pollution at the Kåfjord is very unusual. The artificial pollutions, caused by the usage of smoke screens of the German battleship Tirpitz in 1944 (Fig. 1d), were locally restricted, short-term disturbances in the area around the Kåfjord (Hartl et al. 2019). The Tirpitz (the sister ship of the Bismarck) anchored at the fjord of Alta as a fleet-in-being to disturb allied convoys through the Arctic Ocean (Frère-Cook 1977). Recognized as a military threat, the allies tried to sink the vessel during six raids between April and September 1944. The German battleship protected itself against the aerial attacks with local ignited smoke consisting of sulfur trioxide, zinc/hexachloroethane, sulfuric and hydrochloric acid (Oberkommando der Kriegsmarine 1942). These artificial fogs covered the whole fjord of Alta and were very noxious for the German soldiers (H. Dv. 211/2 1939; H. Dv. 211/1 1940; H. Dv. 211/3 1941, Hampe 1963), furthermore for trees close to the anchor place of the Tirpitz (Hartl et al. 2019).

To analyze the growth responses of the birch trees to disturbances of extreme climate conditions, moth outbreaks and artificial pollution back in time, an analysis of tree rings is an appropriate method due to providing precise detection of annual resolved growth reactions with the opportunity of forest dynamics reconstructions in decadal to even millennial scales (Trotsiuk et al. 2018; Fritts 1976; Fritts and Swetnam 1989). An analysis of birch tree rings is challenging due to the non-straight stem growths combined with a more frequent appearance of missing rings than in many other tree species (Kirchhefer 1996). This growth sensitivity can be traced back to the way B. pubescens cope with (mechanical) disturbances. Mechanical disturbances influence the anatomical structure negatively, so the xylem-specific hydraulic conductivity decreases; the according tree ring will be narrow. Thus, the inter-annual tree-ring width (TRW) variability is depending on the intensity of the disturbance on single tree level (Tumajer and Treml 2019). To consider this behavior, we will investigate the growth response to varying disturbance factors by analyzing growth/climate response using climate data and TRW but also by analyzing growth variability around the years of the different disturbances, i.e. the basal area increment (BAI) of single trees to detect growth decline and trees’ recovery ability (Zang et al. 2014). Finally, we rank the disturbance factors by their harmfulness for B. pubescens at our site.

We assume that (1) the TRW at the latitudinal tree-line of northern Norway will show a correlation to exceptional negative climate conditions that will be visible as a decline of single tree growth and a low growth variability throughout the site in the rings of the current and following years; (2) the most severe moth outbreaks recorded for Alta affected growth of all trees negatively, so that the growth decline for all the tree will be similar, the growth variability will be low. In addition, we expect that the growth decline after moth outbreaks was stronger and the time for recovery took longer than after climate extremes. For the disturbance due to (3) artificial pollution, we hypothesize that there could be a pronounced single tree response in 1944 due to probably different local exposure to the artificial smoke. The results will show a pattern of growth variability in the year of disturbance and in the aftermaths that is not comparable to the two other disturbances due to the local emission of pollution.

Materials and methods

Study site and tree-ring data

In August 2017, 36 downy birch trees were sampled at the Kåfjord close to Alta in northern Norway (69.939° N 23.063° E, Fig. 1). The birches root in a poorly developed soil on a sandy and gravelly glacier moraine at ~ 10 m asl (Geological Survey of Norway 2019). Climatic conditions are characterized by snowy and cold winters, and cool summers (Kottek et al. 2006) including 412 mm precipitation and annual mean temperatures of 1.7 °C.

Each tree was cored twice at breast height (~ 1.30 m) with a 5 mm increment borer. The preparation for TRW measurements included cutting the 64 cores with a microtome (Gärtner and Nievergelt 2010) and staining them with Astra blue to improve the contrast of tree-ring boundaries (Rothe and Hartl-Meier 2014). The samples were measured at a resolution of 10 μm using a LINTAB measurement device and the software TSAPWin (both Rinntech 2011). Missing rings occurred relatively frequently and were set to a minimum thickness of 10 µm to support later data analysis. Cross-dating accuracy was checked visually and statistically using COFECHA (Holmes 1983).

The TRW data were power transformed to reduce heteroscedasticity within the chronology, detrended and standardized using cubic smoothing splines with a 50% frequency cut-off at 30 years to remove age-related growth trends. By this procedure, dimensionless ring width indices (RWI) were calculated in ARSTAN (Cook and Peters 1981; Cook 1983). A mean RWI chronology was produced by averaging the detrended single series using a robust mean (Mosteller and Tukey 1977). As we were also interested in actual growth, we also transformed the raw ring-width series into BAI, the annual ring area increase, based on the ring circumference approximated by a circle (Biondi and Qeadan 2008).

To estimate the internal coherence of the chronology, inter-series correlations (Rbar) were calculated over 30 years windows, lagged by 15 years (Table 1; Wigley et al. 1984). The RWI series were finally truncated at the year 1901 in line with the first year of gridded climate data and after 2000 due to the very high numbers of missing rings hindering an accurate dating (see changing Rbar values in Table 1). Other descriptive statistics for the chronology include the mean series length (MSL), mean sensitivity (year-to-year variability in tree ring width), average growth rate (AGR), and mean BAI, all provided in Table 1.

Table 1 Characteristics of the B. pubescens chronology at the east coast of the Kåfjord

Statistical analysis

The estimation of the influence of the climate elements temperature and precipitation on birch growth at the Kåfjord was calculated with bootstrapped correlations using the treeclim package in R (Zang and Biondi 2015) with the RWI chronology and monthly climate data from the KNMI Climate Explorer (KNMI 2020). We considered previous-year June to current-year December. As the local climate station Alta Lufthavn data covers only the period 1964–2018 but the gridded CRU TS 4.01 data correlates with r = 0.99 for temperature and 0.57 for precipitation, we used the grid point 70.25° N/22.75° E and the period 1901–2000 for further analyses (University of East Anglia Climatic Research Unit, Harris, Jones 2017).

To assess the growth responses of the Norwegian birch trees to different disturbance factors, we considered cold extremes, severe moth outbreaks and pollution as event years (see Fig. 3). Based on the growth/climate response analyses, we derived May/June (MJ) temperatures are most influential for birch growth. Therefore, we selected the five coldest MJ temperatures as cold extremes that are below the 1.96 * standard deviation of the reference period’s (1961–1990) temperature anomalies mean for the twentieth century (see Fig. S1), namely 1902, 1909, 1915, 1951 and 1955. Regional moth outbreaks for northern Norway were recorded in several studies (Fig. 3, green shadings) (e.g. Tenow 1972; Tømmervik et al. 2004; Karlsson et al. 2004; Hagen et al. 2007; Jepsen et al. 2008; Babst et al. 2010a; Karlsen et al. 2013). According to the area around Alta, the most severe moth outbreaks occurred most likely in 1955, 1962, 1975, 1982 and 1989 (Kirchhefer 1996) so that we selected these years (Fig. 1, green lines). For the growth response to pollution, we considered the multiple smoke screen events in 1944 hiding the Tirpitz. To analyze the disturbance factors, we calculated percentage growth changes by dividing the BAI of each single series by the mean BAI of the five preceding years (Hartl et al. 2019). Afterwards, we used the superposed epoch analysis (SEA) (Panofsky and Brier 1958) with the individual BAI series to rank the disturbances. In this study, the time span of 5 years before to 2 years after the particular disturbance were examined to visualize the growth change around the years of disturbance. The variance in BAI were stated as percentual growth changes with respect to the mean of the 5 years preceding disturbance events (− 5 to − 1).

Results

The growth/climate response analysis reveals temperature to be important for tree growth, whereas precipitation correlates insignificantly (Fig. 2). Birch growth is sensitive to early season temperatures including current-year May and June. The seasonal MJ mean temperatures correlate at r1901–2000 = 0.33 (p ≤ 0.05) with the birch chronology.

Fig. 2
figure 2

Climate response analysis of the RWI chronology and the climate parameters a temperature and b precipitation from previous-year June to current-year December. Asterisks (*) mark significant values at p ≤ 0.05

The occurrence of disturbance factors is shown in Fig. 3 together with the RWI chronology and missing ring frequency. Thereby, it is obvious that cold events can appear inconspicuous and do not necessarily cause RWI minima, except for 1902 and 1955. In 1902, there is a growth decline followed by a RWI minimum and a high percentage of absent rings (14%) in the subsequent year 1903. In 1955, the cold temperatures correspond with minimum growth indicated by the RWI chronology.

Fig. 3
figure 3

Characteristics of the birch chronology and temporal distribution of disturbance events. The upper histogram shows the percentage number of missing rings (in grey), the black curve in the middle is the detrended RWI chronology, and the lower grey curve shows the chronology replication, i.e. number of tree-ring series. The blue vertical lines mark the 5 years with the lowest May/June monthly mean temperatures, and the orange line indicates the pollution event with the smoke screen actions of the German Kriegsmarine for hiding the battleship Tirpitz. Green shadings highlight possible moth outbreaks, and the green lines indicate the five most severe outbreaks in the study region

In 1944, the year of the battleship Tirpitz smoke screen incidents, growth rates were relatively low, though not conspicuously compared to other years. However, 8% of the trees already showed an absent ring (Fig. 3, orange line), and these values increased in 1945 to the highest number of missing rings recorded throughout the twentieth century (14 of all 101 absent rings, ~ 22% of the trees). Most of the other growth minima are connected to insect outbreaks, particularly during the more recent decades, including the five most severe outbreaks in 1955, 1962, 1975, 1982 and 1989 (green lines and shadings in Fig. 3) which were all accompanied by missing rings. High numbers of absent rings are also detected in 1975 (~ 16%) and in 1955 (~ 14%).

The single tree growth responses to varying disturbance factors differ substantially (Figs. 4; S2). In 1951, a year of unusual cold MJ temperatures, more than 75% of the trees showed a negative growth change and 50% of the trees showed a growth decline > 25%. The other 3 years of extreme cold MJ temperatures (1909, 1915, 1902) are characterized by a wide spread of growth responses among trees. In 1909, 75% of the trees showed a growth increase compared to the five preceding years. In 1902 as well as in 1915, birch growth is ordinary without any extreme deviations. However, after the year with the coldest MJ temperatures, i.e. in 1903, 86% of the trees show a growth decline followed by a negative growth response of only 38% in the subsequent year.

Fig. 4
figure 4

Percentage growth change of single trees to different disturbance events. The extreme events for temperature (blue) and moth outbreaks (green) are sorted with increasing intensity from upper to lower panels. The growth change for the event year in 1955, a combination of low May/June temperatures and a moth outbreak (green/blue) and the 1944 smoke screen event are shown in the bottom panel. The grey background represents the respective five-year reference period preceding the single events

In 1944, the year of the Tirpitz concealing tasks, the trees reacted differently. Whereas the growth response in the event year is similar to the cold year 1951, with three-quarter of the trees showing a growth decline and 50% of the trees having a growth reduction of more than 25%, in the subsequent year, three-quarter of the trees react again with a growth decline and the median growth reduction is even stronger (~ 45%). In 1946, > 50% of the trees started their recovery while the growth variability was still high (interquartile range of > 89%) and the other trees showed a negative growth response.

The growth response to severe moth outbreaks is less variable, compared to the responses to cold extremes and artificial fog events, as almost every tree reacted with a decline. In 1962, more than 50% of the trees reduced their growth by > 17% in comparison to the 5 years preceding the event. In the years with heavier outbreaks the growth response is even worse, e.g. median growth declines > 54% in 1975 and 1982, and > 57% in 1989. Importantly, the recovery of almost every tree started immediately after the event year while the growth variability between the individual trees increased.

The persistence of growth reductions increases with the intensity of the recorded moth outbreaks, as seen in Fig. 4. After 1962 and 1975, the trees needed only 1 year for (almost) recovery, after 1982 > 50% of the trees and 1989 > 75% of the trees needed 2 years to reach the pre-disturbance growth level. The year 1955 is outstanding, however, as very cold MJ temperatures co-occurred with a grave moth epidemic, resulting in severe growth reductions of all trees by 33% or more, and half of the birches showing decreases > 76% with respect to the five preceding years. In 1956, still 75% of the trees showed a negative growth response and > 50% of the trees did not recover until 1957 2 years after the event.

The SEA juxtaposes the median growth change behavior of the birch trees to different types of disturbance factors (Fig. 5). The years of low MJ temperatures did not cause any noticeable growth change, but growth even increased in the two following years. In contrast, the pollution events in 1944 caused a growth decrease of 28%, but in contrast to the other disturbance factors, the growth decline was even greater in the following year. Finally, there is a distinct growth response in years with severe moth outbreaks including a median growth reduction of 44%. This negative response was not long lasting as the trees showed a 12% growth increase in the year after the mass outbreak. The combination of moth outbreaks and very low MJ temperatures, as seen in 1955, was most harmful, caused a growth decline > 60%, and trees needed 2 years and more to fully recovery.

Fig. 5
figure 5

Percentage growth change to varying disturbance factors combined in a superposed epoch analysis. The percentage growth changes refer to the 5 years preceding the event years represented by the dashed line. Blue indicates the response to cold May/June-temperatures, green to insect outbreaks, blue/green to a combination of low May/June temperatures and insect epidemics (1955) and orange to pollution through smoke screen actions in 1944

Discussion

A detailed analysis of birch growth responses from the Kåfjord in northern Norway revealed varying sensitivities to disturbance factors including unusual cold mean MJ temperatures, foliage feeding moth outbreaks, and artificial pollution. Growth/climate response analysis showed the sensitivity of downy birch to early season temperatures as we found significant positive correlation with mean MJ temperatures. This association suggests a connection of the onset of radial growth with the beginning of budburst, triggered by MJ temperatures as reported by Klaveness and Wiegolaski (1996), Siljamo et al. (2008) and Karlsen et al. (2009). Accordingly, low mean MJ temperatures can cause a late onset of cambial activity or at least slow down the speed of cambial cell productions and xylem expansions triggering reduced growth rates during such “event” years (Ren et al. 2019). A good example is 1951, during which low growth rates of 75% of the trees likely represents a short growing season. However, we also found highly variable growth responses during most cold event years. This reflects the fact that the trees’ foliation and the resulting ability to do photosynthesis are not necessarily affected by low MJ temperatures, the influence of beneficial environmental conditions during their later growth period can balance the disturbance. Thus, the growth responses were heterogenous among trees, regardless of the intensity of the cold years.

The chemical impact of the Tirpitz’ smoke screens in 1944 was more severe. The artificial fogs caused a growth decline for more than 50% of the trees within the growth period whereas the variability of growth responses in the subsequent year 1945 was very high combined with the highest number of missing rings within the studied time span. The high quantity of absent rings are evidences that the downy birch canopies could (partly) be defoliated by the artificial fogs which took place within three months of the growth period (Hartl et al. 2019). Defoliation of birches reduces the photosynthesis rates massively and confines radial growth accordingly (Treshow and Anderson 1991; Hoogesteger and Karlsson 1992). A study by Haselhoff and Lindau (1903) supports this suggestion: They detail deciduous tree’s ability of enfolding new leaves to ease the impact of acid during the same vegetation period, compared to evergreens (see also Hartl et al. 2019). The high growth variability in the subsequent year 1945 suggests that the impact of smoke screens might have varied among the individual trees, i.e. some of the trees were still harmed in 1945 due to high (multiple) defoliation rates whereas others already recovered. In general, these artificial pollutions as a disturbance had the strongest effects on growth of single trees with respect to aftermaths.

Our results also reveal that moth outbreaks have serious effects on the tree growth. In the 5 years of most severe epidemics, more than 75% of the trees showed distinct growth declines. Similar to pollution effects, the growth decline can be explained by defoliation and reduced photosynthesis. The defoliation of birches in consequence of massive moth outbreaks is common in northern Scandinavia (Eckstein et al. 1991; Karlsson et al. 2004; Jepsen et al. 2008; Babst et al. 2010a). Even the high numbers of missing rings, such as in 1955, can be related to a well-known outbreak of E. autumnata in this area (Tenow 1972; Eckstein et al. 1991; Kirchhefer 1996; Karlsson et al. 2004).

Generally, the moth epidemics are the dominating disturbance factor at our study site: The more severe the intensity of moth outbreaks, the greater the growth decline but also the higher the growth synchronization within the studied area. However, the most severe impact on birch growth was not caused by a single disturbance factor, but by the combination of two: In 1955, mean MJ temperatures were low and likely shifted the onset of the growing season towards June, synchronous to the beginning of the egg hatch of birch moths. Tenow (1972) and Beck et al. (2016) concluded that a temporal synchronization of egg hatch and birch foliation supports larvae feeding and sets the foundation for a severe defoliation event at the site-to-regional scale. In severe cases, this can lead to the entire die-off of a birch site. The co-occurring disturbances in 1955 resulted in the strongest growth decline in our studied birch population More than 50% of the trees still suffered in 1956 and it took the trees 2 years to recover. This conclusion is supported by Kirchhefer (1996), who recorded a moth epidemic in 1955 as the most severe outbreak in Alta region. Tenow (1972) even stated that the intensity of this moth outbreak was exceptional for the whole of Fennoscandia within the last century.

Conclusions

In this study, we focused on single tree growth responses of B. pubescens in years of different disturbance factors in the Kåfjord in northern Norway. Birch tree ring growth from 1901 to 2000 is affected by three main factors: low MJ temperature, (artificial) chemical impacts and moth outbreaks. The growth responses do, however, not depend on the severity of low MJ temperatures. We find a close association between the chemical impact of the Tirpitz concealing tasks in 1944 and a growth decline for three quarter of the birch trees with a strong effect on the cambial activity of the following year. This combination of high growth variability and the highest number of missing rings leads to the fact that some trees were recovering fast while other trees were still suffering. In contrast, severe moth outbreaks of E. autumnata and O. brumata have similar effects on the growth of the entire site. The results demonstrate a positive association between the intensity of documented epidemics and strong growth declines. However, strongest growth declines are recorded if both cold MJ temperatures and moth outbreaks co-occur such as in 1955. A synergy of disturbance factors leads to a synchronization of the birch foliation in the early growing season with the moth’s egg hatch. Considering all disturbances, we are able to rank the factors by their growth impact: Cold MJ temperatures have the least influence on birch growth, followed by the artificial smoke in 1944, and moth outbreaks. The most severe disturbance factor is a combination of low monthly MJ temperatures and moth outbreaks.