Reply to the EFSA (2016) on the relevance of recent publications (Hofmann et al. 2014, 2016) on environmental risk assessment and management of Bt-maize events (MON810, Bt11 and 1507)

In this commentary, we respond to a report of the EFSA GMO Panel (EFSA EFSA Supp Publ, 1) that criticises the outcomes of two studies published in this journal (Hofmann et al. Environ Sci Eur 26: 24, 2; Environ Sci Eur 28: 14, 3). Both publications relate to the environmental risk assessment and management of Bt-maize, including maize events MON810, Bt11 and maize 1507. The results of Hofmann et al. (Environ Sci Eur 26: 24, 2), using standardised pollen mass filter deposition measurements, indicated that the EFSA Panel model had underestimated pollen deposition and, hence, exposure of non-target organisms to Bt-maize pollen. The results implied a need for safety buffer distances in the kilometre range for protected nature reserve areas instead of the 20–30 m range recommended by the EFSA Panel. As a result, the EFSA Panel revised their model (EFSA EFSA J 13: 4127, 4), adopting the slope of the empirical data from Hofmann et al. The intercept, however, was substantially reduced to less than 1% at one point by introducing further assumptions based on the estimates of mainly panel members, citing possible ‘uncertainty’. Hofmann et al. (Environ Sci Eur 28: 14, 3) published extensive empirical data regarding pollen deposition on leaves. These results were part of a larger 3-year study involving detailed measurements of pollen release, dispersal and deposition over the maize flowering period. The data collected in situ confirmed the previous predictions of Hofmann et al. (Environ Sci Eur 26: 24, 2). Mean levels and observed variability of pollen deposition on maize and four lepidopteran host plants exceeded the assumptions and disagreed with the conclusions of the EFSA Panel. The EFSA Panel reacted in a report (EFSA EFSA Supp Publ, 1) criticising the methods and outcomes of the two published studies of Hofmann et al. while reaffirming their original recommendations. We respond here point-by-point, showing that the critique is not justified. Based on our results on Urtica leaf pollen density, we confirm the need for specific environmental impact assessments for Bt-maize cultivation with respect to protected habitats within isolation buffer distances in the kilometre range.     

Part I : <Emphasis Type="Italic">a priori</Emphasis> study of erosion and deposition with large eddy simulation of turbulent flow over multiple 2D sandy Gaussian hills

The micro-scale prediction of sand trapping or take-off over hilly terrains is a crucial issue in semi-arid regions for soil depletion. In this context, large eddy simulations around one or several hills are performed in order to provide statistical parameters to characterize the flow at micro-scales and provide data for mesoscale modelling. We focus on the determination of recirculation zones since they play an important role in solid particle erosion or entrapment. A new wall modeling adapted from Huang et al. (J Turbul 17:1–24, 2016) for rough boundary layers is found to improve the prediction of the recirculation zone length downstream of an isolated hill and is used for all the numerical cases presented here. A geometrical parameterization of the recirculation zones is proposed. When the recirculation region is assumed to have an ellipsoidal shape, the total surface of the recirculation can be obtained from this new parameterization and easily extrapolated to more general dune configurations. Numerical results are compared with experiments performed in our laboratory (Simoëns et al. in Procedia IUTAM 17:110–118, 2015) and good agreement is achieved. We explore general aerodynamic cases deduced from the urban canopy scheme of Oke (Energy Build 11:103–113, 1988). In this scheme the momentum and mass exchange between the upper layer and the space between hills is sorted according to the streamwise hill spacing within three basic cases of skimming, wake or isolated flow. The study of the recirculation zones, the mean velocity and Reynolds stress profiles around an isolated or two consecutive hills with different distances shows that the double hill configuration with 3H separation behaves as much as a whole to the upcoming flow. The vortex formed between the crests does not strongly affect the overall evolution of the outer flow. By an a priori prediction of the preferential zones of erosion and accumulation of fictive particles, it is shown that isolated dunes present more deposition and less erosion than two-hill configurations. The results presented in this study will be discussed in the presence of Lagrangian transport of sand particles above 2D Gaussian hills in future work.     

Modeling and experiments of polydisperse particle clouds

