ALBERT

Description: Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric lifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). For the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 Tg CH4 yr−1 (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 Tg CH4 yr−1 or ∼ 60 % is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 Tg CH4 yr−1 or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 Tg CH4 yr−1 larger than our estimate for the previous decade (2000–2009), and 24 Tg CH4 yr−1 larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30 % larger global emissions (737 Tg CH4 yr−1, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼ 65 % of the global budget,

Description: The hydroxyl radical (OH), which is the dominant sink of methane (CH4), plays a key role in closing the global methane budget. Current top-down estimates of the global and regional CH4 budget using 3D models usually apply prescribed OH fields and attribute model–observation mismatches almost exclusively to CH4 emissions, leaving the uncertainties due to prescribed OH fields less quantified. Here, using a variational Bayesian inversion framework and the 3D chemical transport model LMDz, combined with 10 different OH fields derived from chemistry–climate models (Chemistry–Climate Model Initiative, or CCMI, experiment), we evaluate the influence of OH burden, spatial distribution, and temporal variations on the global and regional CH4 budget. The global tropospheric mean CH4-reaction-weighted [OH] ([OH]GM-CH4) ranges 10.3–16.3×105 molec cm−3 across 10 OH fields during the early 2000s, resulting in inversion-based global CH4 emissions between 518 and 757  Tg yr−1. The uncertainties in CH4 inversions induced by the different OH fields are similar to the CH4 emission range estimated by previous bottom-up syntheses and larger than the range reported by the top-down studies. The uncertainties in emissions induced by OH are largest over South America, corresponding to large inter-model differences of [OH] in this region. From the early to the late 2000s, the optimized CH4 emissions increased by 22±6  Tg yr−1 (17–30  Tg yr−1), of which ∼25  % (on average) offsets the 0.7  % (on average) increase in OH burden. If the CCMI models represent the OH trend properly over the 2000s, our results show that a higher increasing trend of CH4 emissions is needed to match the CH4 observations compared to the CH4 emission trend derived using constant OH. This study strengthens the importance of reaching a better representation of OH burden and of OH spatial and temporal distributions to reduce the uncertainties in the global and regional CH4 budgets.

Description: The validation of convective processes in global climate models (GCMs) could benefit from the use of large datasets that provide long-term climatologies of the spatial statistics of convection. To that regard, echo top heights (ETHs), convective areas, and frequencies of mesoscale convective systems (MCSs) from 17 years of data from a C-band polarization (CPOL) radar are analyzed in varying phases of the Madden–Julian Oscillation (MJO) and northern Australian monsoon in order to provide ample validation statistics for GCM validation. The ETHs calculated using velocity texture and reflectivity provide similar results, showing that the ETHs are insensitive to various techniques that can be used. Retrieved ETHs are correlated with those from cloud top heights retrieved by Multifunctional Transport Satellites (MTSATs), showing that the ETHs capture the relative variability in cloud top heights over seasonal scales. Bimodal distributions of ETH, likely attributable to the cumulus congestus clouds and mature stages of convection, are more commonly observed when the active phase of the MJO is over Australia due to greater mid-level moisture during the active phase of the MJO. The presence of a convectively stable layer at around 5 km altitude over Darwin inhibiting convection past this level can explain the position of the modes at around 2–4 km and 7–9 km. Larger cells were observed during break conditions compared to monsoon conditions, but only during the inactive phase of the MJO. The spatial distributions show that Hector, a deep convective system that occurs almost daily during the wet season over the Tiwi Islands, and sea-breeze convergence lines are likely more common in break conditions. Oceanic MCSs are more common during the night over Darwin. Convective areas were generally smaller and MCSs more frequent during active monsoon conditions. In general, the MJO is a greater control on the ETHs in the deep convective mode observed over Darwin, with higher distributions of ETH when the MJO is active over Darwin.

Description: The COnvective Precipitation Experiment (COPE) was designed to investigate the origins of heavy convective precipitation over the southwestern UK, a region that experiences flash flooding due to heavy precipitation from slow-moving convective systems. In this study, the microphysical and dynamical characteristics of developing turrets during 4 days in July and August 2013 are analyzed. In situ cloud microphysical measurements from the University of Wyoming King Air and vertically pointing W-band radar measurements from Wyoming Cloud Radar are examined, together with data from the ground-based NXPol radar. The 4 days presented here cover a range of environmental conditions in terms of wind shear and instability, resulting in a similarly wide variability in observed ice crystal concentrations, both across days as well as between clouds on individual days. The highest concentration of ice was observed on the days in which there was an active warm-rain process supplying precipitation-sized liquid drops. The high ice concentrations observed (〉100 L−1) are consistent with the production of secondary ice particles through the Hallett–Mossop process. Turrets that ascended through remnant cloud layers above the 0 ∘C level had higher ice particle concentrations, suggesting that entrainment of ice particles from older clouds or previous thermals may have acted to aid in the production of secondary ice through the Hallett–Mossop process. Other mechanisms such as the shattering of frozen drops may be more important for producing ice in more isolated clouds.

Description: Laboratory and in-flight evaluations of uncertainties of measurements from a Cloud Droplet Probe (CDP) are presented. A description of a water-droplet-generating device, similar to those used in previous studies, is provided along with validation of droplet sizing and positioning. Seven experiments with droplet diameters of 9, 17, 24, 29, 34, 38, and 46 µm tested sizing and counting performance across a 10 µm resolution grid throughout the sample area of a CDP. Results indicate errors in sizing that depend on both droplet diameter and position within the sample area through which a droplet transited. The CDP undersized 9µm droplets by 1–4 µm. Droplets with diameters of 17 and 24 µm were sized to within 2 µm, which is the nominal CDP bin width for droplets of that size. The majority of droplets larger than 17 µm were oversized by 2–4 µm, while a small percentage were severely undersized, by as much as 30 µm. This combination led to an artificial broadening and skewing of the spectra such that mean diameters from a near-monodisperse distribution compared well (within a few percent), while the median diameters were oversized by 5–15 %. This has implications on how users should calibrate their probes. Errors in higher-order moments were generally less than 10 %. Comparisons of liquid water content (LWC) calculated from the CDP and that measured from a Nevzorov hot-wire probe were conducted for 17 917 1 Hz in-cloud points. Although some differences were noted based on volume-weighted mean diameter and total droplet concentration, the CDP-estimated LWC exceeded that measured by the Nevzorov by approximately 20 %, more than twice the expected difference based on results of the laboratory tests and considerations of Nevzorov collection efficiency.

Description: To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (KDP) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (ZDR) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived KDP and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that ∼ 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign.

Description: Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH4) budget over 2000–2012 (Saunois et al., 2016), we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000–2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000–2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008–2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16–32] Tg CH4 yr−1 higher methane emissions over the period 2008–2012 compared to 2002–2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002–2006 and 2008–2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13CH4. However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13CH4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the "global carbon budget" – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of our imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2007–2016), EFF was 9.4 ± 0.5 GtC yr−1, ELUC 1.3 ± 0.7 GtC yr−1, GATM 4.7 ± 0.1 GtC yr−1, SOCEAN 2.4 ± 0.5 GtC yr−1, and SLAND 3.0 ± 0.8 GtC yr−1, with a budget imbalance BIM of 0.6 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1. Also for 2016, ELUC was 1.3 ± 0.7 GtC yr−1, GATM was 6.1 ± 0.2 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1 and SLAND was 2.7 ± 1.0 GtC yr−1, with a small BIM of −0.3 GtC. GATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the higher fossil emissions and smaller SLAND for that year consistent with El Niño conditions. The global atmospheric CO2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data indicate a renewed growth in EFF of +2.0 % (range of 0.8 % to 3.0 %) based on national emissions projections for China, USA, and India, and projections of Gross Domestic Product corrected for recent changes in the carbon intensity of the economy for the rest of the world. For 2017, initial data indicate an increase in atmospheric CO2 concentration of around 5.3 GtC (2.5 ppm), attributed to a combination of increasing emissions and receding El Niño conditions. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set (Le Quéré et al., 2016; 2015b; 2015a; 2014; 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2017.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use change (ELUC), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2009–2018), EFF was 9.5±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.9±0.02 GtC yr−1 (2.3±0.01 ppm yr−1), SOCEAN 2.5±0.6 GtC yr−1, and SLAND 3.2±0.6 GtC yr−1, with a budget imbalance BIM of 0.4 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in EFF was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC yr−1, reaching 10 GtC yr−1 for the first time in history, ELUC was 1.5±0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5±0.9 GtC yr−1 (42.5±3.3 GtCO2). Also for 2018, GATM was 5.1±0.2 GtC yr−1 (2.4±0.1 ppm yr−1), SOCEAN was 2.6±0.6 GtC yr−1, and SLAND was 3.5±0.7 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in EFF of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in 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 of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated by this work are available at https://doi.org/10.18160/gcp-2019 (Friedlingstein et al., 2019).

