Critical Loads – Is There A Need For A New Concept?

The concept of critical loads has been an important andsuccessful tool for the development of control strategiesfor transboundary air pollution in Europe. The use of theconcept has led us to a situation where very few areas inEurope will have an exceedance of critical loads foracidification in 2010, indicating that the benefits offurther control acidifying substances will be lessuseful. The critical loads concept does not, however,take into account the large benefits of further controlin damaged systems but where critical loads are nolonger exceeded. In this paper we discuss the importanceof widening the critical loads concept to include thesebenefits and we propose an additional effect-relatedmeasure, Dynamic Impact Analysis, to be included infurther control strategies and assessments. With such aconcept the actual situation and its further developmentwill be included in assessments and control strategies.     

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.     

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'.     

European Critical Loads of Cadmium, Lead and Mercury and their Exceedances

Critical loads of cadmium, lead and mercury were computed by 18 countries of the LRTAP Convention. These national data were collated into a single database for the purpose of identifying sensitive areas in Europe. Computing exceedances, i.e. comparing the critical loads to atmospheric deposition, shows that cadmium was not a widespread risk in 2000, that the risk from lead deposition has decreased since 1990 but was still widespread in 2000, and that the risk from mercury remains high without much change from 1990 to 2000 in most of the countries.     

The Role of Mineral Weathering Rate Determinations in Generating Uncertainties in the Calculation of Critical Loads of Acidity and their Exceedance

Current legislation within Europe aimed at limitingecosystem damage resulting from inputs of atmosphericpollution is based on the critical load concept. Mineralweathering rates are central to the calculation ofcritical loads (acceptable levels) of acid deposition.The authors have undertaken a number of studies whichillustrate the complications and limitations inherent inpredicting mineral weathering rates and the implicationswhich these have for critical loads calculations andmapping. Calculated weathering rates and critical loadsfor two acid-sensitive parent materials (greywackes andgranites) are presented and are used to illustrate theimpact that uncertainty can have on critical loadexceedances. The results have obvious implications forportraying the uncertainties of critical loads to policy makers.     

Predicting Freshwater Critical Loads from Catchment Characteristics using National Datasets

Maps of freshwater critical loads are used toguide emission strategies for sulphur and nitrogen bothnationally and internationally. Water chemistry data arerequired to calculate critical loads and the production ofnational maps therefore relies on the existence of extensivechemistry datasets. However, the data required to calculatecritical loads are not readily available for all sites. Thisarticle explores how empirical statistical models mightpotentially be used to predict critical loads using nationallyavailable datasets representing a range of catchmentcharacteristics. Initially a global regression model forexplaining freshwater critical load variation across a broadspectrum of catchment types (from lowland agricultural tomountain lakes) throughout mainland Britain is described. Whenattention is focused on more specific catchment types (i.e.upland and non-arable) it is shown that the global model hasless explanatory power. A regionalisation of Great Britain(based on 100 km grid squares) shows that the global modelcannot necessarily be applied successfully within a narrowerregional context. Separate analyses were undertaken on each ofthe regional subsets using backward selection regression. Thevariables emerging as significant predictors variedsubstantially across the regions, as did the explanatory powerof the models. This was also the case when the analysis wasconfined to upland and non-arable catchments. This approachcould be developed so that critical loads assessments can bemade for populations of standing waters rather than simplythose for which water chemistry is available.     

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.     

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.     

The Importance of Selecting Appropriate Criteria for Calculating Acidity Critical Loads for Terrestrial Ecosystems Using the Simple Mass Balance Equation

The Simple Mass Balance (SMB) equationis commonly used throughout Europe for thecalculation of acidity critical loads for forestsoils. Different criteria can be set in themodel depending on whether the receptor (e.g. treeroots) is more sensitive to the toxic effects ofaluminium or to unfavourable pH conditions. Thispaper examines the effects on critical loadscalculations of using different criteria andcritical limits, and demonstrates the importanceof selecting the most appropriate and justifiablecriteria for the chosen receptor, since they caneffect the critical loads values obtained. Abrief review of the range of different criteriaand limits used throughout Europe is included. In addition, the gibbsite equilibrium constant,used in the SMB equation to represent therelationship between dissolved aluminium andhydrogen ions in soil solution, is discussed. This relationship is not generally described inthe literature as a criterion in the equation,but this work highlights the effects differentgibbsite values have on critical loadcalculations and the importance of applying themost appropriate value for the soil in therooting zone of the receptor.     

Cost-effectiveness Analysis of Reducing the Emission of Nitrogen Oxides in Asia

