Manageable Workflows for Processing Parallel Sequencing Data

ACM Computing Classification System (1998): D.2.11, D.1.3, D.3.1, J.3, C.2.4.

Accurate estimation of isoelectric point of protein and peptide based on amino acid sequences

Motivation: In any macromolecular polyprotic system - for example protein, DNA or RNA - the isoelectric point - commonly referred to as the pI - can be defined as the point of singularity in a titration curve, corresponding to the solution pH value at which the net overall surface charge - and thus the electrophoretic mobility - of the ampholyte sums to zero. Different modern analytical biochemistry and proteomics methods depend on the isoelectric point as a principal feature for protein and peptide characterization. Protein separation by isoelectric point is a critical part of 2-D gel electrophoresis, a key precursor of proteomics, where discrete spots can be digested in-gel, and proteins subsequently identified by analytical mass spectrometry. Peptide fractionation according to their pI is also widely used in current proteomics sample preparation procedures previous to the LC-MS/MS analysis. Therefore accurate theoretical prediction of pI would expedite such analysis. While such pI calculation is widely used, it remains largely untested, motivating our efforts to benchmark pI prediction methods. Results: Using data from the database PIP-DB and one publically available dataset as our reference gold standard, we have undertaken the benchmarking of pI calculation methods. We find that methods vary in their accuracy and are highly sensitive to the choice of basis set. The machine-learning algorithms, especially the SVM-based algorithm, showed a superior performance when studying peptide mixtures. In general, learning-based pI prediction methods (such as Cofactor, SVM and Branca) require a large training dataset and their resulting performance will strongly depend of the quality of that data. In contrast with Iterative methods, machine-learning algorithms have the advantage of being able to add new features to improve the accuracy of prediction. Contact: yperez@ebi.ac.uk Availability and Implementation: The software and data are freely available at https://github.com/ypriverol/pIR. Supplementary information: Supplementary data are available at Bioinformatics online.

Electromagnetic design and loss calculations of a 1.12-MW high-speed permanent-magnet motor for compressor applications

Electromagnetic design of a 1.12-MW, 18 000-r/min high-speed permanent-magnet motor (HSPMM) is carried out based on the analysis of pole number, stator slot number, rotor outer diameter, air-gap length, permanent magnet material, thickness, and pole arc. The no-load and full-load performance of the HSPMM is investigated in this paper by using 2-D finite element method (FEM). In addition, the power losses in the HSPMM including core loss, winding loss, rotor eddy current loss, and air friction loss are predicted. Based on the analysis, a prototype motor is manufactured and experimentally tested to verify the machine design.

Mechanistic studies on the degradation of gasoline oxygenates by advanced oxidation technologies and the reaction of 4-methyl-1,2,4-triazoline-3,5-dione with tetracyclopropylethylene

Gasoline oxygenates (MTBE, methyl tert-butyl ether; DIPE, di-isopropyl ether; ETBE, ethyl tert-butyl ether; TAME, tert-amyl ether) are added to gasoline to boost octane and enhance combustion. The combination of large scale use, high water solubility and only minor biodegradability has now resulted in a significant gasoline oxygenate contamination occurring in surface, ground, and drinking water systems. Combination of hydroxyl radical formation and the pyrolytic environment generated by ultrasonic irradiation (665 kHz) leads to the rapid degradation of MTBE and other gasoline oxygenates in aqueous media. ^ The presence of oxygen promotes the degradation processes by rapid reaction with carbon centered radicals indicating radical processes involving O 2 are significant pathways. A number of the oxidation products were identified. The formation of products (alcohols, ketones, aldehydes, esters, peroxides, etc) could be rationalized by mechanisms which involve hydrogen abstraction by OH radical and/or pyrolysis to form carboncentered radicals which react with oxygen and follow standard oxidation chain processes. ^ The reactions of N-substituted R-triazolinediones (RTAD; R = CH 3 or phenyl) have attracted considerable interest because they exhibit a number of unusual mechanistic characteristics that are analogous to the reactions of singlet oxygen (1O2) and offer an easy way to provide C-N bond(s) formation. The reactions of triazolinedione with olefins have been widely studied and aziridinium imides are generally accepted to be the reactive intermediates. ^ We observed the rapid formation of an unusual intermediate upon mixing tetracyclopropylethylene with 4-methyl-1,2,4-triazoline-3,5-dione in CDCl 3. Detailed characterization by NMR (proton, 13C, 2-D NMRs) indicates the intermediate is 5,5,6,6-tetracyclopropyl-3-methyl-5,6-dihydro-oxazolo[3,2- b][1,2,4]-triazolium-2-olate. Such products are extremely rare and have not been studied. Upon warming the intermediate is converted to 2 + 2 diazetidine (major) and ene product (minor). ^ To further explore the kinetics and dynamics of the reaction activation energies were obtained using Arrhenius plots. Activation energies for the formation of the intermediate from reactants, and 2+2 adduct from the intermediate were determined as 7.48 kcal moll and 19.8 kcal mol−1 with their pre-exponential values of 2.24 × 105 dm 3 mol−1 sec−1 and 2.75 × 108 sec−1, respectively, meaning net slow reactions because of low pre-exponential values caused by steric hindrance. ^

