Hydrometeorological variables predict fecal indicator bacteria densities in freshwater: data-driven methods for variable selection

Statistical models of microbial water quality inform risk management for water recreation. Current research focuses on resource-intensive, location-specific data collection and water quality modeling, but this approach may be cost-prohibitive for risk managers responsible for numerous recreation sites. As an alternative, we tested the ability of two data-driven models, tree regression and random forests with conditional inference trees, to select readily available hydrometeorological variables for use in linear mixed effects (LME) models predicting bacterial density. The study included the Chicago Area Waterway System (CAWS) and Lake Michigan beaches and harbors in Chicago, Illinois, at which Escherichia coli and enterococci were measured seasonally in 2007–2009. Tree regression node variables reduced data dimensionality by >50 %. Variable importance ranks from random forests were used in a forward-step selection based on R ² and root mean squared prediction error (RMSPE). We found two to three variables explained bacteria densities well relative to random forests with all variables. LME models with tree- or forest-selected variables performed reasonably well (0.335?<?R ²?<?0.658). LME models for Lake Michigan had good prediction accuracy with respect to the single sample maximum standard (72–77 %), but limited sensitivity (23–62 %). Results suggest that our alternative approach is feasible and performs similarly to more resource-intensive approaches.     

Effect of diel period and season on seining effort required to detect changes in Lake Erie beach fish assemblages

Using beach seining data from Lake Erie beaches, we examined the effect of diel period and season on two common study objectives considered sensitive to sampling effort: (1) has fish assemblage composition changed between two sampling time periods?, and (2) how many samples are necessary to provide a representative sample of the fish assemblage? Across a range of effect sizes and power levels, sample sizes were estimated for species richness, total fish abundance, round goby (Neogobius melanostomus) abundance and emerald shiner (Notropis atherinoides) abundance. For most sampling periods and effect sizes, several hundred to several thousand seine hauls were estimated to be required for detecting large-scale (e.g. basin-wide) changes in round goby and emerald shiner abundance. Depending on the indicator of interest, night seining (either spring or fall) was interpreted as the most efficient approach. For both seasons, sampling effort beyond 40 to 50 night seine hauls results in few additional species being detected. Both species accumulation curves and estimated sample sizes indicate that less effort is required to detect species richness changes than fish abundance. For most effect sizes and power levels, estimates of the number of seine hauls required to detect changes in fish abundance are unrealistically large for most monitoring programs. More practical monitoring strategies could be achieved by adopting alternative indicators (e.g. guilds), a more liberal significance level (p=0.1), a paired-site sampling design, or including reference sites.     

The occurrence of Shiga toxin-producing Escherichia coli in bathing water of the Sierra de la Ventana region, Buenos Aires Province, Argentina

The region of Sierra de la Ventana is located in the southwest of Buenos Aires Province, Argentina. Traditionally, this area has been devoted to livestock and agriculture, but tourism has had a significant development in recent years. In the region, there are many rivers and streams that are used for swimming and bathing. A survey of the occurrence of Shiga toxin-producing Escherichia coli (STEC) in these waters was conducted, and the microbiological quality of rivers and streams was investigated. No E. coli O157 was recovered by immunomagnetic separation. Nevertheless, the Shiga toxin gene, exclusively stx2 genotype, was detected in four non-O157 E. coli strains. Two STEC strains carried eae factor, but none of them harbored the EHEC-hlyA gene. Three of the STEC isolates belonged to samples obtained in the warm months, and one to the winter sampling. In the sample sites where STEC strains were isolated the counts of E. coli/100 ml exceeded or were close to the limit recommended by the United States Environmental Protection Agency for bathing water. The relationship observed between the rainy season and E. coli counts suggests that among the main causes for the hygienic indicator increase is the runoff of manure deposited on soils that may also induce the entrance of pathogens into the aquatic environment. This research, the first reporting STEC isolation from recreational waters in this area, revealed that streams and rivers from a beef-producing area of Argentina are a reservoir of STEC strains.     

