LIPGENE food-exchange model for alteration of dietary fat quantity and quality in free-living participants from eight European countries

Controlled human intervention trials are required to confirm the hypothesis that dietary fat quality may influence insulin action. The aim was to develop a food-exchange model, suitable for use in free-living volunteers, to investigate the effects of four experimental diets distinct in fat quantity and quality: high SFA (HSFA); high MUFA (HMUFA) and two low-fat (LF) diets, one supplemented with 1.24g EPA and DHA/d (LFn-3). A theoretical food-exchange model was developed. The average quantity of exchangeable fat was calculated as the sum of fat provided by added fats (spreads and oils), milk, cheese, biscuits, cakes, buns and pastries using data from the National Diet and Nutrition Survey of UK adults. Most of the exchangeable fat was replaced by specifically designed study foods. Also critical to the model was the use of carbohydrate exchanges to ensure the diets were isoenergetic. Volunteers from eight centres across Europe completed the dietary intervention. Results indicated that compositional targets were largely achieved with significant differences in fat quantity between the high-fat diets (39.9 (SEM 0.6) and 38.9 (SEM 0.51) percentage energy (%E) from fat for the HSFA and HMUFA diets respectively) and the low-fat diets (29.6 (SEM 0.6) and 29.1 (SEM 0.5) %E from fat for the LF and LFn-3 diets respectively) and fat quality (17.5 (SEM 0.3) and 10.4 (SEM 0.2) %E front SFA and 12.7 (SEM 0.3) and 18.7 (SEM 0.4) %E MUFA for the HSFA and HMUFA diets respectively). In conclusion, a robust, flexible food-exchange model was developed and implemented successfully in the LIPGENE dietary intervention trial.

Successful manipulation of the quality and quantity of fat and carbohydrate consumed by free-living individuals using a food exchange model

Our objective in this study was to develop and implement an effective intervention strategy to manipulate the amount and composition of dietary fat and carbohydrate (CHO) in free-living individuals in the RISCK study. The study was a randomized, controlled dietary intervention study that was conducted in 720 participants identified as higher risk for or with metabolic syndrome. All followed a 4-wk run-in reference diet [high saturated fatty acids (SF)/high glycemic index (GI)]. Volunteers were randomized to continue this diet for a further 24 wk or to I of 4 isoenergetic prescriptions [high monounsaturated fatty acids (MUFA)/high GI; high MUFA/low GI; low fat (LF)/high GI; and LF/low GI]. We developed a food exchange model to implement each diet. Dietary records and plasma phospholipid fatty acids were used to assess the effectiveness of the intervention strategy. Reported fat intake from the LF diets was significantly reduced to 28% of energy (%E) compared with 38% E from the HM and LF diets. SF intake was successfully decreased in the HM and LF diets was similar to 10% E compared with 17% E in the reference diet (P = 0.001). Dietary MUFA in the HIM diets was similar to 17% E, significantly higher than in the reference (12% E) and LF diets (10% E) (P = 0.001). Changes in plasma phospholipid fatty acids provided further evidence for the successful manipulation of fat intake. The GI of the HGI and LGI arms differed by similar to 9 points (P = 0.001). The food exchange model provided an effective dietary strategy for the design and implementation across multiple sites of 5 experimental diets with specific targets for the proportion of fat and CHO. J. Nutr. 139: 1534-1540, 2009.

The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

Development of a food-exchange model to replace saturated fat with MUFAS and n-6 PUFAS in adults at moderate cardiovascular risk

