Der Einsatz von Pflanzentests bei der Sedimentbewertung

Background and aim Aquatic organisms of the three trophic levels (producers, consumers, reducers) have been used for a long time for ecotoxicological assessments of water quality (HABAB 2000; HABAK 1999). In biotest systems that have become established standards, algae, daphnia, and luminescent bacteria as test organisms are often used so far. However, comparative studies have shown that algae are not sufficiently sensitive to all phytotoxic agents and thus cannot be the exclusive indicator organism of any toxic effects on autotrophic life. Moreover, the EU Water Framework Directive (EU-WFD 2000) sets macrobenthos, fish, phytoplankton, and macrophytes equally side by side as indicators of the status of waters. Nevertheless, biotest systems relying on higher plants have been rarely used in assessments of aquatic sediments so far. Against this background, the aquatic duckweed test (DIN EN ISO 20079) was developed and standardized, and a sediment contact test was developed with Myriophyllum aquaticum (Feiler et al. 2004). The latter was subject of a joint research project of the German Federal Ministry of Education and Research (BMBF-Verbundprojekt “SeKT”) that examined and compared several sediment contact tests (Feiler et al. 2005). Compared to tests with algae, a biotest using macrophytes has the advantage to be more representative of higher plants, since target organism and test organism are closer related in terms of evolutionary history. The present study demonstrates that the application of biotests with higher plants in analyses of contaminated sediments yields valuable results that may contribute to a concept for the integrative assessment of water quality.     

Model selection criteria for overdispersed data and their application to the characterization of a host-parasite relationship

In the statistical modeling of a biological or ecological phenomenon, selecting an optimal model among a collection of candidates is a critical issue. To identify an optimal candidate model, a number of model selection criteria have been developed and investigated based on estimating Kullback’s (Information theory and statistics. Dover, Mineola, 1968) directed or symmetric divergence. Criteria that target the directed divergence include the Akaike (2nd international symposium on information theory. Akadémia Kiadó, Budapest, Hungary, pp 267–281, 1973, IEEE Trans Autom Control AC 19:716–723, 1974) information criterion, AIC, and the “corrected” Akaike information criterion (Hurvich and Tsai in Biometrika 76:297–307, 1989), AICc; criteria that target the symmetric divergence include the Kullback information criterion, KIC, and the “corrected” Kullback information criterion, KICc (Cavanaugh in Stat Probab Lett 42:333–343, 1999; Aust N Z J Stat 46:257–274, 2004). For overdispersed count data, simple modifications of AIC and AICc have been increasingly utilized: specifically, the quasi Akaike information criterion, QAIC, and its corrected version, QAICc (Lebreton et al. in Ecol Monogr 62(1):67–118 1992). In this paper, we propose analogues of QAIC and QAICc based on estimating the symmetric as opposed to the directed divergence: QKIC and QKICc. We evaluate the selection performance of AIC, AICc, QAIC, QAICc, KIC, KICc, QKIC, and QKICc in a simulation study, and illustrate their practical utility in an ecological application. In our application, we use the criteria to formulate statistical models of the tick (Dermacentor variabilis) load on a white-footed mouse (Peromyscus leucopus) in northern Missouri.     

Quantifying a significance of sediment particle size to hyporheic sedimentary oxygen demand with a permeable stream bed

A mechanistic model of sedimentary oxygen demand (SOD) for hyporheic flow is presented. The permeable sediment bed, e.g. sand or fine gravel, is considered with hydraulic conductivity in the range $0.1 < K < 20$  cm/s. Hyporheic pore water flow is induced by pressure fluctuations at the sediment/water interface due to near-bed turbulent coherent motions. A 2-D advection–diffusion equation is linked to the pore water flow model to simulate the effect of advection–dispersion driven by interstitial flow on oxygen transfer through the permeable sediment. Microbial oxygen uptake in the sediment is expressed as a function of the microbial growth rate, and is related to the sediment properties, i.e. the grain diameter $(d_{s})$ and porosity $(\phi )$ . The model describes the significance of sediment particle size to oxygen transfer through the sediment and microbial oxygen uptake: With increasing grain diameter $(d_{s})$ , the hydraulic conductivity $(K)$ increases so does the oxygen transfer rate, while particle surface area per volume (the available surface area for colonization by biofilms) decreases reducing the microbial oxygen uptake rate. Simulation results show that SOD increases as the hydraulic conductivity $(K)$ increases before a threshold has been reached. After that, SOD diminishes with the increment of the hydraulic conductivity $(K)$ .     

