Modelling cadmium accumulation at a regional scale in the Netherlands

The cadmium content in soils in the rural environment in the Netherlands may increase towards an exceedance of quality standards due to atmospheric deposition and the use of fertilizers and animal manure. To evaluate this problem, a simple dynamic, process-oriented model SOACAS has been developed which is aimed at predicting the accumulation of heavy metals in the topsoil in a regional context. SOACAS describes the fate of a metal in one completely mixed soil compartment using a Freundlich isotherm and analytical equations to solve the mass balance. We tested if it was possible to reconstruct the soil’s present cadmium contents, using independent estimates of historical cadmium loads (‘hind-cast simulation’). About 2500 recent point observations of cadmium contents in rural areas were available. Before comparison, a map was created to translate the point information obtained from the field data to areal average information (resolution 500 × 500 m2), required for comparison with results of SOACAS. A regression model in combination with a locally-weighted smoother within the framework of Generalized Additive Modelling (GAM) was used for this purpose. A realistic geographical pattern could be obtained with very few a priori assumptions. Comparison of the map obtained by the GAM and the map obtained by hind-cast simulation showed that, despite the large uncertainties about historical cadmium loadings, the current cadmium contents were only slightly underestimated by SOACAS. Moreover, the geographical pattern for the observed and simulated contents compared reasonably well. On this basis of this exercise we believe that the model can be used to predict trends of future metal contents as a function of emission/immission scenario’s. Simulations showed that cadmium contents currently decrease in highly polluted areas around industrial plants in the South-Eastern part of the country, and still increase in arable land.

[1]  R. Borger,et al.  EVOLUTION OF ACIDITY, ORGANIC MATTER CONTENT, AND CEC IN UNCULTIVATED SOILS OF NORTH BELGIUM DURING THE PAST 25 YEARS , 1988 .

[2]  Eco,et al.  Determination of field-based partition coefficients for heavy metals in Dutch soils and the relationships of these coefficients with soil characteristics , 1996 .

[3]  D. Kinniburgh,et al.  Adsorption of alkaline earth transition and heavy metal cations by hydrous oxide gels of iron and aluminum , 1976 .

[4]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[5]  Harald Sverdrup,et al.  Modelling recent and historic soil data from the Rothamsted Experimental Station, UK using SAFE , 1995 .

[6]  Willem Bouten,et al.  Modelling soil water dynamics in a forested ecosystem. III: Model description and evaluation of discretization , 1992 .

[7]  S. Trudgill,et al.  Solute modelling in catchment systems , 1995 .

[8]  Johan Bouma,et al.  Derivation of land qualities to assess environmental problems from soil surveys. , 1986 .

[9]  P. McCullagh.  Quasi-Likelihood Functions , 1983 .

[10]  W. Salomons,et al.  Influence of cattle manure slurry application on the solubility of cadmium, copper, and zinc in a manured acidic, loamy sand soil , 1993 .

[11]  S. Majid Hassanizadeh,et al.  Model definition and model validation , 1994 .

[12]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[13]  Johan Bouma,et al.  Transfer functions and threshold values: from soil characteristics to land qualities. , 1987 .

[14]  Gerard B. M. Heuvelink,et al.  Combining soil maps with interpolations from point observations to predict quantitative soil properties , 1992 .

[15]  A. Boekhold,et al.  Field scale behaviour of cadmium in soil. , 1992 .

[16]  Trevor Hastie,et al.  Statistical Models in S , 1991 .

[17]  A. L. Page,et al.  Sorption of trace quantities of cadmium by soils with different chemical and mineralogical composition , 1978, Water, Air, and Soil Pollution.

[18]  C. Appelo,et al.  Geochemistry, groundwater and pollution , 1993 .

[19]  G. Reinds,et al.  Application of soil acidification models with different degrees of process description (SMART, RESAM and NUCSAM) on an intensively monitored spruce site , 1995 .

[20]  Edzer Pebesma,et al.  Mapping groundwater quality in the Netherlands , 1997 .

[21]  Lennart Folkeson,et al.  Fluxes of Cu, Zn, Pb, Cd, Cr, and Ni in temperate forest ecosystems , 1989 .

[22]  William H. Press,et al.  Numerical recipes , 1990 .

[23]  Lbg,et al.  Landelijke beelden van de diffuse metaalbelasting vande bodem en de metaalgehalten in de bovengrond, alsmede de relatie tussen gehalten en belasting , 1996 .

[24]  H. Selim,et al.  Transport of zinc and cadmium in soils: experimental evidence and modeling approaches , 1994 .

[25]  I. Yamane,et al.  Heavy Metal Pollution in Soils of Japan , 1981 .

[26]  P. McCullagh,et al.  Generalized Linear Models , 1972, Predictive Analytics.

[27]  J. Bouma,et al.  Use of soil-map delineations to improve (Co-)kriging of point data on moisture deficits , 1988 .

[28]  W. Chardon.  Mobiliteit van cadmium in de bodem , 1984 .

[29]  C. Williams,et al.  THE ACCUMULATION IN SOIL OF CADMIUM RESIDUES FROM PHOSPHATE FERTILIZERS AND THEIR EFFECT ON THE CADMIUM CONTENT OF PLANTS , 1976 .

[30]  T. Yee,et al.  Generalized additive models in plant ecology , 1991 .