The recommendation to reduce saturated fatty acid (SFA) consumption to ≤10% of total energy (%TE) is a key public health target aimed at lowering cardiovascular disease (CVD) risk. Replacement of SFA with unsaturated fats may provide greater benefit than replacement with carbohydrates, yet the optimal type of fat is unclear. The aim was to develop a flexible food-exchange model to investigate the effects of substituting SFAs with monounsaturated fatty acids (MUFAs) or n-6 (ω-6) polyunsaturated fatty acids (PUFAs) on CVD risk factors. In this parallel study, UK adults aged 21-60 y with moderate CVD risk (50% greater than the population mean) were identified using a risk assessment tool (n = 195; 56% females). Three 16-wk isoenergetic diets of specific fatty acid (FA) composition (%TE SFA:%TE MUFA:%TE n-6 PUFA) were designed using spreads, oils, dairy products, and snacks as follows: 1) SFA-rich diet (17:11:4; n = 65); 2) MUFA-rich diet (9:19:4; n = 64); and 3) n-6 PUFA-rich diet (9:13:10; n = 66). Each diet provided 36%TE total fat. Dietary targets were broadly met for all intervention groups, reaching 17.6 ± 0.4%TE SFA, 18.5 ± 0.3%TE MUFA, and 10.4 ± 0.3%TE n-6 PUFA in the respective diets, with significant overall diet effects for the changes in SFA, MUFA, and n-6 PUFA between groups (P < 0.001). There were no differences in the changes of total fat, protein, carbohydrate, and alcohol intake or anthropometric measures between groups. Plasma phospholipid FA composition showed changes from baseline in the proportions of total SFA, MUFA, and n-6 PUFA for each diet group, with significant overall diet effects for total SFA and MUFA between groups (P < 0.001). In conclusion, successful implementation of the food-exchange model broadly achieved the dietary target intakes for the exchange of SFA with MUFA or n-6 PUFA with minimal disruption to the overall diet in a free-living population. This trial was registered at clinicaltrials.gov as NCT01478958.

Trinuclear cu(ll) complexes containing peripheral ketonic oxygen bridges and a mu(3)-OH core: Syntheses, crystal structures, spectroscopic and magnetic properties

Four new trinuclear copper(II) complexes, [(CuL1)(3)(mu(3)-OH)](ClO4)(2)center dot H2O (1), [(CuL2)(3)(mu(3)-OH)](CIO4)(2) (2), [(CuL3)(3)-(mu(3)-OH)](ClO4)(4)center dot H2O (3), and [(CuL4)(3)(mu(3)-OH)](ClO4)(2)center dot H2O (4), where HL1 = 8-amino-4,7,7-trimethyl-5-azaoct-3-en-2-one, HL2 = 7-amino-4-methyl-5-azaoct-3-en-2-one, HL3 = 7(ethylamino)-4-methyl-5-azahept-3-en-2-one, and HL4 = 4-methyl-7-(methylamino)-5-azahept-3-en-2-one, have been derived from the four tridentate Schiff bases (HL1, HL2, HL3, and HL4) and structurally characterized by X-ray crystallography. For all compounds, the cationic part is trinuclear with a CU3OH core held by three carbonyl oxygen bridges between each pair of copper(II) atoms. The copper atoms are five-coordinate with a distorted square-pyramidal geometry; the equatorial plane consists of the bridging oxygen atom of the central OH group together with three atoms (N, N, O) from one ligand whereas an oxygen atom of a second ligand occupies the axial position. Magnetic measurements have been performed in the 2-300 K temperature range. The experimental data could be satisfactorily reproduced by using an isotropic exchange model, H = -J(S1S2+S2S3+S1S3) yielding as best-fit parameters: J = -66.7 and g = 2.19 for 1, J = -36.6 and g = 2.20 for 2, J = -24.5 and g = 2.20 for 3, and J = -14.9 and g = 2.05 for 4. EPR spectra at low temperature show the existence of spin frustration in complexes 3 and 4, but it has not been possible to carry out calculations of the antisymmetric exchange parameter, G, from magnetic data. In frozen methanolic solution, at 4 K, hyperfine splitting in all complexes and spin frustration in complex 4 seem to be confirmed. ((c) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2005)

Trinuclear Cu(II) complexes containing peripheral ketonic oxygen bridges and a mu(3)-OH core: steric influence on their structures and existence