Hydrolithon spp. (Rhodophyta, Corallinales) overgrow live corals (Cnidaria, Scleractinia) in Yemen

In Yemen, off the northwestern coast in the Gulf of Aden, the coralline algae Hydrolithon rupestre (Foslie) Penrose 1996 and H. murakoshii Iryu and Matsuda 1996 have been observed to overgrow and kill living Porites lutea Milne-Edwards and Haime, 1860. Similarly, Hydrolithon onkodes (Heydrich) Penrose and Woelkerling 1992 and H. rupestre were observed overgrowing Stylophora pistillata (Esper, 1797). Competitive interactions between P. lutea and H. murakoshii were monitored from 2006 to 2009 at two sites and showed an average linear growth of 8.3 (±1.9 SD) mm?year^?1 over the coral. The small polyps of S. pistillata and P. lutea combined with putative chemical compounds produced by Hydrolithon spp. are likely to allow the coralline overgrowth. Although corallines can locally kill coral tissues, the CCA/coral interactions do not seem to affect the overall live coral cover at the study sites.     

Using surface drifter observations to measure tidal vortices and relative diffusion at Aransas Pass,Texas

Tidal vortices play an important role in the flushing of coastal regions. At the mouth of a tidal inlet, the input of circulation by the ebb tide may force the formation of a starting-jet dipole vortex. The continuous ebb jet current also creates a periodic sequence of secondary vortices shed from the inlet mouth. In each case, these tidal vortices have a shallow aspect ratio, with a lateral extent much greater than the water depth. These shallow vortices affect the transport of passive tracers, such as nutrients and sediment from the estuary to the ocean and vice versa. Field observation of tidal vortices primarily relies on ensemble averaging over several vortex events that are repeatable in space and can be sampled by a fixed Eulerian measurement grid. This paper presents an adaptive approach for locating and measuring within tidal vortices that propagate offshore near inlets and advect along variable trajectories set by the wind-driven currents. A field experiment was conducted at Aransas Pass, Texas to measure these large-scale vortices. Locations of the vortices produced during ebb tide were determined using near real-time updates from surface drifters deployed near or within the inlet during ebb tide, and the paths of towed acoustic Doppler current profiler (ADCP) transects were selected by analysis of the drifter observations. This method allowed ADCP transects to be collected within ebb generated tidal vortices, and the paths of the drifters indicated the presence of both the starting-jet dipole and the secondary vortices of the unstable ebb tidal jet. Drifter trajectories were also used to estimate the size of each observed vortex as well as the statistics of relative diffusion offshore of Aransas Pass. The field data confirmed the starting-jet spin-up time (time until the vortex dipole begins to propagate offshore) measured in the laboratory by Bryant et al. [6] and that the Strouhal condition of \(St=0.2\) predicts the shedding of secondary vortices from the inlet mouth. The size of the rotational core of the vortex is also shown to be approximated physically by the inlet width or by \(0.02UT\) , where U is the maximum velocity through the inlet channel and T is the tidal period, and confirms results found in laboratory experiments by Nicolau del Roure et al. [23]. Additionally, the scale of diffusion was approximately 1–15 km and the apparent diffusivity was between 2–130  \(m^2/s\) following Richardsons law.     

Extension of the skin shear stress Li’s relationship to the flat bed

A proper estimation of the skin shear stress τ _s is necessary for a proper evaluation of sediment flux at the sediment–fluid interface. Several empirical formulas of the skin shear stress have been proposed in the literature for rippled bed as function of the factor form η/λ (η and λ are respectively the height and wavelength of the bedform). These formulas express that in the presence of bedform, τ _s is a partition of the total shear stress τ _b . In contrast, when the bottom is flat, τ _s is exactly equal to τ _b . Based on in situ measurements, Li (J Geophys Res 99:791–799, 1994) has proposed a new formula of τ _s depending on u _*/η (u _* is the friction velocity, ${u_{*}=\sqrt{\tau_{b}/\rho}}$ ), but not as a function of η/λ. This formulation appears physically more realistic, but does not cover the flat bottom range. The purpose of this short note is therefore the extension of the Li’s expression to this type of bottom.     