A model for polydisperse particle clouds has been developed in this study. We extended the monodisperse particle cloud model of Lai et al. (Environ Fluid Mech 13(5):435–463, 2013) to the case of polydisperse particles. The particle cloud is first considered to be a thermal or buoyant vortex ring, with the thermal induced velocity field modeled by an expanding spherical Hill’s vortex. The buoyancy of the composite thermal is assumed to be the sum of buoyancy contributed by the all particles inside the thermal. Individual particles (of different particle properties) in the cloud are then tracked by the particle tracking equation using the computed induced velocity field. The turbulent dispersion effect is also accounted for by using a random walk model. Experiments of polydisperse particle clouds were carried out to validate the model. The agreement between model predictions and experiments was reasonable. We further validate our model by comparing it with the LES study of Wang et al. (J Hydraul Eng ASCE 141(7):06015006, 2014). The limitations of our model are then discussed with reference to the comparison. Overall, although some flow details are not captured by our model, the simplicity and generality of the model makes it useful in engineering applications.     

Response to Kruse-Plass et al. (2017) regarding the risk to non-target lepidopteran larvae exposed to pollen from one or more of three Bt maize events (MON810, Bt11 and 1507)

We respond to the paper of Kruse-Plass et al. (Environ Sci Eur 29:12, 2017), published in this journal, regarding the risk to non-target lepidopteran larvae exposed to pollen from one or more of three Bt maize events (MON810, Bt11 and 1507). We emphasise that what is important for environmental risk assessment is not the number of pollen grains per se, but the degree of exposure of a NT lepidopteran larva to Bt protein contained in maize pollen. The main text of this response deals with general issues which Kruse-Plass et al. have failed to understand; more detailed refutations of each of their claims are given in Additional file 1. Valid environmental risk assessment requires direct measurement of pollen on leaves at varying distances outside a source field(s); such measurements reflect the potential exposure experienced by an individual larva on a host plant. There are no new data in the Kruse-Plass et al. paper, or indeed any data directly quantifying pollen on actual host-plant leaves outside a maize field; only data gathered within or at the edge of maize crops were reported. Values quoted by Kruse-Plass et al. for deposition on host plants outside the field were estimates only. We reiterate the severe methodological criticisms made by EFSA [Relevance of a new scientific publication (Hofmann et al. 2016) for previous environmental risk assessment conclusions and risk management recommendations on the cultivation of Bt-maize events MON810, Bt11 and 1507. EFSA Supp Publ; EN-1070, 2016], which render this estimation procedure unreliable. Furthermore, criticisms of EFSA (EFSA J 2015(13):4127, 2015) and of EFSA [Relevance of a new scientific publication (Hofmann et al. 2016) for previous environmental risk assessment conclusions and risk management recommendations on the cultivation of Bt-maize events MON810, Bt11 and 1507. EFSA Supp Publ; EN-1070, 2016] made by Kruse-Plass et al. are shown in Additional file 1 to be without foundation. We therefore consider that there is no valid evidence presented by Kruse-Plass et al. to justify their conclusions.     

Fitting the truncated negative binomial distribution to count data

Modeling empirical distributions of repeated counts with parametric probability distributions is a frequent problem when studying species abundance. One must choose a family of distributions which is flexible enough to take into account very diverse patterns and possess parameters with clear biological/ecological interpretations. The negative binomial distribution fulfills these criteria and was selected for modeling counts of marine fish and invertebrates. This distribution depends on a vector \(\left( K,\mathfrak {P}\right) \) of parameters, and ranges from the Poisson distribution (when \(K\rightarrow +\infty \)) to Fisher’s log-series, when \(K\rightarrow 0\). Moreover, these parameters have biological/ecological interpretations which are detailed in the literature and in this study. We compared three estimators of K, \(\mathfrak {P}\) and the parameter \(\alpha \) of Fisher’s log-series, following the work of Rao CR (Statistical ecology. Pennsylvania State University Press, University Park, 1971) on a three-parameter unstandardized variant of the negative binomial distribution. We further investigated the coherence underlying parameter values resulting from the different estimators, using both real count data collected in the Mauritanian Exclusive Economic Zone (MEEZ) during the period 1987–2010 and realistic simulations of these data. In the case of the MEEZ, we first built homogeneous lists of counts (replicates), by gathering observations of each species with respect to “typical environments” obtained by clustering the sampled stations. The best estimation of \(\left( K,\mathfrak {P}\right) \) was generally obtained by penalized minimum Hellinger distance estimation. Interestingly, the parameters of most of the correctly sampled species seem compatible with the classical birth-and-dead model of population growth with immigration by Kendall (Biometrika 35:6–15, 1948).     

