Passive sampling of acidic trace gases Trace gases were monitored using the passive sampling technology as applied in the International Geosphere Atmosphere Chemistry–Deposition of Biogeochemically Important Trace Species (IGAC–DEBITS–Africa/IDAF) network.37These measurements are based on the molecular diffusion of gases onto species-specific impregnated filters.38,39 Passive samplers have been tested and validated since the beginning of 1998 in the IDAF network at six African stations40,41 This technique has also been tested in different tropical and subtropical regions.42Precision of the passive samplers, expressed as mean percentage difference between duplicates, has been found to be in the range 5% – 10% for NO2 and 10% – 25% for SO2.43,44 Monitoring sites were selected at integral one-degree grid intersections, with some additional sites positioned at half-degree grid intervals. To avoid the impacts from local air pollution sources, care was taken to position the sites away from urban centres, industrial point sources or other local disturbances. Samplers were also located away from major roads. Sites were mostly in rural farmland, forestry and natural reserves, sparsely inhabited and often non-electrified locations. In instances where the grid point was close to potential interferences, sites were shifted from the original fixed grid point by distances of not more than 15 km. Passive samplers were exposed in the field for 30-day intervals, corresponding to calendar months, for 2 years (two annual cycles running from September 2005 to August 2007). The monitoring sites were jointly managed with field volunteers who serviced them on a monthly basis. Samples were prepared in the laboratory and couriered or mailed to the site operators. Quality assurance and quality control Quality assurance and control for passive sampling has been described previously28 and so only a short overview is given here. An accredited atmospheric chemistry laboratory with its own quality assurance and control measures was used. Throughout the sampling campaign, measures were undertaken for field quality assurance. These ranged from duplicate sampling and seasonal site visits to monitors and procedural inspections of the monitoring campaign. For internal quality control purposes, comparisons were made between the two annual cycles of the monthly, seasonal and annual means. In addition, data analysis included checking for extreme concentrations that prompted the scientific flagging of outliers. Quality assurance and control for the wet chemistry was undertaken by a licensed industrial laboratory. Ion chromatography chemical analyses performed in the laboratory were annually evaluated by the World Meteorological Organization (WMO) and the results compared to reference values.33Dry gaseous deposition calculation Dry deposition rates were calculated from the measured ambient concentration and dry deposition velocity for the species derived from an inferential model.45,46 The dry deposition rate, F is determined as follows: F = -Vd x C [Eqn 1] where C is the atmospheric concentration and Vd is the deposition velocity (downward flux). For an inferential model a number of meteorological parameters are required. Those parameters were measured and included in an inferential model in a deposition calculation presented by Mphepya.33His work, based on a set of 3 years of monitored data (1996 – 1998), was used to calculate the deposition velocities for two sites – Elandsfontein and Palmer. A mix of 66% grass and 34% maize was considered to be representative of the Elandsfontein site, while a mix of 85% grass and 15% forest was considered for Palmer.33Such surface vegetation mixes were found to be representative for the sites of this study and thus both sets of data could be used for this study. Differentiated deposition velocities for day and night were given for each season33in each year of observation as well as for the 3-year means. The 3-year means (per day/night and season) for both Elandsfontein and Palmer were then averaged and the Elandsfontein–Palmer averages applied to each site in the dry deposition calculations of this study. The Vd differed more in terms of day–night and for different seasons for different sites than for different sites under the same daytime and seasons, indicating that surface cover is not the determining factor of Vd (Table 1). Nevertheless, gaseous deposition velocities remain a major uncertainty as no direct deposition velocities for SO2 and NO2 were measured locally. However, the values shown in Table 1 were considered to be the most suitable ones for this study. Another important step taken in this calculation was to record periods of light and dark for each monitoring site for each month of sampling, based on the geoposition of each monitoring site.47Wet acidic deposition In most cases, the long-term spatial variation in wet deposition is determined more by variations in the precipitation and less by variations in concentrations of acidic species in the rain. The wet deposition flux is determined by calculating the measured concentration of the species of interest in the rain and the precipitation rate.33For this study, measured precipitations were available from the meteorological station network of the South African Weather Service.48To compile the study-specific estimate of the cumulative acidic deposition, wet deposition was estimated by using a 10-year volume weighted mean precipitation chemistry record.34Ion concentrations for the compounds of interest were extracted (Table 2). This included ammonium, which on deposition enters into reactions involving plants, thus contributing to net acidification.49The wet deposition calculation followed an inferential process where the measured compound of interest (acidic ion concentration) was multiplied by the measured (per site) precipitation for mean annual volume:F = C x P [Eqn 2]In this equation F is wet flux (mol/m2 per year), C is measured concentration in precipitation (mol/m3) and P is precipitation rate (mm/year). For this study, sulphate, nitrate and ammonium were included to obtain the overall inorganic – sulphur and nitrogen – deposition loads. Dry and wet deposition results were combined in order to derive a study-specific total acidic deposition estimate from sulphur and nitrogen (excluding organic acids) species. Although no aerosols were calculated separately they were indirectly calculated and included from the modelled soil dust deposition and the measured rainfall chemistry composition data. Critical loads exceedance This section is based on the methodology published by Kuylenstierna et al.20However, several modifications were made for this study, which tested the possibilities based on other reports. Base cation deposition estimation Dry base cation deposition: Kuylenstierna et al.20,21 set the carbonate content in wind-blown dust at 10%.50 Elsewhere, percentages have been cited within the range 3% – 20% .50,51,52 Because of the relatively wide range and high uncertainty of the soil deposition data, it was decided to test an upper and a lower limit, within which the acidic loads could be evaluated against the soil or ecosystems sensitivity. These range limits were the lowest and highest percentages plausible. So the base cation (BC) dry deposition rates were estimated on the one hand to be 20% of the wind-blown and deposited soil dust, and 3% on the other hand. Two concentration ranges of soil BC deposition, 10 meq/m2 – 25 meq/m2 per year and 25 meq/m2 – 50 meq/m2 per year (Kuylenstierna et al.20), were based on the model by Tegen and Fung53 for the study region (Figure 1). The midpoint of each range was used as a mean estimate of BC deposition from soil dust, representing an estimate of 10% calcium carbonate content. This mean value was then used to calculate two extreme values, 3% and 20%, of soil dust BC estimates. Wet base cation deposition: An estimate for wet deposition of BC, using long-term measurements of rain chemistry base cations in the study region,34 was included to infer a total BC deposition and to improve the net acidic load estimate (Table 3). Having four long-term precipitation chemistry sets for four different geographic regions provided the rain chemistry data needed to estimate rain chemistry concentration for the sites in this network. Matching of long-term rain chemistry sets depended on geographical proximity to the rain chemistry measurement site and proximity to the pollution source region or relative remoteness in the background. The positions of the long-term rain chemistry sites in relation to the passive sampling network of this study are indicated in Figure 2. The BC ion concentrations (excluding NH4+) were summed and then multiplied by the total annual recorded precipitation. The annual results were then averaged and the means used as the wet BC component. This BC component was not apportioned for the industrial and sea-source BC. On this basis, two cumulative BC estimates were taken into account when the net acidic deposition was calculated. The lower cumulative BC deposition (3% soil dust BC content with wet BC), as the worst-case scenario, was plotted and interpolated into isopleths for a spatial distribution overview (Figure 3). Soil data An international registry of soil information is updated and maintained by the International Soil Reference and Information Centre (ISRIC), Food and Agriculture Organization (FAO) and United Nations Environment Programme (UNEP) under the International Union of Soil Sciences (IUSS) in a compilation named Soil and Terrain Digital Database (SOTER). Each SOTER database is comprised of two main elements: a geographical component and an attribute component. For this study, the regional SOTER base for southern Africa was used, SOTERSAF (SOTER for Southern Africa, release 1.0, FAO and ISRIC 2003).35 Both geographical and attribute data were then handled by a geographic information system (GIS) programme to create polygons of soil types over the study area. The geographical data were related to an attribute database containing cation exchange capacity and base saturation. These values were then combined to calculate the sensitivity classes using the sensitivity allocation as per the Kuylenstierna et al.20 method. Determining critical loads from soil sensitivity classes The 2-year mean sulphur and nitrogen total deposition rates (total; dry and wet) were converted from deposition rate (kg/ha per year) into units suitable for critical loads assessment, i.e. meq/m2 per year.54 Critical loads are defined in relation to the buffering rate of soils according to the Kuylenstierna method.20 In terms of this definition, critical loads are therefore highly correlated to soil buffering rates, which are in turn related to weathering rates. A categorical classification of soil sensitivity has been adapted from Kuylenstierna et al.20 as an intermediate step in combining acidic deposition rates with soil vulnerability. The ranges of acidic deposition at which each of these five soil classes are susceptible to critical load exceedance are indicated in Table 4. Results Total (dry and wet) deposition loads As expected, total deposition estimates were highest for the site with high dry deposition and relatively high wet deposition. Elandsfontein (Site 17), at the centre of the industrial Highveld, had the highest deposition with an overall total acidic deposition of 23.2 kg/ha per year for the first year and 15.8 kg/ha per year for the second year. Sites downwind of the industrial Highveld (Site 1, Standerton and Site 2, Amersfoort) had the second highest cumulative deposition results for both annual cycles. The lowest total acidic depositions were for Site 25 near Phalaborwa (3.1 kg/ha per year in the first year and 2.9 kg/ha per year in the second year), which had low dry and wet deposition, specifically as a result of low precipitation affecting the wet deposition contribution. When the results for all the sites were compared year on year, substantial differences were noted; for certain sites the differences varied by up to 8.6 kg/ha per year (Table 5). This variation can be directly attributed to the differences in precipitation from one year to the next, which directly affects the wet deposition estimate and thus the cumulative (dry and wet) deposition for certain sites. Deposition of acidic compounds over the study region The 2-year mean results for all sites were interpolated into isopleths through a GIS programme for a spatial distribution view of the total (dry and wet) acidic deposition loads (Figure 4). Both the upper and lower BC deposition results were subtracted from the mean total (dry and wet) acidic deposition values to obtain the net lower and higher acidic deposition values, respectively. The lower net acidic deposition had no anion deposition greater than the cation deposition, except in the industrial Highveld, and represents the best-case scenario. It was not used for further assessment. The higher net acidic deposition result represents the worst-case scenario and was used for further assessment (Table 6). The fourth column in Table 6 (higher net acidic deposition) was plotted and interpolated to create a spatial layer (Figure 5). The levels of acid deposition occur over the central industrialised Highveld, close to the major sources. There seem to be two regions of deposition: the first is an intense ovoid-shaped area over the source regions, pointed towards the south-east as expected, covering a 2° by 2° area (~50 000 km2). This region can be considered the most vulnerable to adverse effects from acid deposition. A secondary zone of deposition of lower intensity extended to the south and south-south-west, over the eastern Free State and Lesotho (Figure 5). Soil sensitivity The next stage of the analysis was to compare the net acidic deposition with calculated critical loads for the soils in the study area. The soil sensitivity polygons are presented in Figure 6. These were then assigned critical loads as described above. Overall, most soils fell into the least sensitive class. However, certain regions are characterised by sensitive soils (Classes 1 and 2, respectively).
Acidification risk assessment in the study region The higher net level acid deposition layer was laid over (subtracted from) the soil sensitivity values map corresponding to the higher (25-50-75-100) soil sensitivity estimate as per Table 4. The output is shown as the exceedance of critical loads of acidity (meq/m2 per year) for the soils (Figure 7). Exceedance of buffering capacity: exceedance of critical loads The exceedance map shown in Figure 7 indicates the impacted areas and the degree to which net acidic deposition exceeds critical loads. These areas indicate ecosystems at risk, where acidification and consequent impacts to terrestrial and aquatic ecosystems have occurred or may occur. In contrast, the grey areas indicate areas not at risk from critical load exceedance. The exceedance map shows that only the central industrial Highveld source area has the potential for high levels of exceedance. Two smaller areas have estimated exceedances of the soil buffering capacity in the intermediate range (51 meq/m2 per year – 75 meq/m2 per year) – one just east of Amersfoort and the other north of Standerton (Figure 7). A relatively large zone of potential exceedance in the lowest range (1 meq/m2 per year – 25 meq/m2 per year) occurred in the KwaZulu-Natal Midlands (in a triangle, south-east from Newcastle, north-east from Escourt and north-west from Kranskop). Several scattered small areas of sensitive soils are identified as intermediate and low levels of potential exceedance. Discussion and conclusions This study investigated the risk of acidification from both sulphur and nitrogen non-organic acidifying species, anthropogenic and natural, through the method of critical loads mapping. It follows axiomatically that changes in the ecosystem composition and function could occur as deposition continues. The results presented indicate that the upper limit estimate of acidic deposition may exceed the critical loads capacity of the local soils in several areas. The level of exceedance is comparable to the level of acidic deposition exceedance in other regions of the world with high emission and deposition rates.49,55,56 A high critical loads (for surface deposition) exceedance was determined in the western and central Highveld industrial region and adjacent area to the north and not, as anticipated, downwind towards the escarpment and the major forestry areas. Several smaller areas located downwind from major sources showed exceedance, although the levels of exceedance were low. Unsurprisingly, these areas of exceedance are on the Highveld plateau. Fieldwork is suggested to test soil acidity and check the effects of acidification in the ecosystems of the areas that showed exceedance. Critical loads are based on a steady-state concept. Many ecosystems are not in equilibrium with present or projected depositions, because there are processes which delay the reaching of steady state for years or decades or longer. Soil acidification is also a natural process and it is difficult to quantify the extent to which this process is being accelerated by atmospheric acidic deposition. It is also extremely difficult to demonstrate in the field a clear connection between critical load exceedance, declining base cation to aluminium ratios in the soil solution and biological damage resulting from such changes. Therefore dynamic modelling is suggested for estimation of durations involved in attaining a certain chemical state in response to deposition scenarios.57 Of course, the method applied in this study is not the only method available. The suitability of other methods, particularly direct methods such as the determination of an index of acid neutralising capacity,58 in local circumstances should be explored. Our conclusion challenges the current understanding of the status of acid deposition from Highveld industrial emissions. However, uncertainties in the BC deposition data, wet acidic deposition estimates and soil database attributes remain. In this respect, measurements of soil dust deposition, with specific attention to BC content, is needed. There is a further need for direct wet deposition measurements on an extensive spatial and temporal scale. Such additional measurements would provide a greater level of confidence and would be needed to validate and implement costly control and mitigation policies on trace gas emissions, or to substantiate a decision that such mitigation is not justified on ecological grounds. In either case, the location of sensitive soils with respect to the deposition pattern is a step in the right direction towards identifying and monitoring ecosystems at risk. The study was based on a joint project of the Department of Geography, Environmental Management and Energy Studies, University of Johannesburg; the Climatology Research Group, University of the Witwatersrand, Johannesburg; and the Atmospheric Chemistry Research Group, North West University, Potchefstroom. Thanks are extended to all those involved. The study was funded initially by Anglo American Corporation, and then jointly by Eskom and Sasol, with a matching contribution from the Department of Trade and Industry THRIP programme. M. Josipovic acknowledges the support of a National Research Foundation (NRF) Grantholder’s Bursary, an Oppenheimer Memorial Trust Fellowship (16068) and a Sasol Fund supplementary bursary through the Faculty of Science, University of Johannesburg. The work was supported in part by an NRF Focus Area Project grant to HJA - Sustainability Studies using GIS & Remote Sensing FA2005040600018. 1.McCormickJAcid earth: The politics of acid pollution. 3rd ed2.HeldGScheifingerHSnymanGTosenGZunckelMThe climatology and meteorology of the Highveld. In: Held G, Gore BJ, Surridge AD, Tosen GR, Turner CR, Walmsley RD, editors. Air pollution and its impacts on the South African Highveld.3.GallowayJNAcid deposition: Perspectives in time and space.1995851524doi:10.1007/BF00483685 4.SivertsenBCMataleCPereiraLRMSulphur emissions and trans-frontier air pollution in southern Africa. Report 355.RorichRPGalpinJSAir quality in the Mpumalanga Highveld region1998941091146.WellsRBLloydSMTurnerCRNational air pollution source inventory. In: Held G, Gore BJ, Surridge AD, Tosen GR, Turner CR, Walmsley RD, editors. Air pollution and its impacts on the South African Highveld.7.United Nations Economic Commission for Europe. Conclusions and draft recommendations of the workshops on critical levels for forests, crops and materials and on critical loads for sulphur and nitrogen.8.United Nations Economic Commission for Europe Convention on Long-range Transboundary Air Pollution [homepage on the Internet].Available from: www.unece.org/env/lrtap 9.SkeffingtonRAThe use of critical loads in environmental policy making: A critical appraisal.19993324510.SverdrupHde VriesWHenriksenAMapping critical loads: A guide to the criteria, calculations, data collection and mapping of critical loads. UN-ECE and Nordic Council of Ministers, Report Miljörapport. 11.TroedssonTNykvistNMarklära och markvård [Soil science and soil tending]. 12.RodheHCowlingEGalballyIGallowayJHerreraRAcidification and regional air pollution in the tropics. In: Rodhe H, Herrera R, editors.13.KuylenstiernaJCICambridgeHCinderbySChadwickMJTerrestrial ecosystem sensitivity to acidic deposition in developing countries.19958523192324doi:10.1007/BF01186180 14.Environment Canada15.OmernikJMTotal alkalinity of surface waters, 1:7,500,00016.ZhaoDXiongJMaoJ.Mapping sensitivity, critical loads and exceedance for China.17.WadaHIwasaYArimitsuKMap for assessing susceptibility of Japanese soils to acid precipitation. Scale 1:2,000,000.18.ChartresCJSmithCJMcDonaldWSNobleADPotential effects of acid deposition on Australian soils. Workshop presented at: Validation of Sensitivity Maps to Acidic Deposition; 1996; Bangkok, China.19.HettelinghJ-PSverdrupHZhaoDDeriving critical loads for Asia.19958525652570 doi:10.1007/BF01186220 20.KuylenstiernaJCIRodheHCinderbySHicksKAcidification in developing countries: Ecosystem sensitivity and the critical load approach on a global scale.2001302028doi:10.1579/0044-7447-30.1.20PMid:11351789 21.KuylenstiernaJCIRodheHCinderbySHicksKGlobal assessment of terrestrial ecosystem sensitivity to acidic deposition [homepage on the Internet].Available from: http://www.sei.se/dload/seiy/globalassess.pdf 22.KuylenstiernaJCIChadwickMJThe relative sensitivity of ecosystems in Europe to the indirect effects of acidic depositions in regional acidification models. In: Kämäri J, Brakke DF, Jenkins A, Norton SA, Wright RF, editors. 23.HornungMSkeffingtonRAeditorsCritical loads: Concepts and applications. Proceedings of the 28th ITE Symposium; 1992; London, United Kingdom.24.DowningRJHettelinghJ-Pde SmetPAMeditorsCalculation and mapping of critical loads in Europe: Status report 1993.25.PoschMde SmetPAMHettelinghJ-PDowningRJeditorsCalculation and mapping of critical thresholds in Europe.26.OlbrichKASkoroszewskiRTaljaardJZunckelMAn acid deposition risk advisory system: A test study of Region F to define actual and critical loads.27.Van TienhovenAMOlbrichKASkoroszewskiRTaljaardJZunckelMApplication of the critical loads approach in South Africa.19958525772582doi:10.1007/BF01186222 28.JosipovicMAnnegarnHJKneenMAPienaarJJPikethSJConcentrations, distributions and critical levels exceedance assessment of SO2, NO2 and O3 in South Africa. 2010171181196doi:10.1007/s10661-009-1270-5,PMid:20024614 29.World Meteorological Organization.http://www.wmo.ch/pages/prog/arep/gaw/gaw-reports.html 30.United Nations Environment Programme [homepage on the Internet]. 31.National Acid Precipitation Assessment Program [homepage on the Internet].http://ny.cf.er.usgs.gov/napap/index.html 32.United Nations Economic Commission for Europe.http://www.unece.org/env/lrtap/emep/welcome.html 33.MphepyaNJAtmospheric deposition characteristics of sulphur and nitrogen compounds in South Africa. 34.TurnerCRde Beer GHA review of ten years’ rain quality data from the South African interior35.International Soil Reference and Information Centre [homepage on the Internet]. http://www.isric.org/NR/exeres/BC5BB444-DB94-4FE3-B102-2872B0B1A9F1.htm 36.FreimanMTPikethSJAir transport into and out of the industrial Highveld region of South Africa. 2003429941002doi:10.1175/1520-0450(2003)042<0994:ATIAOO>2.0.CO;2 37.PienaarJJArtaxoPGaly-LacauxCChow PengLBalaRPresentation of the DEBITS’s activity: Past, present and future.38.FermMA sensitive diffusional sampler. 39.FermMLindskogASvanbergP-ABoströmC-A[New measurement technique for air pollutants]. Kemisk Tidskrift.199413032Swedish. 40.Al OurabiHLacauxJPMeasurement of the atmospheric concentrations of NH3, NO2, HNO3 in tropical Africa by use of diffusive samplers. 41.LacauxJPDEBITS’s activity in Africa: Atmospheric deposition in northern hemisphere of tropical Africa. 42.FermMRodheHMeasurements of air concentrations of SO2, NO2 and NH3 at rural and remote sites in Asia.1997271729doi:10.1023/A:1005816621522 43. AyersGPKeywoodMDGillettRManinsPCMalfroyHBardsleyTValidation of passive diffusion samplers for SO2 and NO2. 44.CruzLPSCamposVPSilvaAMCTavaresTMA field evaluation of a SO2 passive sampler in tropical industrial and urban air2004383764256429doi:10.1016/j.atmosenv.2004.07.02245.HicksBBBaldochiDDMeyersTPHoskerRPMattDRA preliminary multiple resistance routine for deriving dry deposition velocity from measured quantities.198736311330doi:10.1007/BF0022967546.HicksBBHoskerRPMeyersTPWomackJDDry deposition inferential measurement technique – I: Design and test of a prototype meteorological and chemical system for determining dry deposition. 19911025A2361237047.US Navy Office. Astronomical Information Center, Data Service, Day/Night Across the Earth [homepage on the Internet].http://aa.usno.navy.mil/data/docs/earthview.php 48.South African Weather Service [homepage on the Internet]. c2006 [cited 2008 May 22]. Available from:http://www.weathersa.co.za/ 49.World Meteorological Organization (WMO). 50.GomesLGilletteDAA comparison of characteristics of aerosol from dust storms in central Asia with soil-derived dust from other regions.199327A2539254451.Food and Agriculture Organization of the United Nations Educational, Scientific and Cultural Organization (FAO-UNESCO). Soil map of the world. 52.AvilaAQueraltIGallartFMartin-Vide JAfrican dust over north-east Spain: Mineralogy and source regions. In: Guerzoni S, Chester R, editors. The impact of desert dust across the Mediterranean.53.TegenIFungIContribution to the atmospheric mineral aerosol load from land surface modification.19951001870718726doi:10.1029/95JD0205154.PoschMde SmetPAMHettelinghJ-PDowning RJ, editors. Modelling and mapping of critical thresholds in Europe: Status report 2001 of the Coordination Centre for Effects, National Institute for Public Health and the Environment, RIVM report no. 55.BouwmanAFvan VuurenDPDerwentRGPoschMA global analysis of acidification and eutrophication of terrestrial ecosystems.2002141349382doi:10.1023/A:1021398008726 56.WhelpdaleDMSummersPWSanhuezaEA global overview of atmospheric acid deposition fluxes.199748217247 doi:10.1023/A:1005708821454 57.United Nations Economic Commission for Europe Convention on Long-range Transboundary Air Pollution. Manual on methodologies and criteria for modelling and mapping critical loads & levels and air pollution effects, risks and trends [homepage on the Internet]. c2004 [cited 2005 Sept 29]. Available from:http://www.icpmapping.org 58.DoddsHAFeyMVEvaluation of some systems for classifying soil sensitivity to acid deposition in the South African Highveld.19981494199