Functional ANOVA starting from discrete data: an application to air quality data

A nonparametric functional approach is proposed to compare the mean functions of $k$ k samples of curves. In practice, curves data are usually collected in a discrete form and hence they must be pre-processed to use purely functional techniques. However, in the context of $k$ k -sample tests, the pre-processing step can have effects in terms of power reduction. Hall and Van Keilegom (Stat Sin 17:1511–1531, 2007) proposed a methodology to minimizing these effects in the context of tests for the equality of two distribution functions. Their procedure is here extended to the case of $k$ k -sample hypothesis tests. The asymptotic validity of the procedure is established and its finite sample performance is analyzed through Monte Carlo experiments. As an illustration, the method is applied to air quality data collected from several monitoring stations placed at different geographical locations at the center of Spain.     

Sediment processes and flow reversal in the undular tidal bore of the Garonne River (France)

A tidal bore is a series of waves propagating upstream as the tidal flow turns to rising, and the bore front corresponds to the leading edge of the tidal wave in a funnel shaped estuarine zone with macro-tidal conditions. Some field observations were conducted in the tidal bore of the Garonne River on 7 June 2012 in the Arcins channel, a few weeks after a major flood. The tidal bore was a flat undular bore with a Froude number close to unity: $\hbox {Fr}_{1} = 1.02$ and 1.19 (morning and afternoon respectively). A key feature of the study was the simultaneous recording of the water elevation, instantaneous velocity components and suspended sediment concentration (SSC) estimates, together with a detailed characterisation of the sediment bed materials. The sediment was some silty material ( $\hbox {d}_{50} \approx 13~\upmu \hbox {m}$ ) which exhibited some non-Newtonion thixotropic behaviour. The velocity and SSC estimate were recorded simultaneously at high frequency, enabling a quantitative estimate of the suspended sediment flux at the end of the ebb tide and during the early flood tide. The net sediment flux per unit area was directed upstream after the bore, and its magnitude was much larger than that at end of ebb tide. The field observations highlighted a number of unusual features on the morning of 7 June 2012. These included (a) a slight rise in water elevation starting about 70 s prior to the front, (b) a delayed flow reversal about 50 s after the bore front, (c) some large fluctuations in suspended sediment concentration (SSC) about 100 s after the bore front and (d) a transient water elevation lowering about 10 min after the bore front passage. The measurements of water temperature and salinity showed nearly identical results before and after the tidal bore, with no evidence of saline and thermal front during the study.     

Bootstrap confidence intervals for adaptive cluster sampling design based on Horvitz–Thompson type estimators

Perez and Pontius (J Stat Comput Simul 76:755–764, 2006) introduced several bootstrap methods under adaptive cluster sampling using a Horvitz–Thompson type estimator. Using a simulation study, they showed that their proposed methods provide confidence intervals with highly understated coverage rates. In this article, we first show that their bootstrap methods provide biased bootstrap estimates. We then define two bootstrap methods, based on the method of Gross (Proceeding of the survey research methods section. American Statistical Association, Alexandria, VA, pp 181–184, 1980) and Bootstrap With Replacement, that provide unbiased bootstrap estimates of the population mean with bootstrap variances matching the corresponding unbiased variance estimator. Using a simulation study, we show that the bootstrap confidence intervals based on our proposed methods have better performance than those based on available bootstrap methods, in the sense of having coverage proportion closer to the nominal coverage level. We also compare the proposed intervals to empirical likelihood based intervals in small samples.     

n-Octyl esters of long-chain fatty acids are not anthropogenic pollution markers

