Gilthead seabream (Sparus aurata) carrying antibiotic resistant enterococci. A potential bioindicator of marine contamination?

Antibiotic resistance in bacteria is a growing problem that is not only restricted to the clinical setting but also to other environments such as marine species that harbor antibiotic resistant bacteria and therefore may serve as reservoirs for antibiotic-resistance genetic determinants. The aim of this study was to evaluate antibiotic resistance phenotypes in enterococci isolated from fecal samples of gilthead seabream and the associated mechanisms of resistance. A collection of 118 samples were analyzed and 73 enterococci were recovered. The strains showed high percentages of resistance to erythromycin and tetracycline (58.9% and 17.8%, respectively). Lower level of resistance (<13%) was detected for quinupristin-dalfopristin, ampicillin, high-level-gentamicin, high-level-streptomycin, high-level-kanamycin, ciprofloxacin and chloramphenicol. The erm(B), tet(L) or tet(M), aac(6′)-aph(2″) and aph(3′)-IIIa genes were shown in isolates resistant to erythromycin, tetracycline, high-level gentamicin and high-level kanamycin, respectively. Antibiotic resistance in natural microbiota is becoming a concern of human and environmental health.     

Spatial distribution and ecology of the Recent Ostracoda from Tangra Yumco and adjacent waters on the southern Tibetan Plateau: A key to palaeoenvironmental reconstruction

We elucidate the ecology of Recent Ostracoda from a deep brackish lake, Tangra Yumco (30°45′—31°22′N and 86°23′—86°49′E, 4595 m a.s.l.) and adjacent waters on the southern Tibetan Plateau. Ostracod associations (living and empty valves) in sixty-six sediment samples collected from diverse aquatic habitats (lakes, estuary-like water and lagoon-like water waters, rivers, ponds and springs) were quantitatively assessed.Eleven Recent Ostracoda were found (nine living and two as empty valves only). Cluster analysis established two significant (p < 0.05) habitat specific associations; (i) Leucocytherella sinensis, Limnocythere inopinata, Leucocythere? dorsotuberosa, Fabaeformiscandona gyirongensis and Candona xizangensis are lacustrine fauna. (ii) Tonnacypris gyirongensis, Candona candida, Ilyocypris sp., Heterocypris incongruens and Heterocypris salina are temporary water species.Ostracod distribution and abundance are significantly (p < 0.05) correlated to physico-chemical variables. The first two axes of a canonical correspondence analysis (CCA) explain 30.9% of the variation in the species abundance data. Conductivity and habitat types are the most influential ecological factors explaining the presence and abundance of ostracods. Spearman correlation analysis reveals that: (i) Two species, L.? dorsotuberosa (r = 0.25) and L. inopinata (r = 0.36) have a significant positive correlation with conductivity while one species, T. gyirongensis (r = −0.68) displays a significant negative correlation with conductivity. Limnocythere inopinata correlates significantly positive (r = 0.37) with alkalinity. Fabaeformiscandona gyirongensis correlates significantly positive (r = 0.28) with water depth.Key indicator living assemblages are: (i) L. sinensis dominates Ca-depleted brackish waters although ubiquitously distributed; (ii) L.? dorsotuberosa dwells in fresh to brackish waters; (iii) L. inopinata predominates in mesohaline to polyhaline waters; (iv) F. gyirongensis inhabits exclusively brackish-lacustrine deeper waters; (v) C. candida populates freshwaters; (vi) T. gyirongensis and Ilyocypris sp. are restricted to shallow temporary waters; (vii) H. incongruens occurs in ponds.Water depth indicators are F. gyirongensis and L.? dorsotuberosa, useful in ostracod assemblages for palaeo-water depth reconstruction.Our results expand the knowledge of the ecological significance of Recent Tibetan Ostracoda ecology. This is a new insight on habitat chacteristics of both living assemblages and sub-Recent associations of ostracods in mountain aquatic ecosystems. The new modern ostracod dataset can be used for the quantitative reconstruction of past environmental variables (e.g., conductivity) and types of water environment. The key indicator ostracods are relevant in palaeolimnological and climate research on the Tibetan Plateau.     

Molecular characterizations of chloramphenicol- and oxytetracycline-resistant bacteria and resistance genes in mariculture waters of China

