The soil compartment is an important interface between the atmosphere and the subsurface hydrosphere. In this paper a conceptual approach for regional hydrologic soil modelling (RHSM) is presented, which provides two important qualities for modelling. First, the soil compartment is directly coupled to the atmosphere via the land surface and to the aquifers. Second, extremely fine (5 cm vertical) resolutions of the soil system can be realized at regional scales (several hundreds of km2). This high-resolution modelling could be achieved by parallel computation techniques. The RHSM approach is applied to the Beerze-Reusel drainage basin, which belongs to the Meuse River basin. Moisture transport in the soil system was calculated with extremely high vertical resolution at a regional scale based on rainfall-evaporation data for the year 2000. As a result, highly resolved regional groundwater recharge pattern addressing the heterogeneity of soil systems could be determined. A real case application of concept of regional hydrologic soil modelling is presented.

Human activity is changing air quality and climate in the US Pacific Northwest. In a first application of non-metric multidimensional scaling to a large-scale, framework dataset, we modeled lichen community response to air quality and climate gradients at 1416 forested 0.4 ha plots. Model development balanced polluted plots across elevation, forest type and precipitation ranges to isolate pollution response. Air and climate scores were fitted for remaining plots, classed by lichen bioeffects, and mapped. Projected 2040 temperatures would create climate zones with no current analogue. Worst air scores occurred in urban-industrial and agricultural valleys and represented 24% of the landscape. They were correlated with: absence of sensitive lichens, enhancement of nitrophilous lichens, mean wet deposition of ammonium >0.06 mg l-1, lichen nitrogen and sulfur concentrations >0.6% and 0.07%, and SO2 levels harmful to sensitive lichens. The model can detect changes in air quality and climate by scoring re-measurements. Lichen-based air quality and climate gradients in western Oregon and Washington are responsive to regionally increasing nitrogen availability and to temperature changes predicted by climate models.

A Biotic Ligand Model was developed predicting the effect of cobalt on root growth of barley (Hordeum vulgare) in nutrient solutions. The extent to which Ca2+, Mg2+, Na+, K+ ions and pH independently affect cobalt toxicity to barley was studied. With increasing activities of Mg2+, and to a lesser extent also K+, the 4-d EC50Co2+ increased linearly, while Ca2+, Na+ and H+ activities did not affect Co2+ toxicity. Stability constants for the binding of Co2+, Mg2+ and K+ to the biotic ligand were obtained: log KCoBL = 5.14, log KMgBL = 3.86 and log KKBL = 2.50. Limited validation of the model with one standard artificial soil and one standard field soil showed that the 4-d EC50Co2+ could only be predicted within a factor of four from the observed values, indicating further refinement of the BLM is needed. Biotic Ligand Models are not only a useful tool to assess metal toxicity in aquatic systems but can also be used for terrestrial plants.

Residence time-dependent distribution patterns of hexachlorobenzene (HCB) and dichlorodiphenyltrichloroethane (DDT) among different soil organic matter fractions of three Chinese soils were investigated. Soil organic matter (SOM) was fractionated into fulvic acid (FA), humic acid (HA), bound-humic acid (BHA), lipid, and insoluble residue (IR) fractions using methyl isobutyl ketone (MIBK) method. Results revealed that as the residence time prolonged, the amounts of HCB and DDT in the FA, HA and BHA fractions decreased, while those in the lipid and IR fractions increased. One- and two-compartment first order, and one- and two-parameter pore-diffusion kinetic models were used to describe the mobility of HCB and DDT from the FA, HA and BHA fractions. The results suggest that excellent agreements were achieved between the experimental data and fits to the two-compartment first order kinetic model (R2 > 0.97). The transfer rates of HCB and DDT followed the order FA > HA > BHA. HCB and DDT tend to transfer from FA, HA and BHA fractions to lipid and IR fractions with increasing residence time.

The United States and Canada currently use exposure-based metrics to protect vegetation from O3. Using 5 years (1999-2003) of co-measured O3, meteorology and growth response, we have developed exposure-based regression models that predict Populus tremuloides growth change within the North American ambient air quality context. The models comprised growing season fourth-highest daily maximum 8-h average O3 concentration, growing degree days, and wind speed. They had high statistical significance, high goodness of fit, include 95% confidence intervals for tree growth change, and are simple to use. Averaged across a wide range of clonal sensitivity, historical 2001-2003 growth change over most of the 26 M ha P. tremuloides distribution was estimated to have ranged from no impact (0%) to strong negative impacts (-31%). With four aspen clones responding negatively (one responded positively) to O3, the growing season fourth-highest daily maximum 8-h average O3 concentration performed much better than growing season SUM06, AOT40 or maximum 1 h average O3 concentration metrics as a single indicator of aspen stem cross-sectional area growth. A new exposure-based metric approach to predict O3 risk to North American aspen forests has been developed.

