iMarNet: an ocean biogeochemistry model intercomparison project within a common physical ocean modelling framework

Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.

Nitrogen budget for the Oldman River Basin, southern Alberta, Canada

The Oldman River Basin (OMRB), located in southern Alberta (Canada), with an area of 28,200 km2, is mainly forested in its western part and is used for intensive agriculture in its eastern part. The objective of this paper is to estimate the nitrogen (N) budget for the Oldman River Basin as a whole and its sub-basins, and to discuss differences in the N budget between various sub-basins. Better knowledge of the N budget in this watershed may be also utilized for understanding N dynamics in similar watersheds within semi-arid climatic regions. The model used is a mass balance spreadsheet model that takes into account N inputs and N export through surface water. During the last 120 years, anthropogenic N inputs to the OMRB have increased circa 40 fold. By the end of the 20th century, the OMRB received an annualN input of about 5174 kg N km-2 yr-1, whereas only about 25 kg N km-2 yr-1 were exported via riverine flow. For the sub-basins, annual N inputs ranged from 2516 to 19011 kg N km-2 yr-1, and annual N export via riverine flows varied between 6 and 277 kg N km-2 yr-1. Over 85% of total N inputs to the OMRB are due to anthropogenic activities, including manure (55%), synthetic fertilizer (27%), and N fixation on agricultural lands (4%). Sewage accounted for less than 1%, and N inputs from atmospheric deposition and fixation in forests represented 6 and 8% respectively. Despite increasing anthropogenic N inputs, N export with riverine flow currently accounts for only 1% of the inputs, indicating thatmost of theNinputs are currently retained in the OMRB or are re-emitted into the atmosphere.

A heuristic model for the simulation of the deformation of elastic and spongy material for virtual reality applications

This paper presents a practical algorithm for the simulation of interactive deformation in a 3D polygonal mesh model. The algorithm combines the conventional simulation of deformation using a spring-mass-damping model, solved by explicit numerical integration, with a set of heuristics to describe certain features of the transient behaviour, to increase the speed and stability of solution. In particular, this algorithm was designed to be used in the simulation of synthetic environments where it is necessary to model realistically, in real time, the effect on non-rigid surfaces being touched, pushed, pulled or squashed. Such objects can be solid or hollow, and have plastic, elastic or fabric-like properties. The algorithm is presented in an integrated form including collision detection and adaptive refinement so that it may be used in a self-contained way as part of a simulation loop to include human interface devices that capture data and render a realistic stereoscopic image in real time. The algorithm is designed to be used with polygonal mesh models representing complex topology, such as the human anatomy in a virtual-surgery training simulator. The paper evaluates the model behaviour qualitatively and then concludes with some examples of the use of the algorithm.

Pressure drop of gas–liquid Taylor flow in round micro-capillaries for low to intermediate Reynolds numbers

In this paper, a model is presented that describes the pressure drop of gas-liquid Taylor flow in round capillaries with a channel diameter typically less than 1 mm. The analysis of Bretherton (J Fluid Mech 10:166-188, 1961) for the pressure drop over a single gas bubble for vanishing liquid film thickness is extended to include a non-negligible liquid film thickness using the analysis of Aussillous and Qu,r, (Phys Fluids 12(10):2367-2371, 2000). This result is combined with the Hagen-Poiseuille equation for liquid flow using a mass balance-based Taylor flow model previously developed by the authors (Warnier et al. in Chem Eng J 135S:S153-S158, 2007). The model presented in this paper includes the effect of the liquid slug length on the pressure drop similar to the model of Kreutzer et al. (AIChE J 51(9):2428-2440, 2005). Additionally, the gas bubble velocity is taken into account, thereby increasing the accuracy of the pressure drop predictions compared to those of the model of Kreutzer et al. Experimental data were obtained for nitrogen-water Taylor flow in a round glass channel with an inner diameter of 250 mu m. The capillary number Ca (gl) varied between 2.3 x 10(-3) and 8.8 x 10(-3) and the Reynolds number Re (gl) varied between 41 and 159. The presented model describes the experimental results with an accuracy of +/- 4% of the measured values.

Degradation of 4-fluorobiphenyl in soil investigated by 19F NMR spectroscopy and 14C radiolabelling analysis

