Cooperative Institute for Research in Environmental Sciences

The ongoing revolution in high-throughput sequencing continues to democratize the ability of small groups of investigators to map the microbial component of the biosphere. In particular, the coevolution of new sequencing platforms and new software tools allows data acquisition and analysis on an unprecedented scale. Here we report the next stage in this coevolutionary arms race, using the Illumina GAIIx platform to sequence a diverse array of 25 environmental samples and three known "mock communities" at a depth averaging 3.1 million reads per sample. We demonstrate excellent consistency in taxonomic recovery and recapture diversity patterns that were previously reported on the basis of metaanalysis of many studies from the literature (notably, the saline/nonsaline split in environmental samples and the split between host-associated and free-living communities). We also demonstrate that 2,000 Illumina single-end reads are sufficient to recapture the same relationships among samples that we observe with the full dataset. The results thus open up the possibility of conducting large-scale studies analyzing thousands of samples simultaneously to survey microbial communities at an unprecedented spatial and temporal resolution.

DNA sequencing continues to decrease in cost with the Illumina HiSeq2000 generating up to 600 Gb of paired-end 100 base reads in a ten-day run. Here we present a protocol for community amplicon sequencing on the HiSeq2000 and MiSeq Illumina platforms, and apply that protocol to sequence 24 microbial communities from host-associated and free-living environments. A critical question as more sequencing platforms become available is whether biological conclusions derived on one platform are consistent with what would be derived on a different platform. We show that the protocol developed for these instruments successfully recaptures known biological results, and additionally that biological conclusions are consistent across sequencing platforms (the HiSeq2000 versus the MiSeq) and across the sequenced regions of amplicons.

Abstract Black carbon aerosol plays a unique and important role in Earth's climate system. Black carbon is a type of carbonaceous material with a unique combination of physical properties. This assessment provides an evaluation of black‐carbon climate forcing that is comprehensive in its inclusion of all known and relevant processes and that is quantitative in providing best estimates and uncertainties of the main forcing terms: direct solar absorption; influence on liquid, mixed phase, and ice clouds; and deposition on snow and ice. These effects are calculated with climate models, but when possible, they are evaluated with both microphysical measurements and field observations. Predominant sources are combustion related, namely, fossil fuels for transportation, solid fuels for industrial and residential uses, and open burning of biomass. Total global emissions of black carbon using bottom‐up inventory methods are 7500 Gg yr −1 in the year 2000 with an uncertainty range of 2000 to 29000. However, global atmospheric absorption attributable to black carbon is too low in many models and should be increased by a factor of almost 3. After this scaling, the best estimate for the industrial‐era (1750 to 2005) direct radiative forcing of atmospheric black carbon is +0.71 W m −2 with 90% uncertainty bounds of (+0.08, +1.27) W m −2 . Total direct forcing by all black carbon sources, without subtracting the preindustrial background, is estimated as +0.88 (+0.17, +1.48) W m −2 . Direct radiative forcing alone does not capture important rapid adjustment mechanisms. A framework is described and used for quantifying climate forcings, including rapid adjustments. The best estimate of industrial‐era climate forcing of black carbon through all forcing mechanisms, including clouds and cryosphere forcing, is +1.1 W m −2 with 90% uncertainty bounds of +0.17 to +2.1 W m −2 . Thus, there is a very high probability that black carbon emissions, independent of co‐emitted species, have a positive forcing and warm the climate. We estimate that black carbon, with a total climate forcing of +1.1 W m −2 , is the second most important human emission in terms of its climate forcing in the present‐day atmosphere; only carbon dioxide is estimated to have a greater forcing. Sources that emit black carbon also emit other short‐lived species that may either cool or warm climate. Climate forcings from co‐emitted species are estimated and used in the framework described herein. When the principal effects of short‐lived co‐emissions, including cooling agents such as sulfur dioxide, are included in net forcing, energy‐related sources (fossil fuel and biofuel) have an industrial‐era climate forcing of +0.22 (−0.50 to +1.08) W m −2 during the first year after emission. For a few of these sources, such as diesel engines and possibly residential biofuels, warming is strong enough that eliminating all short‐lived emissions from these sources would reduce net climate forcing (i.e., produce cooling). When open burning emissions, which emit high levels of organic matter, are included in the total, the best estimate of net industrial‐era climate forcing by all short‐lived species from black‐carbon‐rich sources becomes slightly negative (−0.06 W m −2 with 90% uncertainty bounds of −1.45 to +1.29 W m −2 ). The uncertainties in net climate forcing from black‐carbon‐rich sources are substantial, largely due to lack of knowledge about cloud interactions with both black carbon and co‐emitted organic carbon. In prioritizing potential black‐carbon mitigation actions, non‐science factors, such as technical feasibility, costs, policy design, and implementation feasibility play important roles. The major sources of black carbon are presently in different stages with regard to the feasibility for near‐term mitigation. This assessment, by evaluating the large number and complexity of the associated physical and radiative processes in black‐carbon climate forcing, sets a baseline from which to improve future climate forcing estimates.