On conditional skewness with applications to environmental data

The statistical literature contains many univariate and multivariate skewness measures that allow two datasets to be compared, some of which are defined in terms of quantile values. In most situations, the comparison between two random vectors focuses on univariate comparisons of conditional random variables truncated in quantiles; this kind of comparison is of particular interest in the environmental sciences. In this work, we describe a new approach to comparing skewness in terms of the univariate convex transform ordering proposed by van Zwet (Convex transformations of random variables. Mathematical Centre Tracts, Amsterdam, 1964), associated with skewness as well as concentration. The key to these comparisons is the underlying dependence structure of the random vectors. Below we describe graphical tools and use several examples to illustrate these comparisons.     

Buffer sensitivity of photosynthetic carbon utilisation in eight tropical seagrasses

Some of the mechanisms involved in inorganic carbon (Ci) acquisition by tropical seagrasses from the western Indian Ocean were described by Björk et al. (Mar Biol 129:363–366, 1997). However, since then, it has been found that an additional, buffer-sensitive, system of Ci utilisation may operate in some temperate seagrasses (Hellblom et al. in Aquat Bot 69:55–62, 2001, Hellblom and Axelsson in Photos Res 77:173–191, 2003); this buffer sensitivity indicates a mechanism in which electrogenic H⁺ extrusion may form acidic diffusion boundary layers, in which either HCO ₃ ^? –H⁺ is co-transported into the cells, or where HCO ₃ ^? is converted to CO₂ (as catalysed by carbonic anhydrase) prior to uptake of the latter Ci form. Because a buffer was used in the 1997 study, we found it important to reinvestigate those same eight species, taking into account the direct effect of buffers on this potential mode of Ci acquisition in these plants. In doing so, it was found that all seagrass species investigated except Cymodocea serrulata were sensitive to 50 mM TRIS buffer of the same pH as the natural seawater in which they grew (pH 8.0). Especially sensitive were Halophila ovalis, Halodule wrightii and Cymodocea rotundata, which grow high up in the intertidal zone (only ca. 50–65% of the net photosynthetic activity remained after the buffer additions), followed by the submerged Enhalus acoroides and Syringodium isoetifolium (ca. 75% activity remaining), while Thalassia hemprichii and Thalassodendron ciliatum, which grow in-between the two zones, were less sensitive to buffer additions (ca. 80–85% activity remaining). In addition to buffer sensitivity, all species were also sensitive to acetazolamide (AZ, an inhibitor of extracellular carbonic anhydrase activity) such that ca. 45–80% (but 90% for H. ovalis) of the net photosynthetic activity remained after adding this inhibitor. Raising the pH to 8.8 (in the presence of AZ) drastically reduced net photosynthetic rates (0–14% remaining in all species); it is assumed that this reduction in rates was due to the decreased CO₂ concentration at the higher pH. These results indicate that part of the 1997 results for the same species were due to a buffer effect on net photosynthesis. Based on the present results, it is concluded that (1) photosynthetic Ci acquisition in six of the eight investigated species is based on carbonic anhydrase catalysed HCO ₃ ^? to CO₂ conversions within an acidified diffusion boundary layer, (2) C. serrulata appears to support its photosynthesis by extracellular carbonic anhydrase catalysed CO₂ formation from HCO ₃ ^? without the need for acidic zones, (3) H. ovalis features a system in which H⁺ extrusion may be followed by HCO ₃ ^? –H⁺ co-transport into the cells, and (4) direct, non-H⁺-mediated, uptake of HCO ₃ ^? is improbable for any of the species.     

Ecotoxicological risks of the abandoned F–Ba–Pb–Zn mining area of Osor (Spain)

