Estimating daily methane production in individual cattle with irregular feed intake patterns from short-term methane emission measurements

Spot measurements of methane emission rate (n = 18 700) by 24 Angus steers fed mixed rations from GrowSafe feeders were made over 3- to 6-min periods by a GreenFeed emission monitoring (GEM) unit. The data were analysed to estimate daily methane production (DMP; g/day) and derived methane yield (MY; g/kg dry matter intake (DMI)). A one-compartment dose model of spot emission rate v. time since the preceding meal was compared with the models of Wood (1967) and Dijkstra et al. (1997) and the average of spot measures. Fitted values for DMP were calculated from the area under the curves. Two methods of relating methane and feed intakes were then studied: the classical calculation of MY as DMP/DMI (kg/day); and a novel method of estimating DMP from time and size of preceding meals using either the data for only the two meals preceding a spot measurement, or all meals for 3 days prior. Two approaches were also used to estimate DMP from spot measurements: fitting of splines on a 'per-animal per-day' basis and an alternate approach of modelling DMP after each feed event by least squares (using Solver), summing (for each animal) the contributions from each feed event by best-fitting a one-compartment model. Time since the preceding meal was of limited value in estimating DMP. Even when the meal sizes and time intervals between a spot measurement and all feeding events in the previous 72 h were assessed, only 16.9% of the variance in spot emission rate measured by GEM was explained by this feeding information. While using the preceding meal alone gave a biased (underestimate) of DMP, allowing for a longer feed history removed this bias. A power analysis taking into account the sources of variation in DMP indicated that to obtain an estimate of DMP with a 95% confidence interval within 5% of the observed 64 days mean of spot measures would require 40 animals measured over 45 days (two spot measurements per day) or 30 animals measured over 55 days. These numbers suggest that spot measurements could be made in association with feed efficiency tests made over 70 days. Spot measurements of enteric emissions can be used to define DMP but the number of animals and samples are larger than are needed when day-long measures are made.

Field measurement of beef pen manure methane and nitrous oxide reveals a surprise for inventory calculations

Few data exist on direct greenhouse gas emissions from pen manure at beef feedlots. However, emission inventories attempt to account for these emissions. This study used a large chamber to isolate N2O and CH4 emissions from pen manure at two Australian commercial beef feedlots (stocking densities, 13-27 m(2) head) and related these emissions to a range of potential emission control factors, including masses and concentrations of volatile solids, NO3-, total N, NH4+, and organic C (OC), and additional factors such as total manure mass, cattle numbers, manure pack depth and density, temperature, and moisture content. Mean measured pen N2O emissions were 0.428 kg ha(-1) d(-1) (95% confidence interval [CI], 0.252-0.691) and 0.00405 kg ha(-1) d(-1) (95% CI, 0.00114-0.0110) for the northern and southern feedlots, respectively. Mean measured CH4 emission was 0.236 kg ha(-1) d(-1) (95% CI, 0.163-0.332) for the northern feedlot and 3.93 kg ha(-1) d(-1) (95% CI, 2.58-5.81) for the southern feedlot. Nitrous oxide emission increased with density, pH, temperature, and manure mass, whereas negative relationships were evident with moisture and OC. Strong relationships were not evident between N2O emission and masses or concentrations of NO3- or total N in the manure. This is significant because many standard inventory calculation protocols predict N2O emissions using the mass of N excreted by the animal.

A universal equation to predict methane production of forage-fed cattle in Australia

The methods for estimating methane emissions from cattle as used in the Australian national inventory are based on older data that have now been superseded by a large amount of more recent data. Recent data suggested that the current inventory emissions estimates can be improved. To address this issue, a total of 1034 individual animal records of daily methane production (MP) was used to reassess the relationship between MP and each of dry matter intake (DMI) and gross energy intake (GEI). Data were restricted to trials conducted in the past 10 years using open-circuit respiration chambers, with cattle fed forage-based diets (forage >70%). Results from diets considered to inhibit methanogenesis were omitted from the dataset. Records were obtained from dairy cattle fed temperate forages (220 records), beef cattle fed temperate forages (680 records) and beef cattle fed tropical forages (133 records). Relationships were very similar for all three production categories and single relationships for MP on a DMI or GEI basis were proposed for national inventory purposes. These relationships were MP (g/day) = 20.7 (±0.28) × DMI (kg/day) (R2 = 0.92, P < 0.001) and MP (MJ/day) = 0.063 (±0.008) × GEI (MJ/day) (R2 = 0.93, P < 0.001). If the revised MP (g/day) approach is used to calculate Australia’s national inventory, it will reduce estimates of emissions of forage-fed cattle by 24%. Assuming a global warming potential of 25 for methane, this represents a 12.6 Mt CO2-e reduction in calculated annual emissions from Australian cattle.

