Sizeable net export of base cations from a Carpathian flysch catchment indicates their geogenic origin while the 26Mg/24Mg, 44Ca/40Ca and 87Sr/86Sr isotope ratios in runoff are indistinguishable from atmospheric input

Nutrient imbalances may negatively affect the health status of forests exposed to multiple stress factors, including drought and bark beetle calamities. We studied the origin of base cations in runoff from a small Carpathian catchment underlain by base-poor flysch turbidites using magnesium (Mg), calcium (Ca) and strontium (Sr) isotope composition of 10 ecosystem compartments. Our objective was to constrain conclusions drawn from long-term hydrochemical monitoring of inputs and outputs. Annual export of Mg, Ca and Sr exceeds 5-to-15 times their atmospheric input. Mass budgets per se thus indicate sizeable net leaching of Mg, Ca and Sr from bedrock sandstones and claystones. Surprisingly, δ26Mg, δ44Ca and 87Sr/86Sr isotope ratios of runoff were practically identical to those of atmospheric deposition and soil water but significantly different from bedrock isotope ratios. We did not find any carbonates in the studied area as a hypothetical, easily dissolvable source of base cations whose isotope composition might corroborate the predominance of geogenic base cations in the runoff. Marine carbonates typically have lower δ26 Mg and 87Sr/86Sr ratios, and silicate sediments often have higher δ26Mg and 87Sr/86Sr ratios than runoff at the study site. Mixing of these two sources, if confirmed, could reconcile the flux and isotope data. Supplementary Information The online version contains supplementary material available at 10.1007/s11356-024-32866-1.