In order to gain an understanding of the diversity and distribution of antimicrobial-resistant bacteria and their resistance genes in maricultural environments, multidrug-resistant bacteria were screened for the rearing waters from a mariculture farm of China. Both abalone Haliotis discushannai and turbot Scophthalmus maximus rearing waters were populated with abundant chloramphenicol-resistant bacteria. These bacteria were also multidrug resistant, with Vibriosplendidus and Vibriotasmaniensis being the most predominant species. The chloramphenicol-resistance gene cat II, cat IV or floR could be detected in most of the multidrug-resistant isolates, and the oxytetracycline-resistance gene tet(B), tet(D), tet(E) or tet(M) could also be detected for most of the isolates. Coexistence of chloramphenicol- and oxytetracycline-resistance genes partially explains the molecular mechanism of multidrug resistance in the studied maricultural environments. Comparative studies with different antimicrobial agents as the starting isolation reagents may help detect a wider diversity of the antimicrobial-resistant bacteria and their resistance genes.     

The use of Bayesian methods for fitting rating curves,with case studies

Using data from two very large watersheds and five smaller, this paper explores the use of Bayesian methods for fitting rating curves. Posterior distribution of rating-curve parameters were calculated using Markov Chain Monte Carlo (MCMC) methods, and 95% credible intervals were calculated for predicted discharges, given stage. Expected discharge was related to stage using a link function. For the five smaller watersheds, the assumptions were (a) that the distribution of discharge Q, given stage h, is Normal, with variance proportional to h; (b) that a log link function relates μ_Q, the mean of Q∣h, to a function of stage, of the form μ_Q = β(h + α)^λ. For the two large watersheds, however, a better fit was obtained by taking the distribution of Q to be log-Normal, and the link function as ln μ_Q = β₀ + β₁h. For the two large watersheds, priors for all three parameters were taken as uninformative; for the five smaller, the prior for parameter λ was taken as Normally distributed, N(2, 0.5). Acceptable ratings were obtained for all seven sites. It is argued that distributions of derived variables (such as annual maximum discharge) can be derived directly from (a) the posterior distribution of rating-curve parameters, and (b) the stage record, without recourse to additional assumptions. Estimates thus obtained for the T-year event will incorporate rating-curve uncertainty. It is argued that Bayesian methods are appropriate for rating-curve calculation because their inherent flexibility (a) allows the incorporation of prior information about the nature of a rating curve; (b) yields credible intervals for predicted discharges and quantities derived from them; (c) can be extended to allow for uncertainty in stage measurements.     

Importance of reversible attachment in predicting E. coli transport in saturated aquifers from column experiments

Drinking water wells indiscriminatingly placed adjacent to fecal contaminated surface water represents a significant but difficult to quantify health risk. Here we seek to understand mechanisms that limit the contamination extent by scaling up bacterial transport results from the laboratory to the field in a well constrained setting. Three pulses of Escherichia coli originating during the early monsoon from a freshly excavated pond receiving latrine effluent in Bangladesh were monitored in 6 wells and modeled with a two-dimensional (2-D) flow and transport model conditioned with measured hydraulic heads. The modeling was performed assuming three different modes of interaction of E. coli with aquifer sands: (1) irreversible attachment only (best-fit k_i = 7.6 day⁻¹); (2) reversible attachment only (k_a = 10.5 and k_d = 0.2 day⁻¹); and (3) a combination of reversible and irreversible modes of attachment (k_a = 60, k_d = 7.6, k_i = 5.2 day⁻¹). Only the third approach adequately reproduced the observed temporal and spatial distribution of E. coli, including a 4-log₁₀ lateral removal distance of ∼9 m. In saturated column experiments, carried out using aquifer sand from the field site, a combination of reversible and irreversible attachment was also required to reproduce the observed breakthrough curves and E. coli retention profiles within the laboratory columns. Applying the laboratory-measured kinetic parameters to the 2-D calibrated flow model of the field site underestimates the observed 4-log₁₀ lateral removal distance by less than a factor of two. This is promising for predicting field scale transport from laboratory experiments.     

Seasonal and diurnal variation of electron and iron concentrations at the meteor heights above Arecibo