The purpose of this study is to evaluate cost-effective reduction strategies for nitrogen oxides (NO _x ) in the Asian region. The source-receptor relationships of the Lagrangian “puff” model of long-range transportation, ATMOS-N, were used to calculate the wet/dry deposition of the nitrogen (N) in Asia. Critical loads of N deposition in Asia were calculated from the relationships between the critical load of sulfur (S) and balance of N in and out using the data of S critical load of RAINS-ASIA. The cost functions of N reduction of Asian countries were derived by the regression analysis with the data of cost functions of European countries used in RAINS. In order to assess the environmental impact, the gaps between N deposition and critical load of N were calculated. The emission of NO _x was reduced in some cases of this model, and the changes of gaps between N deposition and critical load were observed as well as the changes of the reduction cost. It is shown that a uniform reduction of NO _x emissions by countries in Asia is not cost-effective strategy.     

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.     

An Analysis of the Structure of the Simple Mass Balance Equation: Implications for Testing National Critical Loads Maps

The critical loads concept is used by the UN-ECEConvention on Long Range Transboundary Air Pollution(CLRTAP) for setting pollution reduction targets.Increasing numbers of countries are adopting the SimpleMass Balance equation (SMB) to calculate critical loads ofacidifying S and N for forest soils. The equation is madeup of a series of mass balances each of which is used tocalculate a leaching flux. The assumptions in the SMBequation were investigated by testing its ability topredict current sulphur load and by comparing each of thecalculated leaching fluxes to measured rates. It was notpossible to predict current sulphur load at our sites usingthe SMB equation. The leaching tests demonstrated that,primarily due to its steady state assumptions, the SMBequation generates critical loads that are theoretical longterm estimates of risk, and are untestable. Thesimplifying assumptions sometimes lead to illogicalresults. Some of these can be improved by adding a final,simple but dynamic, calculation step to determine theexpected time until effects are observed. Theacceptability of combining annual average data, which bestapproximates steady state, with a biological indicator isquestionable. It is not possible to test critical loadscalculated using the SMB equation when applied with all ofits assumptions but it is possible to test its fundamentalapproach using non steady state data.     

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.     

Aluminium: The Need for a Re-Evaluation of Its Toxicity and Solubility in Mature Forest Stands

Aluminium (Al) is a key element in critical loadcalculations for forest. Here, we argue for re-evaluating theimportance of Al. Effects of two levels of enhanced Alconcentrations and lowered Ca:Al ratios in the soil solutionin a field manipulation experiment in a mature spruce stand(1996–1999) on tree vitality parameters were tested. Inaddition, Al solubility controls were tested. Various loads ofAl were added to forest plots by means of an irrigationsystem. Potentially toxic Al concentrations and criticalratios of Ca to inorganic Al were established. The ratio of Cato total Al was not a suitable indicator for unfavourableconditions for plant growth. No significant effects on crowncondition, tree growth and fine root production were observedafter three years of treatment. In 1999, foliar Mg content inthe highest Al addition treatment had declined significantly.This agreed with the known response to Al stress of seedlingsin nutrient solution experiments. No support was found forusing the chemical criterion Ca:Al ratio in soil solution,foliar and root tissue as an indicator for forest damage dueto acidification. Al solubility was considerably lower thanimplied by the assumption of equilibrium with gibbsite,particularly in the root zone. The gibbsite equilibrium iscommonly used in critical load models. Substitution of thegibbsite equilibrium with an Al-organic matter complexationmodel to describe Al solubility in soil water may have largeconsequences for calculation of critical loads. The resultsindicate that critical load maps for forests should bereconsidered.     

Comparison of Critical Load Exceedance and Its Uncertainty Based on National and Site-specific Data

Critical loads have been used to develop international agreements on acidifying air pollution abatement, and within the UK and other countries, to develop national policies for pollution abatement. The Environment Agency (England and Wales) has regulatory obligations to protect sites of high conservation value from the threat of acidification, and hence requires a practical methodology for acidification assessments at the site-specific scale. The Environment Agency has therefore posed the question: Are the national critical load exceedance models sufficiently robust to form the basis for methods to assess harm to individual sites or are they only useful for national policy development? In order to provide one measure of the appropriateness of applying the models at the site-specific scale we incorporated estimates of uncertainty in both national and site-specific data into the calculation of critical load exceedance for individual sites. The exceedance calculations use data from a wide range of sources and the accuracy of the exceedance will be influenced by the accuracy of the input data sets. Using Monte Carlo methods to incorporate the uncertainty in the input data sets into the calculation a distribution of critical load exceedance values is generated rather than a single point estimate. This paper compares uncertainty analyses for coniferous forested sites in England and Wales using both national scale and site-specific data sets and uncertainty ranges.     

The Applicability of National Critical Loads Data in Assessing Designated Sites

