Quarter century of mercury litterfall at a coniferous forest responding to climate change, Central Europe.

This work evaluated the 25-year-long trends (1994-2018) in mercury (Hg) concentrations and fluxes in spruce litterfall at a forest research plot Nacetín (NAC) recovering from acidic deposition in the Ore Mountains, Czech Republic. The mean litterfall Hg deposition averaged 51 ± 18 µg m-2 year-1, which has been the highest litterfall Hg deposition reported up to date on the European continent. In contrast, the wet deposition (2017-2019) was an order of magnitude lower averaging at 2.5 ± 1.5 µg m-2 year-1. All the spruce litterfall components bark, twigs, needles, cones, and a mixture of unidentified fragments had elevated mean Hg concentrations relative to background sites averaging 256 ± 77, 234 ± 62, 119 ± 23, 95 ± 14, and 44 ± 15 µg kg-1, respectively. Elevated litterfall Hg deposition and concentrations were attributed to the nearby local Hg emission source-coal-fired power plants. Temporally, the decrease of Czech Hg emissions since the 1990s was reflected by the decreasing trend of Hg concentrations in litterfall bark, cones, and twigs, while in needles and other material, Hg increased but insignificantly. Total litterfall ratios of Hg/C, Hg/N, and Hg/S were lower than those in soil O horizons averaging at 0.23 ± 0.04, 9.5 ± 2.0, and 170 ± 37 µg g-1, respectively. Since the beginning of monitoring, total litterfall Hg/C exhibited no trend, Hg/N decreased, and Hg/S increased. The litterfall biomass deposition averaging at 469 ± 176 g m-2 year-1 increased through time resulting in an increased Hg litterfall deposition at NAC by 1.1 µg m-2 year-1 despite the decreases in Czech Hg emissions. Peaks of annual litterfall Hg deposition up to 96 µg m-2 year-1 at NAC during the 25 years of monitoring resulted from weather extremes such as rime-snow accumulation, wind gusts, droughts, and insect infestation, which all significantly affected the annual biomass deposition. Based on our observations, further increases in biomass and litterfall Hg deposition rates can be expected due to the onset of bark beetle infestation and the increasing number of droughts caused by climate change.

How do root fungi of Alnus nepalensis and Schima wallichii recover during succession of abandoned land?

Alnus nepalensis and Schima wallichii are native tree species accompanying succession in abandoned agricultural land in the middle mountainous region of central Nepal. To understand how root fungi recover during spontaneous succession, we analyzed the diversity and composition of arbuscular mycorrhizal (AM), ectomycorrhizal (ECM), and total fungi in tree fine roots from three land use types, short-term abandoned land (SA), long-term abandoned land (LA), and regenerated forest (RF) as a reference. Additionally, ECM morphotypes were examined. The results showed different speeds of succession in the studied fungal groups. While the change in the AM fungal community appears to be rapid and LA resembles the composition of RF, the total fungi in the abandoned land types are similar to each other but differed significantly from RF. Interestingly, the relative abundance of Archaeosporaceae followed a trend differing between the tree species (SA < LA in A. nepalensis, but SA > LA in S. wallichii). Unlike AM and total fungi, there was no significant difference in the ECM community of A. nepalensis between land use types, probably due to their low species diversity (9 ECM morphotypes, 31 ECM operational taxonomic units). However, Cortinarius sp. was significantly more abundant in RF than in the other land use types, whereas Alnicola, Tomentella, and Russula preferred young stages. Our results suggest that for both studied tree species the AM fungal succession could reach the stage of regenerated forest relatively fast. In the case of total fungi, because of hyperdiversity and composed of species specialized to a variety of environments and substrates, the transition was expected to be delayed in abandoned land where the vegetation was still developing and the ecosystem was not as complex as that found in mature forests.

Toward spatially polarized human pressure? A dynamic factor analysis of ecological stability and the role of territorial gradients in Czech Republic.

