Processes controlling DOC in pore water during simulated drought cycles in six different UK peats

The effect of episodic drought on dissolved organic carbon (DOC) dynamics in peatlands has been the subject of considerable debate, as decomposition and DOC production is thought to increase under aerobic conditions, yet decreased DOC concentrations have been observed during drought periods. Decreased DOC solubility due to drought-induced acidification driven by sulphur (S) redox reactions has been proposed as a causal mechanism; however evidence is based on a limited number of studies carried out at a few sites. To test this hypothesis on a range of different peats, we carried out controlled drought simulation experiments on peat cores collected from six sites across Great Britain. Our data show a concurrent increase in sulphate (SO4) and a decrease in DOC across all sites during simulated water table draw-down, although the magnitude of the relationship between SO4 and DOC differed between sites. Instead, we found a consistent relationship across all sites between DOC decrease and acidification measured by the pore water acid neutralising capacity (ANC). ANC provided a more consistent measure of drought-induced acidification than SO4 alone because it accounts for differences in base cation and acid anions concentrations between sites. Rewetting resulted in rapid DOC increases without a concurrent increase in soil respiration, suggesting DOC changes were primarily controlled by soil acidity not soil biota. These results highlight the need for an integrated analysis of hydrologically driven chemical and biological processes in peatlands to improve our understanding and ability to predict the interaction between atmospheric pollution and changing climatic conditions from plot to regional and global scales.

Multimodel estimate of the global terrestrial water balance: setup and first results

Six land surface models and five global hydrological models participate in a model intercomparison project (WaterMIP), which for the first time compares simulation results of these different classes of models in a consistent way. In this paper the simulation setup is described and aspects of the multi-model global terrestrial water balance are presented. All models were run at 0.5 degree spatial resolution for the global land areas for a 15-year period (1985-1999) using a newly-developed global meteorological dataset. Simulated global terrestrial evapotranspiration, excluding Greenland and Antarctica, ranges from 415 to 586 mm year-1 (60,000 to 85,000 km3 year-1) and simulated runoff ranges from 290 to 457 mm year-1 (42,000 to 66,000 km3 year-1). Both the mean and median runoff fractions for the land surface models are lower than those of the global hydrological models, although the range is wider. Significant simulation differences between land surface and global hydrological models are found to be caused by the snow scheme employed. The physically-based energy balance approach used by land surface models generally results in lower snow water equivalent values than the conceptual degree-day approach used by global hydrological models. Some differences in simulated runoff and evapotranspiration are explained by model parameterizations, although the processes included and parameterizations used are not distinct to either land surface models or global hydrological models. The results show that differences between model are major sources of uncertainty. Climate change impact studies thus need to use not only multiple climate models, but also some other measure of uncertainty, (e.g. multiple impact models).

Adsorption from aqueous solutions on opened carbon nanotubes: organic compounds speed up delivery of water from inside

We report the results of first systematic studies of organic adsorption from aqueous solutions onto relatively long single walled carbon nanotubes (four tubes, in initial and oxidised forms). Using molecular dynamics simulations (GROMACS package) we discuss the behaviour of tube-water as well as tube-adsorbate systems, for three different adsorbates (benzene, phenol and paracetamol).

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Evaluation of water and energy budgets in regional climate models applied over Europe

This study presents a model intercomparison of four regional climate models (RCMs) and one variable resolution atmospheric general circulation model (AGCM) applied over Europe with special focus on the hydrological cycle and the surface energy budget. The models simulated the 15 years from 1979 to 1993 by using quasi-observed boundary conditions derived from ECMWF re-analyses (ERA). The model intercomparison focuses on two large atchments representing two different climate conditions covering two areas of major research interest within Europe. The first is the Danube catchment which represents a continental climate dominated by advection from the surrounding land areas. It is used to analyse the common model error of a too dry and too warm simulation of the summertime climate of southeastern Europe. This summer warming and drying problem is seen in many RCMs, and to a less extent in GCMs. The second area is the Baltic Sea catchment which represents maritime climate dominated by advection from the ocean and from the Baltic Sea. This catchment is a research area of many studies within Europe and also covered by the BALTEX program. The observed data used are monthly mean surface air temperature, precipitation and river discharge. For all models, these are used to estimate mean monthly biases of all components of the hydrological cycle over land. In addition, the mean monthly deviations of the surface energy fluxes from ERA data are computed. Atmospheric moisture fluxes from ERA are compared with those of one model to provide an independent estimate of the convergence bias derived from the observed data. These help to add weight to some of the inferred estimates and explain some of the discrepancies between them. An evaluation of these biases and deviations suggests possible sources of error in each of the models. For the Danube catchment, systematic errors in the dynamics cause the prominent summer drying problem for three of the RCMs, while for the fourth RCM this is related to deficiencies in the land surface parametrization. The AGCM does not show this drying problem. For the Baltic Sea catchment, all models similarily overestimate the precipitation throughout the year except during the summer. This model deficit is probably caused by the internal model parametrizations, such as the large-scale condensation and the convection schemes.