Critical loads have been successfully used within Europe in the development of effects-based policies for pollution abatement, including the Second Sulphur Protocol and the Protocol to abate acidification, eutrophication and ground-level ozone (CLRTAP, 1979). This success has encouraged the UK Environment Agency and Conservation Agencies to use the national critical load maps as a screening tool in assessing the threats from acidification and eutrophication to designated (Natura 2000) sites. The UK maps of critical loads are based on national-scale data sets appropriate for national-scale assessments, and were never intended for use at the site-specific level. Site-based assessments are often targeted at Special Areas of Conservation, a sub-set of the UK Natura 2000 sites. The spatial data available includes the boundaries of the sites but not the location of the designated features. Ancillary data is variable from one site to another; habitat types may be described in detail with cross-reference to classes of the National Vegetation Classification (NVC: Rodwell, 1991 et seq), but information available on soils and geology is generalised and has not been related to the habitats or species being protected. Hence it can be difficult to relate the individual sites to the national maps, even where appropriate to do so. This paper examines the underlying uncertainties in the national critical load maps showing how the maps could give misleading results if used for site-specific assessments. It also includes advice on how to determine when the national data may be appropriate as a policy-tool at the site-level.     

Critical Loads Exceedances in Germany and their Dependence on the Scale of Input Data

Critical load exceedances have been used as an effects-related parameter for guiding international air emission control negotiations. High-resolution critical load data are combined with low-resolution deposition data.This article shows that doing so systematicallyunderestimates `true' critical load exceedances as obtainedfrom combining critical load and deposition data of identicalhigh spatial resolution. 95th percentile critical loadexceedances in EMEP grids based on high resolution depositiondata are 60 and 150% higher (mean values for nutrientnitrogen and acidity, respectively) than critical loadexceedances based on the low resolution EMEP depositionmodel. The latter are used in international negotiations. Differences in individual EMEP grid squares vary betweeninsignificantly different from zero and 340%, depending onregional deposition and critical load characteristics andcritical load types (nutrient nitrogen versus acidity).Exceedances based on high-resolution deposition values arealso compared to EMEP grid averages of these values forforests only. This comparison excludes the effect ofsystematically higher depositions to forests. Still, thescale difference of (averaged, low-resolution) deposition and(high-resolution) critical loads data yields underestimatesof the 95th percentiles by on average ca. 20%.These systematic errors due to the scale dependence should beborne in mind when interpreting effects of internationalemission control measures.     

Critical Loads of Acidity for Forest Soils: Tentative Modifications

We reviewed the current methods for calculatingcritical loads of acidity for forest soils. The consequencesof four sets of assumptions concerning the soil modelstructure, parameter values and the critical loads criterionwere explored by comparing the values of the averageaccumulated exceedance (AAE) calculated for Finland withdeposition values for the year 1995. The AAE index is given inthe unit of deposition and is a measure of how far a region isfrom being protected in terms of fulfilling a certaincriterion, taking into account the size of the ecosystem areas.Using a critical limit for the molar ratio of theconcentrations of base cations to aluminium in soil solutiongave the lowest average accumulated exceedance. Assumingorgano-aluminium complexes and leaching of organic anions gaveAAE = 4 eq ha^-1 a^-1, which was close to the valueobtained with the standard approach used in Finland, assuminggibbsite equilibrium and no leaching of organic anions,yielding AAE = 5 eq ha^-1 a^-1. With a critical basesaturation limit, instead of the concentrations criterion, theAAE index was 17 eq ha^-1 a^-1. The highest averageaccumulated exceedance (AAE = 25 eq ha^-1 a^-1),corresponding to the lowest critical load, was obtained whenthe effects-based criterion (critical concentration or criticalbase saturation) was substituted with one restricting thedeterioration of the neutralizing capacity of the soil, ANC _le(crit) = 0. These tests illustrate the variabilityof the critical load values for acidity that can be introducedby changing the criterion or by varying the calculation method,without, however, representing the extreme values of criticalloads that could be derived.     

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.     

Concentrations,Cumulative Exposure and Critical Levels of Ozone in Ireland

Concentrations of tropospheric ozone(O₃) and exceedance of critical levels to vegetationhave been investigated and mapped for Ireland. Hourlyozone concentration data (1995–1997) at 7 sevenmonitoring stations and the CORINE landcover database,supported by a Geographical Information System, wereused. AOT40 (Accumulated exposure Over a Threshold of 40ppb) was calculated for daylight hours for each station,and mapped using surface interpolation. Average O₃concentrations vary from year to year, and were estimatedto be 28 ppb, 26 ppb and 24 ppb for 1995, 1996 and 1997respectively. Ozone concentrations show a large diurnalvariation, with a maximum in the afternoon and a minimumat night-time. The critical level for crops and (semi-)natural vegetation was exceeded in all years examined.The highest exceedance occurred in 1995, where thecritical level was exceeded for almost 35% of the mappedarea. Approximately 15% and 1% of the mapped area wasexposed to exceedance levels during 1996 and 1997respectively. The maximum cumulative exposures (AOT40)were approximately 5000, 3890 and 3230 ppbh in 1995, 1996and 1997 respectively. The critical level for forests wasnot exceeded during the period of investigation.