In light of global change, research on ecosystem dynamics and the related environmental policies are increasingly required to face with the inherent polarization in areas with low and high human pressure. Differential levels of human pressure are hypothesized to reflect development paths toward ecological stability of local systems vis à vis socioeconomic resilience. To delineate the latent nexus between socioeconomic development paths and ecological stability of local systems, we proposed a multidimensional, diachronic analysis of 28 indicators of territorial disparities, and ecological stability in 206 homogeneous administrative units of Czech Republic over almost 30 years (1990-2018). Mixing time-invariant factors with time-varying socio-environmental attributes, a dynamic factor analysis investigated the latent relationship between ecosystem functions, environmental pressures, and the background socioeconomic characteristics of the selected spatial units. We identified four geographical gradients in Czech Republic (namely elevation, economic agglomeration, demographic structure, and soil imperviousness) at the base of territorial divides associated with the increased polarization in areas with low and high human pressure. The role of urbanization, agriculture, and loss of natural habitats reflective of rising human pressure was illustrated along the selected gradients. Finally, policy implications of the (changing) geography of ecological disturbances and local development paths in Czech Republic were briefly discussed.

Delineating the Intrinsic, Long-Term Path of Land Degradation: A Spatially Explicit Transition Matrix for Italy, 1960-2010.

Vulnerability to land degradation in southern Europe has increased substantially in the last decades because of climate and land-use change, soil deterioration, and rising human pressure. The present work focuses on a quantitative evaluation of changes over time in the level of vulnerability to land degradation of a Mediterranean country (Italy) using a composite indicator, the environmentally sensitive area index (ESAI), which is the final outcome of a complex model conceived to assess land vulnerability on the basis of climate, soil, vegetation, and human pressure. Considering four different levels of vulnerability to land degradation (not affected, potentially affected, fragile, and critical), the main trajectories of this index were highlighted in a long-time perspective (1960-2010), discriminating dynamics over two sub-periods (1960-1990 and 1990-2010). The empirical results at a very detailed spatial scale (1 km2 grid) reflect spatial consolidation of degradation hot-spots over time. However, aggregated trajectories of change indicate an overall improvement in the environmental conditions between 1990 and 2010 compared with what is observed during the first period (1960-1990). Worse environmental conditions concerned southern Italian regions with a dry climate and poor soil conditions in the first time interval, large parts of northern Italy, traditionally recognized as a wet and affluent agricultural region, experienced increasing levels of land vulnerability in the second time interval. Being classified as an unaffected region according with the Italian national action plan (NAP), the expansion of (originally sparse) degradation hot-spots in northern Italy, reflective of an overall increase in critical areas, suggests a substantial re-thinking of the Italian NAP. This may lead to a redesign of individual regional action plans (RAPs) implementing place-specific approaches and comprehensive measures to be adopted to mitigate land degradation.

Investigating metropolitan change through mathematical morphology and a dynamic factor analysis of structural and functional land-use indicators.

We presented an operational rationale grounded on complex system thinking to quantify structural and functional landscape transformations along three stages representative of post-war metropolitan development in Rome, Italy (urbanisation with population/settlement densification, 1949-1974; suburbanisation with medium-density settlement expansion, 1974-1999; counter-urbanisation with settlement sprawl, 1999-2016). A mathematical morphology approach assessing the geometric form of land patches and a multi-way factor analysis (MFA) of landscape metrics were used to investigate the joint evolution of urban form and land-use functions over time. The empirical results of the MFA delineated the multivariate relationship between nine land-use classes (with distinctive socioeconomic functions) and seven morphological types (reflecting different landscape structures) according to four observation times (1949, 1974, 1999, 2016). Taken as an intrinsic attribute of complex landscape systems experiencing intense transformations, an estimation of the 'rapidity-of-change' in the form-functions relationship at a given development stage was derived from MFA outcomes separately for urbanisation, suburbanisation, and counter-urbanisation. A simplified form-functions relationship, reflecting the spatial polarisation in compact settlements and rural (low-density) landscapes, was observed with compact urbanisation. By stimulating urban sprawl into fringe farmland, suburbanisation resulted in patchy and heterogeneous rural landscapes. Counter-urbanization was associated with the fragmentation of built-up settlements leading to a chaotic mosaic of land structures that mixes urban and rural traits. Rapidity-of-change in form-function relationships was greater during suburbanisation than urbanisation and counter-urbanisation. It reflects the intrinsic pressure of economic growth in contemporary cities.