Two tridentate Schiff bases, HL1(6-amino-3-methyl-1-phenyl-4-azahex-2-en-1-one), and HL2 (6-atnino-3,6-dimethyl-1-phenyl-4-azahex-2-en-1-one) on reaction with Cu(II) perchlorate in the presence of triethyl amine yielded two new trinuclear copper(II) complexes, [(CuL1)(3)(mu(3)-OH)](ClO4)(2) (1) and [(CuL2)(3)(mu(3)-OH)](ClO4)(2) center dot 0.75H(2)O (2), whereas another tridentate ligand HL3 (7-amino-3-methyl-1-phenyl-4-azahept-2-en-1-one) underwent hydrolysis under the same reaction conditions to result in the formation of a mononuclear complex, [Cu(bn)(pn)ClO4] (3) [where bn = 1-benzoylacetonate and pn = 1,3-propanediamine]. All three complexes have been characterized by X-ray crystallography. For both 1 and 2 the cationic part is trinuclear with a [Cu3OH] core held by three carbonyl oxygen bridges between each pair of copper(II) atoms. The structure of 3 is a monomer with a chelating 1,3-propanediamine and a benzoyl acetone moiety. Magnetic measurements of I and 2 have been performed in the 2-300 K temperature range. The experimental data could be satisfactorily reproduced by using an isotropic exchange model, H = -J(S1S2 + S2S3 + S1S3), yielding as best fit parameters: J = -25.6 cm(-1), g = 2.21 for 1 and J = 11.2 cm(-1), g = 2.10 for 2. The EPR spectra at low temperature could be indicative of spin frustration in complex 1. (C) 2006 Elsevier Ltd. All rights reserved.

Use of a reduced Schiff-Base ligand to prepare novel chloro-bridged chains of rare Cu(II) trinuclear complexes with mixed azido/oxo and chloro/oxo bridges

Two mixed bridged one-dimensional (1D) polynuclear complexes, [Cu3L2(mu(1,1)-N-3)(2)(mu-Cl)Cl](n) (1) and {[Cu3L2(mu-Cl)(3)Cl]center dot 0.46CH(3)OH}(n), (2), have been synthesized using the tridentate reduced Schiff-base ligand HL (2-[(2-dimethylamino-ethylamino)-methyl]-phenol). The complexes have been characterized by X-ray structural analyses and variable-temperature magnetic susceptibility measurements. In both complexes the basic trinuclear angular units are joined together by weak chloro bridges to form a 1D chain. The trinuclear structure of 1 is composed of two terminal square planar [Cu(L)(mu(1,1)-N-3)] units connected by a central Cu(II) atom through bridging nitrogen atoms of end-on azido ligands and the phenoxo oxygen atom of the tridentate ligand. These four coordinating atoms along with a chloride ion form a distorted trigonal bipyramidal geometry around the central Cu(II). The structure of 2 is similar; the only difference being a Cl bridge replacing the mu(1,1)-N-3 bridge in the trinuclear unit. The magnetic properties of both trinuclear complexes can be very well reproduced with a simple linear symmetrical trimer model (H = JS(i)S(i+1)) with only one intracluster exchange coupling (J) including a weak intertrimer interaction (.j) reproduced with the molecular field approximation. This model provides very satisfactory fits for both complexes in the whole temperature range with the following parameters: g = 2.136(3), J = 93.9(3) cm(-1) and zj= -0.90(3) cm(-1) (z = 2) for 1 and g = 2.073(7), J = -44.9(4) cm(-1) and zJ = -1.26(6) cm(-1) (z = 2) for 2.

SATgene dietary model to implement diets of differing fat composition in prospectively genotyped groups (apoE) using commercially available foods