n-Octyl esters of higher fatty acids have been reported as markers of urban wastewater in sediments, polychaeta, fish, crabs and oysters (Chaler et al. in J Chromatogr A 1046:203–210, 2004)). However, up to date, there were no subsequent studies to confirm this claim. Likewise, Chaler and co-workers did not consider that the mentioned compounds might occur naturally in the environment. Here we found seven n-octyl straight-chain alkanoates, from C-20 to C-26, in the wild-growing plant Heracleum sphondylium L., Apiaceae, from 14 locations. Those plant metabolites were unambiguously identified by gas chromatographic co-injection of the synthesized esters with the inflorescence washings. All identified octyl esters represent new natural compounds, except for octyl docosanoate. Since we have demonstrated that n-octyl esters occur naturally, and in abundance, they cannot be recognized any longer as wastewater markers. Additionally, here we provide evidences that the compounds identified by Chaler et al. (2004) are in fact 2-ethylhexyl esters, mistakenly identified as n-octyl esters.     

Physical and biological controls of algal blooms in the Río de la Plata

Coupled three-dimensional hydrodynamic and ecological numerical simulations were used to investigate the role of transport, stagnation zones and dispersion on inter-annual blooms of the diatom Aulacoseira sp. in the vicinity of the drinking water intakes of the Buenos Aires city (Argentina) in the upper Río de la Plata. Three different summer events were analyzed. First, a mild biomass bloom year (2006–2007), second, a high biomass bloom year (2007–2008) and third, a “normal” no bloom year (2009–2010). Simulated water height, water temperature, suspended solids and chlorophyll \(a\) concentrations patterns compared well with field data. Results revealed that the advection of phytoplankton cells via inflows to the Río de la Plata triggered Aulacoseira sp. blooms in the domain. In addition, excessive growth observed near the drinking water intakes, along the Argentinean margin, were associated with long retention times (stagnant region) and weak horizontal dispersion. Increased concentrations of suspended solids in the water column, in response to re-suspension events, did not prevent the blooms, however, were found to also play a key role in controlling the rate of phytoplankton growth. Finally, a non-dimensional parameter, R, that considers phytoplankton patch size, e-folding growth and dispersion time scales is shown to determine the potential bloom occurrences, as well as bloom intensity; R values higher than 5.7 suggest intense phytoplankton growth. For the mild biomass bloom year, \(R = 7.5\) , for the high biomass bloom year, \(R = 11\) and for the “normal” no bloom year \(R= 0.4\) .     

Defining behavioural syndromes and the role of ‘syndrome deviation’ in understanding their evolution

This commentary highlights multivariate tools that have been used by evolutionary biologists in the study of syndromes and their evolution and discusses the insights that these methods provide into evolutionary processes relative to the metric ‘syndrome deviation’ that has recently been proposed by Herczeg and Garamszegi (Behav Ecol Sociobiol 66:161–169, 2012). We clarify that non-zero phenotypic correlations arise from the joint influences of within- and between-individual correlations, whereas only non-zero between-individual correlations represent behavioural syndromes, and discuss how acknowledgement of this subtle difference between phenotypic and between-individual correlations affects the applicability of syndrome deviation for the study of behavioural syndromes.     

Statistical modeling to assess drinking water crisis

We modify the beta distribution by adding a parameter ρ ≥ 0 to capture the urgency level with which nations address drinking water crises. In essence, the urgency parameter relates to region’s environmental or development levels. The modified version is named an urgency biased beta distribution (UBBD). Its statistical properties are considered and illustrated by analysis of the drinking water data in the World Health Organization (2011) report. The methodology of this article is applicable to similar data sets arising in ecological, engineering, medical, business, marketing, and communication, etc. studies.     

Analysis of periodic patterns of daily precipitation through simultaneous modeling of its serially observed occurrence and amount