We report observations of seasonal and local time variation of the averaged electron and iron concentrations, as well as simultaneous measurements of the two species, above the Arecibo Observatory (18.35°N, 66.75°N), Puerto Rico. The average Fe profile between 21:00 and 24:00 LT has a single peak at about 85 km with the exception of the summer when an additional peak exists at about 95 km. The higher Fe peak in the summer is correlated with higher electron concentrations in this season. The three nights of simultaneous measurements of electron and iron concentrations show that narrow layers of Fe and electrons are well correlated. Comparison of the climatological and simultaneous Fe and electron data suggests that recombination of Fe⁺ plays an important role in determining the Fe profile in the upper part of the Fe layer. Above 93 km, the Fe concentration appears to increase after sunset if the electron concentration exceeds about 4000 electrons cm⁻³. The average rate of Fe production is about 0.1 atom cm⁻³ s⁻¹ for all seasons at 100 km in the early evening hours. A chemical model reveals that the concentration of Fe⁺ must be 50–80% of the total ionization over Arecibo for typical equinox conditions to explain the observed rate of Fe production. These high relative Fe⁺ concentrations are consistent with in situ observations that Fe⁺ is usually the dominant ion in sporadic E layers in the nighttime lower E region. This suggests that the source of Fe⁺ is provided by sporadic E layers descending over Arecibo after sunset. The Fe density between 80 and 85 km decreases during the night, for all seasons. This is attributed to the formation of stable molecular Fe species, such as FeOH, due to the increase in O₃ and decrease in atomic O and H during the night at these altitudes.     

The delayed occurrence of equatorial ionospheric F2 layer post-sunset height decreases following auroral-zone substorm onsets

Nighttime F₂ layer height decreases have been examined for post-sunset intervals at the equatorial station, Huancayo. The analyses involved mainly R_z max years, (1957–1960), although R_z min years (1974–1977) were also used. The data were obtained from tabulations of the ionogram parameters h'F, f_○F₂ and f_○E_s. The height reductions are delayed by more than 7 h following geomagnetic substorm onsets for locations at longitudes to the east. The reduced occurrence of spread-F and f_○E_s enhancements is found to be associated. These enhancements are also recorded at the European station, Dourbes about 4.5 h before the Huancayo enhancements. It is proposed that an LS-TID which propagates in the 80 km sound channel may be involved. Also, the experimental evidence suggests that the westward propagation of the LS-TIDs allows equatorial disturbances to occur at local times (before midnight) which are similar to the times when the LS-TIDs are generated.     

Coda-Q and its lapse time dependence analysis in the interaction zone of the Dinarides,the Alps and the Pannonian basin

The single backscattering model was used to estimate total attenuation of coda waves (Q_c) of local earthquakes recorded on eight seismological stations in the complex area of the western continental Croatia. We estimated Q₀ and n, parameters of the frequency dependent coda-Q using the relation Q_c = Q₀fⁿ. Lapse time dependence of these parameters was studied using a constant 30 s long time window that was slid along the coda of seismograms. Obtained Q_c were distributed into classes according to their lapse time, t_L. For t_L = 20–50 s we estimated Q₀ = 45–184 and n = 0.49–0.94, and for t_L = 60–100 s we obtained Q₀ = 119–316 and n = 0.37–0.82. There is a tendency of decrease of parameter n with increasing Q₀, and vice versa. The rates of change of both Q₀ and n seem to decrease for lapse times larger than 50–80 s, indicating an alteration in rock properties controlling coda attenuation at depths of about 100–160 km. A very good correlation was found between the frequency dependence parameter n and the Moho depths for lapse times of 50, 60 and 70 s.     

Evaluation and modeling of runoff and sediment yield for different land covers under simulated rain in a semiarid region of Brazil

This paper quantifies the runoff and sediment yield for four different land covers in a semiarid region of Brazil. The WESP model, a distributed, event-oriented runoff-erosion model, was applied and its physical parameters, N_s and K_R, were adjusted based on observed runoff and sediment yield data using simulated rainfall with an average intensity of 53 mm h^-1. The sediment yield obtained was 53.02 kg ha^-1 (caatinga vegetation), 231.96 kg ha^-1 (bare soil), 309.75 kg ha^-1 (beans), and 847.38 kg ha^-1(corn). The results showed that caatinga cover yields the lowest erosion and runoff when compared to the other treatments. The results also show that the sediment yield and runoff values simulated with N_s, K_I, and K_R parameters were well calibrated, within acceptable deviations. The caatinga vegetation was more effective in protecting the soil, when compared to the other types of coverage. The beans and corn covers had the highest values of runoff and sediment yield, even higher than those observed for bare soil.     