esponse to dietary fat manipulation is highly heterogeneous, yet generic population-based recommendations aimed at reducing the burden of CVD are given. The APOE epsilon genotype has been proposed to be an important determinant of this response. The present study reports on the dietary strategy employed in the SATgenɛ (SATurated fat and gene APOE) study, to assess the impact of altered fat content and composition on the blood lipid profile according to the APOE genotype. A flexible dietary exchange model was developed to implement three isoenergetic diets: a low-fat (LF) diet (target composition: 24 % of energy (%E) as fat, 8 %E SFA and 59 %E carbohydrate), a high-saturated fat (HSF) diet (38 %E fat, 18 %E SFA and 45 %E carbohydrate) and a HSF-DHA diet (HSF diet with 3 g DHA/d). Free-living participants (n 88; n 44 E3/E3 and n 44 E3/E4) followed the diets in a sequential design for 8 weeks, each using commercially available spreads, oils and snacks with specific fatty acid profiles. Dietary compositional targets were broadly met with significantly higher total fat (42·8 %E and 41·0 %E v. 25·1 %E, P ≤ 0·0011) and SFA (19·3 %E and 18·6 %E v. 8·33 %E, P ≤ 0·0011) intakes during the HSF and HSF-DHA diets compared with the LF diet, in addition to significantly higher DHA intake during the HSF-DHA diet (P ≤ 0·0011). Plasma phospholipid fatty acid analysis revealed a 2-fold increase in the proportion of DHA after consumption of the HSF-DHA diet for 8 weeks, which was independent of the APOE genotype. In summary, the dietary strategy was successfully implemented in a free-living population resulting in well-tolerated diets which broadly met the dietary targets set.

Crop yield reduction in the tropics under climate change: Processes and uncertainties

Many modelling studies examine the impacts of climate change on crop yield, but few explore either the underlying bio-physical processes, or the uncertainty inherent in the parameterisation of crop growth and development. We used a perturbed-parameter crop modelling method together with a regional climate model (PRECIS) driven by the 2071-2100 SRES A2 emissions scenario in order to examine processes and uncertainties in yield simulation. Crop simulations used the groundnut (i.e. peanut; Arachis hypogaea L.) version of the General Large-Area Model for annual crops (GLAM). Two sets of GLAM simulations were carried out: control simulations and fixed-duration simulations, where the impact of mean temperature on crop development rate was removed. Model results were compared to sensitivity tests using two other crop models of differing levels of complexity: CROPGRO, and the groundnut model of Hammer et al. [Hammer, G.L., Sinclair, T.R., Boote, K.J., Wright, G.C., Meinke, H., and Bell, M.J., 1995, A peanut simulation model: I. Model development and testing. Agron. J. 87, 1085-1093]. GLAM simulations were particularly sensitive to two processes. First, elevated vapour pressure deficit (VPD) consistently reduced yield. The same result was seen in some simulations using both other crop models. Second, GLAM crop duration was longer, and yield greater, when the optimal temperature for the rate of development was exceeded. Yield increases were also seen in one other crop model. Overall, the models differed in their response to super-optimal temperatures, and that difference increased with mean temperature; percentage changes in yield between current and future climates were as diverse as -50% and over +30% for the same input data. The first process has been observed in many crop experiments, whilst the second has not. Thus, we conclude that there is a need for: (i) more process-based modelling studies of the impact of VPD on assimilation, and (ii) more experimental studies at super-optimal temperatures. Using the GLAM results, central values and uncertainty ranges were projected for mean 2071-2100 crop yields in India. In the fixed-duration simulations, ensemble mean yields mostly rose by 10-30%. The full ensemble range was greater than this mean change (20-60% over most of India). In the control simulations, yield stimulation by elevated CO2 was more than offset by other processes-principally accelerated crop development rates at elevated, but sub-optimal, mean temperatures. Hence, the quantification of uncertainty can facilitate relatively robust indications of the likely sign of crop yield changes in future climates. (C) 2007 Elsevier B.V. All rights reserved.

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.

Mechanisms of inverse agonist action at D-2 dopamine receptors