Due to its potential toxic properties, metal mobilization is of major concern in areas surrounding Pb–Zn mines. In the present study, metal contents and toxicity of soils, aqueous extracts from soils and mine drainage waters from an abandoned F–Ba–Pb–Zn mining area in Osor (Girona, NE Spain) were evaluated through chemical extractions and ecotoxicity bioassays. Toxicity assessment in the terrestrial compartment included lethal and sublethal endpoints on earthworms Eisenia fetida, arthropods Folsomia candida and several plant species, whereas aquatic tests involved bacteria Vibrio fischeri, microalgae Raphidocelis subcapitata and crustaceans Daphnia magna. Total concentrations of Ba (250–5110 mg kg^?1), Pb (940 to >5000 mg kg^?1) and Zn (2370–11,300 mg kg^?1) in soils exceeded intervention values to protect human health. Risks for the aquatic compartment were identified in the release of drainage waters and in the potential leaching and runoff of metals from contaminated soils, with Cd (1.98–9.15 µg L^?1), Pb (2.11–326 µg L^?1) and Zn (280–2900 µg L^?1) concentrations in filtered water samples surpassing US EPA Water Quality Criteria (2016a, b). Terrestrial ecotoxicity tests were in accordance with metal quantifications and identified the most polluted soil as the most toxic. Avoidance and reproduction tests with earthworms showed the highest sensitivity to metal contamination. Aquatic bioassays performed in aqueous extracts from soils confirmed the results from terrestrial tests and also detected toxic effects caused by the mine drainage waters. Algal growth inhibition was the most sensitive aquatic endpoint. In view of the results, the application of a containment or remediative procedure in the area is encouraged.     

Microbial quality and physical–chemical characteristics of thermal springs

Microbial quality and physical–chemical properties of recreational spas were surveyed to investigate the health aspect of the spas’ water. A total of 195 samples were collected from pools and springs of the spas in five sites from Ardebil Province of Iran. The effects of an independent factor defined as ‘condition’ and its component sub-factors (i.e., sampling point, location, and sampling date) on microbial quality and physical–chemical properties of the spas were studied by applying path analysis. The influence of physical–chemical properties on microbial quality was also considered. The percentage of samples exceeding the ISIRI (Swimming pool water microbiological specifications (vol 9412), Institute of Standards and Industrial Research of Iran, Tehran, 2007) limits for Staphylococcus (spp.) was up to 55.8 in the springs and 87.8 in the pools, 58.1 and 99.2 for HPC, 90.7 and 97.8 for total coliform and fecal coliform, and 9.3 and 34.4 for Pseudomonas aeruginosa, respectively. There were significant differences between the pools and springs for both physical–chemical properties and microbial quality. From the path analysis, sampling point was the most effective sub-factor of ‘condition’ on both the physical–chemical properties and microbial quality. Among the physical–chemical properties, water color had the most enhancing or additive influence on microbial pollution, while EC indicated a reducing or subtractive effect.     

Four genes,morphology and ecology: distinguishing a new species of <Emphasis Type="Italic">Acesta</Emphasis> (Mollusca; Bivalvia) from the Gulf of Mexico

There is currently an active scientific debate about the most correct and efficient way to identify species. To date, few studies in the marine realm have combined the available taxonomic methods. In this study we have used morphology, ecology and molecular analyses to identify a new species within the bivalve genus Acesta. All four genes studied (12S, 16S, Cytb, COI) suggested that a common cold-seep species in the Gulf of Mexico, A. bullisi, should be divided into two distinct species. This conclusion was supported by morphological traits and by observations of ecological distribution. A. oophaga Järnegren, Schander and Young n. sp. is described here, and A. bullisi Vokes (Tulane Stud Geol 2:75–92, 1963) is re-described. This study shows that DNA barcoding in combination with traditional morphological and ecological analyses may be an important tools to identify hidden biodiversity among deep-water organisms such as bivalves.     

Evaluation of existing equations for temporal scour depth around circular bridge piers