Description and evaluation of a sampling system for monitoring hydrocarbons in sediments

A composite random sampling design was used to estimate the concentrations of hydrocarbons in sediments from two near-shore areas of Scotland (Firth of Clyde and Firth of Forth). The aim of this work was to estimate a mean value for each parameter in these areas, and to determine whether this can be done with more thorough coverage (better representation), better precision and less variance at lower analytical cost through a composite random sampling scheme rather than a simple random sampling scheme, and thereby contribute to the re-design of the UK National Marine Monitoring Programme (NMMP), re-named the UK Clean Seas Environmental Monitoring Programme (CSEMP) in 2006. Samples were collected using a simple random sampling design during 2005. All sediment samples were analysed for their particle size distribution and total organic carbon (TOC). All sediments were analysed for polycyclic aromatic hydrocarbons (PAHs) and n-alkanes. The concentrations of PAHs and n-alkanes in the study areas are described, and sources of PAHs were investigated through the PAH distributions and n-alkane profiles. Individual sediment samples from each area were combined to give a series of composite sub-samples, each comprised of 5 individual sediment samples. These composite samples were re-analysed for the same parameters as the individual samples. Mean total PAH (2- to 6-ring parent and branched) concentrations, based on the individual original sediment samples collected through simple random sampling, were 1858 microg kg(-1) dry weight (SE = 196 microg kg(-1) dry weight, n = 25) and 532.4 microg kg(-1) dry weight (SE = 59 microg kg(-1) dry weight, n = 25) in the Clyde and Forth, respectively. Mean total PAH concentrations of the composite samples were 1745 microg kg(-1) dry weight (SE = 121.0 microg kg(-1) dry weight, n = 5) in the Clyde and 511.6 microg kg(-1) dry weight (SE = 37.4 microg kg(-1) dry weight, n = 5) in the Forth. No significant differences were found between the mean PAH concentrations from the two sampling designs. This study demonstrated that the composite random sampling design gave a mean value with less variance than the simple random sampling design, at significantly reduced analytical effort (and cost).     

Efficient search procedures for extreme pollutant values

Extreme pollutant values are of great interest in water quality monitoring because of their frequent toxicological significance. The principal barrier to the detection of these values, however, is the cost of extensive and comprehensive monitoring. This paper demonstrates an efficient method to determine the maximum sample measurement from a finite set of sequential samples without explicitly testing them all. It is assumed that the process of sample measurement is distinct from collection and has higher costs. It is further assumed that the measurements have high positive autocorrelation.A methodology is presented based on a common industrial testing procedure referred to as composite sampling—the physical pooling or compositing of a set of sequential samples before measurement. A method known as primary first order compositing (PFOC) was found to be superior to the traditional technique of random sampling, particularly if small composite sizes are utilized.The authors are, respectively, operations research analyst, Vancouver, B. C.; Associate Professor and Chairman, Policy Analysis Division, Faculty of Commerce and Business Administration, University of British Columbia; and Associate Professor, Management Science Division, Faculty of Commerce, U.B.C.     

Evaluating performance of stormwater sampling approaches using a dynamic watershed model

Accurate quantification of stormwater pollutant levels is essential for estimating overall contaminant discharge to receiving waters. Numerous sampling approaches exist that attempt to balance accuracy against the costs associated with the sampling method. This study employs a novel and practical approach of evaluating the accuracy of different stormwater monitoring methodologies using stormflows and constituent concentrations produced by a fully validated continuous simulation watershed model. A major advantage of using a watershed model to simulate pollutant concentrations is that a large number of storms representing a broad range of conditions can be applied in testing the various sampling approaches. Seventy-eight distinct methodologies were evaluated by "virtual samplings" of 166 simulated storms of varying size, intensity and duration, representing 14 years of storms in Ballona Creek near Los Angeles, California. The 78 methods can be grouped into four general strategies: volume-paced compositing, time-paced compositing, pollutograph sampling, and microsampling. The performances of each sampling strategy was evaluated by comparing the (1) median relative error between the virtually sampled and the true modeled event mean concentration (EMC) of each storm (accuracy), (2) median absolute deviation about the median or "MAD" of the relative error or (precision), and (3) the percentage of storms where sampling methods were within 10% of the true EMC (combined measures of accuracy and precision). Finally, costs associated with site setup, sampling, and laboratory analysis were estimated for each method. Pollutograph sampling consistently outperformed the other three methods both in terms of accuracy and precision, but was the most costly method evaluated. Time-paced sampling consistently underestimated while volume-paced sampling over estimated the storm EMCs. Microsampling performance approached that of pollutograph sampling at a substantial cost savings. The most efficient method for routine stormwater monitoring in terms of a balance between performance and cost was volume-paced microsampling, with variable sample pacing to ensure that the entirety of the storm was captured. Pollutograph sampling is recommended if the data are to be used for detailed analysis of runoff dynamics.     

