A model for estimating the potential biomagnification of chemicals in a generic food web: Preliminary development

GOAL, SCOPE AND BACKGROUND: Bioaccumulation and biomagnification of organic pollutants have been increasingly assessed and modeled during the last years. Due to the complexity of these processes and the large variability of food webs, setting generic assessments for these parameters is really difficult. Equilibrium models, based on a compound's lipophylicity, are the main tool in regulatory proposals, such as for identifying Persistent, Bioaccumulative and Toxic Substances (PBTs), although a refinement has been claimed by the scientific community. Toxicokinetic studies offer an alternative for these estimations, where biomagnification is modeled as a succession of bioaccumulation processes, each one regulated by toxicokinetic parameters. METHODS: A review of kinetic models covering species belonging to different trophic levels and with different ecological behavior has been conducted. The results were employed for setting a conceptual model for estimating the biomagnification potential in a generic food web, which was mathematically implemented through system dynamic models developed under data sheet software. Crystal Ball was then employed for allowing Monte Carlo based probabilistic calculations. Bioaccumulation laboratory assays have been performed to estimate toxicokinetic parameters in mussels (Mytilus edulis) with two PAHs (chrysene and benzo[a]pyrene). The contamination was delivered via food. The exposure period lasted more than one month followed then by a depuration phase. The contaminant content was determined on an individual basis on five replicates. RESULTS AND DISCUSSION:. The reviewed information suggested the development of a tiered conceptual biomagnification model, starting with a simplified food chain which can be refined to more realistic and complex models in successive levels. CONCLUSIONS: The mathematical implementation of the conceptual model offers tools for estimating the potential for bioaccumulation and biomagnification of chemicals under very different conditions. The versatility of the model can be used for both comparative estimations and for validating the model. RECOMMENDATIONS AND PERSPECTIVES: Since bioaccumulation and biomagnification processes are crucial elements for a proper risk assessment of chemicals, their estimation by mathematical models has been widely tested. However, inregulatory assessments, too simplistic models are still being used quite often. The biomagnification model presented in this study should be amore accurate alternative to these models. In comparison to other previously published biomagnification models, the present one covers the time variation of bioaccumulation using just a few toxicokinetic parameters.     

Retooling the ethanol industry: thermophilic anaerobic digestion of thin stillage for methane production and pollution prevention.

Anaerobic digestion of corn ethanol thin stillage was tested at thermophilic temperature (55 degrees C) with two completely stirred tank reactors. The thin stillage wastestream was organically concentrated with 100 g/L total chemical oxygen demand and 60 g/L volatiles solids and a low pH of approximately 4.0. Steady-state was achieved at 30-, 20-, and 15-day hydraulic retention times (HRTs) and digester failure at a 12-day HRT. Significant reduction of volatile solids was achieved, with a maximum reduction (89.8%) at the 20-day HRT. Methane yield ranged from 0.6 to 0.7 L methane/g volatile solids removed during steady-state operation. Effluent volatile fatty acids below 200 mg/L as acetic acid were achieved at 20- and 30-day HRTs. Ultrasonic pretreatment was used for one digester, although no significant improvement was observed. Ethanol plant natural gas consumption could be reduced 43 to 59% with the methane produced, while saving an estimated $7 to $17 million ($10 million likely) for a facility producing 360 million L ethanol/y.     

Laboratory longitudinal diffusion tests: 2. Parameter estimation by inverse analysis

This study focuses on the verification of test interpretations for different state analyses of diffusion experiments. Part 1 of this study identified that steady, quasi-steady and equilibrium state analyses for the through- and in-diffusion tests with solution reservoirs are generally feasible where the tracer is not highly sorptive. In Part 2 we investigate parameter identifiability in transient-state analysis of reservoir concentration variation using a numerical approach. For increased generality, the analytical models, objective functions and Jacobian matrix necessary for inverse analysis of transient-state data are reformulated using unified dimensionless parameters. In these dimensionless forms, the number of unknown parameters is reduced and a single dimensionless parameter represents the sorption property. The dimensionless objective functions are evaluated for individual test methods and parameter identifiability is discussed in relation to the sorption property. The effects of multiple minima and measurement error on parameter identifiability are also investigated. The main findings are that inverse problems for inlet and outlet reservoir concentration analyses are generally unstable and well-posed, respectively. Where the tracer is sorptive, the inverse problem for the inlet reservoir concentration analysis may have multiple minima. When insufficient measurement data is collected, multiple solutions may result and this should be taken into consideration when inversely analyzing data including that of inlet reservoir concentration. Verification of test interpretation by cross-checking different state analyses is feasible where the tracer is not highly sorptive. In an actual experiment, test interpretation validity is demonstrated through consistency between theory and practice for different state analyses.     