Scour is defined as the processes of removal of sediment particles from water stream bed by the erosive action of activated water, and also carries sediment away from the hydraulic structures. Scour is the main cause of pier failure. Numerous equations are available for estimating temporal and equilibrium scour depth. The present study describes the phenomenon of temporal scour depth variation at bridge piers and deals with the methods for its estimation. The accuracy of six temporal scour depth equations are also checked in this study. After graphical and statistical analysis, it was found that the relationship proposed by Oliveto and Hager (J Hydraul Eng (ASCE) 128(9):811–520, 2002) predicts temporal scour depth better than other equations. Three equations of equilibrium time of scour are also used for computing equilibrium time. Equilibrium time equation proposed by Choi and Choi (Water Environ J 30(1–2):14–21, 2016) gives better agreements with observed values.     

Effects of temperature and salinity acclimation of adults on larval survival,physiology, and early development of Lytechinus variegatus (Echinodermata: Echinoidea)

Larval survival and developmental rates of Lytechinus variegatus (Lamarck) were determined as a function of temperature and salinity in two experiments by: (1) directly transferring fertilized eggs to 35, 30, 27.5, 25, 20, 15, and 10S seawater at 18 and 23°C, and (2) acclimation of adult sea urchins to the conditions described above for 1 to 4 wk prior to spawning. Developmental rates and percent survival of larvae prior to metamorphosis decreased at salinities below 35 (Q₁₀ values for metamorphosis=0.380 to 0.384). Temperature and salinity significantly (P<0.05) affected metabolic rates of L. variegatus plutei. These results show that L. variegatus larvae are stenohaline when compared to larvae of other echinoderm species. LC₅₀ values (S), developmental rates, and survival to metamorphosis indicate that acclimation of adult sea urchins to lower salinity prior to spawing and fertilization does not enhance development or survival of embryos exposed to low salinity.     

Optimal model averaging by data perturbation for spatial noisy data

Geostatistical models play an important role in spatial data analysis, in which model selection is inevitable. Model selection methods, such as AIC and BIC, are popular for selecting appropriate models. In recent years, some model averaging methods, such as smoothed AIC and smoothed BIC, are also applied to spatial data models. However, the corresponding averaging estimators are outperformed by optimal model averaging estimators (Hansen in Econometrica 75:1175–1189, 2007) for the ordinary linear models. Therefore, this paper focuses on the optimal model averaging method for geostatistical models. We propose a weight choice criterion for the model averaging estimator on the basis of the generalized degrees of freedom and data perturbation technique. We further theoretically prove the resultant estimator is asymptotically optimal in terms of the mean squared error, and numerically demonstrate its satisfactory performance. Finally, the proposed method is applied to a mercury data set.     

Bayesian estimation of species relative abundances and habitat preferences using opportunistic data

We develop a new statistical procedure to monitor relative species abundances and their respective preferences for different habitat types, using opportunistic data. Following Giraud et al. (Biometrics 72(2):649–658, 2015), we combine the opportunistic data with some standardized data in order to correct the bias inherent to the opportunistic data collection. Species observations are modeled by Poisson distributions whose parameters quantify species abundances and habitat preferences, and are estimated using Bayesian computations. Our main contributions are (i) to tackle the bias induced by habitat selection behaviors, (ii) to handle data where the habitat type associated to each observation is unknown, (iii) to estimate probabilities of selection of habitat for the species. As an illustration, we estimate common bird species habitat preferences and abundances in the region of Aquitaine (France).     

Geochemical distribution of trace element concentrations in the vicinity of Boroo gold mine,Selenge Province,Mongolia

The environmental impacts of Boroo gold mine project in Mongolia was evaluated by chemical characterization of trace element concentrations in water, soils and tailing dam sediment samples. The results showed that concentrations of B, Cd, Ni and Se in the water samples were within the accepted levels of the Mongolia water quality standard (MNS4586: 1998). However, the concentrations of Al, As, Cu, Mn, Fe, Pb, U and Zn were higher than the maximum allowable concentration especially in the monitoring and heap leach wells. The average concentrations of As, Cd, Cu, Ni, Pb and Zn in the tailing dam sediment were 4419, 58.5, 56.0, 4.8, 20.6 and 25.7 mg/kg, respectively. Generally, arsenic and heavy metals in the soil samples were within the acceptable concentrations of the soil standard of Mongolia (MNS 5850: 2008). The chemical characterization of As solid phase in tailing dam sediment showed that the majority of As were found in the residual fraction comprising about 74% of total As. Assessing the potential risk to humans, simple bioavailability extraction test was used to estimate bioavailability of arsenic and heavy metals, and the concentrations extracted from tailing dam sediment were; 288.2 mg/kg As, 7.2 mg/kg Cd, 41.1 mg/kg Cu, 13.5 mg/kg Pb, 4.7 mg/kg Ni and 23.5 mg/kg Zn, respectively. From these results, the Boroo gold mine project has presently not significantly impacted the environment, but there is a high probability that it may act as a source of future contamination.     