Primary and secondary host selection by Ips typographus depends on Norway spruce crown characteristics and phenolic-based defenses.

Climate change is expected to intensify bark beetle population outbreaks in forests globally, affecting biodiversity and trajectories of change. Aspects of individual tree resistance remain poorly quantified, particularly with regard to the role of phenolic compounds, hindering robust predictions of forest response to future conditions. In 2003, we conducted a mechanical wounding experiment in a Norway spruce forest that coincided with an outbreak of the bark beetle, Ips typographus. We collected phloem samples from 97 trees and monitored tree survival for 5 months. Using high-performance liquid chromatography, we quantified induced changes in the concentrations of phenolics. Classification and regression tools were used to evaluate relationships between phenolic production and bark beetle resistance, in the context of other survival factors. The proximity of beetle source populations was a principal determinant of survival. Proxy measures of tree vigor, such as crown defoliation, mediated tree resistance. Controlling for these factors, synthesis of catechin was found to exponentially increase tree survival probability. However, even resistant trees were susceptible in late season due to high insect population growth. Our results show that incorporating trait-mediated effects improves predictions of survival. Using an integrated analytical approach, we demonstrate that phenolics play a direct role in tree defense to herbivory.

Toward a spatially explicit analysis of land vulnerability to degradation: a country-level approach supporting policy strategies.

Vulnerability to land degradation in Mediterranean Europe increased substantially in the last decades because of the latent interplay of climate and land-use change, progressive soil deterioration, and rising human pressure. The present study provides a quantitative evaluation of the intrinsic change over time in the level of vulnerability to land degradation over a representative Mediterranean area (Italy) using a normative indicator, the percentage of land classified as 'critical' in total area. This indicator derives from a spatially explicit elaboration of the ESA (Environmental Sensitive Area) Index (ESAI), a standard methodology of land classification considering different levels of vulnerability to degradation at a particularly refined spatial scale (1 km2). This indicator was calculated over a relatively long time interval (1960-2010) and aggregated at the geographical scale of administrative regions in Italy, a relevant domain in the implementation of the National Action Plan (NAP) to combat desertification and the adoption of individual Regional Action Plans (RAP). A significant - but spatially heterogeneous - increase in 'critical' land was observed in Italy, leading to distinctive dynamics in northern/central regions and southern regions. Climate aridity and anthropogenic pressure leveraged the sudden vulnerability in some marginal land of Northern Italy - a region classified as unexposed to desertification risk - paralleling the levels observed in some districts of Southern Italy, an 'affected' region to desertification risk. These results suggest a re-thinking of mitigation policies proposed in the Italian NAP and a redesign of the RAPs toward place-specific adaptation measures, especially in the 'less exposed' Northern Italian region.

Desertification risk fuels spatial polarization in 'affected' and 'unaffected' landscapes in Italy.

Southern Europe is a hotspot for desertification risk because of the intimate impact of soil deterioration, landscape transformations, rising human pressure, and climate change. In this context, large-scale empirical analyses linking landscape fragmentation with desertification risk assume that increasing levels of land vulnerability to degradation are associated with significant changes in landscape structure. Using a traditional approach of landscape ecology, this study evaluates the spatial structure of a simulated landscape based on different levels of vulnerability to land degradation using 15 metrics calculated at three time points (early-1960s, early-1990s, early-2010s) in Italy. While the (average) level of land vulnerability increased over time almost in all Italian regions, vulnerable landscapes demonstrated to be increasingly fragmented, as far as the number of homogeneous patches and mean patch size are concerned. The spatial balance in affected and unaffected areas-typically observed in the 1960s-was progressively replaced with an intrinsically disordered landscape, and this process was more intense in regions exposed to higher (and increasing) levels of land degradation. The spread of larger land patches exposed to intrinsic degradation brings to important consequences since (1) the rising number of hotspots may increase the probability of local-scale degradation processes, and (2) the buffering effect of neighbouring (unaffected) land can be less effective on bigger hotspots, promoting a downward spiral toward desertification.