Archaea, Bacteria, and methane production along environmental gradients in fens and bogs

Metanogeenit ovat hapettomissa oloissa eläviä arkkien pääryhmään kuuluvia mikrobeja, joiden ainutlaatuisen aineenvaihdunnan seurauksena syntyy metaania. Ilmakehässä metaani on voimakas kasvihuonekaasu. Yksi suurimmista luonnon metaanilähteistä ovat kosteikot. Pohjoisten soiden metaanipäästöt vaihtelevat voimakkaasti eri soiden välillä ja yhden suon sisälläkin, riippuen muun muassa vuodenajasta, suotyypistä ja kasvillisuudesta. Väitöskirjatyössä tutkittiin metaanipäästöjen vaihtelun mikrobiologista taustaa. Tutkimuksessa selvitettiin suotyypin, vuodenajan, tuhkalannoituksen ja turvesyvyyden vaikutusta metanogeeniyhteisöihin sekä metaanintuottoon kolmella suomalaisella suolla. Lisäksi tutkittiin ei-metanogeenisia arkkeja ja bakteereita, koska ne muodostavat metaanin tuoton lähtöaineet osana hapetonta hajotusta. Mikrobiyhteisöt analysoitiin DNA- ja RNA-lähtöisillä, polymeraasiketjureaktioon (PCR) perustuvilla menetelmillä. Merkkigeeneinä käytettiin metaanin tuottoon liittyvää mcrA-geeniä sekä arkkien ja bakteerien ribosomaalista 16S RNA-geeniä. Metanogeeniyhteisöt ja metaanintuotto erosivat huomattavasti happaman ja vähäravinteisen rahkasuon sekä ravinteikkaampien sarasoiden välillä. Rahkasuolta löytyi lähes yksinomaan Methanomicrobiales-lahkon metanogeeneja, jotka tuottavat metaania vedystä ja hiilidioksidista. Sarasoiden metanogeeniyhteisöt olivat monimuotoisempia, ja niillä esiintyi myös asetaattia käyttäviä metanogeeneja. Vuodenaika vaikutti merkittävästi metaanintuottoon. Talvella havaittiin odottamattoman suuri metaanintuottopotentiaali sekä viitteitä aktiivisista metanogeeneista. Arkkiyhteisön koostumus sen sijaan vaihteli vain vähän. Tuhkalannoitus, jonka tarkoituksena on edistää puiden kasvua ojitetuilla soilla, ei merkittävästi vaikuttanut metaanintuottoon tai -tuottajiin. Ojitetun suon yhteisöt kuitenkin muuttuivat turvesyvyyden mukaan. Vertailtaessa erilaisia PCR-menetelmiä todettiin, että kolmella mcrA-geeniin kohdistuvalla alukeparilla havaittiin pääosin samat ojitetun suon metanogeenit, mutta lajien runsaussuhteet riippuvat käytetyistä alukkeista. Soilla havaitut bakteerit kuuluivat pääjaksoihin Deltaproteobacteria, Acidobacteria ja Verrucomicrobia. Lisäksi löydettiin Crenarchaeota-pääjakson ryhmiin 1.1c ja 1.3 kuuluvia ei-metanogeenisia arkkeja. Tulokset ryhmien esiintymisestä hapettomassa turpeessa antavat lähtökohdan selvittää niiden mahdollisia vuorovaikutuksia metanogeenien kanssa. Tutkimuksen tulokset osoittivat, että metanogeeniyhteisön koostumus heijastaa metaanintuottoon vaikuttavia kemiallisia tai kasvillisuuden vaihteluita kuten suotyyppiä. Soiden metanogeenien ja niiden fysiologian parempi tuntemus voi auttaa ennustamaan ympäristömuutosten vaikutusta soiden metaanipäästöihin.

Oxidation of bis-(2-hydroxy-1-naphthyl)methane with 2,3-dichloro-5,6-dicyano-1,4-benzoquinone. X-Ray crystal structure of a novel product

A novel compound obtained by the oxidation of the title compound with 2,3-dichloro-5,6-dicyano-1,4-benzoquinone has been assigned structure (5) on the basis of spectral data and X-ray crystal structure analysis.

