Modeling carbon sequestration under zero tillage at the regional scale. I. The effect of soil erosion

Zero tillage is recognized as a potential measure to sequester carbon dioxide in soils and to reduce CO₂ emissions from arable lands. An up-scaling approach of the output of the Environmental Policy Integrated Climate (EPIC) model with the information system SLISYS-BW has been used to estimate the CO₂-mitigation potential in the state of Baden-Württemberg (SW-Germany). The state territory of 35,742 km² is subdivided into eight agro-ecological zones (AEZ), which have been further subdivided into a total of 3976 spatial response units. Annual CO₂-mitigation rates where estimated from the changes in soil organic carbon content comparing 30 years simulations under conventional and zero tillage. Special attention was given to the influence of tillage practices on the losses of organic carbon through soil erosion, and consequently on the calculation of CO₂-mitigation rates. Under conventional tillage, mean carbon losses through erosion in the AEZ were estimated to be up to 0.45 Mg C ha⁻¹ a⁻¹. The apparent CO₂-mitigation rate for the conversion from conventional to zero tillage ranges from 0.08 to 1.82 Mg C ha⁻¹ a⁻¹ in the eight AEZ, if the carbon losses through soil erosion are included in the calculations. However, the higher carbon losses under conventional tillage compared to zero tillage are composed of both, losses through enhanced CO₂ emissions, and losses through intensified soil erosion. The adjusted net CO₂-mitigation rates of zero tillage, subtracting the reduced carbon losses through soil erosion, are between 0.07 and 1.27 Mg C ha⁻¹ a⁻¹ and the estimated net mitigation rate for the entire state amounts to 285 Gg C a⁻¹. This equals to 1045 Gg CO₂-equivalents per year with the cropping patterns in the reference year 2000. The results call attention to the necessity to revise those estimation methods for CO₂-mitigation which are exclusively or predominantly based on the measurements of differential changes in total soil organic carbon without taking into account the tillage effects on carbon losses through soil erosion.     

Modeled carbon sequestration variation in a linked erosion–deposition system

Unraveling the consequences of hydrologic transport on carbon (C) storage will help identify feedbacks between land management alternatives, climate change, and soil-vegetation-atmospheric-transfers (SVATs) of C. There is a need for theoretically driven models of erosion and deposition that includes transport induced mineralization to better understand the controls on SVATs of C. Here we present a model developed using a systems-dynamic approach that coupled C-SVATs at a 2-day resolution with a discrete event erosion–deposition model occurring with a prescribed return interval. Five possible mass-balance transformations of C occurring between the two patches were explicitly modeled: net primary production (NPP), decomposition, erosion, transport induced mineralization, and deposition. The net C-SVAT, NPP minus decomposition, exhibited three stable points of no net C flux. Starting with arbitrary initial C pool in each patch above the bifurcation point, the model approached a quasi-steady state, which included both the short-term and longer term consequences of erosion; in the baseline simulation 5080 g C m⁻² was stored prior to erosion and 100 years of low intensity erosion 4840 g C m⁻² SOC remained. Low intensity erosion also generated spatial heterogeneity; from an initial homogeneous distribution to 40% of the C stored in the eroded patch and 60% of the C stored in the deposition patch. Erosion reduction resulted in a corresponding increase in total soil C content that was positively related to the magnitude of erosion reduction. In conjunction with providing a modeling framework for reducing the uncertainty in C-SVAT, this model is a prototype of a growing theory of ecosystem processes within spatially explicit landscapes, a meta-ecosystem model.     

Dynamics of soil organic matter in primary and secondary forest succession on sandy soils in The Netherlands: An application of the ROMUL model