Management After Windstorm Affects the Composition of Ectomycorrhizal Symbionts of Regenerating Trees but Not Their Mycorrhizal Networks.

Due to ongoing climate change, forests are expected to face significant disturbances more frequently than in the past. Appropriate management is intended to facilitate forest regeneration. Because European temperate forests mostly consist of trees associated with ectomycorrhizal (ECM) fungi, understanding their role in these disturbances is important to develop strategies to minimize their consequences and effectively restore forests. Our aim was to determine how traditional (EXT) and nonintervention (NEX) management in originally Norway spruce (Picea abies) forests with an admixture of European larch (Larix decidua) affect ECM fungal communities and the potential to interconnect different tree species via ECM networks 15 years after a windstorm. Ten plots in NEX and 10 plots in EXT with the co-occurrences of Norway spruce, European larch, and silver birch (Betula pendula) were selected, and a total of 57 ECM taxa were identified using ITS sequencing from ECM root tips. In both treatments, five ECM species associated with all the studied tree species dominated, with a total abundance of approximately 50% in the examined root samples. Because there were no significant differences between treatments in the number of ECM species associated with different tree species combinations in individual plots, we concluded that the management type did not have a significant effect on networking. However, management significantly affected the compositions of ECM symbionts of Norway spruce and European larch but not those of silver birch. Although this result is explained by the occurrence of seedlings and ECM propagules that were present in the original forest, the consequences are difficult to assess without knowledge of the ecology of different ECM symbionts.

Decreasing litterfall mercury deposition in central European coniferous forests and effects of bark beetle infestation.

We evaluated a 14-year trend (2003-2017) in mercury (Hg) concentrations and fluxes in six litterfall categories (needles, bark, twigs, cones, lichen, and a mixture of unidentified fragments) at six research plots situated in two central European unmanaged mountain forest stands, dominated by mature Norway spruce. One of the stands (catchment of Plesné Lake, PL) was infested by bark beetle and all mature spruces died at three of four research plots during the study. One PL plot and two plots in the second stand (catchment of Certovo Lake, CT) were intact and used as a control. At the intact plots, the litterfall Hg deposition averaged 45 and 32 µg/m2/year in the PL and CT catchments, respectively, while bulk precipitation Hg deposition was an order of magnitude lower (2.6 µg/m2/year). In the individual litter categories, Hg concentrations averaged 223, 195, 126, 81, 79 and 30 µg/kg in lichen, unidentified fragments, bark, twigs, needles and cones, respectively. As a result of decreasing European Hg emissions, Hg concentrations in most litter categories decreased from 2003 to 2017. Consequently, the litter-associated Hg flux to the forest floor decreased from 66 to 23 µg/m2/year during 2003-2017. The litterfall mass fluxes exhibited no trends at the intact plots. In contrast, the litter-associated Hg flux increased 5-fold after tree dieback due to elevated litterfall, averaging 218 µg/m2/year Hg at the PL infested plots during 2004-2009. The relative contribution of individual litter categories to the total Hg flux shifted from needle to bark and twig dominance. Starting in 2010, Hg flux decreased to pre-disturbance levels for the following decade. The tree mortality in the PL catchment provided a unique opportunity to evaluate changes in litter-associated Hg fluxes to the forest floor during and after natural tree dieback.

Management regime is the most important factor influencing ectomycorrhizal species community in Norway spruce forests after windthrow.

