Use of the Dynamic Model `MAGIC' to Predict Recovery Following Implementation of the Gothenburg Protocol

Reduction of sulphur deposition causesrecovery of acidified surface waters. Processes in thecatchment delay recovery. The acidification model MAGICwas applied to the Vikedal and Tovdal rivers in southernNorway. Response in water chemsitry is delayed by 10–20 yr. The delay is due to release of old sulphate atVikedal and cation exchange at Tovdal. Assuming that theGothenburg protocol is fully implemented by the year2010, much of the predicted increase of ANC will occur inthe next 10 yr with a levelling off by about 2040. Ifnitrogen leaching increases in the future, however,recovery of ANC will not be as rapid, nor as complete.Critical load for acidity calculated by steady-statemodels is confirmed by the MAGIC predictions. Futurerequirement for mitigation measures such as liming willdecrease in the future as acid deposition decreases. Bythe year 2046 the liming requirement will be reduced byabout 45% at Vikedal and 65% at Tovdal. One of the mainpurposes of the Norwegian national monitoring programmeis to provide documentation of changes in environmentalquality due to long-range transported air pollutants.Modelling applications such as this clearly show that thedata fill this purpose.     

A Novel Method for Mapping Critical Loads Across a River Network: Application to the River Dart,Southwest England

To date, estimates of freshwater critical loads have beenbased on a single sample site within a given area, in theUK the `most sensitive' surface water in each 10 km gridsquare. The critical loads obtained are thus highlydependent on the sites chosen, and at a relatively coarsespatial resolution. To produce a higher resolutioncritical load assessment, the PEARLS (Prediction ofAcidification and Recovery on a Landscape Scale)procedure has been used to estimate critical loads acrossa large (248 km²), partially acid-sensitivecatchment in Southwest England. PEARLS utilises availablesoils and land-use databases, and sampled streamchemistry data, to derive characteristic runoffcompositions for a set of landscape types. Mixingequations are then used to calculate runoff chemistry,and subsequently critical loads, throughout the streamnetwork. Results show major spatial variability, withcritical loads lowest in streams draining peat-moorlandheadwaters, and generally increasing downstream asagricultural land contributes an increasing proportion ofrunoff. The 5th percentile freshwater critical loadfor the catchment is estimated at 0.29 keq H⁺ ha^-1yr^-1,and critical loads are exceeded for around 40% of totalstream length. The PEARLS methodology provides a novelopportunity to assess the spatial variability infreshwater critical loads, and to provide estimates ofexceedance at whole catchment scale. It has potentialapplication in the assessment of surface watersensitivity to acidification across wider areas in the UKand elsewhere.     

The Link Between the Exceedance of Acidity Critical Loads for Freshwaters,Current Chemical Status and Biological Damage: A Re-Interpretation

Critical loads have for several years been employed bypolicymakers to aid in the development of strategies for aciddeposition abatement. They provide an effects-based approachwhereby an acid deposition flux greater than the critical load(known as critical load exceedance) implies that long-termharmful effects on a selected target organism will occur.Implicit in this approach are two assumptions: first, theexceedance of a critical load will harm the target organism,and second, the severity of biological impact is related to themagnitude of exceedance. However, static models give noindication of when the predicted damage might occur. One suchmodel, the Steady-State Water Chemistry (SSWC) model, employs aseries of empirical relationships to derive the pre-industrial,baseline leaching rate of base cations from measured waterchemistry using the so-called `F-factor'. The SSWC model setsthe critical load relative to pre-industrial base cationleaching (a permanent buffer of acid deposition) and a selectedacid neutralizing capacity (ANC) value which corresponds with aknown likelihood of damage to a biological target organism.Here we interpret the meaning of critical load exceedance as aprediction of steady-state ANC, and explore the relationshipbetween exceedance of the critical load and current chemistry. We demonstrate that a critical loadexceedance with the SSWC model does not necessarily indicatethat the critical chemical threshold (zero ANC) has alreadybeen crossed, and there may be no correlation betweenexceedance and biological status. A reformulation of the SSWCmodel is proposed which provides a direct link between currentdeposition and current chemical conditions, and is thereforemore likely to indicate current biological damage. Thereformulation illustrates the discrepancy between currentchemical status and that predicted by the SSWC model atsteady-state, which is a function of the `F-factor'.     

