“The Optimal Steady-State in Groundwater Management,” by Keith C. Knapp and Eli Feinerman2

Knapp and Feinerman (1985) pose and solve a problem of steady-state allocation of ground water based directly on the underlying dynamic problem. Their dynamic steady-state formulation incorporates both the equations of transient ground. water flow and the discount rate. We wish to discuss two aspects of their analysis. First, we question their assertion that the computational advantages of the dynamic steady.state formulation will justify its substitution for the full transient problem. In fact, the dynamic steady-state problem will often be harder to solve than the properly formulated transient problem. Second, we argue that the dynamic steady-state is a concept that has limited applicability in ground-water management. In cases where the optimal steady state is indeed useful, the dynamic solution is often identical to the static solution.     

The evaluation of the equilibrium partitioning method using sensitivity distributions of species in water and soil

The equilibrium partitioning method (EqP-method) can be used to calculate soil quality standards (expressed in mg/kg) from aquatic quality standards (expressed in microg/l) using a partitioning coefficient. The validity of this application of the EqP-method was studied comparing aquatic with terrestrial toxicity data. The data set collected for deriving environmental quality standards in the Netherlands, was used for this study. For 10 organic substances (chlorpyrifos, atrazine, carbofuran, pentachlorophenol, chlordane, aldrin, trichlorobenzene, heptachlor, trichlorophenol and trichloroethene) and for 8 metals, sufficient data were available. The aquatic toxicity data were multiplied by the partitioning coefficient in order to obtain aquatic data expressed in mg/kg. For some compounds the terrestrial toxicity data were significantly higher than the aquatic data but for other compounds it was the other way around. These differences indicate that the EqP-method can give significant over-or underestimations, due to inaccurate partitioning coefficients or differences in species sensitivities. These over- or underestimations can have an impact on the setting of environmental quality standards which are based on the hazardous concentration 5% (HC5) values. The uncertainty in the calculation of HC5 values attributed to the use of the EqP-method, was quantified. The HC5 values derived using the EqP-method were in 5% of the cases more than 20 times higher than the corresponding HC5 values that were derived directly from soil toxicity tests. Despite of this uncertainty the use of the EqP-method can still be advocated for setting soil quality guidelines when only a very limited number of terrestrial toxicity data are available.     

Vapour pressures, aqueous solubilities, Henry's Law constants, and octanol/water partition coefficients of a series of mixed halogenated dimethyl bipyrroles

Basic physical-chemical properties of five bromine and chlorine containing mixed halogenated dimethyl bipyrroles (HDBPs) were determined using established methods. Subcooled liquid vapour pressures (P(o)(L,25)), aqueous solubilities (S(w,25)), and octanol/water partition coefficients (K(ow)) were determined using the gas chromatography-retention time, generator column, and slow-stirring methods, respectively. Henry's Law constants (H25) were estimated using experimentally-derived P(o)(L) and S(w,25) data. Values of all four properties were generally similar to those reported for other polyhalogenated aromatic compounds [P(o)(L,25) = (7.55-191) x 10(-6) Pa; S(w,25) = (1.0-1.9) x 10(-5) g/l; log K(ow) = 6.4-6.7; H25 = 0.0020-0.14 Pa m3/mol]. The effect of replacing a chlorine with a bromine atom significantly decreased P(o)(L,25) (log P(o)(L,25) = -0.4197 (# bromine atoms) - 2.643, p<0.01) and H25 (log H25 = -0.508 (# bromine atoms) + 0.394, p<0.02). There were no significant effects of bromine/chlorine substitution on S(w,25) or K(ow). A simple Level I equilibrium partitioning model predicted the environmental behaviour of HDBPs to be similar to a tetrabrominated diphenyl ether. Only slight differences in behaviour amongst HDBP congeners were predicted since substitution of a bromine for a chlorine (Cl/Br substitution) atom had less effect than H/Cl or H/Br substitution on P(o)(L,25), S(w,25), H25, and K(ow).     

An assessment of the emissions inventory processing systems EMS-2001 and SMOKE in grid-based air quality models