Ectomycorrhizal (ECM) fungi, as symbionts of many tree species in temperate forests, are thought to play an important role in forest regeneration processes after large disturbances. Their reaction to different disturbance and management regimes was studied in spruce forests (Lariceto-Piceetum) 10 years after a severe windthrow in the Tatra National Park (Slovak Republic). ECM community structure was compared between different "management types″-cleared area (EXT), area affected by wildfire (FIRE), uncleared area left for natural development (NEX), and mature forest as a control (REF). Based on Illumina sequencing of soil samples, we determined that the percentage of sequences assigned to ECM fungi decreased with increasing disturbance and management intensity (REF → NEX → EXT → FIRE). Similarly, the total number of ECM species per each of ten sampling points per plot (100 ha) differed between managed (EXT-11 species, FIRE-9) and unmanaged (NEX-16, REF-14) treatments. On the other hand, the percentage of sequences belonging to ericoid mycorrhizal fungi increased. Management type significantly influenced the composition of the ECM community, while vegetation and soil characteristics explained less data variation. The ECM species assemblage of the unmanaged site (NEX) was the most similar to the mature forest, while that of the burnt site was the most different. Thelephora terrestris dominated in all treatments affected by windthrow, accompanied by Tylospora fibrillosa (NEX) and Tylospora asterophora (EXT and FIRE). Management regime was also the most important factor affecting ECM species composition on the roots of spruce seedlings assessed by Sanger sequencing.

Discrimination between acute and chronic decline of Central European forests using map algebra of the growth condition and forest biomass fuzzy sets: A case study.

Forest decline is either caused by damage or else by vulnerability due to unfavourable growth conditions or due to unnatural silvicultural systems. Here, we assess forest decline in the Czech Republic (Central Europe) using fuzzy functions, fuzzy sets and fuzzy rating of ecosystem properties over a 1×1km grid. The model was divided into fuzzy functions of the abiotic predictors of growth conditions (Fpred including temperature, precipitation, acid deposition, soil data and relative site insolation) and forest biomass receptors (Frec including remote sensing data, density and volume of aboveground biomass, and surface humus chemical data). Fuzzy functions were designed at the limits of unfavourable, undetermined or favourable effects on the forest ecosystem health status. Fuzzy sets were distinguished through similarity in a particular membership of the properties at the limits of the forest status margins. Fuzzy rating was obtained from the least difference of Fpred-Frec. Unfavourable Fpred within unfavourable Frec indicated chronic damage, favourable Fpred within unfavourable Frec indicated acute damage, and unfavourable Fpred within favourable Frec indicated vulnerability. The model in the 1×1km grid was validated through spatial intersection with a point field of uniform forest stands. Favourable status was characterised by soil base saturation (BS)>50%, BCC/Al>1, Corg>1%, MgO>6g/kg, and nitrogen deposition<1200mol(H+)/ha·year. Vulnerable forests had BShumus 46-60%, BCC/Al 9-20 and NDVI≈0.42. Chronic forest damage occurs in areas with low temperatures, high nitrogen deposition, and low soil BS and Corg levels. In the Czech Republic, 10% of forests were considered non-damaged and 77% vulnerable, with damage considered acute in 7% of forests and chronic in 5%. The fuzzy model used suggests that improvement in forest health will depend on decreasing environmental load and restoration concordance between growth conditions and tree species composition.

Pollution control enhanced spruce growth in the "Black Triangle" near the Czech-Polish border.

Norway spruce (Picea abies (L.) Karst.) stands in certain areas of Central Europe have experienced substantial dieback since the 1970s. Understanding the reasons for this decline and reexamining the response of forests to acid deposition reduction remains challenging because of a lack of long and well-replicated tree-ring width chronologies. Here, spruce from a subalpine area heavily affected by acid deposition (from both sulfur and nitrogen compounds) is evaluated. Tree-ring width measurements from 98 trees between 1000 and 1350m above sea level (a.s.l.) reflected significant May-July temperature signals. Since the 1970s, acid deposition has reduced the growth-climate relationship. Efficient pollution control together with a warmer but not drier climate most likely caused the increased growth of spruce stands in this region, the so-called "Black Triangle," in the 1990s.