Dietary exposure of three passerine species to PCDD/DFs from the Chippewa, Tittabawassee, and Saginaw River floodplains, Midland, Michigan, USA

Dietary exposure of house wrens (Troglodytes aedon), tree swallows (Tachycineta bicolor), and eastern bluebirds (Sialia sialis) to polychlorinated dibenzofurans (PCDFs) and polychlorinated dibenzo-p-dioxins (PCDDs) near Midland, Michigan (USA) was evaluated based on site-specific data, including concentrations of residues in bolus samples and individual invertebrate orders and dietary compositions by study species. Site-specific dietary compositions for the three species were similar to those reported in the literature, but differed in their relative proportions of some dietary items. Oligocheata (non-depurated) and Brachycera (Diptera) contained the greatest average concentrations of ??PCDD/DFs of the major site-specific dietary items collected via food web-based sampling. Average ingestion values of ??PCDD/DFs from site-specific bolus-based and food web-based dietary concentrations for nestlings at study areas (SAs) were 6- to 20-fold and 2- to 9-fold greater than at proximally located reference areas (RAs), respectively. Average ingestion values of total 2,3,7,8-tetrachlorodibenzo-p-dioxin equivalents (TEQ_WHO???Avian) from site-specific bolus-based and food web-based dietary concentrations for nestlings at SAs were 31- to 121-fold and 9- to 64-fold greater than at proximally located RAs, respectively. Estimates of ??PCDD/DFs and TEQ_WHO???Avian tissue concentrations based on nestling dietary exposures were greater than those measured. Plausible explanations include nestling metabolism of 2,3,7,8-tetrachlorodibenzofuran and assimilation rates of less than the 70% assumed to occur over the nestling growth period. Profiles of the relative concentrations of individual PCDD/DF congeners in samples of invertebrates and bolus at SAs on the Tittabawassee River downstream of the source of contamination were dominated by 1,2,3,4,6,7,8,9-octachlorodibenzo-p-dioxin (22% to 44%) and 2,3,7,8-tetrachlorodibenzofuran (18% to 50%).     

Using variance components to estimate power in a hierarchically nested sampling design

We used variance components to assess allocation of sampling effort in a hierarchically nested sampling design for ongoing monitoring of early life history stages of the federally endangered Devils Hole pupfish (DHP) (Cyprinodon diabolis). Sampling design for larval DHP included surveys (5 days each spring 2007–2009), events, and plots. Each survey was comprised of three counting events, where DHP larvae on nine plots were counted plot by plot. Statistical analysis of larval abundance included three components: (1) evaluation of power from various sample size combinations, (2) comparison of power in fixed and random plot designs, and (3) assessment of yearly differences in the power of the survey. Results indicated that increasing the sample size at the lowest level of sampling represented the most realistic option to increase the survey’s power, fixed plot designs had greater power than random plot designs, and the power of the larval survey varied by year. This study provides an example of how monitoring efforts may benefit from coupling variance components estimation with power analysis to assess sampling design.     

Estimating means and variances: The comparative efficiency of composite and grab samples