Although researchers have begun cataloging the incredible diversity of bacteria found in soil, we are largely unable to interpret this information in an ecological context, including which groups of bacteria are most abundant in different soils and why. With this study, we examined how the abundances of major soil bacterial phyla correspond to the biotic and abiotic characteristics of the soil environment to determine if they can be divided into ecologically meaningful categories. To do this, we collected 71 unique soil samples from a wide range of ecosystems across North America and looked for relationships between soil properties and the relative abundances of six dominant bacterial phyla (Acidobacteria, Bacteroidetes, Firmicutes, Actinobacteria, alpha-Proteobacteria, and the beta-Proteobacteria). Of the soil properties measured, net carbon (C) mineralization rate (an index of C availability) was the best predictor of phylum-level abundances. There was a negative correlation between Acidobacteria abundance and C mineralization rates (r2 = 0.26, P < 0.001), while the abundances of beta-Proteobacteria and Bacteroidetes were positively correlated with C mineralization rates (r2 = 0.35, P < 0.001 and r2 = 0.34, P < 0.001, respectively). These patterns were explored further using both experimental and meta-analytical approaches. We amended soil cores from a specific site with varying levels of sucrose over a 12-month period to maintain a gradient of elevated C availabilities. This experiment confirmed our survey results: there was a negative relationship between C amendment level and the abundance of Acidobacteria (r2 = 0.42, P < 0.01) and a positive relationship for both Bacteroidetes and beta-Proteobacteria (r2 = 0.38 and 0.70, respectively; P < 0.01 for each). Further support for a relationship between the relative abundances of these bacterial phyla and C availability was garnered from an analysis of published bacterial clone libraries from bulk and rhizosphere soils. Together our survey, experimental, and meta-analytical results suggest that certain bacterial phyla can be differentiated into copiotrophic and oligotrophic categories that correspond to the r- and K-selected categories used to describe the ecological attributes of plants and animals. By applying the copiotroph-oligotroph concept to soil microorganisms we can make specific predictions about the ecological attributes of various bacterial taxa and better understand the structure and function of soil bacterial communities.

Organic aerosol (OA) particles affect climate forcing and human health, but their sources and evolution remain poorly characterized. We present a unifying model framework describing the atmospheric evolution of OA that is constrained by high-time-resolution measurements of its composition, volatility, and oxidation state. OA and OA precursor gases evolve by becoming increasingly oxidized, less volatile, and more hygroscopic, leading to the formation of oxygenated organic aerosol (OOA), with concentrations comparable to those of sulfate aerosol throughout the Northern Hemisphere. Our model framework captures the dynamic aging behavior observed in both the atmosphere and laboratory: It can serve as a basis for improving parameterizations in regional and global models.

Upon delivery, the neonate is exposed for the first time to a wide array of microbes from a variety of sources, including maternal bacteria. Although prior studies have suggested that delivery mode shapes the microbiota's establishment and, subsequently, its role in child health, most researchers have focused on specific bacterial taxa or on a single body habitat, the gut. Thus, the initiation stage of human microbiome development remains obscure. The goal of the present study was to obtain a community-wide perspective on the influence of delivery mode and body habitat on the neonate's first microbiota. We used multiplexed 16S rRNA gene pyrosequencing to characterize bacterial communities from mothers and their newborn babies, four born vaginally and six born via Cesarean section. Mothers' skin, oral mucosa, and vagina were sampled 1 h before delivery, and neonates' skin, oral mucosa, and nasopharyngeal aspirate were sampled <5 min, and meconium <24 h, after delivery. We found that in direct contrast to the highly differentiated communities of their mothers, neonates harbored bacterial communities that were undifferentiated across multiple body habitats, regardless of delivery mode. Our results also show that vaginally delivered infants acquired bacterial communities resembling their own mother's vaginal microbiota, dominated by Lactobacillus, Prevotella, or Sneathia spp., and C-section infants harbored bacterial communities similar to those found on the skin surface, dominated by Staphylococcus, Corynebacterium, and Propionibacterium spp. These findings establish an important baseline for studies tracking the human microbiome's successional development in different body habitats following different delivery modes, and their associated effects on infant health.