Variability in fecal water genotoxicity, determined using the Comet assay, is independent of endogenous N-nitroso compound formation attributed to red meat consumption

Red meat consumption causes a dose-dependent increase in fecal apparent total N-nitroso compounds (ATNC). The genotoxic effects of these ATNCs were investigated using two different Comet assay protocols to determine the genotoxicity of fecal water samples. Fecal water samples were obtained from two studies of a total of 21 individuals fed diets containing different amounts of red meat, protein, heme, and iron. The first protocol incubated the samples with HT-29 cells for 5 min at 4 degrees C, whereas the second protocol used a longer exposure time of 30 min and a higher incubation temperature of 37 degrees C. DNA strand breaks were quantified by the tail moment (DNA in the comet tail multiplied by the comet tail length). The results of the two Comet assay protocols were significantly correlated (r = 0.35, P = 0.003), however, only the second protocol resulted in detectable levels of DNA damage. Inter-individual effects were variable and there was no effect on fecal water genotoxicity by diet (P > 0.20), mean transit time (P = 0.588), or weight (P = 0.705). However, there was a highly significant effect of age (P = 0.019). There was no significant correlation between concentrations of ATNCs in fecal homogenates and fecal water genotoxicity (r = 0.04, P = 0.74). ATNC levels were lower in fecal water samples (272 microg/kg) compared to that of fecal homogenate samples (895 microg/kg) (P < 0.0001). Failure to find dietary effects on fecal water genotoxicity may therefore be attributed to individual variability and low levels of ATNCs in fecal water samples.

Antigenotoxicity of probiotics and prebiotics on faecal water-induced DNA damage in human colon adenocarcinoma cells

Six strains of lactic acid producing bacteria (LAB) were incubated (1 x 10(8)cfu/ml) with genotoxic faecal water from a human subject. HT29 human adenocarcinoma cells were then challenged with the resultant samples and DNA damage measured using the single cell gel electrophoresis (comet) assay. The LAB strains investigated were Bifidobacterium sp. 420, Bifidobacterium Bb12, Lactobacillus plantarum, Streptococcus thermophilus, Lactobacillus bulgaricus and Enterococcus faecium. DNA damage was significantly decreased by all bacteria used with the exception of Strep. thermophilus. Bif. Bb12 and Lact. plantarum showed the greatest protective effect against DNA damage. Incubation of faecal water with different concentrations of Bif. Bb12 and Lact. plantarum revealed that the decrease in genotoxicity was related to cell density. Non-viable (heat treated) probiotic cells had no effect on faecal water genotoxicity. In a second study, HT29 cells were cultured in the presence of supernatants of incubations of probiotics with various carbohydrates including known prebiotics; the HT29 cells were then exposed to faecal water. Overall, incubations involving Lact. plantarum with the fructooligosaccharide (FOS)-based prebiotics Inulin, Raftiline, Raftilose and Actilight were the most effective in increasing the cellular resistance to faecal water genotoxicity, whereas fermentations with Elixor (a galactooligosaccharide) and Fibersol (a maltodextrin) were less effective. Substantial reductions in faecal water-induced DNA damage were also seen with supernatants from incubation of prebiotics with Bif. Bb12. The supernatant of fermentations involving Ent. faecium and Bif. sp. 420 generally had less potent effects on genotoxicity although some reductions with Raftiline and Elixor fermentations were apparent.

The Aqua-Planet Experiment (APE): CONTROL SST Simulation