This paper compares the efficiencies of two sampling techniques for estimating a population mean and variance. One procedure, called grab sampling, consists of collecting and analyzing one sample per period. The second procedure, called composite sampling, collectsn samples per period which are then pooled and analyzed as a single sample. We review the well known fact that composite sampling provides a superior estimate of the mean. However, it is somewhat surprising that composite sampling does not always generate a more efficient estimate of the variance. For populations with platykurtic distributions, grab sampling gives a more efficient estimate of the variance, whereas composite sampling is better for leptokurtic distributions. These conditions on kurtosis can be related to peakedness and skewness. For example, a necessary condition for composite sampling to provide a more efficient estimate of the variance is that the population density function evaluated at the mean (i.e.f()) be greater than . If , then a grab sample is more efficient. In spite of this result, however, composite sampling does provide a smaller estimate of standard error than does grab sampling in the context of estimating population means.     

Distributions, sources and risk assessment of polycyclic aromatic hydrocarbons (PAHs) in topsoil at Ji’nan city, China

Surface soil (0-5 cm) samples from 17 sampling sites including different functional areas at Ji'nan city in Shandong Province of China were collected and analyzed for 16 EPA priority polycyclic aromatic hydrocarbons (PAHs). The total PAH concentrations were in the range from 1.31 mg kg(-1) to 254.08 mg kg(-1) (dry weight), and the average level of total PAHs was 23.25 mg kg(-1). The highest total PAHs concentrations were found in steel and iron plant at industrial areas. The total PAHs concentrations in industrial areas were markedly higher than those in other different functional areas. According to comparing total PAHs concentration in Ji'nan city to that of other urban areas, it was found that total PAHs concentrations were 6 to 137 times higher than other areas because of some specific sampling sites such as steel and iron plant and one main roadside. The results showed that PAHs in topsoil of Ji'nan city were suffered from strong pyrogenic influence, especially in industrial areas. However about 52.9% soil samples were mainly originated from both pyrogenic and petrogenic mixed sources based on Flu/Pyr ratios and Phe/Ant ratios. Furthermore, It was found that all individual PAHs except Fle were significantly correlated (P < 0.01) with LMW, HMW, total PAHs and SOM, and individual PAHs except Fle in soils were significantly correlated (P < 0.01) with each other. The nemerow composite index to assess the environmental quality showed that the soil sample of steel and iron plant in industrial areas and one main roadside were heavy pollution of PAHs, and about 47% soil sampling sites were safety, about 53% soil sampling sites were got different grades of PAHs pollution.     

Statistical control of hygienic quality of bathing water

In Denmark the hygienic quality of the bathing water has been controlled, based on general guidelines since 1978. Today more than 1100 control sites in marine areas have been established to ensure safe bathing water quality. According to EEC directives and Danish tradition, the control is usually performed by measuring the content of the indicator bacteria Escherichia coli (E. coli) in 5 to 20 test samples each bathing season.In Denmark, control programmes and data are evaluated using basic statistical quality control principles. This paper presents general guidelines for bacterial control, their statistical background and practical application. Furthermore, the evaluation and application of a specific programme for control of Staphylococcus aureus (S. aureus) are presented. This programme was used in a Danish bay where the authorities prescribe direct control of these potentially harmful bacteria.     

Sampling Benthic Macroinvertebrates in a Large Flood-Plain River: Considerations of Study Design, Sample Size, and Cost