Carbon dioxide and methane exchange between a boreal pristine lake and the atmosphere

Lakes serve as sites for terrestrially fixed carbon to be remineralized and transferred back to the atmosphere. Their role in regional carbon cycling is especially important in the Boreal Zone, where lakes can cover up to 20% of the land area. Boreal lakes are often characterized by the presence of a brown water colour, which implies high levels of dissolved organic carbon from the surrounding terrestrial ecosystem, but the load of inorganic carbon from the catchment is largely unknown. Organic carbon is transformed to methane (CH4) and carbon dioxide (CO2) in biological processes that result in lake water gas concentrations that increase above atmospheric equilibrium, thus making boreal lakes as sources of these important greenhouse gases. However, flux estimates are often based on sporadic sampling and modelling and actual flux measurements are scarce. Thus, the detailed temporal flux dynamics of greenhouse gases are still largely unknown. ----- One aim here was to reveal the natural dynamics of CH4 and CO2 concentrations and fluxes in a small boreal lake. The other aim was to test the applicability of a measuring technique for CO2 flux, i.e. the eddy covariance (EC) technique, and a computational method for estimation of primary production and community respiration, both commonly used in terrestrial research, in this lake. Continuous surface water CO2 concentration measurements, also needed in free-water applications to estimate primary production and community respiration, were used over two open water periods in a study of CO2 concentration dynamics. Traditional methods were also used to measure gas concentration and fluxes. The study lake, Valkea-Kotinen, is a small, humic, headwater lake within an old-growth forest catchment with no local anthropogenic disturbance and thus possible changes in gas dynamics reflect the natural variability in lake ecosystems. CH4 accumulated under the ice and in the hypolimnion during summer stratification. The surface water CH4 concentration was always above atmospheric equilibrium and thus the lake was a continuous source of CH4 to the atmosphere. However, the annual CH4 fluxes were small, i.e. 0.11 mol m-2 yr-1, and the timing of fluxes differed from that of other published estimates. The highest fluxes are usually measured in spring after ice melt but in Lake Valkea-Kotinen CH4 was effectively oxidised in spring and highest effluxes occurred in autumn after summer stratification period. CO2 also accumulated under the ice and the hypolimnetic CO2 concentration increased steadily during stratification period. The surface water CO2 concentration was highest in spring and in autumn, whereas during the stable stratification it was sometimes under atmospheric equilibrium. It showed diel, daily and seasonal variation; the diel cycle was clearly driven by light and thus reflected the metabolism of the lacustrine ecosystem. However, the diel cycle was sometimes blurred by injection of hypolimnetic water rich in CO2 and the surface water CO2 concentration was thus controlled by stratification dynamics. The highest CO2 fluxes were measured in spring, autumn and during those hypolimnetic injections causing bursts of CO2 comparable with the spring and autumn fluxes. The annual fluxes averaged 77 (±11 SD) g C m-2 yr-1. In estimating the importance of the lake in recycling terrestrial carbon, the flux was normalized to the catchment area and this normalized flux was compared with net ecosystem production estimates of -50 to 200 g C m-2 yr-1 from unmanaged forests in corresponding temperature and precipitation regimes in the literature. Within this range the flux of Lake Valkea-Kotinen yielded from the increase in source of the surrounding forest by 20% to decrease in sink by 5%. The free water approach gave primary production and community respiration estimates of 5- and 16-fold, respectively, compared with traditional bottle incubations during a 5-day testing period in autumn. The results are in parallel with findings in the literature. Both methods adopted from the terrestrial community also proved useful in lake studies. A large percentage of the EC data was rejected, due to the unfulfilled prerequisites of the method. However, the amount of data accepted remained large compared with what would be feasible with traditional methods. Use of the EC method revealed underestimation of the widely used gas exchange model and suggests simultaneous measurements of actual turbulence at the water surface with comparison of the different gas flux methods to revise the parameterization of the gas transfer velocity used in the models.

Temperature and concentration dependence of adsorption properties of methane in NaY: a molecular dynamics study