Seizure detection in subdural recordings using nonlinear decision functions

This dissertation proposed a new approach to seizure detection in intracranial EEG recordings using nonlinear decision functions. It implemented well-established features that were designed to deal with complex signals such as brain recordings, and proposed a 2-D domain of analysis. Since the features considered assume both the time and frequency domains, the analysis was carried out both temporally and as a function of different frequency ranges in order to ascertain those measures that were most suitable for seizure detection. In retrospect, this study established a generalized approach to seizure detection that works across several features and across patients. ^ Clinical experiments involved 8 patients with intractable seizures that were evaluated for potential surgical interventions. A total of 35 iEEG data files collected were used in a training phase to ascertain the reliability of the formulated features. The remaining 69 iEEG data files were then used in the testing phase. ^ The testing phase revealed that the correlation sum is the feature that performed best across all patients with a sensitivity of 92% and an accuracy of 99%. The second best feature was the gamma power with a sensitivity of 92% and an accuracy of 96%. In the frequency domain, all of the 5 other spectral bands considered, revealed mixed results in terms of low sensitivity in some frequency bands and low accuracy in other frequency bands, which is expected given that the dominant frequencies in iEEG are those of the gamma band. In the time domain, other features which included mobility, complexity, and activity, all performed very well with an average a sensitivity of 80.3% and an accuracy of 95%. ^ The computational requirement needed for these nonlinear decision functions to be generated in the training phase was extremely long. It was determined that when the duration dimension was rescaled, the results improved and the convergence rates of the nonlinear decision functions were reduced dramatically by more than a 100 fold. Through this rescaling, the sensitivity of the correlation sum improved to 100% and the sensitivity of the gamma power to 97%, which meant that there were even less false negatives and false positives detected. ^

Hydrogeophysical characterization of anisotropy in the Biscayne Aquifer using geophysical methods

The anisotropy of the Biscayne Aquifer which serves as the source of potable water for Miami-Dade County was investigated by applying geophysical methods. Electrical resistivity imaging, self potential and ground penetration radar techniques were employed in both regional and site specific studies. In the regional study, electrical anisotropy and resistivity variation with depth were investigated with azimuthal square array measurements at 13 sites. The observed coefficient of electrical anisotropy ranged from 1.01 to 1.36. The general direction of measured anisotropy is uniform for most sites and trends W-E or SE-NW irrespective of depth. Measured electrical properties were used to estimate anisotropic component of the secondary porosity and hydraulic anisotropy which ranged from 1 to 11% and 1.18 to 2.83 respectively. 1-D sounding analysis was used to models the variation of formation resistivity with depth. Resistivities decreased from NW (close to the margins of the everglades) to SE on the shores of Biscayne Bay. Porosity calculated from Archie's law, ranged from 18 to 61% with higher values found along the ridge. Higher anisotropy, porosities and hydraulic conductivities were on the Atlantic Coastal Ridge and lower values at low lying areas west of the ridge. The cause of higher anisotropy and porosity is attributed to higher dissolution rates of the oolitic facies of the Miami Formation composing the ridge. The direction of minimum resistivity from this study is similar to the predevelopment groundwater flow direction indicated in published modeling studies. Detailed investigations were carried out to evaluate higher anisotropy at West Perrine Park located on the ridge and Snapper Creek Municipal well field where the anisotropy trend changes with depth. The higher anisotropy is attributed to the presence of solution cavities oriented in the E-SE direction on the ridge. Similarly, the change in hydraulic anisotropy at the well field might be related to solution cavities, the surface canal and groundwater extraction wells.^