Abstract. Secondary organic aerosol (SOA) accounts for a significant fraction of ambient tropospheric aerosol and a detailed knowledge of the formation, properties and transformation of SOA is therefore required to evaluate its impact on atmospheric processes, climate and human health. The chemical and physical processes associated with SOA formation are complex and varied, and, despite considerable progress in recent years, a quantitative and predictive understanding of SOA formation does not exist and therefore represents a major research challenge in atmospheric science. This review begins with an update on the current state of knowledge on the global SOA budget and is followed by an overview of the atmospheric degradation mechanisms for SOA precursors, gas-particle partitioning theory and the analytical techniques used to determine the chemical composition of SOA. A survey of recent laboratory, field and modeling studies is also presented. The following topical and emerging issues are highlighted and discussed in detail: molecular characterization of biogenic SOA constituents, condensed phase reactions and oligomerization, the interaction of atmospheric organic components with sulfuric acid, the chemical and photochemical processing of organics in the atmospheric aqueous phase, aerosol formation from real plant emissions, interaction of atmospheric organic components with water, thermodynamics and mixtures in atmospheric models. Finally, the major challenges ahead in laboratory, field and modeling studies of SOA are discussed and recommendations for future research directions are proposed.

Soils collected across a long-term liming experiment (pH 4.0-8.3), in which variation in factors other than pH have been minimized, were used to investigate the direct influence of pH on the abundance and composition of the two major soil microbial taxa, fungi and bacteria. We hypothesized that bacterial communities would be more strongly influenced by pH than fungal communities. To determine the relative abundance of bacteria and fungi, we used quantitative PCR (qPCR), and to analyze the composition and diversity of the bacterial and fungal communities, we used a bar-coded pyrosequencing technique. Both the relative abundance and diversity of bacteria were positively related to pH, the latter nearly doubling between pH 4 and 8. In contrast, the relative abundance of fungi was unaffected by pH and fungal diversity was only weakly related with pH. The composition of the bacterial communities was closely defined by soil pH; there was as much variability in bacterial community composition across the 180-m distance of this liming experiment as across soils collected from a wide range of biomes in North and South America, emphasizing the dominance of pH in structuring bacterial communities. The apparent direct influence of pH on bacterial community composition is probably due to the narrow pH ranges for optimal growth of bacteria. Fungal community composition was less strongly affected by pH, which is consistent with pure culture studies, demonstrating that fungi generally exhibit wider pH ranges for optimal growth.

Soils harbor enormously diverse bacterial populations, and soil bacterial communities can vary greatly in composition across space. However, our understanding of the specific changes in soil bacterial community structure that occur across larger spatial scales is limited because most previous work has focused on either surveying a relatively small number of soils in detail or analyzing a larger number of soils with techniques that provide little detail about the phylogenetic structure of the bacterial communities. Here we used a bar-coded pyrosequencing technique to characterize bacterial communities in 88 soils from across North and South America, obtaining an average of 1,501 sequences per soil. We found that overall bacterial community composition, as measured by pairwise UniFrac distances, was significantly correlated with differences in soil pH (r = 0.79), largely driven by changes in the relative abundances of Acidobacteria, Actinobacteria, and Bacteroidetes across the range of soil pHs. In addition, soil pH explains a significant portion of the variability associated with observed changes in the phylogenetic structure within each dominant lineage. The overall phylogenetic diversity of the bacterial communities was also correlated with soil pH (R(2) = 0.50), with peak diversity in soils with near-neutral pHs. Together, these results suggest that the structure of soil bacterial communities is predictable, to some degree, across larger spatial scales, and the effect of soil pH on bacterial community composition is evident at even relatively coarse levels of taxonomic resolution.