The influence of climate change on stomatal ozone flux to a mountain Norway spruce forest.

Daily stomatal ozone flux to a mountain Norway spruce forest stand at the Bily Kriz experimental site in the Beskydy Mts. (Czech Republic) was modelled using a multiplicative model during the 2009 growing season. The multiplicative model was run with meteorological data for the growing season 2009 and ALADIN-CLIMATE/CZ model data for the 2030 growing season. The exceedance of the flux-based critical level of O(3) (Phytotoxic Ozone Dose) might be lower for Norway spruce at the Bily Kriz experimental site in a future climate (around 2030), due to increased stomatal closure induced by climate change, even when taking into account increased tropospheric background O(3) concentration. In contrast, exceedance of the concentration-based critical level (AOT40) of O(3) will increase with the projected increase in background O(3) concentration. Ozone concentration and stomatal flux of ozone significantly decreased NEP under both present and future climatic conditions, especially under high intensities of solar radiation.

Behaviour of arsenic in forested catchments following a high-pollution period.

Due to high availability of adsorption sites, forested catchments could be net sinks for pollutant arsenic both during the period of increasing and decreasing pollution. We tested this hypothesis along a north-south pollution gradient in spruce die-back affected areas of Central Europe. For two water years (2007-2008), we monitored As fluxes via spruce-canopy throughfall, open-area precipitation, and runoff in four headwater catchments (Czech Republic). Since 1980, atmospheric As inputs decreased 26 times in the north, and 13 times in the south. Arsenic export by runoff was similar to atmospheric inputs at three sites, resulting in a near-zero As mass balance. One site exhibited a net export of As (2.2 g ha(-1) yr(-1)). In contrast, the preceding period (1995-2006) showed much higher As fluxes, and higher As export. Czech catchments do not serve as net sinks of atmospheric As. A considerable proportion of old industrial arsenic is flushed out of the soil.

Ozone flux over a Norway spruce forest and correlation with net ecosystem production.

Daily ozone deposition flux to a Norway spruce forest in Czech Republic was measured using the gradient method in July and August 2008. Results were in good agreement with a deposition flux model. The mean daily stomatal uptake of ozone was around 47% of total deposition. Average deposition velocity was 0.39 cm s(-1) and 0.36 cm s(-1) by the gradient method and the deposition model, respectively. Measured and modelled non-stomatal uptake was around 0.2 cm s(-1). In addition, net ecosystem production (NEP) was measured by using Eddy Covariance and correlations with O3 concentrations at 15 m a.g.l., total deposition and stomatal uptake were tested. Total deposition and stomatal uptake of ozone significantly decreased NEP, especially by high intensities of solar radiation.

Culturable bacterial populations associated with ectomycorrhizae of Norway spruce stands with different degrees of decline in the Czech Republic.

The aim of this study was to determine which species of culturable bacteria are associated with ectomycorrhizae (ECM) of Norway spruce (Picea abies (L.) Karst) in the Sudety Mountains, exposed for years to atmospheric pollutants, acid rain, and climatic stress, and to identify particular species that have adapted to those conditions. Biolog identification was performed on bacterial species from ECM of adult spruce trees and seedlings of stands with low, intermediate, and high forest decline. Bacterial diversity in ECM associated with adult spruce trees, seedlings, and seedlings grown on monoliths was calculated; although the expected values appeared to vary widely, no significant differences among sites were observed. Dendrograms based on the identified bacterial species showed that stands with low forest decline clustered separately from the others. Principal component analysis of the normalized data for ECM-associated species showed a clear separation between stands with high forest decline and stands with low forest decline for seedlings and a less evident separation for adult spruce trees. In conclusion, shifts in ECM-associated culturable bacterial populations seem to be associated with forest decline in Norway spruce stands. Some bacterial species were preferentially associated with mycorrhizal roots depending on the degree of forest decline; this was more evident in seedlings where the species Burkholderia cepacia and Pseudomonas fluorescens were associated with, respectively, ECM of the most damaged stands and those with low forest decline.