Polybrominated diphenyl ethers and HBCD in bird eggs of South Africa

In this paper, the first data on brominated flame retardants (BFRs), in particular polybrominated diphenyl ethers (PBDEs) and hexabromocyclododecane (HBCD) in eggs of bird species from South Africa are described (N=43). Concentrations of PBDEs were detected in all the studied species and in all the geographic areas. Highest concentrations of PBDEs were measured in bird eggs from the Vaal River, which is situated downstream of the most industrialized area in South Africa. Sum-PBDE concentrations were highest in eggs of one African sacred ibis (396 ng g(-1) lipid weight (lw)), possibly due to foraging on dumping sites. Lowest mean level of sum-PBDEs (2.3 ng g(-1) lw) was measured in cattle egrets (N=11) from Barberspan Sanctuary, a Ramsar site. The PBDE congener pattern showed large differences, reflecting different trophic levels, migratory behavior, distance to the source, and, exposure to different PBDE mixtures, among others. HBCD was detected only in four species, and highest levels were measured in one egg of African sacred ibis (71 ng g(-1) lw). In most species, levels of PBDEs were one to several orders of magnitude lower than levels of DDTs and PCBs. PBDEs correlated strongly with DDTs, PCBs and some other organochlorines (OCs), indicating the same source. The widespread occurrence of PBDEs in the South African avian species showed a strong need for further investigations of PBDEs in the Southern African environment.     

Trends in snowpack chemistry and comparison to National Atmospheric Deposition Program results for the Rocky Mountains, US, 1993–2004

Seasonal snowpack chemistry data from the Rocky Mountain region of the US was examined to identify long-term trends in concentration and chemical deposition in snow and in snow-water equivalent. For the period 1993–2004, comparisons of trends were made between 54 Rocky Mountain Snowpack sites and 16 National Atmospheric Deposition Program wetfall sites located nearby in the region. The region was divided into three subregions: Northern, Central, and Southern. A non-parametric correlation method known as the Regional Kendall Test was used. This technique collectively computed the slope, direction, and probability of trend for several sites at once in each of the Northern, Central, and Southern Rockies subregions. Seasonal Kendall tests were used to evaluate trends at individual sites.Significant trends occurred during the period in wetfall and snowpack concentrations and deposition, and in precipitation. For the comparison, trends in concentrations of ammonium, nitrate, and sulfate for the two networks were in fair agreement. In several cases, increases in ammonium and nitrate concentrations, and decreases in sulfate concentrations for both wetfall and snowpack were consistent in the three subregions. However, deposition patterns between wetfall and snowpack more often were opposite, particularly for ammonium and nitrate. Decreases in ammonium and nitrate deposition in wetfall in the central and southern rockies subregions mostly were moderately significant (p<0.11) in constrast to highly significant increases in snowpack (p<0.02). These opposite trends likely are explained by different rates of declining precipitation during the recent drought (1999–2004) and increasing concentration. Furthermore, dry deposition was an important factor in total deposition of nitrogen in the region. Sulfate deposition decreased with moderate to high significance in all three subregions in both wetfall and snowpack. Precipitation trends consistently were downward and significant for wetfall, snowpack, and snow-telemetry data for the central and southern rockies subregions (p<0.03), while no trends were noted for the Northern Rockies subregion.     

Implementation of a windblown dust parameterization into MODELS-3/CMAQ: Application to episodic PM events in the US/Mexico border