Climate simulations by 16 atmospheric general circulation models (AGCMs) are compared on an aqua-planet, a water-covered Earth with prescribed sea surface temperature varying only in latitude. The idealised configuration is designed to expose differences in the circulation simulated by different models. Basic features of the aqua-planet climate are characterised by comparison with Earth. The models display a wide range of behaviour. The balanced component of the tropospheric mean flow, and mid-latitude eddy covariances subject to budget constraints, vary relatively little among the models. In contrast, differences in damping in the dynamical core strongly influence transient eddy amplitudes. Historical uncertainty in modelled lower stratospheric temperatures persists in APE. Aspects of the circulation generated more directly by interactions between the resolved fluid dynamics and parameterized moist processes vary greatly. The tropical Hadley circulation forms either a single or double inter-tropical convergence zone (ITCZ) at the equator, with large variations in mean precipitation. The equatorial wave spectrum shows a wide range of precipitation intensity and propagation characteristics. Kelvin mode-like eastward propagation with remarkably constant phase speed dominates in most models. Westward propagation, less dispersive than the equatorial Rossby modes, dominates in a few models or occurs within an eastward propagating envelope in others. The mean structure of the ITCZ is related to precipitation variability, consistent with previous studies. The aqua-planet global energy balance is unknown but the models produce a surprisingly large range of top of atmosphere global net flux, dominated by differences in shortwave reflection by clouds. A number of newly developed models, not optimised for Earth climate, contribute to this. Possible reasons for differences in the optimised models are discussed. The aqua-planet configuration is intended as one component of an experimental hierarchy used to evaluate AGCMs. This comparison does suggest that the range of model behaviour could be better understood and reduced in conjunction with Earth climate simulations. Controlled experimentation is required to explore individual model behaviour and investigate convergence of the aqua-planet climate with increasing resolution.

Novel polyvinylpyrrolidones to improve delivery of poorly-water soluble drugs; from design to synthesis and evaluation

Polyvinylpyrrolidone is a widely used in tablet formulations with the linear form acting as a wetting agent and disintegrant whereas the cross-linked form is a super-disintegrant. We have previously reported that simply mixing the commercial cross-linked polymer with ibuprofen disrupted drug crystallinity with consequent improvements in drug dissolution behavior. In this study, we have designed and synthesized novel cross-linking agents containing a range of oligoether moieties which have then be polymerized with vinylpyrrolidone to generate a suite of novel excipients with enhanced hydrogen-bonding capabilities. The polymers have a porous surface and swell in most common solvents and in water; properties which suggest their value as disintegrants. The polymers were evaluated in simple physical mixtures with ibuprofen as a model poorly-water soluble drug. The results show that the novel PVPs induce the drug to become “X-ray amorphous”, which increased dissolution to a greater extent than that seen with commercial cross-linked PVP. The polymers stabilize the amorphous drug with no evidence for recrystallization seen after 20 weeks storage.

Assessment of potential climate change impacts on peatland dissolved organic carbon release and drinking water treatment from laboratory experiments

Catchments draining peat soils provide the majority of drinking water in the UK. Over the past decades, concentrations of dissolved organic carbon (DOC) have increased in surface waters. Residual DOC can cause harmful carcinogenic disinfection by-products to form during water treatment processes. Increased frequency and severity of droughts combined with and increased temperatures expected as the climate changes, have potentials to change water quality. We used a novel approach to investigate links between climate change, DOC release and subsequent effects on drinking water treatment. We designed a climate manipulation experiment to simulate projected climate changes and monitored releases from peat soil and litter, then simulated coagulation used in water treatment. We showed that the ‘drought’ simulation was the dominant factor altering DOC release and affected the ability to remove DOC. Our results imply that future short-term drought events could have a greater impact than increased temperature on DOC treatability.

Cloud droplet size and liquid water path retrievals from zenith radiance measurements: examples from the Atmospheric Radiation Measurement Program and the Aerosol Robotic Network

The ground-based Atmospheric Radiation Measurement Program (ARM) and NASA Aerosol Robotic Net- work (AERONET) routinely monitor clouds using zenith ra- diances at visible and near-infrared wavelengths. Using the transmittance calculated from such measurements, we have developed a new retrieval method for cloud effective droplet size and conducted extensive tests for non-precipitating liquid water clouds. The underlying principle is to combine a liquid-water-absorbing wavelength (i.e., 1640 nm) with a non-water-absorbing wavelength for acquiring information on cloud droplet size and optical depth. For simulated stratocumulus clouds with liquid water path less than 300 g m−2 and horizontal resolution of 201 m, the retrieval method underestimates the mean effective radius by 0.8μm, with a root-mean-squared error of 1.7 μm and a relative deviation of 13%. For actual observations with a liquid water path less than 450 g m−2 at the ARM Oklahoma site during 2007– 2008, our 1.5-min-averaged retrievals are generally larger by around 1 μm than those from combined ground-based cloud radar and microwave radiometer at a 5-min temporal resolution. We also compared our retrievals to those from combined shortwave flux and microwave observations for relatively homogeneous clouds, showing that the bias between these two retrieval sets is negligible, but the error of 2.6 μm and the relative deviation of 22 % are larger than those found in our simulation case. Finally, the transmittance-based cloud effective droplet radii agree to better than 11 % with satellite observations and have a negative bias of 1 μm. Overall, the retrieval method provides reasonable cloud effective radius estimates, which can enhance the cloud products of both ARM and AERONET.