Belowground ectomycorrhizal communities in three Norway spruce stands with different degrees of decline in the Czech Republic.

We investigated the ectomycorrhizal communities on the roots of adult trees and seedlings associated with three Norway spruce stands in the Czech Republic using morphological and molecular tools. The stands had different degrees of forest decline due to air pollution. The aims of the study were to obtain information about the belowground ectomycorrhizal community in a heavily damaged spruce forest and to assess whether missing ectomycorrhizal fungal partners could be one of the reasons for the observed lack of regeneration. The ectomycorrhizal species richness on the roots of adult trees was significantly lower in the heavily damaged site Mumlavska hora than in the other two, but less drastically so than that found in a fruitbody survey. The roots of adult trees and seedlings were fully mycorrhizal at this site although they were less species-rich. The most abundant ectomycorrhizal species on the root system of adult trees in all three forest stands was Tylospora fibrillosa, a member of the athelioid clade. It made up over 60% of root tips in Mumlavska hora and its proportion was at least twice that in the other two sites. This species was also an efficient colonizer of roots from seedlings, in particular, in the most damaged site. The different soil properties in this site may have caused the observed differences in the ectomycorrhizal species richness and composition. For example, cation exchange capacity and soil base saturation were lower and the soil more often saturated. However, the number of living trees and their defoliation status may well directly impact the ectomycorrhizal species composition by presumably affecting the amount of carbon delivered to the symbiotic fungal partners. Athelioids and thelephoroids are an important component of the belowground ectomycorrhizal community in most temperate and boreal forests, but the role they play might even be more crucial in stressed forest ecosystems. Based on our results, we suggest that factors other than missing ectomycorrhizal inoculum constrain natural regeneration in the heavily damaged site Mumlavska hora.

The effect of reduced atmospheric deposition on soil and soil solution chemistry at a site subjected to long-term acidification, Nacetín, Czech Republic.

During the 1990s the emissions of SO(2) fell dramatically by about 90% in the Czech Republic; the measured throughfall deposition of sulphur to a spruce forest at Nacetín in the Ore Mts. decreased from almost 50 kg ha(-1) in 1994 to 15 kg ha(-1) in 2005. The throughfall flux of Ca decreased from 17 kg ha(-1) in 1994 to 9 kg ha(-1) in 2005; no change was observed for Mg. The deposition of nitrogen ranged between 15 and 30 kg ha(-1) with no statistically significant trend in the period 1994-2005. The desorption of previously stored sulphur and the decrease of Ca deposition are the main factors controlling the recovery of soil solution. The pH of the soil solution at a depth of 30 cm remains unchanged, and the Al concentration decreased from 320 micromol l(-1) in 1997 to 140 micromol l(-1) in 2005. The enhanced leaching of base cations relative to no acidified conditions has continued, although the Ca concentration decreased from 110 microeq l(-1) in 1997 to 25 microeq l(-1) in 2005 in the mineral soil solution at 30 cm depth. This dramatic change was not observed for Mg concentration in soil solution, because its deposition remained stable during the observed period. Similar patterns were observed in the deeper soil solution at 90 cm. The reduction in Ca availability resulted in lower uptake by tree assimilatory tissues, measured as concentration in needles. Since 2005, the leaching of nitrate observed in soil solution at 30 cm depth has disappeared. By 2003 a similar situation occurred at 90 cm. Higher incorporation into the trees after 1997 could be an important factor. With respect to the formerly high sulphur deposition and consequently released aluminium, which could have negatively influenced the biotic immobilization driven by microbes and fungi, the recovery may have positively impacted and therefore improved retention in the ecosystem during recent years. The delay in the successful retention of nitrogen in the ecosystem was probably caused by the high mineralization of organic matter after improvement of chemical parameters in the organic horizon (increase in pH and decrease in Al concentration). It seems that high mineralization of stored organic matter after decades of high acidic deposition could be an important factor affecting the high losses of nitrogen in spruce forest ecosystems.