Mercury emissions and stable isotopic compositions at Vulcano Island (Italy)

Sampling and analyses methods for determining the stable isotopic compositions of Hg in an active volcanic system were tested and optimized at the volcanic complex of Vulcano (Aeolian Islands, Italy). Condensed gaseous fumarole Hg_(fum)^T, plume gaseous elemental Hg_(g)⁰ and plume particulate Hg_(p)^II were obtained at fumaroles F0, F5, F11, and FA. The average total Hg emissions, based on Hg^T/SO₂ in condensed fumarolic gases and plumes, range from 2.5 to 10.1 kg y^? 1, in agreement with published values [Ferrara, R., Mazzolai, B., Lanzillotta, E., Nucaro, E., Pirrone, N., 2000. Volcanoes as emission sources of atmospheric mercury in the Mediterranean Basin. Sci. Total Environ. 259(1–3), 115–121; Aiuppa, A., Bagnato, E., Witt, M.L.I., Mather, T.A., Parello, F., Pyle, D.M., Martin, R.S., 2007. Real-time simultaneous detection of volcanic Hg and SO₂ at La Fossa Crater, Vulcano (Aeolian Islands, Sicily). Geophys. Res. Lett. 34(L21307).]. Plume Hg_(p)^II increases with distance from the fumarole vent, at the expense of Hg_(g)⁰ and indicates significant in-plume oxidation and condensation of fumarole Hg_(fum)^T.Relative to the NIST SRM 3133 Hg standard, the stable isotopic compositions of Hg are δ²⁰²Hg_(fum)^T = ? 0.74‰ ± 0.18 (2SD, n = 4) for condensed gaseous fumarole Hg_(fum)^T, δ²⁰²Hg_(g)⁰ = ? 1.74‰ ± 0.36 (2SD, n = 1) for plume gaseous elemental Hg_(g)⁰ at the F0 fumarole, and δ²⁰²Hg_(p)^II = ? 0.11‰ ± 0.18 (2SD, n = 4) for plume particulate Hg_(p)^II. The enrichment of Hg_(p)^II in the heavy isotopes and Hg_(g)⁰ in the light isotopes relative to the total condensed fumarolic Hg_(fum)^T gas complements the speciation data and demonstrates a gas-particle fractionation occurring after the gas expulsion in ambient T° atmosphere. A first order Rayleigh equilibrium condensation isotope fractionation model yields a fractionation factor α_cond-gas of 1.00135 ± 0.00058.     

Indicators of microbial beach water quality: Preliminary findings from Teluk Kemang beach,Port Dickson (Malaysia)

This study aims to determine the concentrations of total coliforms and Escherichia coli (E. coli) in beach water, Teluk Kemang beach. This study was also aimed to determine relationship between total coliforms, E. coli and physicochemical parameters. As perceived health symptoms among beach visitors are rarely incorporated in beach water studies, this element was also assessed in this study. A total of eight water sampling points were selected randomly along Teluk Kemang beach. Total coliforms concentrations were found between 20 and 1940 cfu/100 ml. E. coli concentrations were between 0 and 90 cfu/100 ml. Significant correlations were found between total coliforms and E. coli with pH, temperature and oxidation reduction potential. Skin and eyes symptoms were the highest reported though in small numbers. Microbiological water quality in Teluk Kemang public beach was generally safe for recreational activities except sampling location near with sewage outfall.     

Potential use of mangroves as constructed wetland for municipal sewage treatment in Futian,Shenzhen, China