1 Mechanisms of inverse agonist action at the D-2(short) dopamine receptor have been examined. 2 Discrimination of G-protein-coupled and -uncoupled forms of the receptor by inverse agonists was examined in competition ligand-binding studies versus the agonist [H-3]NPA at a concentration labelling both G-protein-coupled and -uncoupled receptors. 3 Competition of inverse agonists versus [H-3] NPA gave data that were fitted best by a two-binding site model in the absence of GTP but by a one-binding site model in the presence of GTP. K-i values were derived from the competition data for binding of the inverse agonists to G-protein-uncoupled and -coupled receptors. K-coupled and K-uncoupled were statistically different for the set of compounds tested ( ANOVA) but the individual values were different in a post hoc test only for (+)-butaclamol. 4 These observations were supported by simulations of these competition experiments according to the extended ternary complex model. 5 Inverse agonist efficacy of the ligands was assessed from their ability to reduce agonist-independent [S-35]GTPγ S binding to varying degrees in concentration-response curves. Inverse agonism by (+)-butaclamol and spiperone occurred at higher potency when GDP was added to assays, whereas the potency of (-)-sulpiride was unaffected. 6 These data show that some inverse agonists ((+)-butaclamol, spiperone) achieve inverse agonism by stabilising the uncoupled form of the receptor at the expense of the coupled form. For other compounds tested, we were unable to define the mechanism.

Investigation of the mechanism of agonist and inverse agonist action at D-2 dopamine receptors

This study investigated, for the D-2 dopamine receptor, the relation between the ability of agonists and inverse agonists to stabilise different states of the receptor and their relative efficacies. K-i values for agonists were determined in competition, versus the binding of the antagonist [H-3]spiperone. Competition data were fitted best by a two-binding site model (with the exception of bromocriptine, for which a one-binding site model provided the best fit) and agonist affinities for the higher (K-h) (G protein-coupled) and lower affinity (K-l) (G protein-uncoupled) sites determined. Ki values for agonists were also determined in competition versus the binding of the agonist [H-3]N-propylnorapomorphine (NPA) to provide a second estimate of K-h,. Maximal agonist effects (E-max) and their potencies (EC50) were determined from concentration-response curves for agonist stimulation of guanosine-5'-O-(3-[S-32] thiotriphosphate) ([S-35]GTPgammaS) binding. The ability of agonists to stabilise the G protein-coupled state of the receptor (K-l/K-h, determined from ligand-binding assays) did not correlate with either of two measures of relative efficacy (relative E-max, Kl/EC50) of agonists determined in [S-35]GTPgammaS-binding assays, when the data for all of the compounds tested were analysed For a subset of compounds, however, there was a relation between K-l/K-h and E-max.. Competition-binding data versus [H-3]spiperone and [H-3]NPA for a range of inverse agonists were fitted best by a one-binding site model. K-i values for the inverse agonists tested were slightly lower in competition versus [H-3]NPA compared to [H-3]spiperone. These data do not provide support for the idea that inverse agonists act by binding preferentially to the ground state of the receptor. (C) 2004 Elsevier Inc. All rights reserved.

Analysis of molecular determinants of affinity and relative efficacy of a series of R- and S-2-(dipropylamino)tetralins at the 5-HT1A serotonin receptor