We applied the simulation model ROMUL of soil organic matter dynamics in order to analyse and predict forest soil organic matter (SOM) changes following stand growth and also to identify gaps of data and modelling problems. SOM build-up was analysed (a) from bare sand to forest soil during a primary succession in Scots pine forest and (b) on mature forest soil under Douglas fir plantations as an example of secondary succession in The Netherlands. As some of the experimental data were unreliable we compiled a set of various scenarios with different soil moisture regime, initial SOM pools and amount and quality of above and below ground litter input. This allowed us to find the scenarios that reflect the SOM dynamics more realistically. In the Scots pine forest, total litter input was estimated as 0.50 kg m⁻² year⁻¹. Two scenarios were defined for the test runs: (a) forest floor moisture regimes—‘dry, mesic and hydric’ and (b) augmenting a root litter pool with three ratios of needles and branches to roots: 1:1, 1:1.5 and 1:2.0. The scenario finally compiled had the following characteristics: (a) climate for dry site with summer drought and high winter moisture of forest floor; (b) a litter input of 0.25 kg m⁻² year⁻¹ above ground and 0.50 kg m⁻² year⁻¹ below ground; (c) a low nitrogen and ash content in all litter fall fractions. The test runs for the estimation of the initial SOM pools and the amount and proportion of above and below ground litter fall were also performed in the Douglas fir plantation. The inputs of above ground litter tested in various combinations were 0.30 and 0.60 kg m⁻² year⁻¹, and below ground litter 0.30, 0.60 and 0.90 kg m⁻² year⁻¹. The scenario that fitted the experimental data had an SOM pool of 20–25 kg m⁻², an aboveground litter input of 0.6 kg m⁻² year⁻¹and a below ground litter input of 0.9 kg m⁻² year⁻¹. The long-term simulation corresponded well with the observed patterns of soil organic matter accumulation associated with the forest soil development in primary and secondary succession. During primary succession in Scots pine forest on dry sand there is a consistent accumulation of a raw humus forest floor. The soil dynamics in the Douglas fir plantation also coincide with the observed patterns of SOM changes during the secondary succession, with SOM decreasing significantly under young forest, and SOM being restored in the older stands.     

A two-dimensional simulation model of phosphorus uptake including crop growth and P-response

Modelling nutrient uptake by crops implies considering and integrating the processes controlling the soil nutrient supply, the uptake by the root system and relationships between the crop growth response and the amount of nutrient absorbed. We developed a model that integrates both dynamics of maize growth and phosphorus (P) uptake. The crop part of the model was derived from Monteith's model. A complete regulation of P-uptake by the roots according to crop P-demand and soil P-supply was assumed. The soil P-supply to the roots was calculated using a diffusion equation and assuming that roots behave as zero-sinks. The actual P-uptake and crop growth were calculated at each time step by comparing phosphate and carbohydrate supply–demand ratios. Model calculations for P-uptake and crop growth were compared to field measurements on a long term P-fertilization trial. Three P-fertilization regimes (no P-fertilization, 42.8 kg P ha⁻¹ year⁻¹ and 94.3 kg P ha⁻¹ year⁻¹) have led to a range of P-supply. Our model correctly simulated both the crop development and growth for all P-treatments. P-uptake was correctly predicted for the two non-limiting P-treatments. Nevertheless, for the limiting P-treatment, P-uptake was correctly predicted during the early period of growth but it was underestimated at the last sampling date (61 day after sowing). Several arguments for under-prediction were considered. However, most of them cannot explain the observed magnitude in discrepancy. The most likely reason might be the fact that biomass allocation between shoot and root must be modelled more precisely. Despite this mismatch, the model appears to provide realistic simulations of the soil–plant dynamic of P in field conditions.     

Estimating California ecosystem carbon change using process model and land cover disturbance data: 1951-2000