The incubation of the model pollutant [U-14C]'-4-fluorobiphenyl (4FBP) in soil, in the presence and absence of biphenyl (a co-substrate), was carried out in order to study the qualitative disposition and fate of the compound using 14C-HPLC and 19F NMR spectroscopy. Components accounted for using the radiolabel were volatilization, CO2 evolution, organic solvent extractable and bound residue. Quantitative analysis of these data gave a complete mass balance. After sample preparation. 14C-HPLC was used to establish the number of 4FBP related components present in the organic solvent extract. 19F NMR was also used to quantify the organic extracts and to identify the components of the extract. Both approaches showed that the composition of the solvent extractable fractions comprised only parent compound with no metabolites present. As the 14C radiolabel was found to be incorporated into the soil organic matter this indicates that metabolites were being generated, but were highly transitory as incorporation into the SOM was rapid. The inclusion of the co-substrate biphenyl was to increase the overall rate of degradation of 4FBP in soil. The kinetics of disappearance of parent from the soil using the data obtained were investigated from both techniques. This is the first report describing the degradation of a fluorinated biphenyl in soil.

Evaluation of thermochemical biomass conversion in fluidized bed

Dado o aumento acelerado dos preços dos combustíveis fósseis e as incertezas quanto à sua disponibilidade futura, tem surgido um novo interesse nas tecnologias da biomassa aplicadas à produção de calor, eletricidade ou combustíveis sintéticos. Não obstante, para a conversão termoquímica de uma partícula de biomassa sólida concorrem fenómenos bastante complexos que levam, em primeiro lugar, à secagem do combustível, depois à pirólise e finalmente à combustão ou gasificação propriamente ditas. Uma descrição relativamente incompleta de alguns desses estágios de conversão constitui ainda um obstáculo ao desenvolvimento das tecnologias que importa ultrapassar. Em particular, a presença de elevados conteúdos de matéria volátil na biomassa põe em evidência o interesse prático do estudo da pirólise. A importância da pirólise durante a combustão de biomassa foi evidenciada neste trabalho através de ensaios realizados num reator piloto de leito fluidizado borbulhante. Verificou-se que o processo ocorre em grande parte à superfície do leito com chamas de difusão devido à libertação de voláteis, o que dificulta o controlo da temperatura do reator acima do leito. No caso da gasificação de biomassa a pirólise pode inclusivamente determinar a eficiência química do processo. Isso foi mostrado neste trabalho durante ensaios de gasificação num reator de leito fluidizado de 2MWth, onde um novo método de medição permitiu fechar o balanço de massa ao gasificador e monitorizar o grau de conversão da biomassa. A partir destes resultados tornou-se clara a necessidade de descrever adequadamente a pirólise de biomassa com vista ao projeto e controlo dos processos. Em aplicações de engenharia há particular interesse na estequiometria e propriedades dos principais produtos pirolíticos. Neste trabalho procurou-se responder a esta necessidade, inicialmente através da estruturação de dados bibliográficos sobre rendimentos de carbonizado, líquidos pirolíticos e gases, assim como composições elementares e poderes caloríficos. O resultado traduziu-se num conjunto de parâmetros empíricos de interesse prático que permitiram elucidar o comportamento geral da pirólise de biomassa numa gama ampla de condições operatórias. Para além disso, propôs-se um modelo empírico para a composição dos voláteis que pode ser integrado em modelos compreensivos de reatores desde que os parâmetros usados sejam adequados ao combustível ensaiado. Esta abordagem despoletou um conjunto de ensaios de pirólise com várias biomassas, lenhina e celulose, e temperaturas entre os 600 e 975ºC. Elevadas taxas de aquecimento do combustível foram alcançadas em reatores laboratoriais de leito fluidizado borbulhante e leito fixo, ao passo que um sistema termo-gravimétrico permitiu estudar o efeito de taxas de aquecimento mais baixas. Os resultados mostram que, em condições típicas de processos de combustão e gasificação, a quantidade de voláteis libertada da biomassa é pouco influenciada pela temperatura do reator mas varia bastante entre combustíveis. Uma análise mais aprofundada deste assunto permitiu mostrar que o rendimento de carbonizado está intimamente relacionado com o rácio O/C do combustível original, sendo proposto um modelo simples para descrever esta relação. Embora a quantidade total de voláteis libertada seja estabelecida pela composição da biomassa, a respetiva composição química depende bastante da temperatura do reator. Rendimentos de espécies condensáveis (água e espécies orgânicas), CO2 e hidrocarbonetos leves descrevem um máximo relativamente à temperatura para dar lugar a CO e H2 às temperaturas mais altas. Não obstante, em certas gamas de temperatura, os rendimentos de algumas das principais espécies gasosas (e.g. CO, H2, CH4) estão bem correlacionados entre si, o que permitiu desenvolver modelos empíricos que minimizam o efeito das condições operatórias e, ao mesmo tempo, realçam o efeito do combustível na composição do gás. Em suma, os ensaios de pirólise realizados neste trabalho permitiram constatar que a estequiometria da pirólise de biomassa se relaciona de várias formas com a composição elementar do combustível original o que levanta várias possibilidades para a avaliação e projeto de processos de combustão e gasificação de biomassa.