Application of biomarkers and compound specific stable isotopes for the assessment of hydrology as a driver of organic matter dynamics in the Everglades ecosystem

The Everglades is a sub-tropical coastal wetland characterized among others by its hydrological features and deposits of peat. Formation and preservation of organic matter in soils and sediments in this wetland ecosystem is critical for its sustainability and hydrological processes are important divers in the origin, transport and fate of organic matter. With this in mind, organic matter dynamics in the greater Florida Everglades was studied though various organic geochemistry techniques, especially biomarkers, bulk and compound specific δ13C and δD isotope analysis. The main objectives were focused on how different hydrological regimes in this ecosystem control organic matter dynamics, such as the mobilization of particulate organic matter (POM) in freshwater marshes and estuaries, and how organic geochemistry techniques can be applied to reconstruct Everglades paleo-hydrology. For this purpose organic matter in typical vegetation, floc, surface soils, soil cores, and estuarine suspended particulates were characterized in samples selected along hydrological gradients in the Water Conservation Area 3, Shark River Slough and Taylor Slough. ^ This research focused on three general themes: (1) Assessment of the environmental dynamics and source-specific particulate organic carbon export in a mangrove-dominated estuary. (2) Assessment of the origin, transport and fate of organic matter in freshwater marsh. (3) Assessment of historical changes in hydrological conditions in the Everglades (paleo-hydrology) though biomarkes and compound specific isotope analyses. This study reports the first estimate of particulate organic carbon loss from mangrove ecosystems in the Everglades, provides evidence for particulate organic matter transport with regards to the formation of ridge and slough landscapes in the Everglades, and demonstrates the applicability of the combined biomarker and compound-specific stable isotope approach as a means to generate paleohydrological data in wetlands. The data suggests that: (1) Carbon loss from mangrove estuaries is roughly split 50/50 between dissolved and particulate carbon; (2) hydrological remobilization of particulate organic matter from slough to ridge environments may play an important role in the maintenance of the Everglades freshwater landscape; and (3) Historical changes in hydrology have resulted in significant vegetation shifts from historical slough type vegetation to present ridge type vegetation. ^

A model of school managerial control : the systematic analysis of managerial behaviors, processes, and indicators

The theoretical construct of control has been defined as necessary (Etzioni, 1965), ubiquitous (Vickers, 1967), and on-going (E. Langer, 1983). Empirical measures, however, have not adequately given meaning to this potent construct, especially within complex organizations such as schools. Four stages of theory-development and empirical testing of school building managerial control using principals and teachers working within the nation's fourth largest district are presented in this dissertation as follows: (1) a review and synthesis of social science theories of control across the literatures of organizational theory, political science, sociology, psychology, and philosophy; (2) a systematic analysis of school managerial activities performed at the building level within the context of curricular and instructional tasks; (3) the development of a survey questionnaire to measure school building managerial control; and (4) initial tests of construct validity including inter-item reliability statistics, principal components analyses, and multivariate tests of significance. The social science synthesis provided support of four managerial control processes: standards, information, assessment, and incentives. The systematic analysis of school managerial activities led to further categorization between structural frequency of behaviors and discretionary qualities of behaviors across each of the control processes and the curricular and instructional tasks. Teacher survey responses (N=486) reported a significant difference between these two dimensions of control, structural frequency and discretionary qualities, for standards, information, and assessments, but not for incentives. The descriptive model of school managerial control suggests that (1) teachers perceive structural and discretionary managerial behaviors under information and incentives more clearly than activities representing standards or assessments, (2) standards are primarily structural while assessments are primarily qualitative, (3) teacher satisfaction is most closely related to the equitable distribution of incentives, (4) each of the structural managerial behaviors has a qualitative effect on teachers, and that (5) certain qualities of managerial behaviors are perceived by teachers as distinctly discretionary, apart from school structure. The variables of teacher tenure and school effectiveness reported significant effects on school managerial control processes, while instructional levels (elementary, junior, and senior) and individual school differences were not found to be significant for the construct of school managerial control.

Representações sociais de familiares sobre o atendimento das emergências psiquiátricas