A pilot-scale mangrove wetland was constructed in Futian, Shenzhen for municipal sewage treatment. Three identical belts (length: 33 m, width: 3 m, depth: 0.5 m) were filled with stone (bottom), gravel and mangrove sand (surface). Seedlings of two native mangrove species (Kandelia candel, Aegiceras corniculatum) and one exotic species (Sonneratia caseolaris) were transplanted to the belts with one species for each belt. The hydraulic loading was 5 m³ d⁻¹ and hydraulic retention time 3 d. High levels of removal of COD, BOD₅, TN, TP and NH₃–N were obtained. The treatment efficiency of S. caseolaris and A. corniculatum was higher than that of K. candel. Faster plant growth was obtained for S. caseolaris. The substrate in the S. caseolaris belt also showed higher enzyme activities including dehydrogenase, cellulase, phosphatase, urease and β-glucosidase. The removal rates of organic matter and nutrients were positively correlated with plant growth. The results indicated that mangroves could be used in a constructed wetland for municipal sewage treatment, providing post-treatment to remove coliforms was also included.     

Natural antifouling compounds produced by a novel fungus Aureobasidium pullulans HN isolated from marine biofilm

A fungus, Aureobasidium pullulans, was isolated from marine biofilm and identified. A bioassay-guided fractionation procedure was developed to isolate and purify antifouling compounds from A. pullulans HN. The procedure was: fermentation broth—aeration and addition of sodium thiosulfate–graduated pH and liquid–liquid extraction—SPE purification—GC–MS analysis. Firstly, the fermentation broth was tested for its toxicity. Then it was treated with aeration and addition of sodium thiosulfate, and its toxicity was almost not changed. Lastly, antifouling compounds were extracted at different pH, the extract had high toxicity at pH 2 but almost no toxicity at pH 10, which suggested the toxicants should be fatty acids. The EC₅₀ of the extract against Skeletonema costatum was 90.9 μg ml⁻¹, and its LC₅₀ against Balanus amphitrete larvae was 22.2 μg ml⁻¹. After purified by HLB SPE column, the EC₅₀ of the extract against S. costatum was 49.4 μg ml⁻¹. The myristic and palmitic acids were found as the main toxicants by GC–MS.     

Active segment of the 12 November 2003 Mw 5.6 earthquake at Salsipuedes oceanic basin,Gulf of California. Mexico

We analyzed records of eight seismic stations of the autonomous broadband seismograph network of a joint project between Utrecht University (the Netherlands), California Institute of Technology, and Centro de Investigación Científica y de Estudios Superiores de Ensenada (CICESE). These stations recorded the Mw 5.6 earthquake that occurred on 12 November 2003 at Salsipuedes basin in the middle of the Gulf of California 2 km west of the island Angel de la Guarda. This event was located at 29.16º N and 113.37º W, 30 km northeast of Bahia de los Angeles. A foreshock and hundreds of aftershocks were recorded in the 48 hours after its origin time. With the location of 29 earthquakes we identified the active segment, perpendicular to the main transform fault NW–SE of Canal de Ballenas, representing the transtensional boundary between the Pacific and North American plates. The direction of the active fault described is consistent with the normal fault mechanism reported by the National Earthquake Information Center (strike=39º, dip=34º, slip=–44º).From the duration magnitude of 456 aftershocks, we calculated a b-value of 1.14±0.28; furthermore, we calculated a seismic moment of (3.5 ±3.3) X10¹⁷Nm, a source radius of 3.7 ± 2.63 km, and a static stress drop of 3.94 ± 1.15 MPa (39.4 ± 11.5 bar.).     

A similarity approach to determine response times to steady-state saturation in landscapes

Relationships describing response times for landscape saturation and subsurface flow for idealised hillslopes after a change in water balance are derived in terms of similarity parameters given by their topographic, soil and climatic attributes. The study was carried out under a set of specific assumptions.The work quantitatively describes effects of each of these attributes on travel time. For example, divergent hillslopes can exhibit travel times that are double those of convergent hillslopes, and concave slopes tend to have lower travel times than planar or convex slopes. It is shown that the hillslope travel time T is dependent on hillslope length L, hydraulic conductivity K and slope S such that T = f(L/KS, B, CR, Δq/smd) and if saturation does not occur in the lower reaches of a hillslope, this can be simplified to T = f(L/KS, B). Here, Δq is change in net external flux and smd is soil moisture deficit; CR and B are the convergence ratio and profile factor respectively. The relative importance of these attributes on response time is discussed.It is shown that subsurface flow hydrographs of simple hillslopes of different scale, derived numerically, can be collapsed into a single curve by normalising them by means of a form of response time, the time constant.Results from the analytical derivation were compared with numerical analysis with good agreement. The theory was applied to a natural catchment using observed data from a bushfire event. Comparison of analytically calculated response time with observed response time for the event showed a large discrepancy. The reason is discussed.     