Land use change, natural disturbance, and climate change directly alter ecosystem productivity and carbon stock level. The estimation of ecosystem carbon dynamics depends on the quality of land cover change data and the effectiveness of the ecosystem models that represent the vegetation growth processes and disturbance effects. We used the Integrated Biosphere Simulator (IBIS) and a set of 30- to 60-m resolution fire and land cover change data to examine the carbon changes of California's forests, shrublands, and grasslands. Simulation results indicate that during 1951-2000, the net primary productivity (NPP) increased by 7%, from 72.2 to 77.1 Tg C yr⁻¹ (1 teragram = 10¹² g), mainly due to CO₂ fertilization, since the climate hardly changed during this period. Similarly, heterotrophic respiration increased by 5%, from 69.4 to 73.1 Tg C yr⁻¹, mainly due to increased forest soil carbon and temperature. Net ecosystem production (NEP) was highly variable in the 50-year period but on average equalled 3.0 Tg C yr⁻¹ (total of 149 Tg C). As with NEP, the net biome production (NBP) was also highly variable but averaged −0.55 Tg C yr⁻¹ (total of -27.3 Tg C) because NBP in the 1980s was very low (-5.34 Tg C yr⁻¹). During the study period, a total of 126 Tg carbon were removed by logging and land use change, and 50 Tg carbon were directly removed by wildland fires. For carbon pools, the estimated total living upper canopy (tree) biomass decreased from 928 to 834 Tg C, and the understory (including shrub and grass) biomass increased from 59 to 63 Tg C. Soil carbon and dead biomass carbon increased from 1136 to 1197 Tg C.Our analyses suggest that both natural and human processes have significant influence on the carbon change in California. During 1951-2000, climate interannual variability was the key driving force for the large interannual changes of ecosystem carbon source and sink at the state level, while logging and fire were the dominant driving forces for carbon balances in several specific ecoregions. From a long-term perspective, CO₂ fertilization plays a key role in maintaining higher NPP. However, our study shows that the increase in C sequestration by CO₂ fertilization is largely offset by logging/land use change and wildland fires.     

Nitrogen controls on ecosystem carbon sequestration: a model implementation and application to Saskatchewan,Canada

A plant–soil nitrogen (N) cycling model was developed and incorporated into the Integrated BIosphere Simulator (IBIS) of Foley et al. [Foley, J.A., Prentice, I.C., Ramankutty, N., Levis, S., Pollard, D., Sitch, S., Haxeltine, A., 1996. An integrated biosphere model of land surface process, terrestrial carbon balance and vegetation dynamics. Global Biogeochem. Cycles 10, 603–628]. In the N-model, soil mineral N regulates ecosystem carbon (C) fluxes and ecosystem C:N ratios. Net primary productivity (NPP) is controlled by feedbacks from both leaf C:N and soil mineral N. Leaf C:N determines the foliar and canopy photosynthesis rates, while soil mineral N determines the N availability for plant growth and the efficiency of biomass construction. Nitrogen controls on the decomposition of soil organic matter (SOM) are implemented through N immobilization and mineralization separately. The model allows greater SOM mineralization at lower mineral N, and conversely, allows greater N immobilization at higher mineral N. The model's seasonal and inter-annual behaviours are demonstrated. A regional simulation for Saskatchewan, Canada, was performed for the period 1851–2000 at a 10 km × 10 km resolution. Simulated NPP was compared with high-resolution (1 km × 1 km) NPP estimated from remote sensing data using the boreal ecosystem productivity simulator (BEPS) [Liu, J., Chen, J.M., Cihlar, J., Park, W.M., 1997. A process-based boreal ecosystem productivity simulator using remote sensing inputs. Remote Sens. Environ. 44, 81–87]. The agreement between IBIS and BEPS, particularly in NPP spatial variation, was considerably improved when the N controls were introduced into IBIS.     

PIXGRO: A model for simulating the ecosystem CO2 exchange and growth of spring barley

