Comparison of models of simazine transport and fate in the subsurface environment in a citrus farm

Contamination of groundwater by agrochemicals is now widely recognized as an extremely important environmental problem. Modern agricultural practices involve the combined use of irrigation with the application of large amounts of agrochemicals to maximize crop yield. Due to flood irrigation and natural runoff, agricultural activities might generate soil, surface water and groundwater contamination problems and leaching of pesticides. Modeling of the transport and fate of pesticides, such as simazine, may help understand the long-term potential risk to the subsurface environment. This paper illustrates a comparative study via the use of three different pesticide transport simulation models and the applicability of those models in determining the groundwater vulnerability to pesticides contamination in a citrus orchard located at the Lower Rio Grande Valley (LRGV). The three models used in the study are the pesticide root zone model-3 (PRZM-3), the pesticide analytical model (PESTAN) and integrated pesticide transport modeling (IPTM). The concentration values obtained from all three models are in agreement, and they show a decreasing trend from the surface through the vadose zone. The problem is how to use this information and, specifically, how to combine the testimony of a number of experts into a single useful judgment. With the aid of the fuzzy multiattribute decision making method, PRZM-3 is deemed as the most promising one for such precision farming applications.     

Initial Changes in Refilled Lysimeters Built with Metal Polluted Topsoil and Acidic or Calcareous Subsoils as Indicated by Changes in Drainage Water Composition

Soil translocation for recultivation of soil removed from construction sites and for the preparation of refilled lysimeters inevitably involves disturbance of soil structure, and, if intermediate storage is included, also drying and rewetting of the soil. We report on an experiment with model forest ecosystems, where uncontaminated forest subsoils were covered with non-contaminated or freshly heavy metal (mainly Zn and Cu) contaminated topsoil in large lysimeters. Monitoring of the chemical composition of the drainage water revealed two distinct soil conditioning phases. During an initial phase of about a year strongly elevated nitrate and sulfate concentrations occurred that were attributed to a mineralisation flush caused by the increased accessability of mineralisable nitrogen and sulfur in destroyed aggregates. These effects were significantly larger in lysimeters with calcareous subsoil than in those with acidic subsoil. The second phase was characterised by a gradual decrease in dissolved organic carbon and sulfate concentrations, in particular in the acidic subsoil. This decrease may be attributed to the depletion of pools made accessible during aggregate destruction or the formation of new aggregates. These chemical changes had only little effects on the concentrations of copper and zinc in the drainage water. Based on our results, it can be concluded that large refilled lysimeters can be used for many purposes without risk of compromised results, if a conditioning phase of about 1 year with sufficiently moist soil conditions is respected. Nevertheless, gradual changes in soil chemical characteristics still occur after this initial phase. Implications for the recultivation of sites using relocated soils are discussed.     

Plant and Soil System Responses to Ozone After 3 Years in a Lysimeter Study with Juvenile Beech (Fagus sylvatica L.)

A lysimeter study was performed to monitor effects of elevated ozone on juvenile trees of Fagus sylvatica L. as well as on the plant–soil system. During a fumigation period over almost three growing seasons, parameters related to plant growth, phenological development and physiology as well as soil functions were studied. The data analyses identified elevated ozone to delay leaf phenology at early and to accelerate it at late developmental stages, to reduce growth, some leaf nutrients (Ca, K) as well as some soluble phenolics (hydroxycinnamic acid derivatives, total flavonol glycosides). No or very weak ozone effects were found in mobile carbon pools of leaves (starch, sucrose), and other phenolic compounds (flavans). Altered gene expression related to stress and carbon cycling corresponded well with findings from leaf phenology and chemical composition analyses indicating earlier senescence and oxidative stress in leaves under elevated ozone. Conversely in the soil system, no effects of ozone were detected on soil enzyme activities, rates of litter degradation and lysimeter water balances. Despite the fact that the three reported years 2003–2005 were climatically very contrasting including a hot and dry as well as an extremely wet summer, and also mild as well as cold winters, the influence of ozone on a number of plant parameters is remarkably consistent, further underlining the phytotoxic potential of elevated tropospheric ozone levels.     

Substance Flow Analyses of Organic Pollutants in Stockholm

This paper summarizes substance flow analyses for four organic substances in the City of Stockholm, Sweden: diethylhexyl phthalate (DEHP), alkylphenolethoxylates (APEO), polybrominated diphenylethers (PBDE) and chlorinated paraffins (CP). The results indicate that the stocks of APEO, PBDE and CP all are approximately 200–250 tonnes, whereas the DEHP stock is two orders of magnitude larger. Emissions can be linked to imported consumer goods such as electronics (PBDE) and textiles (APEO), and to construction materials (DEHP, CP). For several of the substances considerable amounts remain in the technosphere for a long time, even after use of the substance in new products has been eliminated. For example, the use of DEHP as plasticizer for PVC plastics in cables and floorings has more or less been phased-out, but still these applications make up a stock of some 20,000 tonnes (85% of the total DEHP stock in Stockholm) and emit 28 tonnes of DEHP annually (93% of overall emissions). Likewise, the use of chlorinated paraffins in sealants has been radically reduced, but there are 170 tonnes of CP in sealants in Stockholm making up 75% of the stock, and causing half of the emissions to water and air. These emissions are likely to continue for decades, and the stocks therefore have to be considered when analysing and managing the impact of urban substance flows on the environment.     

Sustainability, substance flow management and time. Part I Temporal analysis of substance flows