The reformist movements in the field of mental health have pointed battle flags, among which the prioritization of production of mental health care out of the asylum environment should be highlighted, aiming the reduction of psychiatric beds, greater control over the hospitalization, family co-participation and the rescue of the citizenship of the social players involved. With the progressive reduction of asylum beds, associated with a lot of structural problems in the health services, the occurrence of crises outside the hospital environment has been increasingly frequent, thus giving the family an important therapeutic role. In face of this scenario, there is an urgent need to understand the social construction of the care for psychiatric emergencies, identifying the meanings assigned by family members to their constituent aspects. This study seeks to answer the following research question: what are the social representations of family members about the care of psychiatric emergencies in the city of Mossoró, Rio Grande do Norte? Therefore, the aim is to analyze the social representations of family members about the care of psychiatric emergencies in the city of Mossoró, Rio Grande do Norte. This is an exploratory and descriptive study, with a mixed approach, making use of multimethods: for collection, the semi-structured interview and the Technique of Free Association of Words; for data analysis, the Thematic Analysis of Bardin and its steps was used, with the informational support of the softwares ALCESTE (Analyse Lexicale par Contexte d'un Ensemble de Segments de Texte) and Iramuteq (Interface de R pour les Analyses Multidimensionnelles de Textes et de Questionnaires); and the theoretical support of social representations. The study participants totaled 72, and they were selected from the following criteria: older than18 years, with degree of kinship with users suffering from some mental and behavioral disorder, and who have already witnessed a situation of crisis, rescued by the SAMU or other means and taken to the psychiatric hospital or general emergency room. Preliminary results point out: 1.Previous note of the research project with the aim to disseminate it in the scientific community and ensure the intellectual property of the work; 2.The contextual analysis of the care for emergencies in the study place. Reflection about the phenomenon provide a name to the care for the psychiatric emergencies, which is called immediate context; the technical and operational aspects that influence the care, as a specific/ general context; and mental health policies in Brazil are identified as metacontext; 3. The systematic review from randomized clinical trials in the databases PubMed, COCHRANE, LILACS, SciELO and SCIRUS, with the use of the descriptors: ‘Physical restraint’, ‘Psychiatric emergency services’, ‘Restraint’, ‘Physical and Emergency Services’, ‘Psychiatric’. Only one work met the search protocol criteria: a short-term essay that records limited results about the proportion of people who are in restraint and seclusion. It does not show statistically significant results in relation to indications, contraindications and risks of the use of physical restraint; 4. The social representations of the care for psychiatric emergencies. The study results point to the presence of five thematic categories: 1. feeling in the face of the crisis/care; 2. thoughts and perspectives about the crisis/care; 3. centrality of care in the medical- medication-hospitalization triad; 4. the thinking/acting in the face of the use of physical restraint and police force; 5. periodicity of crises. The central core of the representation is in the first category, whilst the peripheral elements are in the third and fifth categories. The contrast zone is in the second and fourth categories. The sadness is the most prominent element of the structure. The social representations about the care for psychiatric crises are at a time of transition between the hegemonic and reformist models, with the traditional aspects being predominant, but already showing peripheral and contrast elements that point to a possible change in the representational field.

Radionuclides measured on 27 water bottle profiles during POLARSTERN cruise ANT-XXIV/3

As part of the GEOTRACES Polarstern expedition ANT XXIV/3 (ZERO and DRAKE) we have measured the vertical distribution of 234Th on sections through the Antarctic Circumpolar Current along the zero meridian and in Drake Passage and on an EW section through the Weddell Sea. Steady state export fluxes of 234Th from the upper 100m, derived from the depletion of 234Th with respect to its parent 238U, ranged from 621±105 dpm/m**2/d to 1773±90 dpm/m**2/d. This 234Th flux was converted into an export flux of organic carbon ranging from 3.1-13.2 mmolC/m**2/d (2.1-9.0 mmolC/m**2/d) using POC/234Th ratio of bulk (respectively >50 µm) suspended particles at the export depth (100 m). Non-steady state fluxes assuming zero flux under ice cover were up to 23% higher. In addition, particulate and dissolved 234Th were measured underway in high resolution in the surface water with a semi-automated procedure. Particulate 234Th in surface waters is inversely correlated with light transmission and pCO2 and positively with fluorescence and optical backscatter and is interpreted as a proxy for algal biomass. High resolution underway mapping of particulate and dissolved 234Th in surface water shows clearly where trace elements are absorbed by plankton and where they are exported to depth. Quantitative determination of the export flux requires the full 234Th profile since surface depletion and export flux become decoupled through changes in wind mixed layer depth and in contribution to export from subsurface layers. In a zone of very low algal abundance (54-58 °S at the zero meridian), confirmed by satellite Chl-a data, the lowest carbon export of the ACC was observed, allowing Fe and Mn to maintain their highest surface concentrations (Klunder et al., this issue, Middag et al., this issue). An ice-edge bloom that had developed in Dec/Jan in the zone 60-65 °S as studied during the previous leg (Strass et al., in prep) had caused a high export flux at 64.5 °S when we visited the area two months later (Feb/March). The ice-edge bloom had then shifted south to 65-69 °S evident from uptake of CO2 and dissolved Fe, Mn and 234Th, without causing export yet. In this way, the parallel analysis of 234Th can help to explain the scavenging behaviour of other trace elements.

Caractérisation tridimensionnelle de l’amplitude articulaire de l’épaule

L’épaule est l’articulation la plus mobile et la plus instable du corps humain dû à la faible quantité de contraintes osseuses et au rôle des tissus mous qui lui confèrent au moins une dizaine de degrés de liberté. La mobilité de l’épaule est un facteur de performance dans plusieurs sports. Mais son instabilité engendre des troubles musculo-squelettiques, dont les déchirures de la coiffe des rotateurs sont fréquentes et les plus handicapantes. L’évaluation de l’amplitude articulaire est un indice commun de la fonction de l’épaule, toutefois elle est souvent limitée à quelques mesures planaires pour lesquelles les degrés de liberté varient indépendamment les uns des autres. Ces valeurs utilisées dans les modèles de simulation musculo-squelettiques peuvent amener à des solutions non physiologiques. L’objectif de cette thèse était de développer des outils pour la caractérisation de la mobilité articulaire tri-dimensionnelle de l’épaule, en passant par i) fournir une méthode et son approche expérimentale pour évaluer l’amplitude articulaire tridimensionnelle de l’épaule incluant des interactions entre les degrés de liberté ; ii) proposer une représentation permettant d’interpréter les données tri-dimensionnelles obtenues; iii) présenter des amplitudes articulaires normalisées, iv) implémenter une amplitude articulaire tridimensionnelle au sein d’un modèle de simulation numérique afin de générer des mouvements sportifs optimaux plus réalistes; v) prédire des amplitudes articulaires sécuritaires et vi) des exercices de rééducation sécuritaires pour des patients ayant subi une réparation de la coiffe des rotateurs. i) Seize sujets ont été réalisé séries de mouvements d’amplitudes maximales actifs avec des combinaisons entre les différents degrés de liberté de l’épaule. Un système d’analyse du mouvement couplé à un modèle cinématique du membre supérieur a été utilisé pour estimer les cinématiques articulaires tridimensionnelles. ii) L’ensemble des orientations définies par une séquence de trois angles a été inclus dans un polyèdre non convexe représentant l’espace de mobilité articulaire prenant en compte les interactions entre les degrés de liberté. La combinaison des séries d’élévation et de rotation est recommandée pour évaluer l’amplitude articulaire complète de l’épaule. iii) Un espace de mobilité normalisé a également été défini en englobant les positions atteintes par au moins 50% des sujets et de volume moyen. iv) Cet espace moyen, définissant la mobilité physiologiques, a été utilisé au sein d’un modèle de simulation cinématique utilisé pour optimiser la technique d’un élément acrobatique de lâcher de barres réalisée par des gymnastes. Avec l’utilisation régulière de limites articulaires planaires pour contraindre la mobilité de l’épaule, seulement 17% des solutions optimales sont physiologiques. En plus, d’assurer le réalisme des solutions, notre contrainte articulaire tridimensionnelle n’a pas affecté le coût de calculs de l’optimisation. v) et vi) Les seize participants ont également réalisé des séries d’amplitudes articulaires passives et des exercices de rééducation passifs. La contrainte dans l’ensemble des muscles de la coiffe des rotateurs au cours de ces mouvements a été estimée à l’aide d’un modèle musculo-squelettique reproduisant différents types et tailles de déchirures. Des seuils de contrainte sécuritaires ont été utilisés pour distinguer les amplitudes de mouvements risquées ou non pour l’intégrité de la réparation chirurgicale. Une taille de déchirure plus grande ainsi que les déchirures affectant plusieurs muscles ont contribué à réduire l’espace de mobilité articulaire sécuritaire. Principalement les élévations gléno-humérales inférieures à 38° et supérieures à 65°, ou réalisées avec le bras maintenu en rotation interne engendrent des contraintes excessives pour la plupart des types et des tailles de blessure lors de mouvements d’abduction, de scaption ou de flexion. Cette thèse a développé une représentation innovante de la mobilité de l’épaule, qui tient compte des interactions entre les degrés de liberté. Grâce à cette représentation, l’évaluation clinique pourra être plus exhaustive et donc élargir les possibilités de diagnostiquer les troubles de l’épaule. La simulation de mouvement peut maintenant être plus réaliste. Finalement, nous avons montré l’importance de personnaliser la rééducation des patients en termes d’amplitude articulaire, puisque des exercices passifs de rééducation précoces peuvent contribuer à une re-déchirure à cause d’une contrainte trop importante qu’ils imposent aux tendons.

