Neurobehavioral Development in Children With Potential Exposure to Pesticides

Children may be at higher risk than adults from pesticide exposure, due to their rapidly developing physiology, unique behavioral patterns, and interactions with the physical environment. This preliminary study conducted in Ecuador examines the association between household and environmental risk factors for pesticide exposure and neurobehavioral development. We collected data over 6 months in the rural highland region of Cayambe, Ecuador (2003–2004). Children age 24–61 months residing in 3 communities were assessed with the Ages and Stages Questionnaire and the Visual Motor Integration Test. We gathered information on maternal health and work characteristics, the home and community environment, and child characteristics. Growth measurements and a hemoglobin finger-prick blood test were obtained. Multiple linear regression analyses were conducted. Current maternal employment in the flower industry was associated with better developmental scores. Longer hours playing outdoors were associated with lower gross and fine motor and problem solving skills. Children who played with irrigation water scored lower on fine motor skills (8% decrease; 95% confidence interval 9.31 to 0.53), problem-solving skills (7% decrease; 8.40 to 0.39), and Visual Motor Integration test scores (3% decrease; 12.00 to 1.08). These results suggest that certain environmental risk factors for exposure to pesticides may affect child development, with contact with irrigation water of particular concern. However, the relationships between these risk factors and social characteristics are complex, as corporate agriculture may increase risk through pesticide exposure and environmental contamination, while indirectly promoting healthy development by providing health care, relatively higher salaries, and daycare options.

Characteristics of Langmuir Turbulence in the Ocean Mixed Layer

This study uses large-eddy simulation (LES) to investigate the characteristics of Langmuir turbulence through the turbulent kinetic energy (TKE) budget. Based on an analysis of the TKE budget a velocity scale for Langmuir turbulence is proposed. The velocity scale depends on both the friction velocity and the surface Stokes drift associated with the wave field. The scaling leads to unique profiles of nondimensional dissipation rate and velocity component variances when the Stokes drift of the wave field is sufficiently large compared to the surface friction velocity. The existence of such a scaling shows that Langmuir turbulence can be considered as a turbulence regime in its own right, rather than a modification of shear-driven turbulence. Comparisons are made between the LES results and observations, but the lack of information concerning the wave field means these are mainly restricted to comparing profile shapes. The shapes of the LES profiles are consistent with observed profiles. The dissipation length scale for Langmuir turbulence is found to be similar to the dissipation length scale in the shear-driven boundary layer. Beyond this it is not possible to test the proposed scaling directly using available data. Entrainment at the base of the mixed layer is shown to be significantly enhanced over that due to normal shear turbulence.

Impact of Resolution on the Tropical Pacific Circulation in a Matrix of Coupled Models

Results are presented from a matrix of coupled model integrations, using atmosphere resolutions of 135 and 90 km, and ocean resolutions of 1° and 1/3°, to study the impact of resolution on simulated climate. The mean state of the tropical Pacific is found to be improved in the models with a higher ocean resolution. Such an improved mean state arises from the development of tropical instability waves, which are poorly resolved at low resolution; these waves reduce the equatorial cold tongue bias. The improved ocean state also allows for a better simulation of the atmospheric Walker circulation. Several sensitivity studies have been performed to further understand the processes involved in the different component models. Significantly decreasing the horizontal momentum dissipation in the coupled model with the lower-resolution ocean has benefits for the mean tropical Pacific climate, but decreases model stability. Increasing the momentum dissipation in the coupled model with the higher-resolution ocean degrades the simulation toward that of the lower-resolution ocean. These results suggest that enhanced ocean model resolution can have important benefits for the climatology of both the atmosphere and ocean components of the coupled model, and that some of these benefits may be achievable at lower ocean resolution, if the model formulation allows.

Health Benefits of Organic Food: effects of the environment