Introduction
Magnesium (Mg) and calcium (Ca) are essential nutrients performing a variety of functions in the biosphere (Marschner 1995).Magnesium serves as a coordinating cation in chlorophyll and activates enzymes needed for the synthesis of organic molecules, while calcium functions as an intracellular messenger and stabilizes cell walls (Sterner 2002).During pedogenesis, the sources of nutrients change.As soils develop on freshly exposed bedrock, geogenic elements are gradually lost (Wiegand et al 2005;van der Heijden et al 2014;Court et al 2021;Bukoski et al 2021;von Blanckenburg et al 2021;Basdedios et al 2023).When the rates of hydrological export of geogenic nutrients, such as Ca, Mg, K, and P, are higher than the rates of dissolution of bedrock minerals, atmospheric inputs gradually replace rock-derived nutrient sources and help to sustain the productivity of forests.This process has been documented for highly weathered soils in humid, warm regions (Chadwick et al 1999), but also inland ecosystems in the temperate climatic zone have been shown to store substantial amounts of atmospheric nutrients.For example, Miller et al (1993) estimated that 50 to 60% of base cations in the organic soil exchangeable pool and biomass in the Adirondack Mts.(NY, USA) have an atmospheric origin.Belanger and Holmden (2010) calculated that 45 to 90% of Ca present in a Canadian catchment originated from atmospheric deposition.Maher and von Blackenburg (2023) showed that increased water export or decreased regolith thickness due to erosion diminishes the total inventory of nutrients, leading to lower rates of recycling via the plant-available soil reservoir and lower plant growth.Carbonate minerals are associated with considerably faster dissolution kinetics than silicates (Chen et al 2023, and references therein).In the long term, phosphorus (P) is believed to be the master regulator of ecosystems, but Ca and Mg are more mobile than P and may be depleted more rapidly (Chadwick et al 1999).Organic charges that dominate the cation exchange capacity in soils are characterized by a higher affinity for Ca than for Mg (Oursin et al 2023), and, therefore, leaching of Mg through the soil profile is faster, compared to Ca (Maguire and Cowan 2002).
Here we present a combined Ca-Mg isotope and massbalance study of a small forested catchment near the border between the Czech Republic and Poland.Our main objective was to constrain the provenance of Ca and Mg in surface runoff from a steep-slope headwater site on flysch bedrock situated in a highly industrialized part of Silesia.Over the past 60 years, this Central European region has experienced numerous types of stress, including acid rain (ca. 1970-1995;Hruška et al 2023), elevated trace metal deposition (since the 1950s; Voldrichova et al 2014), barkbeetle infestations (1992-1996;2015-2023;Soukhovolsky et al 2022) and climate warming (droughts between 2014 and 2019; Oulehle et al 2021).In the Czech Republic alone, spruce die-back due to anthropogenic acidification affected an area of 500 km 2 , while conifer stands disintegration due to the present-day bark beetle calamity has occurred on an area of 800 km 2 (Zahradnik and Zahradnikova 2023).Atmospheric deposition of alkaline earth elements increased in Central Europe after 1950 as a result of elevated dust emissions from numerous Soviet-style thermal power plants.The resulting replenishment of base cations in ecosystems via wet and dry deposition diminished in the mid-1990s due to technological upgrades in the cluster of soft-coal burning power plants located upwind from the studied region (Kopacek et al. 2016).
Abundance ratios of stable isotopes of Mg and Ca can be used not only as tracers of their inputs into catchment reservoirs but also provide insights into mixing history between solutes, secondary mineral formation, and the recycling of atmospheric and lithogenic inputs by vegetation (Bullen et al 2004;de Villiers et al 2005;Tipper et al 2006Tipper et al , 2008Tipper et al , 2012;;Pogge von Strandmann et al 2008;Jacobson and Holmden 2008;Wimpenny et al 2011;Ma et al 2015;Oehlerich et al 2015;Schmitt 2016;Fahad et al 2016).Isotope ratios are reported in δ values as per mil (‰) deviations in the composition of the sample relative to a standard.The commonly used δ 26 Mg values refer to 26 Mg/ 24 Mg ratios, whereas δ 44 Ca values refer to 44 Ca/ 40 Ca ratios.Isotopically heavy samples are characterized by high δ values, while isotopically light samples have low δ values.
Stable isotope systematics of Mg and Ca in the critical zone (defined as the space between the outer extent of vegetation and the lower limits of groundwater; Brantley et al. 2005) exhibit some similarities.High-temperature processes involving both elements are associated with smaller isotope fractionations than those occurring under lower temperatures (Antonelli et al 2023).Bedrock weathering per se does not fractionate Mg and Ca isotopes (Ryu et al 2011;Cobert et al 2011).Present-day sea water is characterized by relatively light Mg compared to bulk crust and mantle (Guo et al 2019).One exception is represented by marine carbonates preferentially accumulating Mg that is isotopically lighter than the sea water.In contrast, a great majority of geological and environmental samples contain isotopically lighter Ca compared to sea water.On the continents, two well-constrained Mg and Ca isotope fractionation processes occur (Wiggenhauser et al 2022;Cai et al 2022).These are associated with the formation of secondary minerals and uptake by plants (Schmitt 2016;Tipper et al 2016;Teng 2017).River waters typically carry lower δ 26 Mg values than the upper crust because clays preferentially bind isotopically heavy Mg (Brewer et al 2018;Ryu et al 2021;Li et al 2021;Fan et al 2023).Magnesium isotope composition of bulk clay, however, is site-specific, depending on the proportion of isotopically heavy structural Mg and isotopically light exchangeable Mg adsorbed on the phyllosilicate surfaces (Opfergelt et al 2012(Opfergelt et al , 2014;;Fries et al 2019;Hindshaw et al 2020;Zhao et al 2022a).The precipitation of secondary carbonates is an important mechanism driving δ 26 Mg of the dissolved load in rivers to higher values (Zhao et al 2019).Some secondary minerals preferentially bind isotopically light Ca, leaving behind soil solutions and runoff water with isotopically heavy Ca (Brazier et al 2019).Plants often prefer isotopically heavy Mg relative to the nutritive substrate (Black et al 2008;Chapela Lara et al 2017;Schuessler et al 2018).From roots and stem wood toward branches and foliage Mg becomes isotopically lighter (Bolou-Bi et al 2012).In contrast, during assimilation, isotopically light Ca preferentially accumulates in the plant tissues, with Ca in the residual growth medium becoming isotopically heavy.From roots upward to leaves or needles, tree organs typically contain isotopically heavier Ca (Schmitt 2016, and references therein).
In the current study, Ca isotope data have been complemented by 87 Sr/ 86 Sr isotope determinations.These two alkaline earth elements have similar ionic radii and charge and, for many years, were assumed to behave similarly in the critical zone (Åberg et al. 1990;Bedel et al 2016;Burger and Lichtscheidl 2019).Strontium isotope ratios could then serve as a proxy of Ca cycling.Indeed, many studies showed that Sr isotope ratios are a useful tracer of hydrological processes (Nuruzzama et al 2020;Demonterova et al 2022).Recent work, however, has indicated that Ca and Sr cycling in biological processes is de-coupled.Strontium isotope data in Ca-nutrient studies should thus be used with caution (Drouet et al 2005;Brenot et al 2008a,b;Ryan et al 2018;Bouchez and von Blackenburg 2021;Nguyen et al 2023).
Cenki-Tok et al ( 2009) pointed out that "determination of the atmospheric contribution of base cations to catchment surface water using mass budget calculations that are based only on annual element fluxes must be taken with caution."In the current study, we focused on a comparison between Mg, Ca and Sr input-output budgets, which represent the traditional "black-box" approach, and stable isotope fingerprinting which additionally reflects within-catchment processes.We hypothesized that the Mg, Ca and Sr isotope signatures of runoff would mainly reflect the isotope composition of bedrock in the case of a large net hydrological export of the studied elements.In the case of net accumulation of atmospheric base cations in the catchment and/or a balanced mass budget, the isotope composition of runoff would indicate a mixture of geogenic and atmospheric sources.We also hypothesized that shallow soil water would mostly carry isotope signatures of present-day atmospheric Mg, Ca and Sr.The isotope signatures of soil water would differ from the isotope signatures of bulk soil and runoff, both of which would, under net Mg, Ca and Sr export, converge to those of bedrock.Collectively, the flux and isotope data on base cations in small catchments can help to evaluate the sustainability of forests in the present era of global change.

Study site
The small catchment Cervik (CER) is located in the Moravian-Silesian Beskydy Mts. in the eastern Czech Republic (Central Europe), close to the borders with Poland and Slovakia (Fig. 1; Table 1; Adamova 1986;Novak et al. 1996Novak et al. , 2005;;Oulehle et al 2017;2021;Vicha 2019).Geologically, the area belongs to the folded Upper Cretaceous flysch complexes of the Silesian Unit (Outer Western Carpathians; Golonka et al. 2019).The bedrock of the studied catchment is formed by sandstone with subordinated amounts of claystone.It is part of the Istebna Formation (Campan-Paleocene; Appendix I; for geological sketch, see Fig. S1 of the Electronic Annex).North and east of CER, the older Godula Formation crops out.The maximum elevation of the catchment is 961 m.Deluvial  (Novak et al 2020a(Novak et al , 2020b(Novak et al , 2020c(Novak et al , 2021(Novak et al , 2023a(Novak et al , 2023b)).Sulfur deposition data by Czech Hydrometeorological Institute sediments, along with cambisol and stagnosol soils, cover nearly 70% of the surface of the catchment (Fig. S1).The soil depth fluctuates between 40 and 80 cm.Mature Norway spruce (Picea abies) stands cover 88% of the catchment's area, with clearings and young deciduous trees encompassing the remaining 12%.Atmospheric deposition of sulfur (S), directly measured in the catchment since 1994, reached 35 kg ha −1 yr −1 in 1996 and decreased to 6 kg ha −1 yr −1 in 2021 (Oulehle et al 2021; Fig. S2).The nearly sixfold decrease in annual S deposition over the 28-year time period followed a peak in nation-wide S industrial emission and atmospheric deposition rates in the late 1980s (Hruška et al 2002;Oulehle et al 2006;Hunova et al 2014Hunova et al , 2018;;Hunova 2020).Acidifying compounds deposited at the study site originated from (i) North Bohemian soft-coal burning power plants located 300 km to the west (Novak et al 2000), (ii) the Silesian industrial area near Ostrava (eastern Czech Republic, stone coal and metallurgy; 40 km to the north) and (iii) Olkusz industrial area (southern Poland, base metals and thermal power plants; 140 km to the northeast; Fig. 1a; Buzek et al 2017Buzek et al , 2023;;Prechova et al 2020Prechova et al , 2023)).After 1994, reactive nitrogen (N r ) deposition at CER fluctuated between 5 and 14 kg ha −1 yr −1 (Oulehle et al 2021; Fig. S2) and lacked clear-cut temporal trends.The mean pH of spruce throughfall was slightly lower than the mean pH of open-area precipitation in most years.During the 1994-2021 time period, pH of both spruce throughfall and open-area precipitation increased almost linearly from 3.9-4.4 to 5.3-5.5.The pH of runoff (mean of 6.8 for the years 1994-2021) was significantly higher than pH of atmospheric deposition.In contrast to atmospheric deposition, pH of runoff did not increase over time.

Sampling
The locations of CER sampling sites are in Fig. 1b.The Mg and Ca input-output flux estimates are based on monthly sampling of open-area precipitation, spruce canopy throughfall, and surface runoff in hydrological years 1994-2021.
Strontium input-output fluxes are based on monthly sampling in hydrological years 2014-2021.Samples of both types of atmospheric deposition and runoff for Mg, Ca, and Sr isotope analysis were collected between November 2020 and October 2021 (Tab.S1 and S2 in the Electronic Annex).Precipitation collectors were described by Fottova and Skorepova (1998).Soil water for δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr determinations was sampled monthly between June and November 2021.The mineral soil solutions were collected in triplicate by Prenart suction lysimeters from a depth of 60 cm and pooled bi-monthly prior to analysis (Tab.S1 and S2).
Solid samples included bulk rock, forest floor (L, HF), soil from the depth intervals 0-20, and 20-40 cm, spruce xylem, class-1 spruce needles, and spruce roots < 2-mm indiameter.Each sample type, except for bulk bedrock, was collected at 10 to 20 sites around soil pits no. 1, 3 and 5 (solid squares in Fig. 1b and Fig. S1) and combined to form three replicates.Bedrock samples (open squares in Fig. 1b and Fig. S1) were collected as follows: Thirty sandstone samples were collected in the area of soil pits no. 1, 3 and 5 and combined to form three replicates for isotope analysis.Ten claystone samples were collected from each of the three small outcrops in the stream and combined to form three replicates for isotope analysis.Conglomerate was identified and sampled only in one outcrop.Bedrock for a thin-section study was sampled at sites I, II, III and IV (solid triangles in Fig. 1b and Fig. S1).Spruce xylem was sampled 1.4 m above ground using a tree-ring corer.Tables S1 and S2 give the numbers of individual sample types analyzed for Mg, Ca, and Sr isotopes.
Pool-size inventory of Mg and Ca in soil was based on 10 soil pits excavated below a 0.5-m 2 area in 2014 (Fig. 1b).Pool sizes of Mg and Ca were calculated per unit area for aboveground and belowground vegetation, forest floor, and four soil depth intervals (0-10, 10-20, 20-40, and 40-80 cm).

Sample processing and analysis
Detailed methodology of sample processing and analysis can be found in Novak et al (2023a).Here, we briefly summarize individual steps.
Homogenization Bedrock fragments (three composite samples of sandstone, three composite samples of claystone, and one sample of conglomerate) were ground in a jaw crusher and milled.The samples were used for Mg, Ca and Sr isotope analysis and whole-rock silicate analysis.The < 2-mm soil fraction was prepared for isotope analysis by sieving and milling to a < 60-μm grain size.Forest floor samples were homogenized without sieving.Stratified soil samples from 10 soil pits (Fig. 1b) were air-dried, sieved (< 2 mm) and homogenized.Exchangeable soil pools of Mg and Ca were determined using BaCl 2 extracts (Oulehle et al 2010).In contrast, total Mg and Ca pools in aboveground vegetation (spruce xylem and needles) and belowground vegetation (spruce roots) were used for the inventory.

Solid sample dissolution
Samples of bulk bedrock and soil were dissolved in HF and HClO 4 and re-dissolved in concentrated HNO 3 .Biomass/organic-rich materials were ashed at 550 °C for 8 h, digested in H 2 O 2 and HNO 3 , and further treated using the same steps as in the case of bedrock and soil.

Sample purification and mass spectrometry
For the chromatographic separation of Mg, we used the methodology by Pogge von Strandmann et al ( 2011) modified to a onecolumn procedure (Pokharel et al 2017).Separation of Ca was carried out using the methodology described by Holmden (2005) and Holmden and Belanger (2010).Strontium processing was performed according to Erban Kochergina et al (2022).
The Mg isotope ratio measurements were conducted using a Neptune MC-ICP-MS (Thermo Scientific).Details are in Pokharel et al (2017).The Ca and Sr isotope ratio measurements were conducted using a Triton Thermal Ionization Mass Spectrometer (TIMS; Thermo Scientific) using the methods described by Holmden and Belanger (2010) and Erban Kochergina et al (2022), respectively.
PXRD whole-rock analysis Powder X-ray diffraction data were collected on a Bruker D8 Advance powder diffractometer in a Bragg-Brentano geometry.A qualitative phase analysis was performed using the DIFFRAC.Eva software (Bruker 2015).The quantitative phase analysis was carried out using the Rietveld method (Topas 5 Program; Bruker 2014).The detection limit of the method was between 0.2 and 1.0 wt.

Analysis of liquid samples
Element concentrations in natural solutions were determined by FAAS (Mg and Ca), Q-ICP-MS (Sr), and ICP OES (Ba and Rb).Aliquots for Mg and Ca isotope analysis were filtered (0.45 μm), evaporated to dryness and treated with concentrated HNO 3 .Rare solutions where solid residues occurred were re-dissolved in concentrated H 2 O 2 and HNO 3 .Water samples for Sr isotope analysis were filtered, acidified with HNO 3, evaporated to dryness, and the residue, if present, was dissolved in a H 2 O 2 -HNO 3 mixture, evaporated to dryness and dissolved in HCl.

Carbon analysis
Carbonate C was released from the sample by reaction with concentrated phosphoric acid.Dried CO 2 was analyzed in an ELTRA CS 500 Analyzer with infrared detection (precision of 6%).Total C and S were analyzed following the thermal decomposition of the sample (precision of 5 and 4%, respectively).

Mg and Ca inventories and input-output budgets
Annual catchment-level budgets of the studied elements were calculated as atmospheric inputs minus runoff outputs.Atmospheric input fluxes were weighed by the percentages of the catchment's surface covered by spruce stands and clearings, respectively.

Results
Mg, Ca and Sr isotope systematics All three studied elements exhibited an analogous isotopic relationship among four types of water samples.For each of the Mg, Ca, and Sr systems, the isotope ratios of open area precipitation, spruce throughfall, soil water and runoff were statistically indistinguishable (p > 0.05; Fig. 2).In contrast, the isotope composition of Mg, Ca and Sr in bedrock was significantly different from the isotope composition of Mg, Ca, and Sr in runoff and soil water (p < 0.05).In Fig. 2a, b, c, the horizontal gray bands highlight the contrast between the isotope composition of the studied base cations in bedrock and in water samples.
The δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr ratios of spruce needles were statistically indistinguishable from δ 26 Mg, δ 44  between those of atmospheric deposition and bulk soil (Fig. 2a).Bulk soil δ 26 Mg values were indistinguishable from bedrock (p > 0.05).Similarly, the δ 44 Ca values of spruce roots and xylem plotted between the δ 44 Ca values of atmospheric deposition and bulk soil (Fig. 2b).The δ 44 Ca values of soil were the lowest in the whole studied system, significantly lower (p < 0.05) than those of bedrock.In the case of the 87 Sr/ 86 Sr ratios, fine roots were indistinguishable from soil water.At the same time, there was no significant difference between the 87 Sr/ 86 Sr ratios of xylem, bulk soil and bedrock.
The time-series in Fig. S3a,b,c show nearly constant δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr ratios of runoff (their temporal variability was close to the uncertainly of mass spectrometric measurements).The variability in δ 26 Mg and δ 44 Ca values of atmospheric deposition was also relatively low (within 0.6 ‰ and 0.4 ‰, respectively).The variability in 87 Sr/ 86 Sr isotope composition of atmospheric deposition was large (0.08).No seasonality in δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr ratios of liquid samples was observed.An 87 Sr/ 86 Sr vs. 1/Sr concentration plot showed a relatively large negative correlation (Fig. S3d; R 2 = 0.72).The two extreme values in Fig. S3d (bottom right) corresponded to the wettest and driest month of the year respectively (April and September).
With an increasing depth, Mg isotope composition of bulk soil did not change outside of the uncertainty of mass spectrometric measurements; δ 44 Ca values decreased, and 87 Sr/ 86 Sr ratios increased (Fig. 3).
Isotope measurements in individual CER samples are in Tab.S1 and S2 and Fig. S4.In Fig. S5, we give two examples of endmember mixing plots according to Capo et al (1998; δ 26 Mg vs. Ca/Mg and 87 Sr/ 86 Sr vs. Ca/Sr).In both cases, runoff plotted inside the triangle is defined by throughfall, soil water and bulk soil.In both models, soil water contributed over 90% to base cations in runoff.
Temporal trends in catchment input-output fluxes For all the studied base cations, annual fluxes decreased in the order: runoff > > throughfall > open-area precipitation (Fig. 4).This relationship was valid throughout the entire 28-year monitoring period for Mg and Sr. Calcium flux via open-area precipitation was slightly higher than Ca flux via throughfall in 2001 and 2021.Runoff Mg flux between 1994 and 2021 per unit area was on average 18 times higher than rainfall Mg input flux.Over the same time period, runoff Ca flux was, on average, six times higher than rainfall Ca flux.The average Sr runoff flux (2014-2021) was 12 times higher than the Sr rainfall flux.Runoff Mg flux was nine times higher than throughfall Mg flux, Ca runoff flux was nine times higher than the throughfall Ca flux, and runoff Sr flux was nine times higher than throughfall Sr flux.
Fig. S6 gives a comparison of average vegetation-type weighted annual input and output fluxes of Mg, Ca and Sr throughout the observation period.Conversion of the measured element fluxes (left bar) from Fig. 4 to total atmospheric deposition (middle bar) using Na concentrations eliminated the effect of scavenging and leaching of base cations by canopy (Kopacek et al. 2016).The correction had only a minor effect on the base cation fluxes.At the catchment scale, throughout the observation period, after the ratio of forested areas and clearings (88:12) was considered and the Na-correction applied, the average annual runoff flux was 15, 5 and 10 times higher for Mg, Ca, and Sr, respectively, compared to the atmospheric input.
A significant decrease in runoff and throughfall fluxes over time was observed for Mg and Ca (p < 0.01; Fig. 4a,b).In the case of open-area precipitation, borderline statistical significance was found (p = 0.043 for Mg, and p = 0.055 for Ca).

Bedrock chemistry and mineralogy
Mean MgO concentration in bedrock sandstone was approximately eight times lower than MgO concentration in claystone (0.21 vs. 1.62 wt.%; Table 2).Mean CaO concentration in bedrock sandstone was approximately three times lower than CaO concentration in claystone (0.15 vs. 0.52 wt.%; Table 2).Concentrations of Fig. 3 Vertical trends in δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr isotope ratios in CER soils.Means and standard errors are given MgO and CaO in conglomerate was close to those in sandstone.The contents of CO 2 indicating the potential presence of carbonate were small in all three bedrock types (Table 2), close to the detection limit (0.01 wt.%).According to the geological sketch in Fig. S1, claystone built up a mere 4.7% of the catchment's area; its occurrence under the deluvial sediments is not known.
PXRD analysis provided semi-quantitative mineralogical composition of CER sandstone and claystone (Table 3).Quartz dominated sandstone (83%), with K-feldspar and albite making up 8 and 6%, respectively.Phyllosilicates constituted less than 4%.Phyllosilicates were the main component of claystone (40%).The applied methodology did not permit to distinguish between mica, illite, and mixed-layered minerals.This group is listed as mica in Table 3. Quartz and kaolinite were the second and third most abundant mineral phases of the studied claystone (< 35 and < 11%).Albite and K-feldspar represented 7 and 4%, respectively.
Thin-section study revealed additional accessory minerals and rock fragments, such as glauconite, ilmenite and zircon in the sandstone (Fig. S7-I), and glauconite and lithic clast of magmatic and metamorphic rocks, including rhyolite (Fig. S7-II, III) in the conglomerate.Fig. S7-IV captured well-preserved silicified relics of benthic foraminifera of the genus Bathysiphon.

Pool size inventory
The aboveground vegetation pool of Mg was only three times larger than the annual runoff flux of Mg (39 vs. 14 kg ha −1 yr −1 ; Fig. S8).The exchangeable Mg pool was slightly lower in deep soil (40-80 cm below the surface) than in all soil horizons above the depth of 40 cm (31 vs. 35 kg ha −1 yr −1 ).The aboveground vegetation pool of Ca was 11 times larger than the annual runoff flux of Ca (377 vs. 33 kg ha −1 yr −1 ; Fig. S9).The exchangeable Ca pool was similar in deep soil (40-80 cm below the surface) and in soil horizons above the depth of 40 cm (149 vs. 150 kg ha −1 yr −1 ).

Apparent discrepancy between the isotope and flux approaches
In light of the 5-15 times larger export of Mg, Ca, and Sr via CER runoff, compared to atmospheric input (Fig. S6), we expected isotope similarity between runoff and bedrock rather than between runoff and atmospheric deposition.The opposite was true (Figs. 2 and 4).Below, we will summarize arguments in favor of the isotope hypothesis ("CER runoff mostly contains atmospheric base cations, flux comparisons are misleading") and arguments in favor of the mass-balance hypothesis ("CER runoff mostly contains geogenic base cations, the isotope fingerprinting is incomplete").Then, we will suggest further steps to reconcile data in Figs. 2 and 4.
CER contrasts with the previously studied Central European catchments that were also underlain by base-poor bedrock (< 2.1 wt.% of MgO and < 1.8 wt.% of CaO) and also exhibited net export of Mg, Ca, and Sr (Fig. 5;Novak et al 2020a, b;Novak et al 2023a,b; for catchment locations see Fig. 1a).Deposition fluxes of the studied base cations at CER were relatively low, closer to deposition fluxes in the rural catchments previously studied in southern and western Czech Republic (LIZ, LYS) than to deposition fluxes in the industrially polluted northern catchments (UDL, UHL).Calcium isotope systematics had been statistically evaluated at four sites.At all these sites, Ca isotope ratios of runoff and bedrock were statistically indistinguishable.Magnesium isotope systematics had been statistically evaluated at two sites.At both, Mg isotope ratios of runoff and bedrock were also statistically indistinguishable.These comparisons highlight the unusual Mg and Ca isotope systematics at CER.

Arguments in favor of the isotope hypothesis
A large body of literature exists on Mg, Ca and Sr isotope composition of marine carbonates (Tang et al 2008;Li et al 2010;Zhao et al 2022a, b).All three elements exhibit a wide range of isotope ratios.The δ 26 Mg values of marine carbonates span more than 5 ‰ (− 5.3 to 0.0‰; Higgins and Schrag 2015; Liu and Li 2019), while the δ 44 Ca values of recent planktic foraminifera alone span nearly 4‰ (− 2.0 to 1.8 ‰; Gussone et al 2016).The 87 Sr/ 86 Sr ratios of marine carbonates depend on the age.For the 160-26 Ma period corresponding to the formation of the Carpathians, the 87 Sr/ 86 Sr ratios are from 0.706 to 0.708 (Nadaskay et al 2019;McArthur et al. 2020).Isotope signatures of marine carbonates thus span the entire range of CER ecosystem values for δ 44 Ca, and nearly the entire range of CER ecosystem values for δ 26 Mg (y-axes in Fig. 2a, b).The 87 Sr/ 86 Sr ratios in the CER ecosystem lie outside the marine carbonate values (y-axis in Fig. 2c).In principle, a missing lithogenic source of base cations, such as sedimentary carbonate layer that was not isotopically analyzed, could shift the combined bedrock isotope signature toward that of runoff.However, it is unlikely that rapid weathering of the unknown carbonate source would simultaneously shift all three investigated isotope ratios (δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr) to values practically identical to local atmospheric deposition.Indeed, the improbability of such simultaneous shifts is the strongest argument for a major role of recent precipitation in CER runoff.
Unlike in previously studied Czech catchments, at CER, we were unable to isotopically analyze Mg, Ca and Sr of individual rock-forming and accessory minerals.In previously studied catchments underlain by granite, orthogneiss and paragneiss, we separated up to seven minerals from ground matrix and analyzed them isotopically.We found largely variable δ 26 Mg (-3.3 to 0.0 ‰), δ 44 Ca (-0.6 to 1.0 ‰), and 87 Sr/ 86 Sr (0.714 to 2.05) ratios (Novak et al 2020a(Novak et al , 2020b(Novak et al , 2020c;;2023a, 2023b).These studies confirmed previous reports of large isotope variability between minerals in crystalline rocks (Nezat et al 2007;Farkas et al 2011;Ryu et al 2011;Drouet et al 2015).At CER, it was not possible to separate individual minerals from ground arkose sandstone due to ubiquitous intergrowths, and it was not possible to mechanically separate individual clay minerals from claystone.In Appendix III, we have listed empirical dissolution rate constants of individual minerals present in CER sandstone and claystone.Even without knowledge of reaction surfaces and kinetic calculations, it is obvious that individual phases may, step by step, release into solution Mg, Ca and Sr carrying distinct, mutually contrasting isotope compositions.Yet, we can here use the same argument as we have done in the case of variable isotope composition of basic cations in hypothetical missing carbonate sediments: It remains improbable that differences in weathering regimes of individual CER minerals would simultaneously push δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr ratios of dissolved geogenic Mg, Ca and Sr in runoff from the whole-rock values to the almost exact atmospheric deposition values, as seen in Fig. 2.
CER is characterized by steep slopes (an average of 30% according to ArcGIS).Previous biogeochemical research has shown that some atmospherically derived nutrients exhibit shorter mean residence time within steep catchments than at mild-sloping sites (Sueker et al 1999;Michel et al 2000;Novak et al 2004).Hydrological control may cause smaller interaction of atmospherically deposited nutrients with the forest ecosystem at steeper sites, favoring their faster export.
Isotope studies have previously shown that runoff at some of the 14 headwater catchments of the GEOMON network contains atmospheric rather than geogenic S, Pb and Zn (Novak et al 2005;Bohdalkova et al 2014;Andronikov et al 2021).These chemical elements are relatively abundant in crustal rocks, but, at the same time, industrial activities have caused their elevated atmospheric fluxes.Again, rather than suggesting an analogy between the hypothetical high export of atmospheric Mg, Ca and Sr and the isotopically documented predominance of pollutants S, Pb and Zn in catchment runoff in industrial areas, we wish to stress that extremely large reserves of elements in the parental rock do not automatically imply their dominance in runoff.

Arguments in favor of the mass-balance hypothesis
Hydrochemical mass balances at the small-catchment level are associated with uncertainties stemming, inter alia, from a limited number of sampling sites along the slope, discontinuous sampling of runoff, difficult quantification of horizontal deposition, and site-specific interactions between rainfall and tree canopy (Probst et al 1990;Miller et al 1993;Moldan and Cerny 1994;Kohler et al 2015;Schwartz et al 2022).At CER, caution must be exercised given the fact that all deposition samplers were installed at the lowest elevation of the catchment (Fig. 1b).With an increasing elevation, precipitation totals generally increase and may cause additional element input that has not been considered.Havel et al. (1996) found a linear 50 mm increase in annual precipitation total for each additional 100 m elevation in the Ore Mts., northern Czech Republic.Kopacek et al. (2012) showed that at elevations over 900-m the linear increase in precipitation totals is better approximated by an exponential regression and inferred an average 105 mm increase per 100 m elevation for the 300-1550 m a.s.l.span.Havel et al. (1996) monitored concentrations of Mg 2+ and  Ca 2+ in precipitation along a 500 m elevational gradient in the northern Czech Republic and concluded that the concentrations at individual elevations were statistically indistinguishable.The 324 m difference between the lowest and highest point within CER (Table 1), along with the above previously published findings regarding changing amounts of precipitation but nearly constant Mg and Ca concentrations in rainfall at different elevations, indicate that the uncertainty in catchment level Mg and Ca atmospheric input fluxes at CER is probably low (cf., Pu et al 2023).
Importantly, Oulehle et al. (2017) reported a longterm chloride (Cl − ) budget for CER.Internal sources of Cl − , mainly weathering of the parental rock, are minor.Atmospheric deposition of Cl − is a mixture of marineaerosol derived and industrial-emission derived Cl − .Chloride is believed to be a conservative tracer in forested catchments, moving rapidly through the ecosystem with small to negligible biogeochemical interactions.At CER, decreases in Cl − deposition and export were significantly related, and the Cl − mass budget over the 21-year period was balanced (Oulehle et al 2017).The construction of input-output mass balances at CER is thus reasonably well-constrained.
The one-month sampling interval for chemical and isotopic analysis of CER runoff introduces an uncertainty, especially during the snowmelt period (February-April).While runoff water flux is monitored continuously, the Mg, Ca and Sr isotope composition between runoff samplings was poorly constrained.The δ 26 Mg, δ 44 Ca, and 87 Sr/ 86 Sr ratios of runoff during the peak snowmelt between samplings can be assumed to converge to those of recent precipitation because the portion of runoff that was derived directly from snowpack had been barely in contact with the soil and vegetation and the rates of most biogeochemical processes decrease with decreasing temperatures.Even without continuous isotope analysis of spring runoff, data in Fig. 2 already coincide with isotope signatures of precipitation.Kram et al. (2012) observed decreasing Mg 2+ and Ca 2+ concentrations at higher runoff water fluxes in three small forested catchments in the western Czech Republic.In principle, unmeasured lower concentrations of base cations in runoff between spring samplings may affect the calculation of annual Mg 2+ and Ca 2+ export fluxes.Such overestimation of export fluxes of the studied elements, however, will likely be minor (cf., Alewell et al 2004).Overall, we conclude that the 5 to 15 times higher export of the studied base cations from CER relative to the atmospheric input fluxes is too high to be explained by the above-discussed uncertainties in budget calculations.Robustness of the input-output flux estimates is the strongest argument for the predominantly geogenic origin of runoff Mg, Ca and Sr.
Two distinct patterns seen in the CER data (Figs. 2 and  4) do not contribute to testing the geogenic vs. atmospheric runoff generation hypotheses: (i) decreasing Mg and Ca runoff fluxes in Fig. 4 are not a result of predominating atmospheric origin of these elements in runoff in an era of easing air pollution.The runoff flux decrease is instead a consequence of retreating acidification.Declining concentrations of strong acids in soil solutions (mostly H 2 SO 4 of atmospheric origin) result in lower leaching of base cations from the soil cation-exchanger (cf., Heliwell et al 2014;Garmo et al. 2014).The input-output timeseries in Fig. 4a,b have different slope.Changes in the exchangeable Mg and Ca soil pool may involve both geogenic and atmospheric base cation sources, and therefore, changing export fluxes per se do not carry information in atmospheric vs. geogenic runoff origin; (ii) the Mg and Ca isotope composition of xylem as the main plant reservoir plots between the isotope signatures of throughfall and soil (dotted lines in Fig. 2a,b).Hence, simple mixing between soil-derived and atmospheric Mg and Ca may provide a sufficient explanation of xylem isotope ratios.At the same time, it is well established that assimilation by trees prefers isotopically heavy Mg and isotopically light Ca (Page et al 2008;Cobert et al 2011;Hindshaw et al 2013;van der Heijden et al 2015;Kimmig et al 2018).The spruce stands at CER are relatively young ( ≤ 70 yrs), and tree biomass continues to build up.The residual isotopically light Mg and isotopically heavy Ca may shift the isotope signatures of runoff in the direction that has actually been observed (Fig. 2).It is therefore impossible to separate mixing and biological fractionation in Fig. 2. Predominance of residual Mg and Ca in runoff following isotope-selective partial uptake of these elements by plants cannot be taken as an exclusive explanation of isotope patterns in runoff (Fig. 2a,b).
Kinetic-limited chemical erosion that likely predominates at CER over supply-limited erosion (Riebe et al 2004) may have a complex effect on the provenance of base cations in runoff: limited accumulation of soils on the steep slopes results in limited capture of atmospheric inputs by surface soil horizons, and relatively fast export of partly weathered regolith rich in minerals stable under environmental conditions exposes unweathered substrate rich also in unstable minerals.The first phenomenon may enrich runoff in atmospheric base cations, the second phenomenon preferentially supplies stream water with base cations from reactive mineral phases.

Search for the unknown Mg, Ca and Sr source
Based on the abundances of individual minerals (Table 3) and solid phase stoichiometry, mica combined with illitegroup phyllosilicates and plagioclase are the main potential Ca sources in CER claystone.Mica combined with illitegroup phyllosilicates and chlorite are the main Mg sources in CER claystone.In alkaline-earth elements poor sandstone, mica, illite-group minerals and plagioclase are the main Ca source, while mica and illite-group minerals are the main Mg source.All these Mg and Ca sources contributed to the whole rock isotope ratios in Fig. 2. If Mg and Ca in silicate analyses in Table 2 are attributed to individual minerals, starting with carbonates expressed as CO 2 , the noncarbonate phases still contain 96.9 to 99.7% of all present Mg and Ca.In other words, we found negligible amounts of an easily soluble (Ca, Mg) CO 3 phase in the sampled rocks.
Flysch sediments, or turbidites, forming CER parental rock resulted from a combination of fluidal and sediment gravity flow.They were formed as underwater avalanches of uncemented clastic sediments that slid down the steep slope of the continental shelf edge into the deep ocean.Their lithology is highly variable but rarely includes carbonate material.Fig. S10 depicts the lithology of two drill cores (453 and 20705 m deep) which are located at Stare Hamry, 4 km east of CER.These profiles illustrate (i) the large spatial variability of the lithology and (ii) absence of carbonate strata.Despite a mere 30 m distance between the two drill cores, rapid changes in rock types in the shorter core (sandstone, shale, claystone) contrast with the monotonous stratigraphy of the longer core (mainly sandstone, only one layer of claystone intercalations).Neither of these drill cores contains carbonate-rich facies, such as marlstones.Mencik et al. (1983), Golonka and Picha (2006) and Stranik et al. (2021) published a detailed lithology of the older Godula and younger Istebna Formation.Within the 2900-3100 m thick Godula Formation, the claystones were carbonate-free, but a limited occurrence of calcareous sandstones was reported.Within the 1000-1200-m-thick Istebna Formation, arkose sandstones did not contain carbonate.In contrast, Jurassic and Devonian carbonate cobbles were occasionally present in conglomerates.According to Adamova (1986), dark claystone in the Silesian Unit may also contain subordinated amounts of carbonates (calcite, dolomite and siderite).While we were unable to identify any easily weathering carbonate rocks at CER, their occurrence in the broader region cannot be ruled out.On a smaller scale, small amounts of remobilized carbonate may have remained in the Silesian Unit after silicification of the carbonate fossils, such as foraminifera in Fig. S7-IV.
We suggest that there may be an analogy between carbonatederived basic cations in runoff from catchments underlain by heterogeneous flysch sediments and the predominance of carbonate-derived basic cations in rivers draining silicate catchments (Moore et al 2013;Fantle and Tipper 2014;Jacobson et al 2015;Xu et al 2022).Small volumes of carbonate veinlets were shown to supply a large proportion of riverine Ca 2+ in granitoid terranes.In the Carpathian Flysch sediments, we do not expect hydrothermal calcite veinlets, but subordinated amounts of carbonate components of conglomerates, marlstones and calcareous sandstones may contribute to reconciling CER data in Fig. 2 and 4.
Fig. S11 outlines a hypothetical isotope scenario of silicate and carbonate dissolution and its possible effect on δ 26 Mg and 87 Sr/ 86 Sr ratios of runoff (cf., Zhao et al 2022a).The scenario assumes a global average δ 26 Mg and 87 Sr/ 86 Sr isotope ratios of silicate and marine carbonate rocks for the studied area.For both Mg and Sr, the silicate mixing endmember is isotopically heavier, and the carbonate mixing endmember is isotopically lighter than Mg and Sr in CER runoff.Such mixing, if empirically confirmed, would corroborate the above-discussed mass balance hypothesis stating the dominance of geogenic sources of base cations in CER runoff.Global mean δ 44 Ca values are similar for silicate and carbonate rocks (Blättler and Higgins 2017), and hence, Ca isotopes cannot be used for a model analogous to Fig. S11.

Relationship between Sr isotope ratios in runoff and vertical soil profiles
Generally, a linear relationship between 1/concentration and isotope ratios in stream water (Fig. S3d) suggests mixing between two sources of an element (e.g., Zhao et al 2022b).Comparison between Figs.S3d and 3c could, in principle, indicate whether Sr in runoff in a wet month converges to a shallow or a deep soil horizon as its potential source.Unfortunately, the two low Sr-concentration, low-87 Sr/ 86 Sr runoff samples come from hydrologically contrasting times of the year (wet April and dry September).Both these runoff samples converge to the topmost soil horizon having the lowest 87 Sr/ 86 Sr ratio, which, however, remains strikingly different from the runoff Sr isotope signature (0.7225 vs. 0.714-0.715).Longer time-series of isotope data would be needed to more rigorously address possible within-soil sources of runoff Sr.

Conclusions
In a small forested catchment underlain by turbidite sediments, we found nearly identical values for atmospheric deposition and runoff in each of the δ 26 Mg, δ 44 Ca and 87 Sr/ 86 Sr isotope systems.In contrast, each of the Mg, Ca and Sr isotope ratios was significantly different between bedrock and runoff.This was counter-intuitive since longterm hydrogeochemical monitoring of catchment fluxes showed a large net export of each of the three studied base cations via runoff.Over the last 28 hydrological years, the average Mg input via atmospheric deposition was 15 times lower than Mg runoff flux, Ca deposition was 5 times lower than Ca runoff, and Sr deposition was 10 times lower, compared to its hydrological export.Bedrock sandstone and claystone contained only traces of easily dissolving carbonates.We explored the possibility of sizeable export of legacy pollutants Mg, Ca, and Sr that would have been deposited in the high-pollution years 1950-1995.Our mass-budget monitoring started in 1994, and already then, atmospheric Mg, Ca and Sr inputs were small relative to their runoff fluxes.Therefore, the remobilization of legacy pollution from power plant dust temporarily stored in the soil and biomass could not explain the large net export of base cations during our study.We discussed possible sources of uncertainty in both isotope and non-isotope flux data.Uncertainties in catchment-level mass budgets caused by placing precipitation collectors only at the foot of the steep catchment, unknown upslope gradients in base cation concentrations in atmospheric deposition, and discontinuous runoff sampling for concentration and isotope analysis were probably too low to minimize the difference between export and deposition flux of Mg, Ca and Sr.An unknown deeper source of carbonate Ca, Mg and Sr that has not been analyzed isotopically may help to reconcile the apparent isotope-flux discrepancy.Still, it appears to be unlikely that the isotope composition of unknown carbonates would shift the isotope compositions of all three studied elements in runoff almost exactly to those of present-day atmospheric deposition.

Appendix 1. Lithology and mineralogy of the Godula and Istebna Formations
The Silesian Unit builds an outer part of a synorogenic accretionary prism of the Carpathians.A wedge-shaped structure was formed during the closure of oceanic Tethys domains and the subsequent collision of the central Carpathian plate with the Brunovistulian foreland.The Mesosoic pre-flysch, syn-rift, and syn-orogenic flysch units of the Outer Western Carpathians were folded, thrusted on the foreland and uplifted to the present-day mountain-range position.The older Godula Formation is composed of rhythmic alteration of sandstones and claystones/siltstones, with conglomerate intercalations.The younger Istebna Formation is composed of slopeapron and submarine-fan sandstones and arkosic sandstones with conglomerate intercalations.The flysch deposits are intercalated by dark marlstones and claystones and dark shales.The bedrock at CER is formed by a brachynticline of flysch sediments of the Istebna Fomation, whereas rocks of the Godula Formation crop out north and east of CER.The lithology of the psammitic members of both formations is similar.Istebna sandstones are mostly coarser-grained in comparison with the sandstones of the lower segment of the Godula Formation.The Istebna Formation sandstones contain more feldspars and conglomerate interbeds.The presence of pebbles and heavy minerals (assemblages dominated by garnet and zircon) indicates that the dominant source of basin turbidites originated in the crystalline basement of the Silesian cordillera.In addition to quartz pebbles and stable silica-rich clast (cherts, muscovite-biotite quartzites), polymict conglomerates contain frequent unstable components, such as orthogneisses, felsic volcanics (tuffs and spherulithic lavas), aplitic granitoids, and large feldspar clast.Recycled glauconitic sandstones are present scarcely.Siderite is the prevailing cementing mineral in psamites with low carbonate content.Godula Formation sandstones are richer in diagenetic glauconite than Istebna Formation.Dark greenish to dark grey/ brown aleuropelites/shales rich in Fe oxides (hematite, limonite) and siderite contain sulfides (pyrite) and carbonates only as accessory phases.Pelitic and aleuritic members of both formations are similar.In addition to dominant clay minerals, micas, chlorite, and glauconite, these members contain a silty siliciclastic admixture (quartz, feldspar).Organic matter contents (1-8%) tend to be higher in aleuropelites and shales of the Istebna Formation.

Appendix 2. Isotope mixing models for Mg, Ca and Sr in forest ecosystems
To constrain the sources and relative contributions of Mg, Ca and Sr into local ecosystem compartments, we calculated theoretical mixing trends for selected 'end-members' (see perimeters of grey shapes in Fig. S5).The mixing trends for δ 26 Mg values, δ 44 Ca values, 87 Sr/ 86 Sr isotope ratios, Ca/Mg and Ca/Sr concentration ratios were constructed considering the following end-members: (i) bulk soil, (ii) throughfall, and (iii) soil water.The presented mixing trends were modelled for CER using the elemental and isotope compositions for the above potential end-members given in Tab.S1 and Tab.S2.
For example, the corresponding two-end member mixing trends for δ 26 Mg and 87 Sr/ 86 Sr were calculated based on an isotope mass balance equation below for the mixing of (i) bulk soil (SOIL) and (ii) throughfall (TF) derived Mg sources in a theoretical mixture: where δMg MIX represents the Mg isotope composition of the 'mixed' (soil + throughfall) reservoir, and the parameters F SOIL , δ SOIL and C SOIL represent, respectively, the fraction (F, from 0 to 1), the isotope signature and the Mg concentration of the bulk soil (SOIL) reservoir.An identical coding but with a subscript 'TF' was used for the Mg sourced from local throughfall.
A similar approach was used for the calculation of Sr isotope composition ( 87 Sr/ 86 Sr) of the mixture, labelled below as Sr|Sr MIX , based on the following equation: Finally, the Ca/Mg ratio of the mixture was calculated using simple elemental mass balance equations (see below), calculated separately for Ca and Sr sources or end-members: where F SOIL and C SOIL represent, respectively, the fraction (F, from 0 to 1) and the concentration of Ca derived from the bulk soil (SOIL) reservoir; the fraction of Ca derived from precipitation can be calculated as F TF = 1 − F SOIL .
The concentration of Sr in the mixture (Sr MIX ) was calculated using an identical approach, and the corresponding Ca/Sr ratio of the mixture was then calculated by dividing Ca MIX by Sr MIX .
The concentration of Mg in the mixture (Mg MIX ) was calculated using on identical approach, and the corresponding Ca/Mg ratio of the mixture was then calculated by dividing Ca MIX by Mg MIX .The mass balance of the weathering is where A is an accumulation of the elements in the rock, D is the atmospheric deposition of the elements, R is the rate of dissolution of the elements from the rock, and Q is the removal of elements by percolating water.Units are mol.m −3 s −1 .
The content of albite (An10) in sandstone and claystone is similar, 5.8 and 7.1 wt.% (Table 3 in the main text).Unknown is the difference in the reacting surface area.This will certainly be much larger in claystone than in sandstone.If we assume a 50-fold difference, the rate of Ca release from claystone is 4 10 -6 to 2.0 10 -5 g m −3 s −1 .
The calculation of the dissolution rate of Mg is even more uncertain.Let us assume that Mg contribution from chlorite is negligible.The main Mg sources are then illite and biotite.In the X-ray analysis, illite is included in muscovite.Biotite was identified in thin sections.The dissolution rate constant for illite is from 1.0 10 -14 to 1.0 10 -13 mol m −2 s −1 .The dissolution rate constant for biotite is from 2.0 10 -13 to 1.0 10 -9 mol m −2 s −1 .The higher rate constants for dissolution of biotite indicate that the main source of Mg at the weathering front is biotite rather than illite.The average MgO content in sandstone is 0.21 wt.%.The average MgO content in claystone is 1.6 wt.% (Table 1 in the main text).Assuming that the active surface of biotite is the same in sandstone and in claystone, the rate of dissolution of Mg from claystone will be faster by a factor of 7.6.

Fig. 1
Fig. 1 Study site location (a), and sampling scheme in the Cervik (CER) catchment (b).Previous Mg, Ca and Sr isotope studies in headwater catchments in panel (a) are marked by open circles (Novak Comparisons of Mg, Ca, and Sr isotope compositions of individual sample types were based on a one-way analysis of variance and Tukey's multiple comparisons method (Tukey 1953).The estimated slopes and p-values of decreasing fluxes of dissolved Mg and Ca were based on a linear model with autoregressive errors over time.ArcGIS, version 10.3, was used to average slopes in the catchment.The slope raster was generated from a digital elevation model (DEM; spatial resolution of 1 m).A model byCapo et al. (1998) was used to estimate the contributions of possible major mixing endmembers to base cations in runoff (Appendix II).

Fig. 4
Fig. 4 Long-term input-output fluxes of Mg, Ca and Sr at CER

Fig. 5
Fig. 5 Comparison of magnitude and long-term trends in input Mg, Ca and Sr fluxes in four headwater catchments situated in the Czech Republic differing in their industrial pollution levels (polluted north and unpolluted south, see Fig. 1a for locations).All sites are underlain by base poor crystalline bedrock Mg MIX = F SOIL SOIL C SOIL + F TF TF C TF ∕ F SOIL C SOIL + F TF C TF Sr|Sr MIX = F SOIL Sr|Sr SOIL C SOIL + F TF Sr|Sr TF C TF ∕ F SOIL C SOIL + F TF C TF Ca MIX = F SOIL C SOIL + F TF C TF ;Mg MIX = F SOIL C SOIL + F TF C TF Appendix 3.Relative rates of dissolution of CER bedrock mineralsRates of dissolution of minerals in the studied sandstones and claystones depend on rate constants (k d ), reaction surfaces of the minerals, and the rate of water percolation through the rock complexes.

Table 1
Study site characteristics * Ca and  87Sr/ 86 Sr ratios of precipitation and throughfall (p > 0.05).The δ 26 Mg values of spruce roots and xylem were plotted

Table 2
Chemical compositions (wt.% ± SD) of the sandstone, claystone and conglomerate bedrock at the study site