Precipitation is of great importance to agriculture, environment and ecosystem as a regular precipitation pattern is usually vital to healthy plants; excessive or insufficient rainfall can be harmful. Periodic patterns of precipitation can be studied based on regularly observed data over time. Since regularly observed precipitation data are generally skewed with many zeros, two common analysis approaches have been proposed recently. One approach investigates precipitation using a two-part model where the occurrence and positive amount of precipitation are analyzed separately (Piantadosi et al. in Environ Model Assess 14:431–438, 2009), whereas the other approach handles occurrence and amount simultaneously using a Tweedie’s compound Poisson model for independent observations (Hasan and Dunn in Int J Climatol 32:1006–1017, 2012). The former approach fails to maintain the regular temporal structure of serially observed precipitation, whereas the latter approach ignores serial dependence. As there is generally substantial serial correlation in the observed sequence of precipitation data over time, we introduce a compound Poisson state-space model with serially correlated random effects for daily precipitation data. This approach characterizes both occurrence and amount of precipitation simultaneously while accounting for the corresponding serial correlation. Our main results depend only on the first- and second-moment assumptions of unobserved random effects. We illustrate our method with the analysis of the daily precipitation data recorded at Mount Washington, NH, USA.     

A robust weighted EM algorithm for use-availability data

Wildlife sampling for habitat selection often combines a random background sample with a random sample of used sites, because the background sample could contain too few used sites to be informative for rare species. This approach is referred to as use-availability sampling. Two variants are considered where there is: (1) a random background sample including used and unused sites augmented with a sample of used sites, and (2) a sample of used sites augmented with a contaminated background sample, i.e. use is not recorded. A weighted estimator first proposed by Manski and Lerman (Econometrica 45(8):1977?C1988, 1977) forms the basis for our suggested approach. The weighted estimator has been shown to perform better than the usual unweighted approach with uncontaminated data and mis-specified logit models (Xie and Manski in Sociol Methods Res 17(3):283?C302, 1989). A weighted EM algorithm is developed for use with contaminated background data. We show that the weighted estimator continues to perform well with contaminated data and maintains its robustness to model mis-specification. The weighted estimator has not been previously used for use-availability sampling due to reliance on the assumption that only the intercept is biased, which is valid for a correct logit model. We show that adjusting the intercept may not eliminate the bias with an incorrect logit model. In this case, the weighted estimator is a relatively simple and effective alternative.     

Finite area smoothing with generalized distance splines

Most conventional spatial smoothers smooth with respect to the Euclidean distance between observations, even though this distance may not be a meaningful measure of spatial proximity, especially when boundary features are present. When domains have complicated boundaries leakage (the inappropriate linking of parts of the domain which are separated by physical barriers) can occur. To overcome this problem, we develop a method of smoothing with respect to generalized distances, such as within domain distances. We obtain the generalized distances between our points and then use multidimensional scaling to find a configuration of our observations in a Euclidean space of 2 or more dimensions, such that the Euclidian distances between points in that space closely approximate the generalized distances between the points. Smoothing is performed over this new point configuration, using a conventional smoother. To mitigate the problems associated with smoothing in high dimensions we use a generalization of thin plate spline smoothers proposed by Duchon (Constructive theory of functions of several variables, pp 85–100, 1977). This general method for smoothing with respect to generalized distances improves on the performance of previous within-domain distance spatial smoothers, and often provides a more natural model than the soap film approach of Wood et al. (J R Stat Soc Ser B Stat Methodol 70(5):931–955, 2008). The smoothers are of the linear basis with quadratic penalty type easily incorporated into a range of statistical models.     

Ground water contamination with 238U, 234U, 235U, 226Ra and 210Pb from past uranium mining: cove wash,Arizona