Towards data warehousing and mining of protein unfolding simulation data

OBJECTIVES: The prediction of protein structure and the precise understanding of protein folding and unfolding processes remains one of the greatest challenges in structural biology and bioinformatics. Computer simulations based on molecular dynamics (MD) are at the forefront of the effort to gain a deeper understanding of these complex processes. Currently, these MD simulations are usually on the order of tens of nanoseconds, generate a large amount of conformational data and are computationally expensive. More and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. METHODS: To adequately organize, manage, and analyze the data generated by unfolding simulation studies, we designed a data warehouse system that is embedded in a grid environment to facilitate the seamless sharing of available computer resources and thus enable many groups to share complex molecular dynamics simulations on a more regular basis. RESULTS: To gain insight into the conformational fluctuations and stability of the monomeric forms of the amyloidogenic protein transthyretin (TTR), molecular dynamics unfolding simulations of the monomer of human TTR have been conducted. Trajectory data and meta-data of the wild-type (WT) protein and the highly amyloidogenic variant L55P-TTR represent the test case for the data warehouse. CONCLUSIONS: Web and grid services, especially pre-defined data mining services that can run on or 'near' the data repository of the data warehouse, are likely to play a pivotal role in the analysis of molecular dynamics unfolding data.

Numerical simulation of the Filchner overflow

The plume of Ice Shelf Water (ISW) flowing into the Weddell Sea over the Filchner sill contributes to the formation of Antarctic Bottom Water. The Filchner overflow is simulated using a hydrostatic, primitive equation three-dimensional ocean model with a 0.5–2 Sv ISW influx above the Filchner sill. The best fit to mooring temperature observations is found with influxes of 0.5 and 1 Sv, below a previous estimate of 1.6 ± 0.5 Sv based on sparse mooring velocities. The plume first moves north over the continental shelf, and then turns west, along slope of the continental shelf break where it breaks up into subplumes and domes, some of which then move downslope. Other subplumes run into the eastern submarine ridge and propagate along the ridge downslope in a chaotic manner. The next, western ridge is crossed by the plume through several paths. Despite a number of discrepancies with observational data, the model reproduces many attributes of the flow. In particular, we argue that the temporal variability shown by the observations can largely be attributed to the unstable structure of the flow, where the temperature fluctuations are determined by the motion of the domes past the moorings. Our sensitivity studies show that while thermobaricity plays a role, its effect is small for the flows considered. Smoothing the ridges out demonstrate that their presence strongly affects the plume shape around the ridges. An increase in the bottom drag or viscosity leads to slowing down, and hence thickening and widening of the plume

Evaluation of water and energy budgets in regional climate models applied over Europe

Abstract This study presents a model intercomparison of four regional climate models (RCMs) and one variable resolution atmospheric general circulation model (AGCM) applied over Europe with special focus on the hydrological cycle and the surface energy budget. The models simulated the 15 years from 1979 to 1993 by using quasi-observed boundary conditions derived from ECMWF re-analyses (ERA). The model intercomparison focuses on two large atchments representing two different climate conditions covering two areas of major research interest within Europe. The first is the Danube catchment which represents a continental climate dominated by advection from the surrounding land areas. It is used to analyse the common model error of a too dry and too warm simulation of the summertime climate of southeastern Europe. This summer warming and drying problem is seen in many RCMs, and to a less extent in GCMs. The second area is the Baltic Sea catchment which represents maritime climate dominated by advection from the ocean and from the Baltic Sea. This catchment is a research area of many studies within Europe and also covered by the BALTEX program. The observed data used are monthly mean surface air temperature, precipitation and river discharge. For all models, these are used to estimate mean monthly biases of all components of the hydrological cycle over land. In addition, the mean monthly deviations of the surface energy fluxes from ERA data are computed. Atmospheric moisture fluxes from ERA are compared with those of one model to provide an independent estimate of the convergence bias derived from the observed data. These help to add weight to some of the inferred estimates and explain some of the discrepancies between them. An evaluation of these biases and deviations suggests possible sources of error in each of the models. For the Danube catchment, systematic errors in the dynamics cause the prominent summer drying problem for three of the RCMs, while for the fourth RCM this is related to deficiencies in the land surface parametrization. The AGCM does not show this drying problem. For the Baltic Sea catchment, all models similarily overestimate the precipitation throughout the year except during the summer. This model deficit is probably caused by the internal model parametrizations, such as the large-scale condensation and the convection schemes.