Cape Verde coastal ecosystem: a study of community structure, trophic interactions and exploitation pattern

Tese de dout., Ciências e Tecnologia das Pescas, Faculdade de Ciências do Mar e do Ambiente, Universidade do Algarve, 2005

Studies on Advanced Oxidation Processes for the Removal of Acetamiprid from Wastewater

Pollution of water with pesticides has become a threat to the man, material and environment. The pesticides released to the environment reach the water bodies through run off. Industrial wastewater from pesticide manufacturing industries contains pesticides at higher concentration and hence a major source of water pollution. Pesticides create a lot of health and environmental hazards which include diseases like cancer, liver and kidney disorders, reproductive disorders, fatal death, birth defects etc. Conventional wastewater treatment plants based on biological treatment are not efficient to remove these compounds to the desired level. Most of the pesticides are phyto-toxic i.e., they kill the microorganism responsible for the degradation and are recalcitrant in nature. Advanced oxidation process (AOP) is a class of oxidation techniques where hydroxyl radicals are employed for oxidation of pollutants. AOPs have the ability to totally mineralise the organic pollutants to CO2 and water. Different methods are employed for the generation of hydroxyl radicals in AOP systems. Acetamiprid is a neonicotinoid insecticide widely used to control sucking type insects on crops such as leafy vegetables, citrus fruits, pome fruits, grapes, cotton, ornamental flowers. It is now recommended as a substitute for organophosphorous pesticides. Since its use is increasing, its presence is increasingly found in the environment. It has high water solubility and is not easily biodegradable. It has the potential to pollute surface and ground waters. Here, the use of AOPs for the removal of acetamiprid from wastewater has been investigated. Five methods were selected for the study based on literature survey and preliminary experiments conducted. Fenton process, UV treatment, UV/ H2O2 process, photo-Fenton and photocatalysis using TiO2 were selected for study. Undoped TiO2 and TiO2 doped with Cu and Fe were prepared by sol-gel method. Characterisation of the prepared catalysts was done by X-ray diffraction, scanning electron microscope, differential thermal analysis and thermogravimetric analysis. Influence of major operating parameters on the removal of acetamiprid has been investigated. All the experiments were designed using central compoiste design (CCD) of response surface methodology (RSM). Model equations were developed for Fenton, UV/ H2O2, photo-Fenton and photocatalysis for predicting acetamiprid removal and total organic carbon (TOC) removal for different operating conditions. Quality of the models were analysed by statistical methods. Experimental validations were also done to confirm the quality of the models. Optimum conditions obtained by experiment were verified with that obtained using response optimiser. Fenton Process is the simplest and oldest AOP where hydrogen peroxide and iron are employed for the generation of hydroxyl radicals. Influence of H2O2 and Fe2+ on the acetamiprid removal and TOC removal by Fenton process were investigated and it was found that removal increases with increase in H2O2 and Fe2+ concentration. At an initial concentration of 50 mg/L acetamiprid, 200 mg/L H2O2 and 20 mg/L Fe2+ at pH 3 was found to be optimum for acetamiprid removal. For UV treatment effect of pH was studied and it was found that pH has not much effect on the removal rate. Addition of H2O2 to UV process increased the removal rate because of the hydroxyl radical formation due to photolyis of H2O2. An H2O2 concentration of 110 mg/L at pH 6 was found to be optimum for acetamiprid removal. With photo-Fenton drastic reduction in the treatment time was observed with 10 times reduction in the amount of reagents required. H2O2 concentration of 20 mg/L and Fe2+ concentration of 2 mg/L was found to be optimum at pH 3. With TiO2 photocatalysis improvement in the removal rate was noticed compared to UV treatment. Effect of Cu and Fe doping on the photocatalytic activity under UV light was studied and it was observed that Cu doping enhanced the removal rate slightly while Fe doping has decreased the removal rate. Maximum acetamiprid removal was observed for an optimum catalyst loading of 1000 mg/L and Cu concentration of 1 wt%. It was noticed that mineralisation efficiency of the processes is low compared to acetamiprid removal efficiency. This may be due to the presence of stable intermediate compounds formed during degradation Kinetic studies were conducted for all the treatment processes and it was found that all processes follow pseudo-first order kinetics. Kinetic constants were found out from the experimental data for all the processes and half lives were calculated. The rate of reaction was in the order, photo- Fenton>UV/ H2O2>Fenton> TiO2 photocatalysis>UV. Operating cost was calculated for the processes and it was found that photo-Fenton removes the acetamiprid at lowest operating cost in lesser time. A kinetic model was developed for photo-Fenton process using the elementary reaction data and mass balance equations for the species involved in the process. Variation of acetamiprid concentration with time for different H2O2 and Fe2+ concentration at pH 3 can be found out using this model. The model was validated by comparing the simulated concentration profiles with that obtained from experiments. This study established the viability of the selected AOPs for the removal of acetamiprid from wastewater. Of the studied AOPs photo- Fenton gives the highest removal efficiency with lowest operating cost within shortest time.