A model, PIXGRO, developed by coupling a canopy flux sub-model (PROXEL_NEE; PROcess-based piXEL Net Ecosystem CO₂ Exchange) to a vegetation structure submodel (CGRO), for simulating both net ecosystem CO₂ exchange (NEE) and growth of spring barley is described. PIXGRO is an extension of the stand-level CO₂ and H₂O-flux model PROXEL_NEE, that simulates the NEE on a process basis, but goes further to include the dry matter production, partitioning, and crop development for spring barley. Dry matter partitioned to the leaf was converted to leaf area index (LAI) using relationships for the specific leaf area (SLA). The canopy flux component, PROXEL_NEE was calibrated using information from the literature on C3 plants and was tested using CO₂ flux data from an eddy-covariance (EC) method in Finland with long-term observations. The growth component (CGRO) was calibrated using data from the literature on spring barley as well as data from the Finland site. It was then validated against field data from two sites in Germany and partly via the use of MODIS remotely sensed LAI from the Finland site.Both the diurnal and the seasonal patterns of gross CO₂ uptake were very well simulated (R² = 0.92). A slight seasonal bias may be attributed to leaf ageing. Crop growth was also well simulated; simulated dry matter agreed with field observed data from Germany (R² = 0.90). For LAI, the agreement between the simulated and observed was good (R² = 0.80), giving an indication that functions describing the conversion of fixed CO₂ to dry matter and the subsequent partitioning leaf dry matter and LAI simulation were robust and provided reliable estimates.The MODIS LAI at a resolution of 1000 m agreed poorly (R² = 0.45) with the PIXGRO simulated LAI and the observed LAI at the Finland site in 2001. We attributed this to the coarse resolution of the image and/or the small size of the barley field (about 17 ha or 0.25 km²) at the Finland site. By deriving a regression relation between the observed LAI and NDVI from a higher resolution MODIS (500 m resolution), the MODIS-recalculated LAI agreed better with the PIXGRO-simulated LAI (R² = 0.86).PIXGRO provides a prototype model bridging the disciplines of plant physiology, crop modeling and remote sensing, for use in a spatial context in evaluating carbon balances and plant growth at stand level, landscape, regional, and with some care, continental scales. Since almost 50% of the European land surface is covered by crops, such a model is needed for the dynamic estimation of LAI and NEE of croplands.     

Estimates of species- and ecosystem-level respiration of woody stems along an elevational gradient in the Luquillo Mountains, Puerto Rico

We measured CO₂ efflux from stems of seven subtropical tree species situated along an elevational gradient in the Luquillo Mountains, Puerto Rico and scaled these measurements up to the landscape level based on modeled and empirical relations. The most important determinants of ecosystem stem respiration were species composition and stem temperature. At a species scale, measured CO₂ efflux per unit bole surface area at a given temperature was highest in the early successional species Cecropia schreberiana and lowest in species that inhabit high elevations such as Micropholis garciniifolia and Tabebuia rigida. Carbon dioxide efflux rates followed a diel pattern that lagged approximately 6 h behind changes in sapwood temperatures. At an ecosystem scale, our simulation model indicates a decreasing trend of stem respiration rates with increasing elevation due to shifts in species composition, lower temperatures and reductions in branch surface area. The highest estimated stem respiration rates were present in the lowland tabonuco forest type and the lowest rates were present in the elfin forest type (mean 7.4 and 2.1 Mg C ha⁻¹ yr⁻¹, respectively). There was slight temperature-induced seasonal variation in simulated stem respiration rates at low elevations, with a maximum difference of 19% between the months of February and July. Our results coincide well with those of Odum and Jordan [Odum, H.T., Jordan, C.F., 1970. Metabolism and evapotranspiration of the lower forest in a giant plastic cylinder. In: Odum, H.T., Pigeon, R.F. (Eds.), A Tropical Rain Forest: A Study of Irradiation and Ecology at El Verde, Puerto Rico. U.S. Atomic Energy Commission, Oak Ridge, TN, pp. I165–I189] for the tabonuco forest type and extend their work by presenting estimates and spatial patterns of woody tissue respiration for the entire mountain rather than for a single forested plot.     

Dynamics of surface albedo of a boreal forest and its simulation