Molecular dynamics calculations on methane sorbed in NaY (Si/Al = 3.0) employing realistic methane-methane and methane-zeolite intermolecular potential functions at different temperatures (50, 150, 220, and 300 K) and concentrations (2, 4, 6, and 8 molecules/cage) are reported. The thermodynamic results are in agreement with the available experimental data. Guest-guest and guest-host radial distribution functions (rdfs), energy distribution functions, distribution of cage occupancy, center-of-cage-center-of-mass (coc-com) rdfs, velocity autocorrelation functions for com and angular motion and the Fourier transformed power spectra, and diffusion coefficients are presented as a function of temperature and concentration. At 50 K, methane is localized near the adsorption site. Site-site migration and essentially free rotational motion are observed at 150 K. Molecules preferentially occupy the region near the inner surface of the alpha-cage. The vibrational frequencies for the com of methane shift toward higher values with decreasing temperature and increasing adsorbate concentration. The observed frequencies for com motion are 36, 53, and 85 cm-1 and for rotational motion at 50 K, 95 and 150 cm-1 in agreement with neutron scattering data. The diffusion coefficients show a type I behavior as a function of loading in agreement with NMR measurements. Cage-to-cage diffusion is found to be always mediated by the surface.

Field intercomparison of four methane gas analysers suitable for eddy covariance flux measurements

Eddy covariance (EC)-flux measurement technique is based on measurement of turbulent motions of air with accurate and fast measurement devices. For instance, in order to measure methane flux a fast methane gas analyser is needed which measures methane concentration at least ten times in a second in addition to a sonic anemometer, which measures the three wind components with the same sampling interval. Previously measurement of methane flux was almost impossible to carry out with EC-technique due to lack of fast enough gas analysers. However during the last decade new instruments have been developed and thus methane EC-flux measurements have become more common. Performance of four methane gas analysers suitable for eddy covariance measurements are assessed in this thesis. The assessment and comparison was performed by analysing EC-data obtained during summer 2010 (1.4.-26.10.) at Siikaneva fen. The four participating methane gas analysers are TGA-100A (Campbell Scientific Inc., USA), RMT-200 (Los Gatos Research, USA), G1301-f (Picarro Inc., USA) and Prototype-7700 (LI-COR Biosciences, USA). RMT-200 functioned most reliably throughout the measurement campaign and the corresponding methane flux data had the smallest random error. In addition, methane fluxes calculated from data obtained from G1301-f and RMT-200 agree remarkably well throughout the measurement campaign. The calculated cospectra and power spectra agree well with corresponding temperature spectra. Prototype-7700 functioned only slightly over one month in the beginning of the measurement campaign and thus its accuracy and long-term performance is difficult to assess.

Sensitivity of wetland methane emissions to model assumptions: application and model testing against site observations

Abstract. Methane emissions from natural wetlands and rice paddies constitute a large proportion of atmospheric methane, but the magnitude and year-to-year variation of these methane sources is still unpredictable. Here we describe and evaluate the integration of a methane biogeochemical model (CLM4Me; Riley et al., 2011) into the Community Land Model 4.0 (CLM4CN) in order to better explain spatial and temporal variations in methane emissions. We test new functions for soil pH and redox potential that impact microbial methane production in soils. We also constrain aerenchyma in plants in always-inundated areas in order to better represent wetland vegetation. Satellite inundated fraction is explicitly prescribed in the model because there are large differences between simulated fractional inundation and satellite observations. A rice paddy module is also incorporated into the model, where the fraction of land used for rice production is explicitly prescribed. The model is evaluated at the site level with vegetation cover and water table prescribed from measurements. Explicit site level evaluations of simulated methane emissions are quite different than evaluating the grid cell averaged emissions against available measurements. Using a baseline set of parameter values, our model-estimated average global wetland emissions for the period 1993–2004 were 256 Tg CH4 yr−1, and rice paddy emissions in the year 2000 were 42 Tg CH4 yr−1. Tropical wetlands contributed 201 Tg CH4 yr−1, or 78 % of the global wetland flux. Northern latitude (>50 N) systems contributed 12 Tg CH4 yr−1. We expect this latter number may be an underestimate due to the low high-latitude inundated area captured by satellites and unrealistically low high-latitude productivity and soil carbon predicted by CLM4. Sensitivity analysis showed a large range (150–346 Tg CH4 yr−1) in predicted global methane emissions. The large range was sensitive to: (1) the amount of methane transported through aerenchyma, (2) soil pH (± 100 Tg CH4 yr−1), and (3) redox inhibition (± 45 Tg CH4 yr−1).

Synthesis and structures of oxyanion encapsulated copper(I)-dppm complexes (dppm = bis(diphenylphosphino) methane)