“Unmixing” Tissue Gene Expression Signatures from Tumor Biopsies

Emergent molecular measurement methods, such as DNA microarray, qRTPCR, and many others, offer tremendous promise for the personalized treatment of cancer. These technologies measure the amount of specific proteins, RNA, DNA or other molecular targets from tumor specimens with the goal of “fingerprinting” individual cancers. Tumor specimens are heterogeneous; an individual specimen typically contains unknown amounts of multiple tissues types. Thus, the measured molecular concentrations result from an unknown mixture of tissue types, and must be normalized to account for the composition of the mixture. For example, a breast tumor biopsy may contain normal, dysplastic and cancerous epithelial cells, as well as stromal components (fatty and connective tissue) and blood and lymphatic vessels. Our diagnostic interest focuses solely on the dysplastic and cancerous epithelial cells. The remaining tissue components serve to “contaminate” the signal of interest. The proportion of each of the tissue components changes as a function of patient characteristics (e.g., age), and varies spatially across the tumor region. Because each of the tissue components produces a different molecular signature, and the amount of each tissue type is specimen dependent, we must estimate the tissue composition of the specimen, and adjust the molecular signal for this composition. Using the idea of a chemical mass balance, we consider the total measured concentrations to be a weighted sum of the individual tissue signatures, where weights are determined by the relative amounts of the different tissue types. We develop a compositional source apportionment model to estimate the relative amounts of tissue components in a tumor specimen. We then use these estimates to infer the tissuespecific concentrations of key molecular targets for sub-typing individual tumors. We anticipate these specific measurements will greatly improve our ability to discriminate between different classes of tumors, and allow more precise matching of each patient to the appropriate treatment

Selection-mutation balance models with epistatic selection

We present an application of birth-and-death processes on configuration spaces to a generalized mutation4 selection balance model. The model describes the aging of population as a process of accumulation of mu5 tations in a genotype. A rigorous treatment demands that mutations correspond to points in abstract spaces. 6 Our model describes an infinite-population, infinite-sites model in continuum. The dynamical equation which 7 describes the system, is of Kimura-Maruyama type. The problem can be posed in terms of evolution of states 8 (differential equation) or, equivalently, represented in terms of Feynman-Kac formula. The questions of interest 9 are the existence of a solution, its asymptotic behavior, and properties of the limiting state. In the non-epistatic 10 case the problem was posed and solved in [Steinsaltz D., Evans S.N., Wachter K.W., Adv. Appl. Math., 2005, 11 35(1)]. In our model we consider a topological space X as the space of positions of mutations and the influence of epistatic potentials

How many facets are needed to represent the surface energy balance of an urban area?