Using Multiple Chemical Criteria for Critical Loads of Acidity in Maritime Regions

Critical loads have become a well-establishedpart of the work programme of the UnitedNations Economic Commission for Europe (UNECE) Conventionon Long-Range Transboundary Air Pollution (CLRTAP). Thelinking of ecosystem response to deposition level is thecentral principal of the critical loads approach. Foreach ecosystem, a biological indicator is chosen, asuitable chemical criterion selected and a criticalchemical limit assigned. The Bc:Al ratio is the mostwidely used chemical criterion for setting criticalloads. However, critical loads based on this criterionare very sensitive to marine deposition. In regions whichreceive high depositions of marine-derived base cations,such as the Republic of Ireland (ROI) and the UnitedKingdom (UK), critical loads based on the Bc:Al criterionwill inevitably be high. Therefore, it is proposed thatcritical loads are estimated using multiple chemicalcriteria with appropriate critical limits to protect thechosen biological indicators. The range of publishedchemical criteria have been applied to the ROI and theUK. The chemical criterion corresponding to the mostsensitive critical load have been mapped and thecontribution of each to the final maps investigated. Thesimulations indicate that the most sensitive criteria forsetting critical loads are based on specifying critical Hor Al concentrations. However, the choice of critical limits andmodel parameters will ultimately effect the criticalloads. Therefore, it is important that appropriate criticallimits are chosen to protect the biological indicator andreceptor ecosystem from long-term damage.     

Critical Loads and Dynamic Modelling to Assess European Areas at Risk of Acidification and Eutrophication

European critical loads and novel dynamic modelling data have been compiled under the LRTAP Convention by the Coordination Centre for Effects. In 2000 9.8% of the pan-European and 20.8% of the EU25 ecosystem area were at risk of acidification. For eutrophication (nutrient N) the areas at risk were 30.1 and 71.2%, respectively. Dynamic modelling results reveal that 95% of the area at risk of acidification could recover by 2030 provided acid deposition is reduced according to present legislation. Insight into the timing of effects of exceedances of critical loads for nutrient N necessitates the further development of dynamic models.     

On the Calculation and Interpretation of Target Load Functions

In this study critical load functions and target load functions of nitrogen and sulphur deposition with respect to acidity and minimum base cation to aluminium ratio were calculated with the SAFE model using three different averaging strategies: (1) averaging based on current forest generation, (2) averaging based on next generation and (3) averaging based on the entire simulation period. From the results it is evident that although target load calculation and indeed critical load calculation is straight forward, there is a problem in translating a predicted recovery according to the target load calculation back to a site-specific condition. We conclude that a policy strategy for emission reductions that ensures recovery, according to calculated target load functions, is likely to be beneficial from an ecosystem point of view. However, such a strategy may not be sufficient to achieve actual non-violation of the chemical criteria throughout the seasonal or rotational variations. To address this issue we propose a method for calculating dynamic critical loads which ensures that the chosen criteria is not violated.     

A Conceptual Model of Spatially Heterogeneous Nitrogen Leaching from a Welsh Moorland Catchment

Soil- and stream-water data from the Plynlimon research area, mid-Wales, have been used to develop a conceptual model of spatial variations in nitrogen (N) leaching within moorland catchments. Extensive peats, in both hilltop and valley locations, are considered near-complete sinks for inorganic N, but leach the most dissolved organic nitrogen (DON). Peaty mineral soils on hillslopes also retain inorganic N within upper organic horizons, but a proportion percolates into mineral horizons as nitrate (NO ₃ ^– ), either through incomplete immobilisation in the organic layer, or in water bypassing the organic soil matrix via macropores. This NO ₃ ^– reaches the stream where mineral soilwaters discharge (via matrix throughflow or pipeflow) directly to the drainage network, or via small N-enriched flush wetlands. NO ₃ ^– in hillslope waters discharging into larger valley wetlands will be removed before reaching the stream. A concept of catchment nitrate leaching zones is proposed, whereby most stream NO ₃ ^– derives from localised areas of mineral soil hillslope draining directly to the stream; the extent of these zones within a catchment may thus determine its overall susceptibility to elevated surface water NO ₃ ^– concentrations.     

