Application of a novel open-source program for measuring the effects of toxicants on the swimming behavior of large groups of unmarked fish

The aim of this study was to demonstrate a novel method for measuring the effects of toxicants on the behavior of groups of unmarked fish. Ctrax is an open source Python-based machine vision program that was originally designed to track fruit flies (Drosophila). The effects of the fish anaesthetic quinaldine sulfate on the behavior of groups of 10 unmarked Arabian killifish (Aphanias dispar) in a small arena were measured. Ctrax was able to identify and track the X-Y coordinates of the individual fish over 1 min recording periods. The data was analyzed using specifically written Matlab (The Mathworks™) toolboxes and the effects of quinaldine sulfate on absolute swimming velocity, forward swimming velocity, rate of change in orientation and distance to wall were calculated. Additionally the effects of quinaldine sulfate on path tortuosity were also measured. Ctrax has significant benefits over other available technologies for tracking the 2-dimensional coordinates of fish. In particular the software is open source and thus freely available, it is accurate, requires only simple easily available equipment and able to track potentially large groups of fish (possibly >50) while maintaining their individual identities.     

profiling areas of ground water contamination by pesticides in california: phase ii – evaluation and modification of a statistical model

A well sampling study was conducted to evaluate anempirical approach to classifying areasof land in California as vulnerable to ground watercontamination by pesticides (Troiano et al., 1994). Wells were sampled from sections of land that had noprevious detections of pesticideresidues. The sections had been classified into vulnerablesoil clusters or into a not-classified groupusing a procedure based on Principal Components Analysis(PCA). Grape, citrus, and olive growingareas of Fresno and Tulare Counties were targeted, areas wherepre-emergence herbicide residues hadbeen detected in well water. Overall, herbicide residues weredetected in 75 of 176 sampled wells, ahigh frequency of detection in relation to results fromprevious targeted well sampling studies. Sinceresidues were also detected in the not-classified group, theclassification procedure was modified usingan approach based on Canonical Variates Analysis (CVA). Moresections were classified intovulnerable soil clusters with the CVA approach than with thePCA method. Data from two otherexplanatory variables, depth to ground water and amount ofpesticide used per section, were includedto illustrate how additional information can be incorporatedinto this approach of identifying vulnerable areas.     

Effects of invertebrate predators and a pesticide on temporary pond microcosms used for aquatic toxicity testing

The effects of increased trophic complexity, through the addition of predatory notonectids (Anisops deanei), on temporary pond microcosms used for aquatic toxicity testing were studied. Replicate microcosms were established using sediment from a dried temporary pond, and treated with one of four concentrations of the organochlorine pesticide endosulfan (0, 1, 10 or 50 microg/L), in the presence or absence of six A. deanei. The tanks were sampled regularly for nine weeks following the addition of the predators and the entire contents of each tank counted after 12 weeks. Analysis using non-metric multidimensional scaling (MDS) and non-parametric MANOVA showed that both Anisops and endosulfan at concentrations >10 microg/L significantly altered community structure. However, an interaction between the effects of Anisops and the effects of endosulfan was not detected. The addition of Anisops did not increase the variability of response and thus did not reduce the sensitivity of the test method.     

Probability distributions and statistical inference for axial data

Observations on axes which lack information on the direction of propagation are referred to as axial data. Such data are often encountered in enviromental sciences, e.g. observations on propagations of cracks or on faults in mining walls. Even though such observations are recorded as angles, circular probability models are inappropriate for such data since the constraint that observations lie only in [0, π) needs to be enforced. Probability models for such axial data are argued here to have a general structure stemming from that of wrapping a circular distribution on a semi-circle. In particular, we consider the most popular circular model, the von Mises or circular normal distribution, and derive the corresponding axial normal distribution. Certain properties of this distribution are established. Maximum likelihood estimation of its parameters are shown to be surprisingly, in contrast to trigonometric moment estimation, numerically quite appealing. Finally we illustrate our results by several real life axial data sets. Received: September 2004/ Revised: December 2004     

Gravity currents with double stratification: a numerical and analytical investigation

We consider high-Reynolds-number Boussinesq gravity current and intrusion systems in which both the ambient and the propagating “current” are linearly stratified. The main focus is on a current of fixed volume released from a rectangular lock; the height ratio of the fluids $H$ , the stratification parameter of the ambient $S$ , and the internal stratification parameter of the current, $\sigma $ , are quite general. We perform two-dimensional Navier–Stokes simulation and compare the results with those of a previously-published one-layer shallow-water model. The results provide insights into the behavior of the system and enhance the confidence in the approximate model while also revealing its limitations. The qualitative predictions of the model are confirmed, in particular: (1) there is an initial “slumping” stage of propagation with constant speed $u_N$ , after which $u_N$ decays with time; (2) for fixed $H$ and $S$ , the increase of $\sigma $ causes a slower propagation of the current; (3) for some combinations of the parameters $H,S, \sigma $ the fluid released from the lock lacks initially (or runs out quickly of) buoyancy “driving power” in the horizontal direction, and does not propagate like a gravity current. There is also a fair quantitative agreement between the predictions of the model and the simulations concerning the spread of the current.     

Modelling gas release event behaviour in hazardous waste tanks

Certain radioactive waste storage tanks at the United States Department of Energy Hanford facilities continuously generate gases as a result of radiolysis and chemical reactions. The congealed sludge in these tanks traps the gases beneath it and causes the level of the waste within the tanks to rise. The waste level continues to rise until the sludge loses its buoyancy and rolls over, changing places with fluid on top. During a rollover, the trapped gases are released, resulting in a sudden drop in the waste level. This is known as a gas release event (GRE). After a GRE, the waste re-congeals and gas again accumulates, leading to another GRE. We are interested in the time between consecutive GREs. Understanding the probabilistic behaviour of the time between consecutive GREs is important because the hydrogen and nitrous oxide gases released during a GRE are flammable and the ammonia that is released is a health risk. From a safety perspective, activity around such waste tanks should be halted when a GRE is imminent. With a credible probability model for the time between consecutive GREs, we can establish time windows in which waste tank research and maintenance activities can be safely performed. We discuss the application of non-linear time series models to this problem.     

Empirical and theoretical evidence concerning the response of the earth's ice and snow cover to a global temperature increase

As a guide to the possible effects of a CO₂-induced warming on the cryosphere, we review the effects of three warm periods in the past, and out theoretical understanding of fluctuations in mountain glaciers, the Greenland and Antarctic ice sheets, ground ice, sea ice and seasonal snow cover. Between 1890 and 1940 A.D. the glaciated area in Switzerland was reduced by over 25%. In the Hypsithermal, at about 6000 BP, ground ice in Eurasia retreated northward by several hundred kilometres. In the interglacial Stage 5e, at about 120 000 BP, global sea-level rose by over 6 m. Fluctuations of mountain glaciers depend on mesoscale “weather” and on their mechanical response to it. Any melting of the Greenland ice sheet is likely to be slow in human terms. The West Antarctic ice sheet (its base below sea-level) is susceptible to an ungrounding, and such an event may have been the cause of the sea-level rise above. The East Antarctic ice sheet is susceptible to mechanical “surges”, which might be triggered by a warming at its margin. Both an ungrounding and a surge might occupy less than 100 yr, and are potetially the most important ice changes in human terms. Modelling studies suggest that a 5°C warming would remove the Arctic pack ice in summer, and this may be the most significant effect for further climatic change.