Abstract The Twentieth Century Reanalysis (20CR) project is an international effort to produce a comprehensive global atmospheric circulation dataset spanning the twentieth century, assimilating only surface pressure reports and using observed monthly sea‐surface temperature and sea‐ice distributions as boundary conditions. It is chiefly motivated by a need to provide an observational dataset with quantified uncertainties for validations of climate model simulations of the twentieth century on all time‐scales, with emphasis on the statistics of daily weather. It uses an Ensemble Kalman Filter data assimilation method with background ‘first guess’ fields supplied by an ensemble of forecasts from a global numerical weather prediction model. This directly yields a global analysis every 6 hours as the most likely state of the atmosphere, and also an uncertainty estimate of that analysis. The 20CR dataset provides the first estimates of global tropospheric variability, and of the dataset's time‐varying quality, from 1871 to the present at 6‐hourly temporal and 2° spatial resolutions. Intercomparisons with independent radiosonde data indicate that the reanalyses are generally of high quality. The quality in the extratropical Northern Hemisphere throughout the century is similar to that of current three‐day operational NWP forecasts. Intercomparisons over the second half‐century of these surface‐based reanalyses with other reanalyses that also make use of upper‐air and satellite data are equally encouraging. It is anticipated that the 20CR dataset will be a valuable resource to the climate research community for both model validations and diagnostic studies. Some surprising results are already evident. For instance, the long‐term trends of indices representing the North Atlantic Oscillation, the tropical Pacific Walker Circulation, and the Pacific–North American pattern are weak or non‐existent over the full period of record. The long‐term trends of zonally averaged precipitation minus evaporation also differ in character from those in climate model simulations of the twentieth century. Copyright © 2011 Royal Meteorological Society and Crown Copyright.

