Evaluating multi-stakeholder perceptions of project impacts: a participatory value-based multi-criteria approach

This paper describes a participatory multi-stakeholder impact assessment approach based on the concept of a value tree (VT), generated using the value focused thinking method (VFT) developed by Keeney (1992 Keeney, R. 1992. Value-focused thinking, Cambridge, MA: Harvard University Press.  [Google Scholar]). This approach allows stakeholders to specify impact parameters to be evaluated, which in the context of the VT, are organized as goals, objectives and alternatives. The approach is implemented in two phases: Phase I is designed to formulate stakeholders' collective VT following VFT concepts and the cognitive mapping method, and Phase II is a participatory valuation approach based on the VT formed in Phase I. The VT and the VFT concepts are adopted as tools to capture different stakeholders' values, goals and perceptions. The resulting VT is structured as a hierarchy between goals and objectives, and a network consisting of relationships, linkages and cross-impacts of the different alternatives and objectives. The hierarchy and network structure enables stakeholders to decompose complex assessment problem into ‘smaller’ units, which makes for easier and clearer assessment context, without ignoring linkages of the units or assessment elements. The second phase allows stakeholders to express their preferences with respect to each assessment element, through a voting system that ultimately leads to measures of importance or relative weights associated with each element. A modified Analytic Hierarchy Process (called Analytic Network Process) was used to distill relative weights from the voting results. Results obtained from a case study in a Zimbabwean community forest show that the proposed approach is easy to implement and can address questions about whether a project can lead to a positive change in attitudes, and whether the changes actually lead to a propensity to adopt alternatives that the project supports (e.g. conservation-oriented alternatives).     

On the front conditions for gravity currents in channels of general cross-section