Public concern over impacts of chemicals in plant and animal production on health and the environment has led to increased demand for organic produce, which is usually promoted and often perceived as containing fewer contaminants, more nutrients, and being positive for the environment. These benefits are difficult to quantify, and potential environmental impacts on such benefits have not been widely studied. This book addresses these key points, examining factors such as the role of certain nutrients in prevention and promotion of chronic disease, potential health benefits of bioactive compounds in plants, the prevalence of food-borne pesticides and pathogens and how both local and global environmental factors may affect any differences between organic and conventionally produced food. This book is an essential resource for researchers and students in human health and nutrition, environmental science, agriculture and organic farming. Main Contents 1. Organic farming and food systems: definitions and key characteristics. 2. The health benefits of n-3 fatty acids and their concentrations in organic and conventional animal-derived foods. 3. Environmental impacts on n-3 content of foods from ruminant animals. 4. Health benefits and selenium content of organic vs conventional foods. 5. Environmental impacts concerning the selenium content of foods. 6. Contaminants in organic and conventional food: the missing link between contaminant levels and health effects. 7. Mycotoxins in organic and conventional foods and effects of the environment. 8. Human pathogens in organic and conventional foods and effects of the environment. 9. What does consumer science tell us about organic foods? 10. The beneficial effects of dietary flavonoids: sources, bioavailability and biological functions. 11. Environmental regulation of flavonoid biosynthesis. 12. Nitrates in the human diet. 13. Impacts of environment and management on nitrate in vegetables and water. 14. Effects of the environment on the nutritional quality and safety of organically produced foods: Round-up and summary.

Measuring and modelling mixture toxicity of imidacloprid and thiacloprid on Caenorhabditis elegans and Eisenia fetida

While the standard models of concentration addition and independent action predict overall toxicity of multicomponent mixtures reasonably, interactions may limit the predictive capability when a few compounds dominate a mixture. This study was conducted to test if statistically significant systematic deviations from concentration addition (i.e. synergism/antagonism, dose ratio- or dose level-dependency) occur when two taxonomically unrelated species, the earthworm Eisenia fetida and the nematode Caenorhabditis elegans were exposed to a full range of mixtures of the similar acting neonicotinoid pesticides imidacloprid and thiacloprid. The effect of the mixtures on C. elegans was described significantly better (p<0.01) by a dose level-dependent deviation from the concentration addition model than by the reference model alone, while the reference model description of the effects on E. fetida could not be significantly improved. These results highlight that deviations from concentration addition are possible even with similar acting compounds, but that the nature of such deviations are species dependent. For improving ecological risk assessment of simple mixtures, this implies that the concentration addition model may need to be used in a probabilistic context, rather than in its traditional deterministic manner. Crown Copyright (C) 2008 Published by Elsevier Inc. All rights reserved.

The analysis of ranked observations of soil structure using indicator geostatistics

Structure is an important physical feature of the soil that is associated with water movement, the soil atmosphere, microorganism activity and nutrient uptake. A soil without any obvious organisation of its components is known as apedal and this state can have marked effects on several soil processes. Accurate maps of topsoil and subsoil structure are desirable for a wide range of models that aim to predict erosion, solute transport, or flow of water through the soil. Also such maps would be useful to precision farmers when deciding how to apply nutrients and pesticides in a site-specific way, and to target subsoiling and soil structure stabilization procedures. Typically, soil structure is inferred from bulk density or penetrometer resistance measurements and more recently from soil resistivity and conductivity surveys. To measure the former is both time-consuming and costly, whereas observations made by the latter methods can be made automatically and swiftly using a vehicle-mounted penetrometer or resistivity and conductivity sensors. The results of each of these methods, however, are affected by other soil properties, in particular moisture content at the time of sampling, texture, and the presence of stones. Traditional methods of observing soil structure identify the type of ped and its degree of development. Methods of ranking such observations from good to poor for different soil textures have been developed. Indicator variograms can be computed for each category or rank of structure and these can be summed to give the sum of indicator variograms (SIV). Observations of the topsoil and subsoil structure were made at four field sites where the soil had developed on different parent materials. The observations were ranked by four methods and indicator and the sum of indicator variograms were computed and modelled for each method of ranking. The individual indicators were then kriged with the parameters of the appropriate indicator variogram model to map the probability of encountering soil with the structure represented by that indicator. The model parameters of the SIVs for each ranking system were used with the data to krige the soil structure classes, and the results are compared with those for the individual indicators. The relations between maps of soil structure and selected wavebands from aerial photographs are examined as basis for planning surveys of soil structure. (C) 2007 Elsevier B.V. All rights reserved.

MEASUREMENT AND MODELING OF THE TOXICITY OF BINARY MIXTURES IN THE NEMATODE CAENORHABDITIS ELEGANS-A TEST OF INDEPENDENT ACTION