The nuclear ribosomal internal transcribed spacer (ITS) region is the formal fungal barcode and in most cases the marker of choice for the exploration of fungal diversity in environmental samples. Two problems are particularly acute in the pursuit of satisfactory taxonomic assignment of newly generated ITS sequences: (i) the lack of an inclusive, reliable public reference data set and (ii) the lack of means to refer to fungal species, for which no Latin name is available in a standardized stable way. Here, we report on progress in these regards through further development of the UNITE database (http://unite.ut.ee) for molecular identification of fungi. All fungal species represented by at least two ITS sequences in the international nucleotide sequence databases are now given a unique, stable name of the accession number type (e.g. Hymenoscyphus pseudoalbidus|GU586904|SH133781.05FU), and their taxonomic and ecological annotations were corrected as far as possible through a distributed, third-party annotation effort. We introduce the term 'species hypothesis' (SH) for the taxa discovered in clustering on different similarity thresholds (97-99%). An automatically or manually designated sequence is chosen to represent each such SH. These reference sequences are released (http://unite.ut.ee/repository.php) for use by the scientific community in, for example, local sequence similarity searches and in the QIIME pipeline. The system and the data will be updated automatically as the number of public fungal ITS sequences grows. We invite everybody in the position to improve the annotation or metadata associated with their particular fungal lineages of expertise to do so through the new Web-based sequence management system in UNITE.

Elucidating the biogeography of bacterial communities on the human body is critical for establishing healthy baselines from which to detect differences associated with diseases. To obtain an integrated view of the spatial and temporal distribution of the human microbiota, we surveyed bacteria from up to 27 sites in seven to nine healthy adults on four occasions. We found that community composition was determined primarily by body habitat. Within habitats, interpersonal variability was high, whereas individuals exhibited minimal temporal variability. Several skin locations harbored more diverse communities than the gut and mouth, and skin locations differed in their community assembly patterns. These results indicate that our microbiota, although personalized, varies systematically across body habitats and time; such trends may ultimately reveal how microbiome changes cause or prevent disease.

Abstract An overview of the Community Earth System Model Version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long preindustrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low‐top” (40 km, with limited chemistry) and “high‐top” (130 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low‐latitude precipitation and shortwave cloud forcing biases; better representation of the Madden‐Julian Oscillation; better El Niño‐Southern Oscillation‐related teleconnections; and a global land carbon accumulation trend that agrees well with observationally based estimates. Most tropospheric and surface features of the low‐ and high‐top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1–5.3 °C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2's transient climate response of 1.9–2.0 °C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and they represent CESM2's contributions to the Coupled Model Intercomparison Project Phase 6.

Exploring large environmental datasets generated by high-throughput DNA sequencing technologies requires new analytical approaches to move beyond the basic inventory descriptions of the composition and diversity of natural microbial communities. In order to investigate potential interactions between microbial taxa, network analysis of significant taxon co-occurrence patterns may help to decipher the structure of complex microbial communities across spatial or temporal gradients. Here, we calculated associations between microbial taxa and applied network analysis approaches to a 16S rRNA gene barcoded pyrosequencing dataset containing >160 000 bacterial and archaeal sequences from 151 soil samples from a broad range of ecosystem types. We described the topology of the resulting network and defined operational taxonomic unit categories based on abundance and occupancy (that is, habitat generalists and habitat specialists). Co-occurrence patterns were readily revealed, including general non-random association, common life history strategies at broad taxonomic levels and unexpected relationships between community members. Overall, we demonstrated the potential of exploring inter-taxa correlations to gain a more integrated understanding of microbial community structure and the ecological rules guiding community assembly.

Organic aerosol (OA) data acquired by the Aerosol Mass Spectrometer (AMS) in 37 field campaigns were deconvolved into hydrocarbon‐like OA (HOA) and several types of oxygenated OA (OOA) components. HOA has been linked to primary combustion emissions (mainly from fossil fuel) and other primary sources such as meat cooking. OOA is ubiquitous in various atmospheric environments, on average accounting for 64%, 83% and 95% of the total OA in urban, urban downwind, and rural/remote sites, respectively. A case study analysis of a rural site shows that the OOA concentration is much greater than the advected HOA, indicating that HOA oxidation is not an important source of OOA, and that OOA increases are mainly due to SOA. Most global models lack an explicit representation of SOA which may lead to significant biases in the magnitude, spatial and temporal distributions of OA, and in aerosol hygroscopic properties.

Our growing awareness of the microbial world's importance and diversity contrasts starkly with our limited understanding of its fundamental structure. Despite recent advances in DNA sequencing, a lack of standardized protocols and common analytical frameworks impedes comparisons among studies, hindering the development of global inferences about microbial life on Earth. Here we present a meta-analysis of microbial community samples collected by hundreds of researchers for the Earth Microbiome Project. Coordinated protocols and new analytical methods, particularly the use of exact sequences instead of clustered operational taxonomic units, enable bacterial and archaeal ribosomal RNA gene sequences to be followed across multiple studies and allow us to explore patterns of diversity at an unprecedented scale. The result is both a reference database giving global context to DNA sequence data and a framework for incorporating data from future studies, fostering increasingly complete characterization of Earth's microbial diversity.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

The development of a new high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) is reported. The high-resolution capabilities of this instrument allow the direct separation of most ions from inorganic and organic species at the same nominal m/z, the quantification of several types of organic fragments (CxHy, CxHyOz, CxHyNp, CxHyOzNp), and the direct identification of organic nitrogen and organosulfur content. This real-time instrument is field-deployable, and its high time resolution (0.5 Hz has been demonstrated) makes it well-suited for studies in which time resolution is critical, such as aircraft studies. The instrument has two ion optical modes: a single-reflection configuration offers higher sensitivity and lower resolving power (up to approximately 2100 at m/z 200), and a two-reflectron configuration yields higher resolving power (up to approximately 4300 at m/z 200) with lower sensitivity. The instrument also allows the determination of the size distributions of all ions. One-minute detection limits for submicrometer aerosol are <0.04 microg m(-3) for all species in the high-sensitivity mode and <0.4 microg m(-3) in the high-resolution mode. Examples of ambient aerosol data are presented from the SOAR-1 study in Riverside, CA, in which the spectra of ambient organic species are dominated by CxHy and CxHyOz fragments, and different organic and inorganic fragments at the same nominal m/z show different size distributions. Data are also presented from the MIRAGE C-130 aircraft study near Mexico City, showing high correlation with independent measurements of surrogate aerosol mass concentration.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate – the “global carbon budget” – is important tobetter understand the global carbon cycle, support the development ofclimate policies, and project future climate change. Here we describe andsynthesize data sets and methodology to quantify the five major componentsof the global carbon budget and their uncertainties. Fossil CO2emissions (EFOS) are based on energy statistics and cement productiondata, while emissions from land-use change (ELUC), mainlydeforestation, are based on land use and land-use change data andbookkeeping models. Atmospheric CO2 concentration is measured directlyand its growth rate (GATM) is computed from the annual changes inconcentration. The ocean CO2 sink (SOCEAN) and terrestrialCO2 sink (SLAND) are estimated with global process modelsconstrained by observations. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the lastdecade available (2010–2019), EFOS was 9.6 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.4 ± 0.5 GtC yr−1 when the cement carbonation sink is included), andELUC was 1.6 ± 0.7 GtC yr−1. For the same decade, GATM was 5.1 ± 0.02 GtC yr−1 (2.4 ± 0.01 ppm yr−1), SOCEAN 2.5 ± 0.6 GtC yr−1, and SLAND 3.4 ± 0.9 GtC yr−1, with a budgetimbalance BIM of −0.1 GtC yr−1 indicating a near balance betweenestimated sources and sinks over the last decade. For the year 2019 alone, thegrowth in EFOS was only about 0.1 % with fossil emissions increasingto 9.9 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.7 ± 0.5 GtC yr−1 when cement carbonation sink is included), and ELUC was 1.8 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5 ± 0.9 GtC yr−1 (42.2 ± 3.3 GtCO2). Also for 2019, GATM was5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEANwas 2.6 ± 0.6 GtC yr−1, and SLAND was 3.1 ± 1.2 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminarydata for 2020, accounting for the COVID-19-induced changes in emissions,suggest a decrease in EFOS relative to 2019 of about −7 % (medianestimate) based on individual estimates from four studies of −6 %, −7 %,−7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in thecomponents of the global carbon budget are consistently estimated over theperiod 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for therepresentation of semi-decadal variability in CO2 fluxes. Comparison ofestimates from diverse approaches and observations shows (1) no consensusin the mean and trend in land-use change emissions over the last decade, (2)a persistent low agreement between the different methods on the magnitude ofthe land CO2 flux in the northern extra-tropics, and (3) an apparentdiscrepancy between the different methods for the ocean sink outside thetropics, particularly in the Southern Ocean. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set (Friedlingstein et al.,2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014,2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2020 (Friedlingstein et al., 2020).

In 1996, version 2.5 of the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk algorithm was published, and it has become one of the most frequently used algorithms in the air–sea interaction community. This paper describes steps taken to improve the algorithm in several ways. The number of iterations to solve for stability has been shortened from 20 to 3, and adjustments have been made to the basic profile stability functions. The scalar transfer coefficients have been redefined in terms of the mixing ratio, which is the fundamentally conserved quantity, rather than the measured water vapor mass concentration. Both the velocity and scalar roughness lengths have been changed. For the velocity roughness, the original fixed value of the Charnock parameter has been replaced by one that increases with wind speeds of between 10 and 18 m s−1. The scalar roughness length parameterization has been simplified to fit both an early set of NOAA/Environmental Technology Laboratory (ETL) experiments and the Humidity Exchange Over the Sea (HEXOS) program. These changes slightly increase the fluxes for wind speeds exceeding 10 m s−1. For interested users, two simple parameterizations of the surface gravity wave influence on fluxes have been added (but not evaluated). This new version of the algorithm (COARE 3.0) was based on published results and 2777 1-h covariance flux measurements in the ETL inventory. To test it, 4439 new values from field experiments between 1997 and 1999 were added, which now dominate the database, especially in the wind speed regime beyond 10 m s−1, where the number of observations increased from 67 to about 800. After applying various quality controls, the database was used to evaluate the algorithm in several ways. For an overall mean, the algorithm agrees with the data to within a few percent for stress and latent heat flux. The agreement is also excellent when the bulk and directly measured fluxes are averaged in bins of 10-m neutral wind speed. For a more stringent test, the average 10-m neutral transfer coefficients were computed for stress and moisture in wind speed bins, using different averaging schemes with fairly similar results. The average (mean and median) model results agreed with the measurements to within about 5% for moisture from 0 to 20 m s−1. For stress, the covariance measurements were about 10% higher than the model at wind speeds over 15 m s−1, while inertial-dissipation measurements agreed closely at all wind speeds. The values for stress are between 8% (for inertial dissipation) and 18% (for covariance) higher at 20 m s−1 than two other classic results. Twenty years ago, bulk flux schemes were considered to be uncertain by about 30%; the authors find COARE 3.0 to be accurate within 5% for wind speeds of 0–10 m s−1 and 10% for wind speeds of between 10 and 20 m s−1.