Synthèse et fonctionnalisation des 4-aminopyrrole-2-carboxylates

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Synthèse et fonctionnalisation des 4-aminopyrrole-2-carboxylates

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Seawater carbonate chemistry and calcification during tropical reef studies at Moorea (French Polynesia), 1996

Community metabolism was investigated using a Lagrangian flow respirometry technique on 2 reef flats at Moorea (French Polynesia) during austral winter and Yonge Reef (Great Barrier Reef) during austral summer. The data were used to estimate related air-sea CO2 disequilibrium. A sine function did not satisfactorily model the diel light curves and overestimated the metabolic parameters. The ranges of community gross primary production and respiration (Pg and R; 9 to 15 g C m-2 d-1) were within the range previously reported for reef flats, and community net calcification (G; 19 to 25 g CaCO3 m-2 d-1) was higher than the 'standard' range. The molar ratio of organic to inorganic carbon uptake was 6:1 for both sites. The reef flat at Moorea displayed a higher rate of organic production and a lower rate of calcification compared to previous measurements carried out during austral summer. The approximate uncertainty of the daily metabolic parameters was estimated using a procedure based on a Monte Carlo simulation. The standard errors of Pg,R and Pg/R expressed as a percentage of the mean are lower than 3% but are comparatively larger for E, the excess production (6 to 78%). The daily air-sea CO2 flux (FCO2) was positive throughout the field experiments, indicating that the reef flats at Moorea and Yonge Reef released CO2 to the atmosphere at the time of measurement. FCO2 decreased as a function of increasing daily irradiance.