In the United States, emission processing models such as Emissions Modeling System-2001 (EMS-2001), Emissions Preprocessor System-Version 2.5 (EPS2.5), and the Sparse Matrix Operator Kernel Emissions (SMOKE) model are currently being used to generate gridded, hourly, speciated emission inputs for urban and regional-scale photochemical models from aggregated pollutant inventories. In this study, two models, EMS-2001 and SMOKE, were applied with their default internal data sets to process a common inventory database for a high ozone (O3) episode over the eastern United States using the Carbon Bond IV (CB4) chemical speciation mechanism. A comparison of the emissions processed by these systems shows differences in all three of the major processing steps performed by the two models (i.e., in temporal allocation, spatial allocation, and chemical speciation). Results from a simulation with a photochemical model using these two sets of emissions indicate differences on the order of +/- 20 ppb in the predicted 1-hr daily maximum O3 concentrations. It is therefore critical to develop and implement more common and synchronized temporal, spatial, and speciation cross-reference systems such that the processes within each emissions model converge toward reasonably similar results. This would also help to increase confidence in the validity of photochemical grid model results by reducing one aspect of modeling uncertainty.     

The OECD Validation Program of the H295R Steroidogenesis Assay for the Identification of In Vitro Inhibitors and Inducers of Testosterone and Estradiol Production. Phase 2: Inter-Laboratory Pre-Validation Studies (8 pp)

Background, Goals and Scope In response to concerns that have been raised about chemical substances that may alter the function of endocrine systems and result in adverse effects on human health, an OECD initiative was undertaken to develop and validate in vitro and in vivo assays to identify chemicals that may interfere with endocrine systems of vertebrates. Here we report on studies that were conducted to develop and standardize a cell-based screening assay using the H295R cell line to prioritize chemicals that may act on steroidogenic processes in humans and wildlife. These studies are currently ongoing as part of the ‘Special Activity on the Testing and Assessment of Endocrine Disruptors’ within the OECD Test Guidelines Program to review, develop, standardize, and validate a number of in vitro and in vivo toxicological assays for testing and assessment of chemicals concerning their potential to interact with the endocrine system of vertebrates. Study Design Six laboratories from five countries participated in the pre-validation studies. Each laboratory tested the effects of three model chemicals on the production of testosterone (T) and estradiol (E2) using the H295R Steroidogenesis Assay. Chemicals tested were well described inducers or inhibitors of steroidogenic pathways (forskolin, prochloraz and fadrozole). All experiments were conducted in 24 well plates following standard protocols. Six different doses per compound were analyzed in triplicate per plate. A quality control (QC) plate was run in conjunction with the chemical exposure plate to account for inter-assay variation. Each chemical exposure was conducted two or three times. Results All laboratories successfully detected increases and/or decreases in hormone production by H295R cells after exposure to the different model compounds and there was good agreement in the pattern of response for all groups. Forskolin increased both T and E2 while fadrozole and prochloraz decreased production of both hormones. All chemicals affected hormone production in a dose-dependent manner with the exception of fadrozole which caused maximum inhibition of E2 at the two least concentrations tested. Some inter-laboratory differences were noted in the alteration of hormone production measured in chemically exposed cells. However, with the exception of the production of T measured at one laboratory in cells exposed to forskolin, the EC₅₀s calculated were comparable (coefficients of variation 34–49%) for all hormones. Discussion and Perspectives The results indicated that the H295R Steroidogenesis Assay protocol was robust, transferable and reproducible among all laboratories. However, in several instances that were primarily related to one laboratory there were unexplained minor uncertainties related to the inter-laboratory hormone production variation. Based on the findings from this Phase 2 prevalidation study, the H295R Steroidogenesis Assay protocol is currently being refined. The next phase of the OECD validation program will test the refined protocol among the same group of laboratories using an extended set of chemicals (∼30) that will include positive and negative chemical controls as well as a broad spectrum of different potential inducers and inhibitors of steroidogenic pathways. Submission Editor: Dr. Carsten Brühl (bruehl@uni-landau.de)     

Metal behaviour in an estuary polluted by acid mine drainage: the role of particulate matter