Uncertainties in Large Scale Assessments of Critical Load Exceedances

In the 1999 Gothenburg protocol to the UN/ECE LRTAP Convention andin the negotiations for an EU acidification strategy the area withexceedances of critical loads has been the preferred measure forenvironmental impacts. The aim of this study has been to assessthe influence of the uncertainty and spatial variation of both thecritical loads and deposition values on the calculated area withexceedances of critical loads. This has been done on a nationalscale for Denmark and on the European scale based on the dataincluded in the RAINS model. It is demonstrated that includinguncertainty and spatial variation in exceedance calculations, ingeneral gives larger exceeded area for the critical load ofacidity, CL(A). The picture for the critical load of nutrientnitrogen, CL_nut(N), is more mixed because of the higherproportion of exceeded areas. A further point of interest is thepossibility of validating relationships between critical loadexceedances and observable damage based on large scale monitoringand model data. It is demonstrated that it will probably not bepossible to use exceedance calculations on European scale as basisfor validation exercises, linking exceedances to observable damage.     

Application of Static Critical Load Models for Acidity to High Mountain Lakes in Europe

Critical load models for acidityprovide a measure of the sensitivity of surfacewaters to acid deposition, and can be used todetermine critical load exceedance and potentiallong-term harmful effects. Three static models,the Steady-State Water Chemistry model, diatommodel and First-order Acidity Balance model, arehere applied to 11 high mountain lakes in Norway,Scotland, the Alps, the Pyrenees and the Tatras.Between five and seven of the lakes show criticalload exceedance, depending on the model used.Nitrogen as well as sulphur deposition isimportant in causing exceedance. Since soil andvegetation cover are generally sparse, geologyand lake retention time appear to be key factorsin the determination of critical load. Retentionof nitrogen is observed, but it is unclearwhether this occurs within the lake or theterrestrial part of the catchment.     

Invasives, Introductions and Acidification: The Dynamics of a Stressed River Fish Community

We describe the development of the fish community in the acidified and limed river Litleåna in southern Norway, and describe how chemical restoration, compensatory introductions of exotics, and accidental invasion of exotics interact to influence the population of the naturally occurring brown trout (Salmo trutta). The river Litleåna is a tributary to the river Kvina in Vest-Agder County, southern Norway. During the years 1996–2004, annual mean pH was 4.9–5.0 and 6.1–6.4 above and below the liming facility, which was installed in 1994. Originally, brown trout was the only fish species in the river, but brook trout (Salvelinus fontinalis) have been intentionally introduced, whereas European minnow (Phoxinus phoxinus) was introduced by accident. Fish densities were recorded by means of electrofishing annually over the ten year period 1995–2004. Although close to extinction before liming was initiated, brown trout fry densities increased from 1995 to 1999, with subsequent varying densities. There has been a simultaneous major increase in the occurrence and density of European minnow since 1997. Our results show that both brown trout and European minnow increase after liming. Minnow densities are negatively affected by low pH episodes in the river. The growth rates of brown trout fry are negatively correlated to minnow densities, indicating competition between the species. Brook trout densities have decreased since liming started, and during the brown trout recovery.     

总氮总磷联合测定方法研究

介绍总氮、总磷联合测定的一种方法。实验中，调整了消解过程的氧化剂用量和酸度，以含有氮和磷的混合标准溶液配制一系列标准样品，经高温、高压消解后，分别进行比色，测定样品中总氮和总磷。该方法与标准方法相比，大大简化了分别测定总氮、总磷的工作量和时间，测定结果令人满意。     

Determination of Nitrogen and Sulfur in Plants and Soils Sampled from the Area Near a Fertilizer Plant