We investigate the question of how many facets are needed to represent the energy balance of an urban area by developing simplified 3-, 2- and 1-facet versions of a 4-facet energy balance model of two-dimensional streets and buildings. The 3-facet model simplifies the 4-facet model by averaging over the canyon orientation, which results in similar net shortwave and longwave balances for both wall facets, but maintains the asymmetry in the heat fluxes within the street canyon. For the 2-facet model, on the assumption that the wall and road temperatures are equal, the road and wall facets can be combined mathematically into a single street-canyon facet with effective values of the heat transfer coefficient, albedo, emissivity and thermodynamic properties, without further approximation. The 1-facet model requires the additional assumption that the roof temperature is also equal to the road and wall temperatures. Idealised simulations show that the geometry and material properties of the walls and road lead to a large heat capacity of the combined street canyon, whereas the roof behaves like a flat surface with low heat capacity. This means that the magnitude of the diurnal temperature variation of the street-canyon facets are broadly similar and much smaller than the diurnal temperature variation of the roof facets. Consequently, the approximation that the street-canyon facets have similar temperatures is sound, and the road and walls can be combined into a single facet. The roof behaves very differently and a separate roof facet is required. Consequently, the 2-facet model performs similarly to the 4-facet model, while the 1-facet model does not. The models are compared with previously published observations collected in Mexico City. Although the 3- and 2-facet models perform better than the 1-facet model, the present models are unable to represent the phase of the sensible heat flux. This result is consistent with previous model comparisons, and we argue that this feature of the data cannot be produced by a single column model. We conclude that a 2-facet model is necessary, and for numerical weather prediction sufficient, to model an urban surface, and that this conclusion is robust and therefore applicable to more general geometries.

Aqueous geochemistry and oxygen isotope compositions of acid mine drainage from the Río Tinto, SW Spain, highlight inconsistencies in current models

The Rio Tinto river in SW Spain is a classic example of acid mine drainage and the focus of an increasing amount of research including environmental geochemistry, extremophile microbiology and Mars-analogue studies. Its 5000-year mining legacy has resulted in a wide range of point inputs including spoil heaps and tunnels draining underground workings. The variety of inputs and importance of the river as a research site make it an ideal location for investigating sulphide oxidation mechanisms at the field scale. Mass balance calculations showed that pyrite oxidation accounts for over 93% of the dissolved sulphate derived from sulphide oxidation in the Rio Tinto point inputs. Oxygen isotopes in water and sulphate were analysed from a variety of drainage sources and displayed delta O-18((SO4-H2O)) values from 3.9 to 13.6 parts per thousand, indicating that different oxidation pathways occurred at different sites within the catchment. The most commonly used approach to interpreting field oxygen isotope data applies water and oxygen fractionation factors derived from laboratory experiments. We demonstrate that this approach cannot explain high delta O-18((SO4-H2O)) values in a manner that is consistent with recent models of pyrite and sulphoxyanion oxidation. In the Rio Tinto, high delta O-18((SO4-H2O)) values (11.2-13.6 parts per thousand) occur in concentrated (Fe = 172-829 mM), low pH (0.88-1.4), ferrous iron (68-91% of total Fe) waters and are most simply explained by a mechanism involving a dissolved sulphite intermediate, sulphite-water oxygen equilibrium exchange and finally sulphite oxidation to sulphate with O-2. In contrast, drainage from large waste blocks of acid volcanic tuff with pyritiferous veins also had low pH (1.7). but had a low delta O-18((SO4-H2O)) value of 4.0 parts per thousand and high concentrations of ferric iron (Fe(III) = 185 mM, total Fe = 186 mM), suggesting a pathway where ferric iron is the primary oxidant, water is the primary source of oxygen in the sulphate and where sulphate is released directly from the pyrite surface. However, problems remain with the sulphite-water oxygen exchange model and recommendations are therefore made for future experiments to refine our understanding of oxygen isotopes in pyrite oxidation. (C) 2009 Elsevier B.V. All rights reserved.

Salt marsh morphodynamics: An investigation of tidal flows and marsh channel equilibrium