The objectives of the study are to present a critical review of the ²³⁸U, ²³⁴U, ²³⁵U, ²²⁶Ra and ²¹⁰Pb levels in water samples from the EPA studies (U.S. EPA in Abandoned uranium mines and the Navajo Nation: Red Valley chapter screening assessment report. Region 9 Superfund Program, San Francisco, 2004, Abandoned uranium mines and the Navajo Nation: Northern aum region screening assessment report. Region 9 Superfund Program, San Francisco, 2006, Health and environmental impacts of uranium contamination, 5-year plan. Region 9 Superfund Program, San Franciso, 2008) and the dose assessment for the population due to ingestion of water containing ²³⁸U and ²³⁴U. The water quality data were taken from Sect. “Data analysis” of the published report, titled Abandoned Uranium Mines Project Arizona, New Mexico, Utah–Navajo Lands 1994–2000, Project Atlas. Total uranium concentration was above the maximum concentration level for drinking water (7.410–1 Bq/L) in 19 % of the water samples, while ²³⁸U and ²³⁴U concentrations were above in 14 and 17 % of the water samples, respectively. ²²⁶Ra and ²¹⁰Pb concentrations in water samples were in the range of 3.7 × 10^?1 to 5.55 × 102 Bq/L and 1.11 to 4.33 × 102 Bq/L, respectively. For only two samples, the ²²⁶Ra concentrations exceeded the MCL for total Ra for drinking water (0.185 Bq/L). However, the ²¹⁰Pb/²²⁶Ra ratios varied from 0.11 to 47.00, and ratios above 1.00 were observed in 71 % of the samples. Secular equilibrium of the natural uranium series was not observed in the data record for most of the water samples. Moreover, the ²³⁵U/^totalU mass ratios ranged from 0.06 to 5.9 %, and the natural mass ratio of ²³⁵U to ^totalU (0.72 %) was observed in only 16 % of the water samples, ratios above or below the natural ratio could not be explained based on data reported by U.S. EPA. In addition, statistical evaluations showed no correlations among the distribution of the radionuclide concentrations in the majority of the water samples, indicating more than one source of contamination could contribute to the sampled sources. The effective doses due to ingestion of the minimum uranium concentrations in water samples exceed the average dose considering inhalation and ingestion of regular diet for other populations around the world (1 μSv/year). The maximum doses due to ingestion of ²³⁸U or ²³⁴U were above the international limit for effective dose for members of the public (1 mSv/year), except for inhabitants of two chapters. The highest effective dose was estimated for inhabitants of Cove, and it was almost 20 times the international limit for members of the public. These results indicate that ingestion of water from some of the sampled sources poses health risks.     

Comment on “Use of local tidal records to identify relative sea level change: accuracy and error for decision makers” by Powell VA, McGlashan DJ, Duck RW (2012) J Coast Conserv

In this Comment we refer to our strong reservations concerning the paper by Powell et al. J Coast Conserv, (2012) recently published on the Online First web site of the Journal of Coastal Conservation. The paper makes a number of comments on data obtained from the Permanent Service for Mean Sea Level (PSMSL) and British Oceanographic Data Centre (BODC) which are incorrect or misleading. In addition, some of their comments on sea-level science in general need to be challenged and corrected.     

A note on a non-stationary point source spatial model

A point source, non-stationary covariance structure model is proposed, having only one additional parameter over a standard, stationary covariance structure, spatial model. Additionally, the proposed model is demonstrated to fit better than the three extra parameter, point source, non-stationary spatial model proposed by Ecker and De Oliveira (Commun Stat Theory Methods 37:2066–2078, 2008). The proposed model is fit from a Bayesian perspective and illustrated using a house sales dataset from Cedar Falls, Iowa.     

Large eddy simulation of dispersion around an isolated cubic building: evaluation of localized dynamic k_\mathrm{SGS}-equation sub-grid scale model

In the present study, the prediction accuracy of a dynamic one-equation sub-grid scale model for the large eddy simulation of dispersion around an isolated cubic building is investigated. For this purpose, the localized dynamic $k_\mathrm{SGS} $ -equation model (LDKM) is employed and the results are compared with the available experimental data and two other classic sub-grid scale models, namely, standard Smagorinsky–Lilly model (SSLM) and dynamic Smagorinsky–Lilly model (DSLM). It is shown that the three SGS models give results in good agreement with experiment. However, near the ground level of the leeward wall, dimensionless time-averaged concentration, $\langle K\rangle $ , profile is not quite similar to the experimental data. It is also demonstrated that the LDKM predicts the values of $\langle K\rangle $ on the roof, leeward and side walls more acceptably than the SSLM and DSLM. Whereas, the streamwise elongation of time-averaged structures of the plume shape is more over-estimated with the LDKM than with the other two SGS models. In terms of numerical difficulty, the LDKM is found to be stable and computationally reasonable. In addition, it does not suffer from a flow dependent constant such as the Smagorinsky coefficient employed in the SSLM model.