Flows of chemical substances need to be managed in a sustainable way. Sustainable development as a whole and the sustainable management of substance flows in particular are both time issues. These include the importance of the dynamics of substance flows and the way these interconnect with the use of resources, the avoidance of environmental pollution, and their effects on health and food production. Another prerequisite for the proper management of substance flows is justice within and between generations. This requires a systematic approach and a systematic analysis of the issues as well as of the actions to be taken. One tool for such a systematic approach is temporal analysis. It brings the temporal aspects of the substances themselves and of their intended use, as well as factors affecting the stakeholders, such as decision makers, producers and consumers, into focus. In the past, timing factors were rarely taken into account. Knowledge of the temporal dynamics of substance flows and their resultant outcomes, as well as of their interaction with ecological, economic and social systems, is a basic requirement for successful substance flow management. The need to include temporal aspects into substance flow management and how to do so is outlined here. Included are not only politicians but also practitioners and scientists who must explicitly take into account adequate time scales, points in time, breaks and other forms of time in planning and acting.     

Modeling the Effects of Potential Salinity Shifts on the Recovery of Striped Bass in the Savannah River Estuary,Georgia–South Carolina,United States

Increased salinity in spawning and nursery grounds in the Savannah River estuary was cited as the primary cause of a 97% decrease in adult striped bass (Morone saxatilis) and a concomitant 96% decrease in striped bass egg production. Restoration efforts focused on environmental remediation and stock enhancement have resulted in restored salinity patterns and increased egg and adult abundances. However, future water needs or harbor development may preclude further recovery by reducing freshwater inflow or increasing salinity intrusion. To assess the effect of potential changes in the salinity regime, we developed models relating discharge, tidal phase, and salinity to striped bass egg and early larval survival and re-cast these in a quantitative Bayesian belief network. The model indicated that a small upstream shift (≤1.67 km) in the salinity regime would have the least impact on striped bass early life history survival, whereas shifts >1.67 km would have progressively larger impacts, with a 8.33-km shift potentially reducing our estimated survival probability by >28%. Such an impact could have cumulative and long-term detrimental effects on the recovery of the Savannah River striped bass population. The available salinity data were collected during average and low flows, so our model represents some typical and some extreme conditions during a striped bass spawning season. Our model is a relatively simplistic, “first-order” attempt at evaluating potential effects of changes in the Savannah River estuarine salinity regime and points to areas of concern and potential future research.     

A framework for valuing the health benefits of improved bathing water quality in the River Irvine catchment

A simple model predicting bathing water concentrations of Escherichia coli from livestock in the Irvine catchment in SW Scotland has been adapted for intestinal enterococci (IE). This has been used to predict risk of bather illness by extrapolation of published data on bather IE exposure vs incidence of gastro-enteritis. Simulated reduction in the risk of illness by reduced faecal loading was multiplied by a willingness to pay for risk reduction to estimate the annual benefits of mitigation. Health benefits of reducing loading by 75% at Irvine Beach were estimated by a willingness to pay method to be about pound 276k pa. Estimated annualised costs of diffuse pollution mitigation measures across the catchment were higher (> pound 1m), and it is very unlikely that 75% mitigation is achievable with current stocking rates. Further work should explore the influence of uncertainty of model parameters, and use emerging epidemiological information on specific zoonotic pathogens such as E. coli O157 and Cryptosporidium. Other components of the value of clean water should also be included to obtain a complete estimate of the cost:benefit of mitigation.     

Identification of key parameters controlling dissolved oxygen migration and attenuation in fractured crystalline rocks

In the crystalline rocks of the Canadian Shield, geochemical conditions are currently reducing at depths of 500-1000 m. However, during future glacial periods, altered hydrologic conditions could potentially result in enhanced recharge of glacial melt water containing a relatively high concentration of dissolved oxygen (O2). It is therefore of interest to investigate the physical and geochemical processes, including naturally-occurring redox reactions, that may control O2 ingress. In this study, the reactive transport code MIN3P is used in combination with 2k factorial analyses to identify the most important parameters controlling oxygen migration and attenuation in fractured crystalline rocks. Scenarios considered are based on simplified conceptual models that include a single vertical fracture, or a fracture zone, contained within a rock matrix that extends from the ground surface to a depth of 500 m. Consistent with field observations, Fe(II)-bearing minerals are present in the fractures (i.e. chlorite) and the rock matrix (biotite and small quantities of pyrite). For the parameter ranges investigated, results indicate that for the single fracture case, the most influential factors controlling dissolved O2 ingress are flow velocity in the fracture, fracture aperture, and the biotite reaction rate in the rock matrix. The most important parameters for the fracture zone simulations are flow velocity in the individual fractures, pO2 in the recharge water, biotite reaction rate, and to a lesser degree the abundance and reactivity of chlorite in the fracture zone, and the fracture zone width. These parameters should therefore receive increased consideration during site characterization, and in the formulation of site-specific models intended to predict O2 behavior in crystalline rocks.     

Downward Migration of Dense Bottom Currents

The focus of this paper is on the dynamics of a dense current flowing along the continental slope, and the frictionally induced downward motion it experiences. In particular, the movement of the lateral boundaries where the isopycnals meet the bottom are considered. The current is taken to be wide compared to the Rossby radius, which is in accordance with observations and makes the dynamics quasigeostrophic. The time development of the plume thickness is studied, using three different parameterisations for the bottom friction.Independently of the choice of parameterisation, the following results are obtained: In the central part of the plume friction acts as a diffusive process to minimise the curvature of the upper surface of the plume. At the upper edge the plume quickly approaches a state with small slope, i.e. small geostrophic velocity and small frictionally induced downward flow. At the lower edge a tongue of dense water shoots out creating a widening layer. The thickness of the migrating layer is approximately equal to the boundary layer depth and its downward speed is comparable to the along-slope geostrophic velocity. The downslope end of the migrating layer may form a steep front which requires some precautions in the numerical procedure.