Copper(I)-dppm complexes encapsulating the oxyanions ClO4-, NO3-, CH3C6H4CO2-, SO42-, and WO42- have been synthesized either by reduction of the corresponding Cu(II) salts and treatment with dppm, or by treating the complex [Cu-2(dppm)(2)(dmcn)(3)](BF4)(2) (1) (dmcn = dimethyl cyanamide) with the respective anion. The isolated complexes [Cu-2(dppm)(2)(dmcn)(2)(ClO4)] (ClO4) (2), [Cu-2(dppm)(2)(dmcn)(2)(NO3)] (NO3) (3), Cu-2(dppm)(2)(NO3)(2) (4), [Cu-2(dppm)(2)(CH3C6H4CO2)(2)]dmcn.2THF (5), Cu-2(dppm)(2)(SO4) (6), and [Cu-3(dppm)(3)(Cl)(WO4)] 0.5H(2)O (7) have been characterized by IR, H-1 and P-31{H-1} NMR, UV-vis, and emission spectroscopy. The solid-state molecular structure of complexes 1, 2, 4, and 7 were determined by single-crystal X-ray diffraction. Pertinent crystal data are as follows: for 1, monoclinic P2(1)/c, a = 11.376(10) Angstrom, b = 42.503(7) Angstrom, c = 13.530(6) Angstrom, beta = 108.08(2)degrees, V = 6219(3) Angstrom(3), Z = 4; for 2, monoclinic P2(1)/c, a = 21.600(3) Angstrom, b = 12.968(3) Angstrom, c = 23.050(3) Angstrom, beta = 115.97(2)degrees, V = 5804(17) Angstrom(3), Z = 4; for 4, triclinic , a = 10.560(4) Angstrom, b = 10.553(3) Angstrom, c = 22.698(3) Angstrom, alpha = 96.08(2)degrees, beta = 96.03(2)degrees, gamma = 108.31(2)degrees, V = 2362(12) Angstrom(3), Z = 2; and for 7, orthorhombic P2(1)2(1)2(1), a = 14.407(4) Angstrom, b = 20.573(7) Angstrom, c = 24.176(6) Angstrom, V = 7166(4) Angstrom(3), Z = 4. Analyses of the crystallographic and spectroscopic data of these complexes reveal the nature of interactions between the Cu-I-dppm core and oxyanion. The anchoring of the oxyanion to the Cu-n(dppm)(n) unit is primarily through coordination to the metal, but the noncovalent C-H ... O interactions between the methylene and phenyl protons of the dppm and oxygen atoms of the oxyanion play a significant role. The solid-state emission spectra for complexes 1-6 are very similar but different from 7. In CDCl3 solution, addition of ClO4- or NO3- (as their tetrabutylammonium salts) to 1 establishes a rapid equilibrium between the anion-complexed and uncomplexed forms. The association constant values for ClO4- and NO3- have been estimated from the P-31{H-1} NMR spectra.

A comparative study of Pt/CeO2 catalysts for catalytic partial oxidation of methane to syngas for application in fuel cell electric vehicles

In the framework of a project aimed at developing a reliable hydrogen generator for mobile polymer electrolyte fuel cells (PEFCs), particular emphasis has been addressed to the analysis of catalysts able to assure high activity and stability in transient operations (frequent start-up and shut-down cycles). In this paper, the catalytic performance of 1 at.% Pt/ceria samples prepared by coprecipitation, impregnation and combustion, has been evaluated in the partial oxidation of methane. Methane conversion and hydrogen selectivity of 96 and 99%, respectively, associated with high stability during 100h of reaction under operative conditions (start-up and shut-down cycles), have been obtained. (C) 2002 Elsevier Science B.V. All rights reserved.

An investigation of carbon nanotubes obtained from the decomposition of methane over reduced Mg1− xM xAl2O4 spinel catalysts

Carbon nanotubes produced by the treatment of Mg1−xMxAl2O4 (M = Fe, Co, or Ni; x = 0.1, 0.2, 0.3, or 0.4) spinels with an H2–CH4 mixture at 1070 °C have been investigated systematically. The grains of the oxide-metal composite particles are uniformly covered by a weblike network of carbon nanotube bundles, several tens of micrometers long, made up of single-wall nanotubes with a diameter close to 4 nm. Only the smallest metal particles (<5 nm) are involved in the formation of the nanotubes. A macroscopic characterization method involving surface area measurements and chemical analysis has been developed in order to compare the different nanotube specimens. An increase in the transition metal content of the catalyst yields more carbon nanotubes (up to a metal content of 10.0 wt% or x = 0.3), but causes a decrease in carbon quality. The best compromise is to use 6.7 wt% of metal (x = 0.2) in the catalyst. Co gives superior results with respect to both the quantity and quality of the nanotubes. In the case of Fe, the quality is notably hampered by the formation of Fe3C particles.