Description: The modeling study presented here aims to estimate how uncertainties in global hydroxyl radical (OH) distributions, variability, and trends may contribute to resolving discrepancies between simulated and observed methane (CH4) changes since 2000. A multi-model ensemble of 14 OH fields was analyzed and aggregated into 64 scenarios to force the offline atmospheric chemistry transport model LMDz (Laboratoire de Meteorologie Dynamique) with a standard CH4 emission scenario over the period 2000–2016. The multi-model simulated global volume-weighted tropospheric mean OH concentration ([OH]) averaged over 2000–2010 ranges between 8.7×105 and 12.8×105 molec cm−3. The inter-model differences in tropospheric OH burden and vertical distributions are mainly determined by the differences in the nitrogen oxide (NO) distributions, while the spatial discrepancies between OH fields are mostly due to differences in natural emissions and volatile organic compound (VOC) chemistry. From 2000 to 2010, most simulated OH fields show an increase of 0.1–0.3×105 molec cm−3 in the tropospheric mean [OH], with year-to-year variations much smaller than during the historical period 1960–2000. Once ingested into the LMDz model, these OH changes translated into a 5 to 15 ppbv reduction in the CH4 mixing ratio in 2010, which represents 7 %–20 % of the model-simulated CH4 increase due to surface emissions. Between 2010 and 2016, the ensemble of simulations showed that OH changes could lead to a CH4 mixing ratio uncertainty of 〉±30 ppbv. Over the full 2000–2016 time period, using a common state-of-the-art but nonoptimized emission scenario, the impact of [OH] changes tested here can explain up to 54 % of the gap between model simulations and observations. This result emphasizes the importance of better representing OH abundance and variations in CH4 forward simulations and emission optimizations performed by atmospheric inversions.