Estimation of benthic macroinvertebrate populations over large spatial scales is difficult due to the high variability in abundance and the cost of sample processing and taxonomic analysis. To determine a cost-effective, statistically powerful sample design, we conducted an exploratory study of the spatial variation of benthic macroinvertebrates in a 37 km reach of the Upper Mississippi River. We sampled benthos at 36 sites within each of two strata, contiguous backwater and channel border. Three standard ponar (525 cm²) grab samples were obtained at each site ('Original Design'). Analysis of variance and sampling cost of strata-wide estimates for abundance of Oligochaeta, Chironomidae, and total invertebrates showed that only one ponar sample per site ('Reduced Design') yielded essentially the same abundance estimates as the Original Design, while reducing the overall cost by 63%. A posteriori statistical power analysis ( = 0.05, = 0.20) on the Reduced Design estimated that at least 18 sites per stratum were needed to detect differences in mean abundance between contiguous backwater and channel border areas for Oligochaeta, Chironomidae, and total invertebrates. Statistical power was nearly identical for the three taxonomic groups. The abundances of several taxa of concern (e.g., Hexagenia mayflies and Musculium fingernail clams) were too spatially variable to estimate power with our method. Resampling simulations indicated that to achieve adequate sampling precision for Oligochaeta, at least 36 sample sites per stratum would be required, whereas a sampling precision of 0.2 would not be attained with any sample size for Hexagenia in channel border areas, or Chironomidae and Musculium in both strata given the variance structure of the original samples. Community-wide diversity indices (Brillouin and 1-Simpsons) increased as sample area per site increased. The backwater area had higher diversity than the channel border area. The number of sampling sites required to sample benthic macroinvertebrates during our sampling period depended on the study objective and ranged from 18 to more than 40 sites per stratum. No single sampling regime would efficiently and adequately sample all components of the macroinvertebrate community.     

Representativity of a Mid-lake Surface Water Chemistry Sample

A mid-lake sample is the conventional sampling strategy used in lake chemistry monitoring programs. Hence, addressing the question of how representative a mid-lake sample is of in-lake conditions is important for interpreting changes in long-term monitoring programs. Our hypotheses were that; (i) the location of a sampling site within a lake will affect its representativity and (ii) a mid-lake sample could be considered representative of average surface water physico-chemistry. We focused on the surface water chemistry in 34 Swedish lakes. To obtain a wide diversity of lake types, the lakes were classified according to mean depth, water residence time, and three main ecoregions. Two sampling seasons were selected, representative of high and low biological activity and stratified and non-stratified conditions, respectively. Ten samples were collected in each lake. The mid-lake samples differed significantly from the nine remaining samples in less than 4% of the total number of two-sample t-tests performed. Samples collected close to lake inlets differed from the other samples more often than samples collected in the central parts of the lakes. We concluded that the location of a sampling site within a lake may affect the sample's representativity of lake water chemistry, and that a mid-lake surface water sample could be considered representative of surface water chemistry.     

Assessing temporal representativeness of water quality monitoring data

The effectiveness of different monitoring methods in detecting temporal changes in water quality depends on the achievable sampling intervals, and how these relate to the extent of temporal variation. However, water quality sampling frequencies are rarely adjusted to the actual variation of the monitoring area. Manual sampling, for example, is often limited by the level of funding and not by the optimal timing to take samples. Restrictions in monitoring methods therefore often determine their ability to estimate the true mean and variance values for a certain time period or season. Consequently, we estimated how different sampling intervals determine the mean and standard deviation in a specific monitoring area by using high frequency data from in situ automated monitoring stations. Raw fluorescence measurements of chlorophyll a for three automated monitoring stations were calibrated by using phycocyanin fluorescence measurements and chlorophyll a analyzed from manual water samples in a laboratory. A moving block bootstrap simulation was then used to estimate the standard errors of the mean and standard deviations for different sample sizes. Our results showed that in a temperate, meso-eutrophic lake, relatively high errors in seasonal statistics can be expected from monthly sampling. Moreover, weekly sampling yielded relatively small accuracy benefits compared to a fortnightly sampling. The presented method for temporal representation analysis can be used as a tool in sampling design by adjusting the sampling interval to suit the actual temporal variation in the monitoring area, in addition to being used for estimating the usefulness of previously collected data.     

Baseline ambient gaseous ammonia concentrations in the Four Corners area and eastern Oklahoma, USA