Surface albedo determines the distribution of solar radiation between the earth's surface and the atmosphere. It affects the global climate directly by altering surface energy balance, and indirectly by controlling ecosystem processes and greenhouse gas exchange. In this study, a land surface albedo model was constructed based on the gap probability approach for ray tracing and the basic optical parameters of ecosystem elements. The model was applied to a boreal deciduous forest and results were compared with field measurements. Results show that seasonal and diurnal albedo dynamics were well simulated by the model. The standard deviation between the simulated and measured reflected radiation was 2.5–5.0 W m⁻² in different seasons. The model also provided an insight into the relationships between surface albedo and radiation components (direct versus diffuse), solar zenith angle, and different wave bands. Model sensitivity analyses show that the surface albedo in winter is very sensitive to the forest wood area index for this boreal aspen forest, suggesting that accurate estimates of wood area index are necessary to improve the accuracy of surface albedo simulation in leafless seasons.     

Ecosystem service value losses from coastal erosion in Europe: historical trends and future projections

Coastal zones experience increased rates of coastal erosion, due to rising sea levels, increased storm surge frequencies, reduced sediment delivery and anthropogenic transformations. Yet, coastal zones host ecosystems that provide associated services which, therefore, may be lost due to coastal erosion. In this paper we assess to what extent past and future coastal erosion patterns lead to losses in land cover types and associated ecosystem service values. Hence, historical (based on CORINE land cover information) and projected (based on Dynamic and Interactive Vulnerability Assessment - DIVA - simulations) coastal erosion patterns are used in combination with a benefits transfer approach. DIVA projections are based on regionalized IPCC scenarios. Relative to the period 1975–2050, a case study is provided for selected European coastal country member states. For historical (1975–2006) coastal erosion trends, we observe territory losses in coastal agricultural, water body and forest & semi-natural areas – total coastal erosion equaling over 4,500 km². Corresponding coastal ecosystem service values decrease from about €22.3 billion per year in 1975 to about €21.6 billion per year in 2006. For future (2006–2050) coastal erosion projections, total territory losses equal between ~3,700 km² and ~5,800 km² – coastal wetland areas being affected most severely. Corresponding coastal ecosystem service values decrease to between €20.1 and €19.4 billion per year by 2050. Hence, we argue that the response strategy of EU member states to deal with coastal erosion and climate change impacts should be based on the economic as well as the ecological importance of their coastal zones.     

Modeling seasonal course of carbon fluxes and evapotranspiration in response to low temperature and moisture in a boreal Scots pine ecosystem

Environmental conditions act above and below ground, and regulate carbon fluxes and evapotranspiration. The productivity of boreal forest ecosystems is strongly governed by low temperature and moisture conditions, but the understanding of various feedbacks between vegetation and environmental conditions is still unclear. In order to quantify the seasonal responses of vegetation to environmental factors, the seasonality of carbon and heat fluxes and the corresponding responses for temperature and moisture in air and soil were simulated by merging a process-based model (CoupModel) with detailed measurements representing various components of a forest ecosystem in Hyytiälä, southern Finland. The uncertainties in parameters, model assumptions, and measurements were identified by generalized likelihood uncertainty estimation (GLUE). Seasonal and diurnal courses of sensible and latent heat fluxes and net ecosystem exchange (NEE) of CO₂ were successfully simulated for two contrasting years. Moreover, systematic increases in efficiency of photosynthesis, water uptake, and decomposition occurred from spring to summer, demonstrating the strong coupling between processes. Evapotranspiration and NEE flux both showed a strong response to soil temperature conditions via different direct and indirect ecosystem mechanisms. The rate of photosynthesis was strongly correlated with the corresponding water uptake response and the light use efficiency. With the present data and model assumptions, it was not possible to precisely distinguish the various regulating ecosystem mechanisms. Our approach proved robust for modeling the seasonal course of carbon fluxes and evapotranspiration by combining different independent measurements. It will be highly interesting to continue using long-term series data and to make additional tests of optional stomatal conductance models in order to improve our understanding of the boreal forest ecosystem in response to climate variability and environmental conditions.     