Although numerous field studies have evaluated flow and transport processes in salt marsh channels, the overall role of channels in delivering and removing material from salt marsh platforms is still poorly characterised. In this paper, we consider this issue based on a numerical hydrodynamic model for a prototype marsh system and on a field survey of the cross-sectional geometry of a marsh channel network. Results of the numerical simulations indicate that the channel transfers approximately three times the volume of water that would be estimated from mass balance considerations alone. Marsh platform roughness exerts a significant influence on the partitioning of discharge between the channel and the marsh platform edge, alters flow patterns on the marsh platform due to its effects on channel-to-platform transfer and also controls the timing of peak discharge relative to marsh-edge overtopping. Although peak channel discharges and velocities are associated with the flood tide and marsh inundation, a larger volume of water is transferred by the channel during ebb flows, a portion of which transfer takes place after the tidal height is below the marsh platform. Detailed surveys of the marsh channels crossing a series of transects at Upper Stiffkey Marsh, north Norfolk, England, show that the total channel cross-sectional area increases linearly with catchment area in the inner part of the marsh, which is consistent with the increase in shoreward tidal prism removed by the channels. Toward the marsh edge, however, a deficit in the total cross-sectional area develops, suggesting that discharge partitioning between the marsh channels and the marsh platform edge may also be expressed in the morphology of marsh channel systems.

An integrative model of amino acid metabolism in the liver of the lactating dairy cow

The objective of this work was to construct a dynamic model of hepatic amino acid metabolism in the lactating dairy cow that could be parameterized using net flow data from in vivo experiments. The model considers 22 amino acids, ammonia, urea, and 13 energetic metabolites, and was parameterized using a steady-state balance model and two in vivo, net flow experiments conducted with mid-lactation dairy cows. Extracellular flows were derived directly from the observed data. An optimization routine was used to derive nine intracellular flows. The resulting dynamic model was found to be stable across a range of inputs suggesting that it can be perturbed and applied to other physiological states. Although nitrogen was generally in balance, leucine was in slight deficit compared to predicted needs for export protein synthesis, suggesting that an alternative source of leucine (e.g. peptides) was utilized. Simulations of varying glucagon concentrations indicated that an additional 5 mol/d of glucose could be synthesized at the reference substrate concentrations and blood flows. The increased glucose production was supported by increased removal from blood of lactate, glutamate, aspartate, alanine, asparagine, and glutamine. As glucose Output increased, ketone body and acetate release increased while CO2 release declined. The pattern of amino acids appearing in hepatic vein blood was affected by changes in amino acid concentration in portal vein blood, portal blood flow rate and glucagon concentration, with methionine and phenylalanine being the most affected of essential amino acids. Experimental evidence is insufficient to determine whether essential amino acids are affected by varying gluconeogenic demands. (C) 2004 Published by Elsevier Ltd.

Influence of drought-induced acidification on the mobility of dissolved organic carbon in peat soils

A strong relationship between dissolved organic carbon (DOC) and sulphate (SO42−) dynamics under drought conditions has been revealed from analysis of a 10-year time series (1993–2002). Soil solution from a blanket peat at 10 cm depth and stream water were collected at biweekly and weekly intervals, respectively, by the Environmental Change Network at Moor House-Upper Teesdale National Nature Reserve in the North Pennine uplands of Britain. DOC concentrations in soil solution and stream water were closely coupled, displaying a strong seasonal cycle with lowest concentrations in early spring and highest in late summer/early autumn. Soil solution DOC correlated strongly with seasonal variations in soil temperature at the same depth 4-weeks prior to sampling. Deviation from this relationship was seen, however, in years with significant water table drawdown (>−25 cm), such that DOC concentrations were up to 60% lower than expected. Periods of drought also resulted in the release of SO42−, because of the oxidation of inorganic/organic sulphur stored in the peat, which was accompanied by a decrease in pH and increase in ionic strength. As both pH and ionic strength are known to control the solubility of DOC, inclusion of a function to account for DOC suppression because of drought-induced acidification accounted for more of the variability of DOC in soil solution (R2=0.81) than temperature alone (R2=0.58). This statistical model of peat soil solution DOC at 10 cm depth was extended to reproduce 74% of the variation in stream DOC over this period. Analysis of annual budgets showed that the soil was the main source of SO42− during droughts, while atmospheric deposition was the main source in other years. Mass balance calculations also showed that most of the DOC originated from the peat. The DOC flux was also lower in the drought years of 1994 and 1995, reflecting low DOC concentrations in soil and stream water. The analysis presented in this paper suggests that lower concentrations of DOC in both soil and stream waters during drought years can be explained in terms of drought-induced acidification. As future climate change scenarios suggest an increase in the magnitude and frequency of drought events, these results imply potential for a related increase in DOC suppression by episodic acidification.