The concentrations of dissolved and suspended particulate Cd, Cu and Zn have been determined in water samples obtained during two axial transects of the Rio Tinto-Huelva Ria system in south-west Spain, which is severely impacted by acid mine drainage. Although the metal concentrations in both phases were elevated, dissolved metals were dominant and, in the upper estuary, constituted > 99% of total metal in the water column. Dissolved metals behaved non-conservatively on each transect, with maximum concentrations in the low salinity region. There was no evidence of metal adsorption within the turbidity maximum zone, despite the high specific surface areas of resuspending particles. Measurements of electrophoretic mobility showed that the suspended particulate matter (SPM) had a positive surface charge in the salinity range 0-4, where the waters had a pH < 3. Desorption experiments were carried out in which SPM from the turbidity maximum zone was resuspended in coastal seawater. The desorption of the metals was monitored for 24 h, using anodic stripping voltammetry (ASV) to detect the variation in total dissolved Cd, Cu and Zn and the species of Cu and Zn. Total dissolved Cd concentrations doubled during the incubation period, whereas the concentration of total dissolved Cu declined and that of Zn remained rather constant. The ASV-labile fraction of dissolved Cu and Zn showed an initial sharp release followed by a slower uptake. However, desorption was shown to be a minor source of dissolved metals and made little contribution to the non-conservative behaviour in the low salinity zone. The results are used to predict the effects of acid mine drainage on estuarine ecology.     

Environmental impact of two successive chemical treatments in a small shallow eutrophied lake: Part II. Case of copper sulfate

The appearance of cyanobacteria ( > 10 colony per ml) was not prevented after alum treatment. In order to prevent cyanobacteria efflorescences in a small shallow polymictic lake (Courtille, France), copper sulfate was applied. Treatment level was 63 microg 1(-1) as Cu2+ from CUSO4, 5 H2O. Cyanobacteria were kept under control during the summer. Microcystis sp. completely disappeared, which allowed swimming in the lake throughout the tourist season. Microcystis only reappeared 2 months after the treatment. Copper content in the water column only returned to its background level 2 months after copper addition. This high residence time of copper in the water might have been caused by complexation and adsorption of copper on natural organic matter, whose level was high in the ecosystem studied. A mechanism of transfer of 'truly' dissolved copper towards particulate copper has been underlined and explains the disappearance of this fraction of copper in the water column.     

Optimal patch time allocation for time-limited foragers

The Charnov Marginal Value Theorem (MVT) predicts the optimal foraging duration of animals exploiting patches of resources. The predictions of this model have been verified for various animal species. However, the model is based on several assumptions that are likely too simplistic. One of these assumptions is that animals are living forever (i.e., infinite horizon). Using a simple dynamic programming model, we tested the importance of this assumption by analysing the optimal strategy for time-limited foragers. We found that, for time-limited foragers, optimal patch residence times should be greater than those predicted from the classic, static MVT, and the deviation should increase when foragers are approaching the end of their life. These predictions were verified for females of the parasitoid Anaphes victus (Hymenoptera: Mymaridae) exploiting egg patches of its host, the carrot weevil Listronotus oregonensis (Coleoptera: Curculionidae). As predicted by the model, females indeed remained for a longer time on host patches when they approached the end of their life. Experimental results were finally analysed with a Cox regression model to identify the patch-leaving decision rules females used to behave according to the model’s predictions.     

The Residence Time of Settling Particles in the Surface Mixed Layer

The transport from the upper mixed layer into the pycnocline of particles with negative buoyancy is considered. Assuming the hydrodynamic parameters to be time- independent, an adjoint model is resorted to that provides a general expression of the residence time in the mixed layer of the constituent under study. It is seen that the residence time decreases as the settling velocity increases or the diffusivity decreases. Furthermore, it is demonstrated that the residence time must be larger than z/w and smaller than h/w, where z, h and w denote the distance to the pycnocline, the thickness of the mixed layer and the sinking velocity. In the vicinity of the pycnocline, the residence time is not necessarily zero; its behaviour critically depends on the eddy diffusivity profile in this region. Closed-form solutions are obtained for constant and quadratic diffusivity profiles, which allows for an analysis of the sensitivity of the residence time to the Peclet number. Finally, an approximate value is suggested of the depth-averaged value of the residence time.