Ambient ammonia monitoring using Ogawa passive samplers was conducted in the Four Corners area and eastern Oklahoma, USA during 2007. The resulting data will be useful in the multipollutant management of ozone, nitrogen oxides, and visibility (atmospheric regional haze) in the Four Corners area, an area with growing oil/gas production and increasing coal-based power plant construction. The passive monitoring data also add new ambient ammonia concentration information for the U.S. and will be useful to scientists involved in present and future visibility modeling exercises. Three week integrated passive ammonia samples were taken at five sites in the Four Corners area and two sites in eastern Oklahoma from December, 2006 through December, 2007 (January, 2008 for two sites). Results show significantly higher regional background ammonia concentrations in eastern Oklahoma (1.8 parts per billion (ppb) arithmetic mean) compared to the Four Corners area (0.2 ppb arithmetic mean). Annual mean ammonia concentrations for all Four Corners area sites for the 2007 study ranged from 0.2 ppb to 1.5 ppb. Peak ambient ammonia concentrations occurred in the spring and summer in both areas. The passive samplers deployed at the Stilwell, Oklahoma site compared favorably with other passive samplers and a continuous ammonia monitoring instrument.     

Strategically designed sample composition for fastest screening of polychlorinated biphenyl congeners in water samples

A new strategy for preparing composite samples of special interest and applicability in environmental screening studies is presented. The use of supersaturated experimental design matrices to conduct the sample composition of water sample specimens in screening studies is demonstrated. In contrast to well known conventional sample composition, this strategic approach provides analytical objects allowing the accurate prediction of analyte concentration levels in the original individual sample specimens while fixing the number of experiments to be carried out down to the very number of sample specimens. This will be of special importance when dealing with analytes that require complicated, labour intensive and expensive analytical processes. To reach this goal, two main conditions must be fulfilled. The first one is the sparsity effect (Pareto principle) which holds for the specimens in the sampling campaign. This means that the number of really anomalous or contaminated specimens, as compared to the total number of specimens to be analysed, is low. In environmental screening studies, frequently this situation can be reasonably assumed. The second condition is to have an effective manner to develop and solve the experimental designs required to build-up the composite samples. The challenging problem of screening PCBs in water samples has been tackled to show the usefulness of this strategic approach by combining chemometrically assisted sample composition and rapid analysis using solid-phase microextraction of the composite samples.     

A simplified sampling procedure for the estimation of methane emission in rice fields

Manual closed chamber methods are widely used for CH₄ measurement from rice paddies. Despite diurnal and seasonal variations in CH₄ emissions, fixed sampling times, usually during the day, are used. Here, we monitored CH₄ emission from rice paddies for one complete rice-growing season. Daytime CH₄ emission increased from 0800 h, and maximal emission was observed at 1200 h. Daily averaged CH₄ flux increased during plant growth or fertilizer application and decreased upon drainage of plants. CH₄ measurement results were linearly interpolated and matched with the daily averaged CH₄ emission calculated from the measured results. The time when daily averaged emission and the interpolated CH₄ curve coincided during the daytime was largely invariant within each of the five distinctive periods. One-hourly sampling during each of these five periods was utilized to estimate the emission during each period, and we found that five one-hourly samples during the season accurately reflected the CH₄ emission calculated based on all 136 hourly samples. This new sampling scheme is simple and more efficient than current sampling practices. Previously reported sampling schemes yielded estimates 9 to 32% higher than the measured CH₄ emission, while our suggested scheme yielded an estimate that was only 5% different from that based on all 136-h samples. The sampling scheme proposed in this study can be used in rice paddy fields in Korea and extended worldwide to countries that use similar farming practices. This sampling scheme will help in producing more accurate global methane budget from rice paddy fields.     

Pollution levels of thiobencarb, propanil, and pretilachlor in rice fields of the muda irrigation scheme, Kedah, Malaysia