Two different indices have been proposed for estimation of the risk caused to forest trees across Europe by ground-level ozone, (i) the concentration based AOT40 index (Accumulated Over a Threshold of 40 ppb) and (ii) the recently developed flux based AFstY index (Accumulated stomatal Flux above a flux threshold Y). This paper compares the AOT40 and AFstY indices for three forest trees species at different locations in Europe. The AFstY index is estimated using the DO3SE (Deposition of Ozone and Stomatal Exchange) model parameterized for Scots pine (Pinus sylvestris), beech (Fagus sylvatica) and holm oak (Quercus ilex). The results show a large difference in the perceived O3 risk when using AOT40 and AFstY indices both between species and regions. The AOT40 index shows a strong north-south gradient across Europe, whereas there is little difference between regions in the modelled values of AFstY. There are significant differences in modelled AFstY between species, which are predominantly determined by differences in the timing and length of the growing season, the periods during which soil moisture deficit limits stomatal conductance, and adaptation to soil moisture stress. This emphasizes the importance of defining species-specific flux response variables to obtain a more accurate quantification of O3 risk. A new flux-based model provides a revised assessment of risks of ozone impacts to European forests.

The level II approach of the critical loads concept adopted by the UNECE aims at a flux based evaluation and takes into account environmental factors governing stomatal conductance. These factors will probably be affected by global change. The flux concept predicts that a decrease in stomatal conductance would protect trees from air pollution effects by decreasing uptake. However, experimental evidence is inconclusive. Numerous results suggest that pollutants and factors subject to global change (drought, CO2) may interact and even exacerbate effects, probably because antioxidative defense systems are involved in both, defense against pollutant effects and protection from natural stress. An effective pollutant dose, which is weighted by physiological defense capacity, would better predict such effects. In this review paper we argue that the flux-based approach is imperfect, because global change effects may also modify the physiological susceptibility to ozone. Instead, a flux concept weighted by defense capacity should be tested. Modeling of ozone effects on plants should include a measure for the plant defense capacity.

The purpose of this study is to develop a mechanistic model to describe a conceptually new “flux-biological response” approach based on biotic ligand model (BLM) and Michaelis-Menten (M-M) kinetics to allow the linkage between valve closure behavior and sodium (Na) transport mechanism in freshwater clam Corbicula fluminea in response to waterborne copper (Cu). We test the proposed model against published data regarding Na uptake kinetics in rainbow trout and Na uptake profile in C. fluminea, confirming that the predictive model is robust. Here, we show that the predicted M-M maximum Cu internalization flux in C. fluminea is 0.369 μmol g-1 h-1 with a half-saturation affinity constant of 7.87 x 10-3 μM. Dynamics of Na uptake and valve closure daily rhythm driven by external Cu can also be predicted simultaneously. We suggest that this “Na transport-valve closure behavior” approach might provide the basis of a future design of biomonitoring tool. A new flux-biological response model can link valve closure and sodium transport mechanisms in freshwater clam in response to copper.

Multimedia mass balance models assume well-mixed homogeneous compartments. Particularly for soils, this does not correspond to reality, which results in potentially large uncertainties in estimates of transport fluxes from soils. A theoretically expected exponential decrease model of chemical concentrations with depth has been proposed, but hardly tested against empirical data. In this paper, we explored the correspondence between theoretically predicted soil concentration profiles and 84 field measured profiles. In most cases, chemical concentrations in soils appear to decline exponentially with depth, and values for the chemical specific soil penetration depth (dp) are predicted within one order of magnitude. Over all, the reliability of multimedia models will improve when they account for depth-dependent soil concentrations, so we recommend to take into account the described theoretical exponential decrease model of chemical concentrations with depth in chemical fate studies. In this model the dp-values should estimated be either based on local conditions or on a fixed dp-value, which we recommend to be 10 cm for chemicals with a log Kow > 3. Multimedia mass model predictions will improve when taking into account depth dependent soil concentrations.

Simple mass balance equations (SMBE) of critical acid loads (CAL) in forest soil were developed to assess potential risks of air pollutants to ecosystems. However, to apply SMBE reliably at large scales, SMBE must be tested for adequacy and uncertainty. Our goal was to provide a detailed analysis of uncertainty in SMBE so that sound strategies for scaling up CAL estimates to the national scale could be developed. Specifically, we wanted to quantify CAL uncertainty under natural variability in 17 model parameters, and determine their relative contributions in predicting CAL. Results indicated that uncertainty in CAL came primarily from components of base cation weathering (BCw; 49%) and acid neutralizing capacity (46%), whereas the most critical parameters were BCw base rate (62%), soil depth (20%), and soil temperature (11%). Thus, improvements in estimates of these factors are crucial to reducing uncertainty and successfully scaling up SMBE for national assessments of CAL.