The purpose of this study was to investigate NO_x and SO₂-inducednitrogen and sulfur levels accumulated over vegetation and soil in the surrounds of the fertilizer plant in Kütahya, Turkey. Therefore, leaf samples collected from the seven locations adjacent to the plant were morphologically observed (pollution-related injuries and color deformation) and chemically analyzed. Besides, nitrogen and sulfur contents were determined in the soil samples collected from each location. A significant correlation was achieved between locations and different vegetation samples and soil samples. The highest nitrogen level present in the trees was 3.07 ± 0.105% (3rd location, Robinia pseudoacacia), with the highest sulfur level being 5.79 ± 0.085% (3rd location, Salix alba). The highest nitrogen and sulfur levels in soils were 0.31 ± 0.045 and 0.57 ± 0.052% (3rd location), respectively. Furthermore, the highest nitrogen and sulfur levels were 4.69 and 2.20%, respectively,in agricultural plants in the surrounds of the fertilizer plant. A higher level of pollutants was observed in the samples collectedfrom the locations in the direction of the prevailing wind. Also, differences between locations as well as between plant specimens were statistically evaluated.     

Atmospheric Levels of Nitrogen Oxides at a Greek Oil Refinery Compared with the Urban Measurements in Athens

Oil refining is among the industrial activities that emit considerable amounts of air pollutants into the atmosphere. Nitrogen oxides are important air pollutants that are emitted by oil refineries as products of combustion processes. The ambient air concentrations of nitrogen oxide (NO) and nitrogen dioxide (NO₂) were monitored continuously at a site close to an oil refinery, near the city of Corinth in Greece, during autumn 1997 together with the main meteorological parameters. The contribution of the oil refinery to the measured atmospheric levels of nitrogen oxides was estimated. The ambient air concentration of nitrogen oxides in the area surrounding the oil refinery were generally lower than the ambient air concentrations in the urban area of Athens in Greece, and the NO₂ levels were always below the existing air quality standards. The influence of the refinery emitted NO_x in the photochemical production of ozone seems to be more important in terms of human and vegetation exposure given the high ozone backgrounds measured in the area.     

Reconstruction of Long-Term Changes in Lake water Chemistry,Zooplankton and Benthos of a Small,Acidified High-Mountain Lake: Magic Modelling and Palaeolimnogical Analysis

Starolesnianske pleso is a small and shallow acid lake in the High Tatra Mountains, situated at 2000 m above sea level, on granitic bedrock, with sparse and thin soil cover. When detailed measurements began in the 1980s Starolesnianske pleso had pH below 5 and only one species of cladoceran zooplankton, the ubiquitous Chydorus sphaericusPalaeolimnological investigations show changes in cladoceran zooplankton and chironomid zoobenthic assemblages since about 1920 and a major change in 1970–1980. The dynamic acidification model MAGIC was used to reconstruct changes in water chemistry over the past 150 years. The results from MAGIC agree well with the sediment record; pH levels gradually decreased from 6.5 in the mid-1800s to about 5.6 in 1920 (first response in biota) to below 5.0 by 1970, and concentrations of inorganic aluminium rose dramatically beginning about 1960. In the 1990s the lake water chemistry showed clear signs of reversal in acidification brought about by a major decline in S deposition.     

Nitrous Oxide Concentrations in the Soil of a Mown Grassland: Comparison of Model Results with Soil Profile Measurements

Nitrous oxide (N₂O) concentrations in the soilprofile of a fertilized grassland on the Swiss plateauwere measured at irregular intervals during one year.Air samples were taken from air-permeable tubesinstalled at depths between 2 and 100 cm belowthe surface. Highest concentrations in thetopsoil were observed following precipitationafter fertilization. In the subsoil,concentrations were highest when the soil watercontent fell below about 92% after a wet periodand mineral nitrogen was available from priorfertilization. N₂O concentrations in thetopsoil were simulated with the process-basedgrassland ecosystem model PaSim 2.5 (PastureSimulation Model) and compared with measured soilN₂O concentrations. The model simulated wellthe concentration peaks after fertilizeradditions and the concentrations in winter. Butthe simulated baseline concentrations during thegrowing season were overestimated. Possiblecauses for this discrepancy are discussed andsuggestions are made to improve the model.