Coupling biochemical and biophysical processes at the leaf level: an equilibrium photosynthesis model for leaves of C3 plants

The paper presents a generic computer model for estimating short-term steady-state fluxes of CO₂, water vapor, and heat from broad leaves and needle-leaved coniferous shoots of C₃ plant species. The model explicitly couples all major processes and feedbacks known to impact leaf biochemistry and biophysics including biochemical reactions, stomatal function, and leaf-boundary layer heat- and mass-transport mechanisms. The ability of the model to successfully predict measured photosynthesis and stomatal-conductance data as well as to simulate a variety of observed leaf responses is demonstrated. A model application investigating physiological and environmental regulation of leaf water-use efficiency (WUE) under steady-state conditions is discussed. Simulation results suggest that leaf physiology has a significant control over the environmental sensitivity of leaf WUE. The implementation of a highly efficient solution technique allows the model to be directly incorporated into plant-canopy and terrestrial ecosystem models.     

Simulating future land use and ecosystem services in Northern Thailand

Enhancing ecosystem services is important as it provides foundation for the wellbeing of people. This paper presents the future land use simulation for enhancing ecosystem services using CLUMondo dynamic spatial model. The land use change was assessed from 1989 to 2013 in Wang Thong watershed of Northern Thailand using GIS and a set of ecosystem services was assessed using the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) model. Future land uses until 2030 were projected for three policy scenarios, namely business-as-usual, integrated land use development, and enhancing environmental services with different levels of emphasis on ecosystem services. In business-as-usual scenario, it was observed that ecosystem services will decline by 2030 from the base year of 2013, whereas in integrated land use development scenario, the ecosystem services will increase by 5% compared to base year due to anticipated effective protection of remaining forests in all existing and proposed protected areas of the study area. In enhancing environmental services scenario, the ecosystem services will increase by 15%. Such results can serve as useful information in policy formulation in developing land use options, which help enhance ecosystem services in future.     

Concentrations of inorganic elements in bottled waters on the Swedish market

This study presents the concentrations of about 50 metals and ions in 33 different brands of bottled waters on the Swedish market. Ten of the brands showed calcium (Ca) concentrations ≤10 mg L⁻¹ and magnesium (Mg) levels <3 mg L⁻¹, implying very soft waters. Three of these waters had in addition low concentrations of sodium (Na; <7 mg L⁻¹), potassium (K; <3 mg L⁻¹) and bicarbonate (HCO₃ ≤31 mg L⁻¹). These brands were collected from barren districts. Nine of the brands were collected from limestone regions. They showed increased Ca-levels exceeding 50 mg L⁻¹ with a maximum of 289 mg L⁻¹. Corresponding Mg-levels were also raised in two brands exceeding 90 mg L⁻¹. Two soft and carbonated waters were supplemented with Na₂CO₃ and NaCl, resulting in high concentrations of Na (644 and 648 mg L⁻¹) and chloride (Cl; 204 and 219 mg L⁻¹). Such waters may make a substantial contribution to the daily intake of NaCl in high water consumers. The storage of carbonated drinking water in aluminum (Al) cans increased the Al-concentration to about 70 μg L⁻¹. Conclusion As there was a large variation in the material as regards concentrations of macro-elements such as Ca, Mg, Na, K and Cl. Supplementation with salts, e.g., Na₂CO₃, K₂ CO₃ and NaCl, can lead to increased concentrations of Na, K and Cl, as well as decreased ratios of Ca/Na and larger ratios of Na/K. Water with high concentrations of e.g., Ca and Mg, may make a substantial contribution to the daily intake of these elements in high water consumers. Al cans are less suited for storage of carbonated waters, as the lowered pH-values may dissolve Al. The levels of potentially toxic metals in the studied brands were generally low.     