1 Factors influencing agonist affinity and relative efficacy have been studied for the 5-HT1A serotonin receptor using membranes of CHO cells expressing the human form of the receptor and a series of R-and S-2-(dipropylamino)tetralins (nonhydroxylated and monohydroxylated (5-OH, 6-OH, 7-OH, 8-OH) species). 2 Ligand binding studies were used to determine dissociation constants for agonist binding to the 5HT(1A) receptor: (a) K-i values for agonists were determined in competition versus the binding of the agonist [H-3]-8-OH DPAT. Competition data were all fitted best by a one-binding site model. (b) K-i values for agonists were also determined in competition versus the binding of the antagonist [H-3]-NAD-199. Competition data were all fitted best by a two-binding site model, and agonist affinities for the higher (K-h) and lower affinity (K-1) sites were determined. 3 The ability of the agonists to activate the 5-HT1A receptor was determined using stimulation of [S-35]-GTPgammaS binding. Maximal effects of agonists (E-max) and their potencies (EC50) were determined from concentration/response curves for stimulation of [S-35]-GTPgammaS binding. 4 K-1/K-h determined from ligand binding assays correlated with the relative efficacy (relative Em) of agonists determined in [S-35]-GTPgammaS binding assays. There was also a correlation between K-1/K-h and K-1/EC50 for agonists determined from ligand binding and [S-35]-GTPgammaS binding assays. 5 Simulations of agonist binding and effect data were performed using the Ternary Complex Model in order to assess the use of K-1/K-h for predicting the relative efficacy of agonists. British Journal of Pharmacology (2003) 138, 1129-1139. doi: 10. 1038/sj.bjp.705085.

Mechanisms of inverse agonist action at D-2 dopamine receptors

1 Mechanisms of inverse agonist action at the D-2(short) dopamine receptor have been examined. 2 Discrimination of G-protein-coupled and -uncoupled forms of the receptor by inverse agonists was examined in competition ligand-binding studies versus the agonist [H-3]NPA at a concentration labelling both G-protein-coupled and -uncoupled receptors. 3 Competition of inverse agonists versus [H-3] NPA gave data that were fitted best by a two-binding site model in the absence of GTP but by a one-binding site model in the presence of GTP. K-i values were derived from the competition data for binding of the inverse agonists to G-protein-uncoupled and -coupled receptors. K-coupled and K-uncoupled were statistically different for the set of compounds tested ( ANOVA) but the individual values were different in a post hoc test only for (+)-butaclamol. 4 These observations were supported by simulations of these competition experiments according to the extended ternary complex model. 5 Inverse agonist efficacy of the ligands was assessed from their ability to reduce agonist-independent [S-35]GTPγ S binding to varying degrees in concentration-response curves. Inverse agonism by (+)-butaclamol and spiperone occurred at higher potency when GDP was added to assays, whereas the potency of (-)-sulpiride was unaffected. 6 These data show that some inverse agonists ((+)-butaclamol, spiperone) achieve inverse agonism by stabilising the uncoupled form of the receptor at the expense of the coupled form. For other compounds tested, we were unable to define the mechanism.

Groundwater residence time distributions in peatlands: implications for peat decomposition and accumulation

Peat soils consist of poorly decomposed plant detritus, preserved by low decay rates, and deep peat deposits are globally significant stores in the carbon cycle. High water tables and low soil temperatures are commonly held to be the primary reasons for low peat decay rates. However, recent studies suggest a thermodynamic limit to peat decay, whereby the slow turnover of peat soil pore water may lead to high concentrations of phenols and dissolved inorganic carbon. In sufficient concentrations, these chemicals may slow or even halt microbial respiration, providing a negative feedback to peat decay. We document the analysis of a simple, one-dimensional theoretical model of peatland pore water residence time distributions (RTDs). The model suggests that broader, thicker peatlands may be more resilient to rapid decay caused by climate change because of slow pore water turnover in deep layers. Even shallow peat deposits may also be resilient to rapid decay if rainfall rates are low. However, the model suggests that even thick peatlands may be vulnerable to rapid decay under prolonged high rainfall rates, which may act to flush pore water with fresh rainwater. We also used the model to illustrate a particular limitation of the diplotelmic (i.e., acrotelm and catotelm) model of peatland structure. Model peatlands of contrasting hydraulic structure exhibited identical water tables but contrasting RTDs. These scenarios would be treated identically by diplotelmic models, although the thermodynamic limit suggests contrasting decay regimes. We therefore conclude that the diplotelmic model be discarded in favor of model schemes that consider continuous variation in peat properties and processes.