Windblown dust is known to impede visibility, deteriorate air quality and modify the radiation budget. Arid and semiarid areas with unpaved and unvegetated land cover are particularly prone to windblown dust, which is often attributed to high particulate matter (PM) pollution in such areas. Yet, windblown dust is poorly represented in existing regulatory air quality models. In a study by the authors on modeling episodic high PM events along the US/Mexico border using the state-of-the-art CMAQ/MM5/SMOKE air quality modeling system [Choi, Y.-J., Hyde, P., Fernando, H.J.S., 2006. Modeling of episodic particulate matter events using a 3D air quality model with fine grid: applications to a pair of cities in the US/Mexico border. Atmospheric Environment 40, 5181–5201], some of the observed PM₁₀ NAAQS exceedances were inferred as due to windblown dust, but the modeling system was incapable of dealing with time-dependent episodic dust entrainment during high wind periods. In this paper, a time-dependent entrainment parameterization for windblown dust is implemented in the CMAQ/MM5/SMOKE modeling system with the hope of improving PM predictions. An approach for realizing windblown dust emission flux for each grid cell over the study domain on an hourly basis, which accounts for the influence of factors such as soil moisture content, atmospheric stability and wind speed, is presented in detail. Comparison of model predictions with observational data taken at a pair of US/Mexico border towns shows a clear improvement of model performance upon implementation of the dust emission flux parameterization.     

Solubility of toluene, benzene and TCE in high-microbial concentration systems

We report measurements of solubility limits for benzene, toluene, and TCE in systems that contain varying levels of biomass up to 0.13 g mL⁻¹ for TCE and 0.25 g mL⁻¹ for benzene and toluene. The solubility limit increased from 21 to 48 mM when biomass (in the form of yeast) was added to aqueous batch systems containing benzene. The toluene solubility limit increased from 4.9 to greater than 20 mM. For TCE, the solubility increased from 8 mM to more than 1000 mM. Solubility for TCE (trichloroethylene) was most heavily impacted by biomass levels, changing by two orders of magnitude as the microbial concentrations approach those in biofilms.     

Coupled infrared extinction spectra and size distribution measurements for several non-clay components of mineral dust aerosol (quartz, calcite, and dolomite)

Simultaneous size distributions and Fourier transform infrared (FTIR) extinction spectra have been measured for several representative components of mineral dust aerosol (quartz, calcite, and dolomite) in the fine particle size mode (D=0.1–1 μm). Optical constants drawn from the published literature have been used in combination with the experimentally determined size distributions to simulate the extinction spectra. In general, Mie theory does not accurately reproduce the peak position or band shape for the prominent IR resonance features in the 800–1600 cm⁻¹ spectral range. The resonance peaks in the Mie simulation are consistently blue shifted relative to the experimental spectra by 20–50 cm⁻¹. Spectral simulations, derived from a simple Rayleigh-based analytic theory for a “continuous distribution of ellipsoids” particle shape model, better reproduce the experimental spectra, despite the fact that the Rayleigh approximation is not strictly satisfied in these experiments. These results differ from our previous studies of particle shape effects in silicate clay mineral dust aerosols where a disk-shaped model for the particles was found to be more appropriate.     

A directional passive air sampler for monitoring polycyclic aromatic hydrocarbons (PAHs) in air mass

A passive air sampler was developed for collecting polycyclic aromatic hydrocarbons (PAHs) in air mass from various directions. The airflow velocity within the sampler was assessed for its responses to ambient wind speed and direction. The sampler was examined for trapped particles, evaluated quantitatively for influence of airflow velocity and temperature on PAH uptake, examined for PAH uptake kinetics, calibrated against active sampling, and finally tested in the field. The airflow volume passing the sampler was linearly proportional to ambient wind speed and sensitive to wind direction. The uptake rate for an individual PAH was a function of airflow velocity, temperature and the octanol-air partitioning coefficient of the PAH. For all PAHs with more than two rings, the passive sampler operated in a linear uptake phase for three weeks. Different PAH concentrations were obtained in air masses from different directions in the field test.     

Long term changes in atmospheric N and S throughfall deposition and effects on soil solution chemistry in a Scots pine forest in the Netherlands

In a Scots pine forest the throughfall deposition and the chemical composition of the soil solution was monitored since 1984. (Inter)national legislation measures led to a reduction of the deposition of nitrogen and sulphur. The deposition of sulphur has decreased by approximately 65%. The total mineral-nitrogen deposition has decreased by ca. 25%, which is mainly due to a reduction in ammonium-N deposition (−40%), since nitrate-N deposition has increased (+50%). The nitrogen concentration in the upper mineral soil solution at 10 cm depth has decreased, leading to an improved nutritional balance, which may result in improved tree vitality. In the drainage water at 90 cm depth the fluxes of NO₃⁻ and SO₄²⁻ have decreased, resulting in a reduced leeching of accompanying base cations, thus preserving nutrients in the ecosystem. It may take still several years, however, before this will meet the prerequisite of a sustainable ecosystem.