Potential impact model to assess agricultural pressure to landscape ecological functions

Conservation and protection of soil and water resources and visual aspects of landscape, as well as the promotion of biodiversity, are some of the central tasks of environmental policy development and social politics in the future. One of the main questions is: ‘which agricultural systems are able to guarantee sustained resource-conserving land use?’ Based on the ecological risk concepts of the 1970s and 1980s, a potential impact model was developed using a universal assessment algorithm derived from fuzzy logic. The model estimated the potential impact of agricultural land use on ecosystem function using a few resource indicators. Intervention intensities of agricultural land-use are set in relation to site conditions and aggregated for each of several defined potential impact categories. The interpretation with respect to risk and the calculation of potential impact values are explained.     

Modelling CO2 and energy exchanges in a northern semiarid grassland using the carbon- and nitrogen-coupled Canadian Land Surface Scheme (C-CLASS)

The development of carbon (C) and nitrogen (N) simulations is one of the ongoing efforts in the land surface schemes of climate models. The C- and N-coupled Canadian Land Surface Scheme (C-CLASS) was recently modified to better represent grassland ecosystems. Improvements include revised plant growth and senescence calculations that are driven by the plant C balance between fixation and respiration, and leaf-out and leaf-fall schemes that are regulated by the seasonal dynamics of C and N reserves. These revisions were developed to better simulate the stress-related senescence and regrowth of perennials. The model was tested with observations of surface carbon and energy fluxes, soil temperature and moisture, and plant growth during 3 years of declining precipitation at a northern semiarid grassland near Lethbridge, Alberta, Canada. The R² and standard deviations between the simulated and observed half-hourly fluxes were 0.95 and 22.5 W m⁻² for net radiation, 0.82 and 42.1 W m⁻² for sensible heat, 0.66 and 29.2 W m⁻² for latent heat, and 0.63 and 0.95 μmol C m⁻² s⁻¹ for net CO₂ exchange. The model and observations both showed a strong impact of declining precipitation on annual carbon budgets in this semi-arid grassland. In a wet year (1998, precipitation = 482 mm), the ecosystem acted as a strong C sink (92 g C m⁻² modelled and 109 g C m⁻² measured from June 20th to December 31st). In a near-normal year (1999, precipitation = 341 mm), a smaller C sink was indicated (24 g C m⁻² modelled and 21 g C m⁻² measured). In a dry year (2000, precipitation = 276 mm), the ecosystem acted as a small C source (−18 g C m⁻² modelled and −17 g C m⁻² measured).     

A process-based model of nitrogen cycling in forest plantations: Part II. Simulating growth and nitrogen mineralisation of Eucalyptus globulus plantations in south-western Australia

We applied a new version of the G’DAY ecosystem model to short-rotation plantations of Eucalyptus globulus growing under a Mediterranean climate in south-western Australia. The new version, that includes modified submodels for biomass production, water balance, litter and soil organic matter (SOM) decomposition, and soil inorganic N balance, was parameterised and applied to three experimental eucalypt sites (Mumballup, Darkan and Northcliffe) of contrasting productivity. With a common base set of parameter values, the model was able to correctly reproduce observed time series of soil water content, canopy leaf area index and stemwood data at the three sites. The model's ability to simulate soil N supply under forest plantations was tested by simulating N mineralisation at each of the three sites over the duration of the experiment (10 years). Simulated annual net N mineralisation in the litter and top 20 cm soil layer ranged from 50 to 170 kg N ha⁻¹ across the sites as a result of differences in rates of litter production, SOM and litter decomposition, and microbial N immobilisation and (re-)mineralisation. Simulations of annual soil N mineralisation were similar to measured rates over a 3-year period, except for an overestimation in 1 year at Mumballup and 2 years at Darkan. Model results indicated the importance of fine root production and turnover for N supply. As plantations age, supply of N to trees increasingly originates from litter decomposition, while the contribution from decomposition of SOM decreases. Although major soil feedbacks associated with litter production, decomposition and N availability are adequately integrated into G’DAY, further work is required in some aspects of the model, including the utility of the C-allocation submodel over a wide range of site conditions and silvicultural treatments.     