Seawater carbonate chemistry and calcification at Anse des Cuivres, NW Mediterranean, 2004

The community metabolism of a shallow infralittoral ecosystem dominated by the calcareous macroalgae Corallina elongata was investigated in Marseilles (NW Mediterranean), by monitoring hourly changes of seawater pH and total alkalinity over 6 d in February 2000. Fair weather conditions prevailed over the study period as indicated by oceanographic (temperature, salinity, and current velocity and direction) and meteorological variables, which validated the standing water hypothesis. This temperate ecosystem exhibited high community gross primary production (GPP = 519 ± 106 mmol C m-2 d-1, n = 6) and also supported high rates of community respiration (R). As a result, the system was slightly autotrophic (net community production, NCP = 20 mmol C m-2 d-1), with a GPP/R ratio of 1.06. NCP exhibited circadian variations with 2- to 3-fold changes in community respiration, both in the light and in the dark. Rates of net community calcification also exhibited circadian variations, with positive rates (up to 24 mmol CaCO3 m-2 h-1) for irradiance values >300 W m-2 (about 1380 µmol photon m-2 s-1). Below this irradiance threshold, net community dissolution prevailed. Daily net calcification (G) was on average 8 mmol CaCO3 m-2 d-1. CO2 fluxes generated by primary production, respiration, and calcification suggest that the study site was a potential atmospheric CO2 sink of 15 mmol CO2 m-2 d-1 at the time of measurement.