The purpose of this study was to investigate the potential risk of pretilachlor, thiobencarb, and propanil pollutants in the water system of the rice fields of the Muda area. The study included two areas that used different irrigation systems namely non-recycled (N-RCL) and recycled (RCL) water. Regular water sampling was carried out at the drainage canals during the weeding period from September to October 2006 in the main season of 2006/2007 and April-May 2007 in off season of 2007. The herbicides were extracted by the solid-phase extraction method and identified using a GC-ECD. Results showed that the procedure for identification of the three herbicides was acceptable based on the recovery test values, which ranged from 84.1% to 96.9%. A wide distribution pattern where more than 79% of the water samples contained the herbicide pollutants was observed at both the areas where N-RCL and RCL water was supplied for the two seasons. During September to October 2006, high weedicide residue concentration was observed at the N-RCL area and it ranged from 0.05 to 1.00 μg/L for pretilachlor and propanil and 10-25 μg/L for thiobencarb. In the case of the area with RCL water, the weedicide residue ranged from 1 to 5 μg/L for pretilachlor and propanil and 10-25 μg/L for thiobencarb. The highest residue level reached was 25-50, 50-100, and 100-200 μg/L for pretilachlor, propanil, and thiobencarb, respectively. During April to May 2007, high residue concentration frequently occurred at the area supplied with N-RCL irrigation water and it ranged from 0.05 to 1.00, 10 to 25, and 25 to 50 μg/L for pretilachlor, propanil, and thiobencarb, respectively. The highest residue level reached was 25-50 μg/L for pretilachlor and 100-200 μg/L for propanil and thiobencarb. There was an accelerated increase in the concentration of the herbicide residues, with the maximum levels reached at the early period of weedicide application, followed by a sharp decrease after the rice fields were completely covered with the rice crop. During the main season of 2006/2007, the concentration of propanil residue gradually rose, although that of the other herbicides declined.     

Assessment of water quality using cluster analysis in coastal region of Mumbai, India

The coastal water quality of Mumbai is deteriorating due to various point and non-point wastewater sources. Hence, it is desirable to monitor coastal water quality for various water-related activities like bathing, contact water sports, recreation, and commercial fishing. The objective of this paper is to assess the seasonal water quality on the basis of seawater standards. Based on water-quality analysis of 17 seafronts and beaches, most of the parameters were exceeding the standards. The statistical cluster analysis was carried out for evaluating impact of wastewater and sewage discharges. The hierarchical cluster analysis resulted into three clustered groups, namely less polluted, moderately polluted, and highly polluted sites with similar characteristics of water quality. Mahim was found to be worst-affected beach due to incoming organic load from the Mithi river in comparison to other seafronts and beaches. Unaccounted sources of sewage and wastewater should be identified and rerouted through sewerage system by improving collection efficiency, treatment, and proper disposal for achieving designated receiving water quality standards.     

Speciation and antimicrobial resistance of Enterococci isolated from recreational beaches in Malaysia

We report the first study on the occurrence of antibiotic-resistant enterococci in coastal bathing waters in Malaysia. One hundred and sixty-five enterococci isolates recovered from two popular recreational beaches in Malaysia were speciated and screened for antibiotic resistance to a total of eight antibiotics. Prevalence of Enterococcus faecalis and Enterococcus faecium was highest in both beaches. E. faecalis/E. faecium ratio was 0.384:1 and 0.375:1, respectively, for isolates from Port Dickson (PD) and Bagan Lalang (BL). Analysis of Fisher’s exact test showed that association of prevalence of E. faecalis and E. faecium with considered locations was not statistically significant (p?<?0.05). Chi-square test revealed significant differences (χ ²?=?82.630, df?=?20, p?<?0.001) in the frequency of occurrence of enterococci isolates from the considered sites. Resistance was highest to nalidixic acid (94.84 %) and least for chloramphenicol (8.38 %). One-way ANOVA using Tukey–Kramer multiple comparison test showed that resistance to ampicillin was higher in PD beach isolates than BL isolates and the difference was extremely statistically significant (p?<?0.0001). Frequency of occurrence of multiple antibiotic resistance (MAR) isolates were higher for PD beach water (64.29 %) as compared to BL beach water (13.51 %), while MAR indices ranged between 0.198 and 0.48. The results suggest that samples from Port Dickson may contain MAR bacteria and that this could be due to high-risk faecal contamination from sewage discharge pipes that drain into the sea water.