Arsenic and heavy metals in native plants at tailings impoundments in Queretaro,Mexico

Ten native plants species that grow in three tailings dams from Ag, Pb, Cu and Zn mine in Queretaro, Mexico were studied. Total concentrations in tailings were 183–14,660 mg/kg As, 45–308 mg/kg Cd, 327–1754 mg/kg Pb, 149–459 mg/kg Cu and 448–505 mg/kg Zn. In the three tailings dams, the solubility of these elements is low. Tailings in dam 1 are acid generating while tailings in dams 2 and 3 are not acid-generating potential. Plants species that accumulate arsenic and heavy metals was identified; Nicotina glauca generally presented the highest concentrations (92 mg/kg As, 106 mg/kg Cd, 189 mg/kg Pb, 95 mg/kg Cu and 1985 mg/kg Zn). Other species that accumulate these elements are Flaveria pubescens, Tecoma stans, Prosopis Sp, Casuarina Sp and Maurandia antirrhiniflora. Two species were found that accumulates a large amount of metals in the root, Cenchrus ciliaris and Opuntia lasiacantha. Concentrations in soils in which plants grow were 488–5990 mg/kg As, 5–129 mg/kg Cd, 169–3638 mg/kg Pb, 159–1254 mg/kg Cu and 1431–13,488 mg/kg Zn. The Accumulation Factor (AF) determined for plants was less than 1, with exception of N. glauca for Cd. The correlation between arsenic and heavy metals found in soils and plants was low. Knowledge of plant characteristics allows it use in planning the reforestation of tailings dams in controlled manner. This will reduce the risk of potentially toxic elements are integrated into the food chain of animal species.     

Evapotranspiration in the Pampean Region using field measurements and satellite data

Evapotranspiration (LE) is an important factor for monitoring crops, water requirements, and water consumption at local and regional scale. In this paper, we applied the semi-empirical model to estimate the daily latent heat flux (LE_d = Rn_d + A − B(Ts − Ta)). LE_d has been estimated using satellite images (Thematic Mapper sensor) and a local dataset (incoming and outgoing short- and long-wave radiation) measured during three years. We first estimated the daily net Radiation (Rn_d) from a linear equation derived from the instantaneous net Radiation (Rn_d = CRn_i + D). Subsequently, coefficients A and B have been estimated for two different cover vegetations (pasture and soybean). For each vegetation cover, an error analysis combining Rn_d, A, B, and surface and air temperatures has been calculated. Results showed that Rn_d had good performance (nonbias and low RMSE). LE_d errors for pasture and soybean were ±28 W m⁻² and ±40 W m⁻² respectively.     

The unique environment of the most acidified permanently meromictic lake in the Czech Republic

Changes in the water properties and biological characteristics of the highly acidic Hromnice Lake (Western Bohemia) were investigated. This 110-year-old lake, formed as a consequence of the mining of pyritic shales, is permanently meromictic. Two chemoclines separate an extremely acidic (pH ～ 2.6) mixolimnion from a metal-rich anoxic monimolimnion. The absence of spring mixolimnetic turnover due to ice melting and very slow heat propagation through the chemocline with a 6-month delay were observed. Extreme mixolimnetic oxygen maxima (up to 31 mg l^?1) in phosphorus-rich lake (PO₄^3? up to 1.6 mg l^?1) well correlated with outbursts of phytoplankton. Phytoplankton consist of several acido-tolerant species of the genera Coccomyxa, Lepocinclis, Chlamydomonas and Chromulina. Surface phytoplankton biomass expressed as chlorophyll-a varies from 2 to 140 μg l^?1. Multicellular zooplankton are almost absent with the exception of Cephalodella acidophila, a small rotifer occurring in low numbers. Large red larvae of the midge Chironomus gr. plumosus were found at the bottom close to the shore, with larvulae in the open water. Developmental stages (protonemata) of a moss, resembling filamentous algae, dwell in the otherwise plant-free littoral zone.     