Estimating the spatial covariance structure using the geoadditive model

In geostatistics, both kriging and smoothing splines are commonly used to generate an interpolated map of a quantity of interest. The geoadditive model proposed by Kammann and Wand (J R Stat Soc: Ser C (Appl Stat) 52(1):1–18, 2003) represents a fusion of kriging and penalized spline additive models. Complex data issues, including non-linear covariate trends, multiple measurements at a location and clustered observations are easily handled using the geoadditive model. We propose a likelihood based estimation procedure that enables the estimation of the range (spatial decay) parameter associated with the penalized splines of the spatial component in the geoadditive model. We present how the spatial covariance structure (covariogram) can be derived from the geoadditive model. In a simulation study, we show that the underlying spatial process and prediction of the spatial map are estimated well using the proposed likelihood based estimation procedure. We present several applications of the proposed methods on real-life data examples.     

The role of black carbon as a catalyst for environmental redox transformation

Black carbon (BC) is an important class of geosorbents that control the fate and transport of organic pollutants in soil and sediment. We previously demonstrated a new role of BC as an electron transfer mediator in the abiotic reduction of nitroaromatic and nitramine compounds by Oh and Chiu (Environ Sci Technol 43:6983–6988, 2009). We proposed that BC can catalyze the reduction of nitro compounds because it contains microscopic graphitic (graphene) domains, which facilitate both sorption and electron transfer. In this study, we assessed the ability of different types of BC—graphite, activated carbon, and diesel soot—to mediate the reduction of 2,4-dinitrotoluene (DNT) and 2,4-dibromophenol (DBP) by H₂S. All three types of BC enhanced DNT and DBP reduction. H₂S supported BC-mediated reduction, as was observed previously with a thiol reductant. The results suggest that BC may influence the fate of organic pollutants in reducing subsurface environments through redox transformation in addition to sorption.     

Geographical expansion of the invader <Emphasis Type="Italic">Caprella scaura</Emphasis> (Crustacea: Amphipoda: Caprellidae) to the East Atlantic coast

Caprella scaura (Templeton in Trans Entomol Soc Lond 1:185–198, 1836) is a native species to the western Indian Ocean. It was first described from Mauritius and later reported from several regions of the world. During the last decade, the species spread out of the Adriatic Sea, and in 2005, C. scaura was reported in Gerona, on the north-eastern coast of Spain. The present study shows the recent rapid expansion of the species to the East Atlantic. During a sampling survey of harbours along the Strait of Gibraltar (Dec 2009–June 2010), we found high densities of C. scaura associated with the bryozoans Zoobotryon verticillatum Della Chiaje, 1822, in Chipiona and Bugula neritina (L.) in Cádiz, Atlantic coast of southern Spain; on the other hand, during May 2009, an important population of C. scaura was also found in the Canary Islands (Tenerife) associated with aquaculture resources. A total of 1,034 individuals of C. scaura were studied, and we confirm that the populations are effectively established, with high densities and reproducing females during the whole year. The population from Cádiz was characterised by smaller specimens but a higher number of eggs per female. These differences could be related to environmental conditions, mainly sea water temperature and to other factors such as competition with other species, or availability of food or substrate. Significant correlation was found between female size and number of eggs in the three populations, but egg size did not vary with female body length. The morphometric relation length/width of the second gnathopod was also calculated and compared between sexes: larger males presented longer gnathopods at the three sites. The most probable introduction vector of this species is shipping; in fact, the Strait of Gibraltar is characterised by an intense shipping traffic. Also, we have observed that the native Caprella equilibra Say, 1818, is being displaced by C. scaura in the harbour of Cádiz, where this species reach densities of 35,945 ind/1,000 g of B. neritina. Further experimental studies are necessary to explore the potential danger of this species, which might spread into natural habitats too.