Influence of the Methane-Zeolite A Interaction Potential on the Concentration Dependence of Self-diffusivity

Studies on the diffusion of methane in a zeolite structure type LTA (as per IZA nomenclature) have indicated that different types of methane zeolite potentials exist in the literature in which methane is treated within the united-atom model. One set of potentials, referred to as model A, has a methane oxygen diameter of 3.14 angstrom, while another set of potential parameters, model B, employs a larger value of 3.46 angstrom. Fritzsche and co-workers (1993) have shown that these two potentials lead to two distinctly different energetic barriers for the passage of methane through the eight-ring window in the cation-free form of zeolite A. Here, we compute the variation of the self-diffusivity (D) with loading (c) for these two types of potentials and show that this slight variation in the diameter changes the concentration dependence qualitatively: thus, D decreases monotonically with c for model A, while D increases and goes through a maximum before finally decreasing for model B. This effect and the surprising congruence of the diffusion coefficients for both models at high loadings is examined in detail at the molecular level. Simulations for different temperatures reveal the Arrhenius behaviour of the self-diffusion coefficient. The apparent activation energy is found to vary with the loading. We conclude that beside the cage-to-cage jumps, which are essential for the migration of the guest molecules, at high concentrations migration within the cage and guest guest interactions with other molecules become increasingly dominant influences on the diffusion coefficient and make the guest zeolite interaction less important for both model A and model B.

Analysis of Parameter Values in the van der Waals and Platteeuw Theory for Methane Hydrates Using Monte Carlo Molecular Simulations

The van der Waals and Platteuw (vdVVP) theory has been successfully used to model the thermodynamics of gas hydrates. However, earlier studies have shown that this could be due to the presence of a large number of adjustable parameters whose values are obtained through regression with experimental data. To test this assertion, we carry out a systematic and rigorous study of the performance of various models of vdWP theory that have been proposed over the years. The hydrate phase equilibrium data used for this study is obtained from Monte Carlo molecular simulations of methane hydrates. The parameters of the vdWP theory are regressed from this equilibrium data and compared with their true values obtained directly from simulations. This comparison reveals that (i) methane-water interactions beyond the first cage and methane-methane interactions make a significant contribution to the partition function and thus cannot be neglected, (ii) the rigorous Monte Carlo integration should be used to evaluate the Langmuir constant instead of the spherical smoothed cell approximation, (iii) the parameter values describing the methane-water interactions cannot be correctly regressed from the equilibrium data using the vdVVP theory in its present form, (iv) the regressed empty hydrate property values closely match their true values irrespective of the level of rigor in the theory, and (v) the flexibility of the water lattice forming the hydrate phase needs to be incorporated in the vdWP theory. Since methane is among the simplest of hydrate forming molecules, the conclusions from this study should also hold true for more complicated hydrate guest molecules.

Methane and carbon dioxide adsorption on edge-functionalized graphene: a comparative DFT study

With a view towards optimizing gas storage and separation in crystalline and disordered nanoporous carbon-based materials, we use ab initio density functional theory calculations to explore the effect of chemical functionalization on gas binding to exposed edges within model carbon nanostructures. We test the geometry, energetics, and charge distribution of in-plane and out-of-plane binding of CO2 and CH4 to model zigzag graphene nanoribbons edge-functionalized with COOH, OH, NH2, H2PO3, NO2, and CH3. Although different choices for the exchange-correlation functional lead to a spread of values for the binding energy, trends across the functional groups are largely preserved for each choice, as are the final orientations of the adsorbed gas molecules. We find binding of CO2 to exceed that of CH4 by roughly a factor of two. However, the two gases follow very similar trends with changes in the attached functional group, despite different molecular symmetries. Our results indicate that the presence of NH2, H2PO3, NO2, and COOH functional groups can significantly enhance gas binding, making the edges potentially viable binding sites in materials with high concentrations of edge carbons. To first order, in-plane binding strength correlates with the larger permanent and induced dipole moments on these groups. Implications for tailoring carbon structures for increased gas uptake and improved CO2/CH4 selectivity are discussed. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4736568]