Variations in heavy metal contamination of stream water and groundwater affected by an abandoned lead–zinc mine in Korea

This study evaluated variations in heavy metal contamination of stream waters and groundwaters affected by an abandoned lead–zinc mine, where a rockfill dam for water storage will be built 11 km downstream. For these purposes, a total of 10 rounds of stream and groundwater samplings and subsequent chemical analyses were performed during 2002–2003. Results of an exploratory investigation of stream waters in 2000 indicated substantial contamination with heavy metals including zinc (Zn), iron (Fe) and arsenic (As) for at least 6 km downstream from the mine. Stream waters near the mine showed metal contamination as high as arsenic (As) 8,923 μg L⁻¹, copper (Cu) 616 μg L⁻¹, cadmium (Cd) 223 μg L⁻¹ and lead (Pb) 10,590 μg L⁻¹, which greatly exceeded the Korean stream water guidelines. Remediation focused on the mine tailing piles largely improved the stream water qualities. However, there have still been quality problems for the waters containing relatively high concentrations of As (6–174 μg L⁻¹), Cd (1–46 μg L⁻¹) and Pb (2–26 μg L⁻¹). Rainfall infiltration into the mine tailing piles resulted in an increase of heavy metals in the stream waters due to direct discharge of waste effluent, while dilution of the contaminated stream waters improved the water quality due to mixing with metal free rain waters. Levels of As, Cu and chromium (Cr) largely decreased after heavy rain but that of Pb was rather elevated. The stream waters were characterized by high concentrations of calcium (Ca) and sulfate (SO₄), which were derived from dissolution and leaching of carbonate and sulfide minerals. It was observed that the proportions of Ca and SO₄ increased while those of bicarbonate (HCO₃) and sodium and potassium (Na+K) decreased after a light rainfall event. Most interestingly, the reverse was generally detected for the groundwaters. The zinc, being the metal mined, was the most dominant heavy metal in the groundwaters (1758–10,550 μg L⁻¹) near the mine, which far exceeded the Korean standard of 1000 μg L⁻¹ for drinking water. The decreases in the heavy metals contents in the groundwaters associated with reduced rainfall were quite different from the increases observed for the stream waters, which is not clearly understood at this time and warrants further investigation.     

Atmospheric Sulfur Deposition onto Different Ecosystems Over China

The regional acid deposition model system (RegADMS) was applied to simulate the air sulfur deposition onto different landuse types over China, in which the dry deposition velocities of SO₂ and sulfate aerosol (SO ₄ ^2– ) were estimated by use of a big leaf resistance analogy model and the wet scavenging coefficients were parameterized in terms of precipitation rate. Investigations show that the annual total sulfur deposition over mainland China is 7.24 mt (1 mt = 10⁶ ton) , in which dry deposition and wet deposition accounts for 56 and 44%, respectively. The sulfur deposition onto agriculture land, grass land, and forest land is 1.09, 3.6 and 1.41 mt, respectively, which sums 6.1 mt and accounts for 84% of the total sulfur deposition. The modeled sulfur deposition was in agreement with the measurement conducted at farmland in Yingtan, a typical read soil region in Jiangxi province of China, during the period of November 1998–October 1999. The total sulfur deposition at the Yingtan site is about 10.3 gm^–2 year^–1 of which 83% is dry deposition. The modeling sulfur deposition at the same site is 8.4 gm^–2 year^–1. Furthermore, the comparison between RegADMS and RAINS-ASIA on modeling regional sulfur deposition shows the consistence of the two models. The correlation coefficient between the simulated sulfur deposition at the medium-large cities reaches 0.72.