Ecological risk assessments must increasingly consider the effects of chemical mixtures on the environment as anthropogenic pollution continues to grow in complexity. Yet testing every possible mixture combination is impractical and unfeasible; thus, there is an urgent need for models that can accurately predict mixture toxicity from single-compound data. Currently, two models are frequently used to predict mixture toxicity from single-compound data: Concentration addition and independent action (IA). The accuracy of the predictions generated by these models is currently debated and needs to be resolved before their use in risk assessments can be fully justified. The present study addresses this issue by determining whether the IA model adequately described the toxicity of binary mixtures of five pesticides and other environmental contaminants (cadmium, chlorpyrifos, diuron, nickel, and prochloraz) each with dissimilar modes of action on the reproduction of the nematode Caenorhabditis elegans. In three out of 10 cases, the IA model failed to describe mixture toxicity adequately with significant or antagonism being observed. In a further three cases, there was an indication of synergy, antagonism, and effect-level-dependent deviations, respectively, but these were not statistically significant. The extent of the significant deviations that were found varied, but all were such that the predicted percentage effect seen on reproductive output would have been wrong by 18 to 35% (i.e., the effect concentration expected to cause a 50% effect led to an 85% effect). The presence of such a high number and variety of deviations has important implications for the use of existing mixture toxicity models for risk assessments, especially where all or part of the deviation is synergistic.

Doppler lidar measurements of oriented planar ice crystals falling from supercooled and glaciated layer clouds

The properties of planar ice crystals settling horizontally have been investigated using a vertically pointing Doppler lidar. Strong specular reflections were observed from their oriented basal facets, identified by comparison with a second lidar pointing 4° from zenith. Analysis of 17 months of continuous high-resolution observations reveals that these pristine crystals are frequently observed in ice falling from mid-level mixed-phase layer clouds (85% of the time for layers at −15 °C). Detailed analysis of a case study indicates that the crystals are nucleated and grow rapidly within the supercooled layer, then fall out, forming well-defined layers of specular reflection. From the lidar alone the fraction of oriented crystals cannot be quantified, but polarimetric radar measurements confirmed that a substantial fraction of the crystal population was well oriented. As the crystals fall into subsaturated air, specular reflection is observed to switch off as the crystal faces become rounded and lose their faceted structure. Specular reflection in ice falling from supercooled layers colder than −22 °C was also observed, but this was much less pronounced than at warmer temperatures: we suggest that in cold clouds it is the small droplets in the distribution that freeze into plates and produce specular reflection, whilst larger droplets freeze into complex polycrystals. The lidar Doppler measurements show that typical fall speeds for the oriented crystals are ≈ 0.3 m s−1, with a weak temperature correlation; the corresponding Reynolds number is Re ∼ 10, in agreement with light-pillar measurements. Coincident Doppler radar observations show no correlation between the specular enhancement and the eddy dissipation rate, indicating that turbulence does not control crystal orientation in these clouds. Copyright © 2010 Royal Meteorological Society

Evaluation of the model representation of the evolution of convective systems using satellite observations of outgoing longwave radiation

We introduce a technique for assessing the diurnal development of convective storm systems based on outgoing longwave radiation fields. Using the size distribution of the storms measured from a series of images, we generate an array in the lengthscale-time domain based on the standard score statistic. It demonstrates succinctly the size evolution of storms as well as the dissipation kinematics. It also provides evidence related to the temperature evolution of the cloud tops. We apply this approach to a test case comparing observations made by the Geostationary Earth Radiation Budget instrument to output from the Met Office Unified Model run at two resolutions. The 12km resolution model produces peak convective activity on all lengthscales significantly earlier in the day than shown by the observations and no evidence for storms growing in size. The 4km resolution model shows realistic timing and growth evolution although the dissipation mechanism still differs from the observed data.

Numerical study of 2D heat transfer in a scraped surface heat exchanger

A numerical study of fluid mechanics and heat transfer in a scraped surface heat exchanger with non-Newtonian power law fluids is undertaken. Numerical results are generated for 2D steady-state conditions using finite element methods. The effect of blade design and material properties, and especially the independent effects of shear thinning and heat thinning on the flow and heat transfer, are studied. The results show that the gaps at the root of the blades, where the blades are connected to the inner cylinder, remove the stagnation points, reduce the net force on the blades and shift the location of the central stagnation point. The shear thinning property of the fluid reduces the local viscous dissipation close to the singularity corners, i.e. near the tip of the blades, and as a result the local fluid temperature is regulated. The heat thinning effect is greatest for Newtonian fluids where the viscous dissipation and the local temperature are highest at the tip of the blades. Where comparison is possible, very good agreement is found between the numerical results and the available data. Aspects of scraped surface heat exchanger design are assessed in the light of the results. (C) 2003 Elsevier Ltd. All rights reserved.

Gene transcription in Daphnia magna: effects of acute exposure to a carbamate insecticide and an acetanilide herbicide.