We consider the propagation of a high-Reynolds-number gravity current in a horizontal channel with general cross-section whose width is \(f(z), 0 \le z\le H\), and the gravity acceleration g acts in \(-z\) direction. (The classical rectangular cross-section is covered by the particular case \(f(z) =\) const.) We assume a two-layer system of homogeneous fluids of constant densities \(\rho _{c}\) (current, of height \(h < H \)) and smaller \(\rho _{a}\) (ambient, filling the remaining part of the channel). We focus attention on the calculation and assessment of the nose Froude-number condition \(Fr = U/(g' h)^{1/2}\); here U is the speed of propagation of the current and \(g' = (\rho _{c}/\rho _{a}-1) g\) is the reduced gravity. We first revisit the steady-state current, and derive compact insightful expressions of Fr and energy dissipation as a function of \(\varphi \) (\(=\) area fraction occupied by the current in the cross-section). We show that the head loss \(\delta _0\) on the stagnation line is formally a degree of freedom in the determination of \(Fr(\varphi )\), and we clarify the strong connections with the head loss \(\delta \) in the ambient fluid, and with the overall rate of dissipation \(\dot{{\mathcal{D}}}\). We demonstrate that the closure \(\delta _0 = 0\) [suggested by Benjamin (J Fluid Mech 31, 209–248, 1968) for the rectangular cross-section] produces in general the smallest Fr for a given \(\varphi \); the results are valid for a significant range \([0, \varphi _{\max }]\), in which the current is dissipative, except for the point \(\varphi _{\max }\) where \(\delta = \dot{{\mathcal{D}}} = 0\). We show that imposing the closure \(\delta = \dot{{\mathcal{D}}} = 0\), which corresponds to an energy-conserving or non-dissipative current, produces in general unacceptable restrictions of the range of validity, and large values of Fr; in particular, deep currents (\(\varphi < 0.3\) say) must be excluded because they are inherently dissipative. On the other hand, the compromise closure \(\delta (\varphi ) =\delta _0(\varphi )\) produces the simple \(Fr(\varphi ) = \sqrt{2}(1 - \varphi )\) formula whose values and dissipation properties are very close, and the range of validity is identical, to these obtained with Benjamin’s closure (moreover, we show that this corresponds to circulation-conservation solutions). The results are illustrated for practical cross-section geometries (rectangle, \(\Delta \) and \(\nabla \) triangle, circle, and the general power-law \(f(z) = b z ^\alpha \) (\(b>0, \alpha \ge 0, 0< z \le H\)). Next, we investigate the connection of the steady-state results with the time-dependent current, and show that in a lock-released current the rate of dissipation of the system is equal to, or larger than, that obtained for Fr corresponding to the conditions at the nose of the current. The results and insights of this study cover a wide range of cross-section geometry and apply to both Boussinesq and non-Boussinesq systems; they reveal a remarkable robustness of Fr as a function of \(\varphi \).     

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).     

Data depth for the uniform distribution

Given a set $X$ of $k$ points and a point $z$ in the $n$ -dimensional euclidean space, the Tukey depth of $z$ with respect to $X$ , is defined as $m/k$ , where $m$ is the minimum integer such that $z$ is not in the convex hull of some set of $k-m$ points of $X$ . If $z$ belongs to the closed region $B$ delimited by an ellipsoid, define the continuous depth of $z$ with respect to $B$ as the quotient $V(z)/\text{ Vol }(B)$ , where $V(z)$ is the minimum volume of the intersection of $B$ with the halfspaces defined by any hyperplane passing through $z$ , and $\text{ Vol }(B)$ is the volume of $B$ . We consider $z$ a random variable and prove that, if $z$ is uniformly distributed in $B$ , the continuous depth of $z$ with respect to $B$ has expected value $1/2^{n+1}$ . This result implies that if $z$ and $X$ are uniformly distributed in $B$ , the expected value of Tukey depth of $z$ with respect to $X$ converges to $1/2^{n+1}$ as the number of points $k$ goes to infinity. These findings have applications in ecology, namely within the niche theory, where it is useful to explore and characterize the distribution of points inside species niche.     

Nitrate uptake in marine phytoplankton: Energy sources and the interaction with carbon fixation

Field studies of whole natural phytoplankton communities from Knight Inlet, B. C., Canada and laboratory cultures of the diatom Skeletonema costatum indicate inorganic carbon fixation may be temporarily suppressed following 10 to 15% enrichment with NO ₃ ^- or NH ₄ ⁺ . (This effect is suggested to be due to competition between inorganic carbon and nitrogen for adenosine triphosphate (ATP), and is reduced when chlorophyll a is increased intracellularly after 6 to 8 h.) Results imply that the source of ATP for nitrate uptake is primarily from Photosystem I (cyclic photophosphorylation) in the presence of light. It would appear that a transient nutrient-adaptive response occurs upon addition of extracellular nitrogen.     

Evolution and decay of gravity wavefields in weak-rotating environments: a laboratory study

Gravity waves are prominent physical features that play a fundamental role in transport processes of stratified aquatic ecosystems. In a two-layer stratified basin, the equations of motion for the first vertical mode are equivalent to the linearised shallow water equations for a homogeneous fluid. We adopted this framework to examine the spatiotemporal structure of gravity wavefields weakly affected by the background rotation of a single-layer system of equivalent thickness \(h_{\ell }\), via laboratory experiments performed in a cylindrical basin mounted on a turntable. The wavefield was generated by the release of a diametral linear tilt of the air–water interface, \(\eta _{\ell }\), inducing a basin-scale perturbation that evolved in response to the horizontal pressure gradient and the rotation-induced acceleration. The basin-scale wave response was controlled by an initial perturbation parameter, \({\mathcal{A}}_{*} = \eta _{0}/h_{\ell }\), where \(\eta _{0}\) was the initial displacement of the air–water interface, and by the strength of the background rotation controlled by the Burger number, \({\mathcal{S}}\). We set the experiments to explore a transitional regime from moderate- to weak-rotational environments, \(0.65\le {\mathcal{S}} \le 2\), for a wide range of initial perturbations, \(0.05\le {\mathcal{A}}_{*}\le 1.0\). The evolution of \(\eta _{\ell }\) was registered over a diametral plane by recording a laser-induced optical fluorescence sheet and using a capacitive sensor located near the lateral boundary. The evolution of the gravity wavefields showed substantial variability as a function of the rotational regimes and the radial position. The results demonstrate that the strength of rotation and nonlinearities control the bulk decay rate of the basin-scale gravity waves. The ratio between the experimentally estimated damping timescale, \(T_{d}\), and the seiche period of the basin, \(T_{g}\), has a median value of \(T_{d}/T_{g}\approx 11\), a maximum value of \(T_{d}/T_{g}\approx 10^{3}\) and a minimum value of \(T_{d}/T_{g}\approx 5\). The results of this study are significant for the understanding the dynamics of gravity waves in waterbodies weakly affected by Coriolis acceleration, such as mid- to small-size lakes.     

Passive scalar roughness lengths for atmospheric boundary layer flow over complex, fractal topographies

When modeling atmospheric boundary layer flow over rough landscapes, surface fluxes of flow quantities (momentum, temperature, etc.) can be described with equilibrium logarithmic law expressions, all of which require specification of a roughness length that is, physically, the elevation at which the flow quantity equals its surface value. In high Reynolds number flows, such as the atmospheric boundary layer, inertial forces associated with turbulent eddy motions are responsible for surface momentum fluxes (form, or pressure drag). Surface scalar fluxes, on the other hand, occur exclusively via diffusion in the immediate vicinity of the topography—the interfacial region—before being advected by turbulent eddy motions into the bulk of the flow. Owing to this difference in surface transfer mechanism, the passive scalar roughness length, $z_{0S}$ , is known to be less than the momentum roughness length, $z_0$ . In this work, classical relations are used to specify $z_{0S}$ during large-eddy simulation of atmospheric boundary layer flow over aerodynamically rough, synthetic, fractal topographies which exhibit power-law height energy spectrum, $E_h (k) \sim k^{\beta _s}$ , where $\beta _s$ is a (predefined) spectral exponent. These topographies are convenient since they resemble natural landscapes and $\beta _s$ can be varied to change the topography’s aerodynamic roughness (the study considers a suite of topographies with $-2.4 \le \beta _s \le -1.2$ , where $-2.4$ and $-1.2$ are the “most smooth” and “most rough” cases, respectively, corresponding with roughness Reynolds number, $Re_0 \approx 10$ and $300$ ). It is often assumed that $z_{0S}/z_{0} \approx 10^{-1}$ for all $Re_0$ . But results from this work show that the roughness length ratio, $z_{0S}/z_{0}$ , depends strongly on $Re_0$ , ranging between $10^{-3}$ and $10^{-1}$ .     

Investigations of carotenoids in some fauna of the Adriatic Sea. II. Anemonia sulcata and Actinia equina (Coelenterata: Anthozoa)

Investigations have been carried out on the carotenoids in Anemonia sulcata (Pennant) and Actinia equina (L.) (Coelenterata: Anthozoa) from the Adriatic Sea. The presence of the carotenoids was determined by means of columnar and thin-layer chromatography. The following carotenoids were found: Anemonia sulcata: β-carotene, lutein-like (epoxy and free) and astaxanthin; Actinia equina: α-carotene, β-carotene, actinioerythrin, lutein-like (?) and astaxanthin (ester and free).     

Gravity currents with tailwaters in Boussinesq and non-Boussinesq systems: two-layer shallow-water dam-break solutions and Navier–Stokes simulations

We consider the dam-break initial stage of propagation of a gravity current of density $\rho _{c}$ released from a lock (reservoir) of height $h_0$ in a channel of height $H$ . The channel contains two-layer stratified fluid. One layer, called the “tailwater,” is of the same density as the current and is of thickness $h_T (< h_0)$ , and the other layer, called the “ambient,” is of different density $\rho _{a}$ . Both Boussinesq ( $\rho _{c}/\rho _{a}\approx 1$ ) and non-Boussinesq systems are investigated. By assuming a large Reynolds number, we can model the flow with the two-layer shallow-water approximation. Due to the presence of the tailwater, the “jump conditions” at the front of the current are different from the classical Benjamin formula, and in some circumstances (clarified in the paper) the interface of the current matches smoothly with the horizontal interface of the tailwater. Using the method of characteristics, analytical solutions are derived for various combinations of the governing parameters. To corroborate the results, two-dimensional direct numerical Navier–Stokes simulations are used, and comparisons for about 80 combinations of parameters in the Boussinesq and non-Boussinesq domains are performed. The agreement of speed and height of the current is very close. We conclude that the model yields self-contained and fairly accurate analytical solutions for the dam-break problem under consideration. The results provide reliable insights into the influence of the tailwater on the propagation of the gravity current, for both heavy-into-light and light-into-heavy motions. This is a significant extension of the classical gravity-current theory from the particular $h_T=0$ point to the $h_T > 0$ domain.     

Initial mixing of inclined dense jet in perpendicular crossflow

A comprehensive experimental investigation for an inclined ( $60^{\circ }$ to vertical) dense jet in perpendicular crossflow—with a three-dimensional trajectory—is reported. The detailed tracer concentration field in the vertical cross-section of the bent-over jet is measured by the laser-induced fluorescence technique for a wide range of jet densimetric Froude number $Fr$ and ambient to jet velocity ratios $U_r$ . The jet trajectory and dilution determined from a large number of cross-sectional scalar fields are interpreted by the Lagrangian model over the entire range of jet-dominated to crossflow-dominated regimes. The mixing during the ascent phase of the dense jet resembles that of an advected jet or line puff and changes to a negatively buoyant thermal on descent. It is found that the mixing behavior is governed by a crossflow Froude number $\mathbf{F} = U_r Fr$ . For $\mathbf{F} < 0.8$ , the mixing is jet-dominated and governed by shear entrainment; significant detrainment occurs and the maximum height of rise $Z_{max}$ is under-predicted as in the case of a dense jet in stagnant fluid. While the jet trajectory in the horizontal momentum plane is well-predicted, the measurements indicate a greater rise and slower descent. For $\mathbf{F} \ge 0.8$ the dense jet becomes significantly bent-over during its ascent phase; the jet mixing is dominated by vortex entrainment. For $\mathbf{F} \ge 2$ , the detrainment ceases to have any effect on the jet behavior. The jet trajectory in both the horizontal momentum and buoyancy planes are well predicted by the model. Despite the under-prediction of terminal rise, the jet dilution at a large number of cross-sections covering the ascent and descent of the dense jet are well-predicted. Both the terminal rise and the initial dilution for the inclined jet in perpendicular crossflow are smaller than those of a corresponding vertical jet. Both the maximum terminal rise $Z_{max}$ and horizontal lateral penetration $Y_{max}$ follow a $\mathbf{F}^{-1/2}$ dependence in the crossflow-dominated regime. The initial dilution at terminal rise follows a $S \sim \mathbf{F}^{1/3}$ dependence.     

The role of metabolic nitrogen in coral calcification

When pieces of the staghorn coral Acropora acuminata are incubated with ¹⁴C-urea, the label is incorporated into skeletal carbonate. Incorporation of this label differs from that of H¹⁴CO ₃ ^- , suggesting urea is not immediately hydrolysed to provide a further source of HCO ₃ ^- . The effects of certain organic substrates upon calcification suggest the ornithine cycle is involved. Citrulline, an ornithine cycle intermediate, is found in high concentrations in the tissues of hermatypic corals. Urea, allantoins, NH₃ and arginine are also present. These compounds are barely detectable in zooxanthellae or an ahermatypic coral. The allantoins may be present as calcium salts. It is suggested that allantoins are the medium by which Ca²⁺ and CO₂ are transported to sites of calcification. Hydrolysis of urea, formed by breakdown of allantoins, yields CO₂ and NH₃. The NH₃ may neutralise protons formed during precipitation of CaCO₃ and bring about their removal from sites of calcification. As well as providing urea, the ornithine cycle may also be involved in the removal of NH₃ from sites of calcification.     

Reducing the environmental impact of coal-washing practice in Australia

All Australian coal for export is washed — that is, the low-grade high-ash material is separated out and disposed of at the washery site. The coarse rejects are generally placed in embankments and the tailings are pumped into settling ponds. These methods of waste disposal can create environmental problems. In addition, they represent a significant energy loss, because about one-third of the waste material is combustible and is thus a potentially valuable source of energy. Laboratory studies at the CSIRO Division of Fossil Fuels at North Ryde, New South Wales, have shown that energy can be successfully recovered from washery wastes by the process of fluidized-bed combustion. A pilot fluidized-bed combustor, with a capacity of 2 tonnes/hour, is undergoing prolonged trial at the Clutha washery near Camden, New South Wales, in a joint project between the Joint Coal Board and CSIRO. The results from the pilot plant tests have provided the basis for a conceptual design and feasibility study for a full-scale tailings treatment plant. This study has indicated that fluidized-bed combustion:

1. offers an environmentally attractive and economically competitive alternative to the disposal of tailings in settling ponds, and
2. can be used to generate large quantities of energy from the wastes, reducing the energy lost in coal washing from around 16% to 6% of the coal mined.

     

Hierarchical modeling of abundance in closed population capture–recapture models under heterogeneity

Hierarchical modeling of abundance in space or time using closed-population mark-recapture under heterogeneity (model \(\hbox {M}_{\text {h}}\) ) presents two challenges: (i) finding a flexible likelihood in which abundance appears as an explicit parameter and (ii) fitting the hierarchical model for abundance. The first challenge arises because abundance not only indexes the population size, it also determines the dimension of the capture probabilities in heterogeneity models. A common approach is to use data augmentation to include these capture probabilities directly into the likelihood and fit the model using Bayesian inference via Markov chain Monte Carlo (MCMC). Two such examples of this approach are (i) explicit trans-dimensional MCMC, and (ii) superpopulation data augmentation. The superpopulation approach has the advantage of simple specification that is easily implemented in BUGS and related software. However, it reparameterizes the model so that abundance is no longer included, except as a derived quantity. This is a drawback when hierarchical models for abundance, or related parameters, are desired. Here, we analytically compare the two approaches and show that they are more closely related than might appear superficially. We exploit this relationship to specify the model in a way that allows us to include abundance as a parameter and that facilitates hierarchical modeling using readily available software such as BUGS. We use this approach to model trends in grizzly bear abundance in Yellowstone National Park from 1986 to 1998.     

Waves and turbulence in katabatic winds

The measurements taken during the Vertical Transport and Mixing Experiment (VTMX, October, 2000) on a northeastern slope of Salt Lake Valley, Utah, were used to calculate the statistics of velocity fluctuations in a katabatic gravity current in the absence of synoptic forcing. The data from ultrasonic anemometer-thermometers placed at elevations 4.5 and 13.9 m were used. The contributions of small-scale turbulence and waves were isolated by applying a high-pass digital (Elliptical) filter, whereupon the filtered quantities were identified as small-scale turbulence and the rest as internal gravity waves. Internal waves were found to play a role not only at canonical large gradient Richardson numbers $(\overline{\hbox {Ri}_\mathrm{g} } >1)$ , but sometimes at smaller values $(0.1 < \overline{\hbox {Ri}_\mathrm{g}}<1)$ , in contrast to typical observations in flat-terrain stable boundary layers. This may be attributed, at least partly, to (critical) internal waves on the slope, identified by Princevac et al. [1], which degenerate into turbulence and help maintain an active internal wave field. The applicability of both Monin-Obukhov (MO) similarity theory and local scaling to filtered and unfiltered data was tested by analyzing rms velocity fluctuations as a function of the stability parameter z/L, where L is the Obukhov length and z the height above the ground. For weaker stabilities, $\hbox {z/L}<1$ , the MO similarity and local scaling were valid for both filtered and unfiltered data. Conversely, when $\hbox {z/L}>1$ , the use of both scaling types is questionable, although filtered data showed a tendency to follow local scaling. A relationship between z/L and $\overline{\hbox {Ri}_\mathrm{g} }$ was identified. Eddy diffusivities of momentum $\hbox {K}_\mathrm{M}$ and heat $\hbox {K}_\mathrm{H}$ were dependent on wave activities, notably when $\overline{\hbox {Ri}_\mathrm{g} } > 1$ . The ratio $\hbox {K}_{\mathrm{H}}/\hbox {K}_{\mathrm{M}}$ dropped well below unity at high $\overline{\hbox {Ri}_\mathrm{g} }$ , in consonance with previous laboratory stratified shear layer measurements as well as other field observations.     

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.     

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.