Balancing the global oceanic neodymium budget: Evaluating the role of groundwater

The distinctly different, ε_Nd(0) values of the Atlantic, Indian, and Pacific Oceans requires that the residence time of Nd in the ocean (i.e., τ_Nd) be on the order of, or less than, the ocean mixing time of ∼ 500–1500 yr. However, estimates of τ_Nd, based on river influxes, range from 4000 to 15,000 yr, thus exceeding the ocean mixing time. In order to reconcile the oceanic Nd budget and lower the residence time by roughly a factor of 10, an additional, as yet unidentified, and hence “missing Nd flux” to the ocean is necessary. Dissolution of materials deposited on continental margins has previously been proposed as a source of the missing flux. In this contribution, submarine groundwater discharge (SGD) is examined as a possible source of the missing Nd flux. Neodymium concentrations (n = 730) and ε_Nd(0) values (n = 58) for groundwaters were obtained from the literature in order to establish representative groundwater values. Mean groundwater Nd concentrations and ε_Nd(0) values were used along with recent estimates of the terrestrial (freshwater) component of SGD (6% of river discharge on a global basis) to test whether groundwater discharge to the coastal oceans could account for the missing flux. Employing mean Nd concentrations of the compiled data base (i.e., 31.8 nmol/kg for all 730 analyses and 11.3 nmol/kg for 141 groundwater samples from a coastal aquifer), the global, terrestrial-derived SGD flux of Nd is estimated to range between 2.9 × 10⁷ and 8.1 × 10⁷ mol/yr. These estimates are of the same order of magnitude, and within a factor of 2, of the missing Nd flux (i.e., 5.4 × 10⁷ mol/yr). Applying the SGD Nd flux estimates, the global average ε_Nd(0) of SGD is predicted to be − 9.1, which is similar to our estimate for the missing Nd flux (− 9.2), and in agreement with the mean (± S.D.) ε_Nd(0) measured in groundwaters (i.e., ε_Nd(0) = −8.9 ± 4.2). The similarities in the estimated SGD Nd flux and corresponding ε_Nd(0) values to the magnitude and isotope composition of the missing Nd flux are compelling, and suggest that discharge of groundwater to the oceans could account for the missing Nd flux. Future investigations should focus on quantifying the Nd concentrations and isotope compositions of groundwater from coastal aquifers from a variety of coastal settings, as well as the important geochemical reactions that effect Nd concentrations in subterranean estuaries in order to better constrain contributions of SGD to the oceanic Nd budget.     

Toxicities of antifouling biocide Irgarol 1051 and its major degraded product to marine primary producers

Irgarol 1051 (2-methythiol-4-tert-butylamino-6-cyclopropylamino-s-triazine) is an algaecide commonly used in antifouling paints. It undergoes photodegradation which yields M1 (2-methylthio-4-tert-butylamino-6-amino-s-triazine) as its major and most stable degradant. Elevated levels of both Irgarol and M1 have been detected in coastal waters worldwide; however, ecotoxicity effects of M1 to various marine autotrophs such as cyanobacteria are still largely unknown. This study firstly examined and compared the 96 h toxicities of Irgarol and M1 to the cyanobacterium Chroococcus minor and two marine diatom species, Skeletonema costatum and Thalassiosira pseudonana. Our results suggested that Irgarol was consistently more toxic to all of the three species than M1 (96 h EC50 values: C. minor, 7.71 μg L⁻¹ Irgarol vs. >200 μg L⁻¹ M1; S. costatum, 0.29 μg L⁻¹ Irgarol vs. 11.32 μg L⁻¹ M1; and T. pseudonana, 0.41 μg L⁻¹ Irgarol vs. 16.50 μg L⁻¹ M1). Secondly, we conducted a meta-analysis of currently available data on toxicities of Irgarol and M1 to both freshwater and marine primary producers based on species sensitivity distributions (SSDs). Interestingly, freshwater autotrophs are more sensitive to Irgarol than their marine counterparts. For marine autotrophs, microalgae are generally more sensitive to Irgarol than macroalgae and cyanobacteria. With very limited available data on M1 (i.e. five species), M1 might be less toxic than Irgarol; nonetheless this finding warrants further confirmation with additional data on other autotrophic species.