Daphnia magna is a key invertebrate in the freshwater environment and is used widely as a model in ecotoxicological measurements and risk assessment. Understanding the genomic responses of D. magna to chemical challenges will be of value to regulatory authorities worldwide. Here we exposed D. magna to the insecticide methomyl and the herbicide propanil to compare phenotypic effects with changes in mRNA expression levels. Both pesticides are found in drainage ditches and surface water bodies standing adjacent to crops. Methomyl, a carbamate insecticide widely used in agriculture, inhibits acetylcholinesterase, a key enzyme in nerve transmission. Propanil, an acetanilide herbicide, is used to control grass and broad-leaf weeds. The phenotypic effects of single doses of each chemical were evaluated using a standard immobilisation assay. Immobilisation was linked to global mRNA expression levels using the previously estimated 48h-EC(1)s, followed by hybridization to a cDNA microarray with more than 13,000 redundant cDNA clones representing >5000 unique genes. Following exposure to methomyl and propanil, differential expression was found for 624 and 551 cDNAs, respectively (one-way ANOVA with Bonferroni correction, Ppesticides promoted transcriptional changes in energy metabolism (e.g., mitochondrial proteins, ATP synthesis-related proteins), moulting (e.g., chitin-binding proteins, cuticular proteins) and protein biosynthesis (e.g., ribosomal proteins, transcription factors). Methomyl induced the transcription of genes involved in specific processes such as ion homeostasis and xenobiotic metabolism. Propanil highly promoted haemoglobin synthesis and up-regulated genes specifically related to defence mechanisms (e.g., innate immunity response systems) and neuronal pathways. Pesticide-specific toxic responses were found but there is little evidence for transcriptional responses purely restricted to genes associated with the pesticide target site or mechanism of toxicity.

A method for estimating the turbulent kinetic energy dissipation rate from a vertically pointing Doppler lidar, and independent evaluation from balloon-borne in situ measurements

A method of estimating dissipation rates from a vertically pointing Doppler lidar with high temporal and spatial resolution has been evaluated by comparison with independent measurements derived from a balloon-borne sonic anemometer. This method utilizes the variance of the mean Doppler velocity from a number of sequential samples and requires an estimate of the horizontal wind speed. The noise contribution to the variance can be estimated from the observed signal-to-noise ratio and removed where appropriate. The relative size of the noise variance to the observed variance provides a measure of the confidence in the retrieval. Comparison with in situ dissipation rates derived from the balloon-borne sonic anemometer reveal that this particular Doppler lidar is capable of retrieving dissipation rates over a range of at least three orders of magnitude. This method is most suitable for retrieval of dissipation rates within the convective well-mixed boundary layer where the scales of motion that the Doppler lidar probes remain well within the inertial subrange. Caution must be applied when estimating dissipation rates in more quiescent conditions. For the particular Doppler lidar described here, the selection of suitably short integration times will permit this method to be applicable in such situations but at the expense of accuracy in the Doppler velocity estimates. The two case studies presented here suggest that, with profiles every 4 s, reliable estimates of ϵ can be derived to within at least an order of magnitude throughout almost all of the lowest 2 km and, in the convective boundary layer, to within 50%. Increasing the integration time for individual profiles to 30 s can improve the accuracy substantially but potentially confines retrievals to within the convective boundary layer. Therefore, optimization of certain instrument parameters may be required for specific implementations.

On the structure of Langmuir turbulence

The Stokes drift induced by surface waves distorts turbulence in the wind-driven mixed layer of the ocean, leading to the development of streamwise vortices, or Langmuir circulations, on a wide range of scales. We investigate the structure of the resulting Langmuir turbulence, and contrast it with the structure of shear turbulence, using rapid distortion theory (RDT) and kinematic simulation of turbulence. Firstly, these linear models show clearly why elongated streamwise vortices are produced in Langmuir turbulence, when Stokes drift tilts and stretches vertical vorticity into horizontal vorticity, whereas elongated streaky structures in streamwise velocity fluctuations (u) are produced in shear turbulence, because there is a cancellation in the streamwise vorticity equation and instead it is vertical vorticity that is amplified. Secondly, we develop scaling arguments, illustrated by analysing data from LES, that indicate that Langmuir turbulence is generated when the deformation of the turbulence by mean shear is much weaker than the deformation by the Stokes drift. These scalings motivate a quantitative RDT model of Langmuir turbulence that accounts for deformation of turbulence by Stokes drift and blocking by the air–sea interface that is shown to yield profiles of the velocity variances in good agreement with LES. The physical picture that emerges, at least in the LES, is as follows. Early in the life cycle of a Langmuir eddy initial turbulent disturbances of vertical vorticity are amplified algebraically by the Stokes drift into elongated streamwise vortices, the Langmuir eddies. The turbulence is thus in a near two-component state, with suppressed and . Near the surface, over a depth of order the integral length scale of the turbulence, the vertical velocity (w) is brought to zero by blocking of the air–sea interface. Since the turbulence is nearly two-component, this vertical energy is transferred into the spanwise fluctuations, considerably enhancing at the interface. After a time of order half the eddy decorrelation time the nonlinear processes, such as distortion by the strain field of the surrounding eddies, arrest the deformation and the Langmuir eddy decays. Presumably, Langmuir turbulence then consists of a statistically steady state of such Langmuir eddies. The analysis then provides a dynamical connection between the flow structures in LES of Langmuir turbulence and the dominant balance between Stokes production and dissipation in the turbulent kinetic energy budget, found by previous authors.

A statistical mechanical approach for the computation of the climatic response to general forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. In this paper we show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied with great success to analyze the climatic response to general forcings. The crucial value of the Ruelle theory lies in the fact that it allows to compute the response of the system in terms of expectation values of explicit and computable functions of the phase space averaged over the invariant measure of the unperturbed state. We choose as test bed a classical version of the Lorenz 96 model, which, in spite of its simplicity, has a well-recognized prototypical value as it is a spatially extended one-dimensional model and presents the basic ingredients, such as dissipation, advection and the presence of an external forcing, of the actual atmosphere. We recapitulate the main aspects of the general response theory and propose some new general results. We then analyze the frequency dependence of the response of both local and global observables to perturbations having localized as well as global spatial patterns. We derive analytically several properties of the corresponding susceptibilities, such as asymptotic behavior, validity of Kramers-Kronig relations, and sum rules, whose main ingredient is the causality principle. We show that all the coefficients of the leading asymptotic expansions as well as the integral constraints can be written as linear function of parameters that describe the unperturbed properties of the system, such as its average energy. Some newly obtained empirical closure equations for such parameters allow to define such properties as an explicit function of the unperturbed forcing parameter alone for a general class of chaotic Lorenz 96 models. We then verify the theoretical predictions from the outputs of the simulations up to a high degree of precision. The theory is used to explain differences in the response of local and global observables, to define the intensive properties of the system, which do not depend on the spatial resolution of the Lorenz 96 model, and to generalize the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable for finite as well as infinite time. Finally, we propose a simple yet general methodology to study general Climate Change problems on virtually any time scale by resorting to only well selected simulations, and by taking full advantage of ensemble methods. The specific case of globally averaged surface temperature response to a general pattern of change of the CO2 concentration is discussed. We believe that the proposed approach may constitute a mathematically rigorous and practically very effective way to approach the problem of climate sensitivity, climate prediction, and climate change from a radically new perspective.

Suitability of pesticide risk indicators for less developed countries: a comparison

Pesticide risk indicators provide simple support in the assessment of environmental and health risks from pesticide use, and can therefore inform policies to foster a sustainable interaction of agriculture with the environment. For their relative simplicity, indicators may be particularly useful under conditions of limited data availability and resources, such as in Less Developed Countries (LDCs). However, indicator complexity can vary significantly, in particular between those that rely on an exposure–toxicity ratio (ETR) and those that do not. In addition, pesticide risk indicators are usually developed for Western contexts, which might cause incorrect estimation in LDCs. This study investigated the appropriateness of seven pesticide risk indicators for use in LDCs, with reference to smallholding agriculture in Colombia. Seven farm-level indicators, among which 3 relied on an ETR (POCER, EPRIP, PIRI) and 4 on a non-ETR approach (EIQ, PestScreen, OHRI, Dosemeci et al., 2002), were calculated and then compared by means of the Spearman rank correlation test. Indicators were also compared with respect to key indicator characteristics, i.e. user friendliness and ability to represent the system under study. The comparison of the indicators in terms of the total environmental risk suggests that the indicators not relying on an ETR approach cannot be used as a reliable proxy for more complex, i.e. ETR, indicators. ETR indicators, when user-friendly, show a comparative advantage over non-ETR in best combining the need for a relatively simple tool to be used in contexts of limited data availability and resources, and for a reliable estimation of environmental risk. Non-ETR indicators remain useful and accessible tools to discriminate between different pesticides prior to application. Concerning the human health risk, simple algorithms seem more appropriate for assessing human health risk in LDCs. However, further research on health risk indicators and their validation under LDC conditions is needed.