Search

Small lakes represent an important source of atmospheric CH4 from northern wetlands. However, spatiotemporal variations in flux magnitudes and the lack of knowledge about their main environmental controls contribute large uncertainty into the global CH4 budget. In this study, we measured methane fluxes from small lakes using chambers and bubble traps. Field investigations were carried out in July–August 2014 within the West Siberian middle and southern taiga zones. The average and median of measured methane chamber fluxes were 0.32 and 0.30 mgCH4 m−2 h−1 for middle taiga lakes and 8.6 and 4.1 mgCH4 m−2 h−1 for southern taiga lakes, respectively. Pronounced flux variability was found during measurements on individual lakes, between individual lakes and between zones. To analyze these differences and the influences of environmental controls, we developed a new dynamic process-based model. It shows good performance with emission rates from the southern taiga lakes and poor performance for individual lakes in the middle taiga region. The model shows that, in addition to well-known controls such as temperature, pH and lake depth, there are significant variations in the maximal methane production potential between these climatic zones. In addition, the model shows that variations in gas-filled pore space in lake sediments are capable of controlling the total methane emissions from individual lakes. The CH4 emissions exhibited distinct zonal differences not only in absolute values but also in their probability density functions: the middle taiga lake fluxes were best described by a lognormal distribution while the southern taiga lakes followed a power-law distribution. The latter suggests applicability of self-organized criticality theory for methane emissions from the southern taiga zone, which could help to explain the strong variability within individual lakes.

This study investigates the capacity of a prognostic biosphere model to simulate global variability in atmospheric CO2 concentrations and vegetation carbon dynamics under current environmental conditions. Global data sets of atmospheric CO2 concentrations, above-ground biomass (AGB), and net primary productivity (NPP) in terrestrial vegetation were assimilated into the biosphere model using an inverse modeling method combined with an atmospheric transport model. In this process, the optimal physiological parameters of the biosphere model were estimated by minimizing the misfit between observed and modeled values, and parameters were generated to characterize various biome types. Results obtained using the model with the optimized parameters correspond to the observed seasonal variations in CO2 concentration and their annual amplitudes in both the Northern and Southern Hemispheres. In simulating the mean annual AGB and NPP, the model shows improvements in estimating the mean magnitudes and probability distributions for each biome, as compared with results obtained using prior simulation parameters. However, the model is less efficient in its simulation of AGB for forest type biomes. This misfit suggests that more accurate values of input parameters, specifically, grid mean AGB values and seasonal variabilities in physiological parameters, are required to improve the performance of the simulation model.

We have developed an improved version of the National Institute for Environmental Studies (NIES) three-dimensional chemical transport model (TM) designed for accurate tracer transport simulations in the stratosphere, using a hybrid sigma-isentropic (σ-θ) vertical coordinate that employs both terrain-following and isentropic parts switched smoothly around the tropopause. The air-ascending rate was derived from the effective heating rate and was used to simulate vertical motion in the isentropic part of the grid (above level 350 K), which was adjusted to fit to the observed age of the air in the stratosphere. Multi-annual simulations were conducted using the NIES TM to evaluate vertical profiles and dry-air column-averaged mole fractions of CO2 and CH4. Comparisons with balloon-borne observations over Sanriku (Japan) in 2000–2007 revealed that the tracer transport simulations in the upper troposphere and lower stratosphere are performed with accuracies of ~5% for CH4 and SF6, and ~1% for CO2 compared with the observed volume-mixing ratios. The simulated column-averaged dry air mole fractions of atmospheric carbon dioxide (XCO2) and methane (XCH4) were evaluated against daily ground-based high-resolution Fourier Transform Spectrometer (FTS) observations measured at twelve sites of the Total Carbon Column Observing Network (TCCON) (Bialystok, Bremen, Darwin, Garmisch, Izaña, Lamont, Lauder, Orleans, Park Falls, Sodankylä, Tsukuba, and Wollongong) between January 2009 and January 2011. The comparison shows the model's ability to reproduce the site-dependent seasonal cycles as observed by TCCON, with correlation coefficients typically on the order 0.8–0.9 and 0.4–0.8 for XCO2 and XCH4, respectively, and mean model biases of ±0.2% and ±0.5%, excluding Sodankylä, where strong biases are found. The ability of the model to capture the tracer total column mole fractions is strongly dependent on the model's ability to reproduce seasonal variations in tracer concentrations in the planetary boundary layer (PBL). We found a marked difference in the model's ability to reproduce near-surface concentrations at sites located some distance from multiple emission sources and where high emissions play a notable role in the tracer's budget. Comparisons with aircraft observations over Surgut (West Siberia), in an area with high emissions of methane from wetlands, show contrasting model performance in the PBL and in the free troposphere. Thus, the PBL is another critical region for simulating the tracer total column mole fractions.

We designed a method to simulate atmospheric CO2 concentrations at several continuous observation sites around the globe using surface fluxes at a very high spatial resolution. The simulations presented in this study were performed using the Global Eulerian-Lagrangian Coupled Atmospheric model (GELCA), comprising a Lagrangian particle dispersion model coupled to a global atmospheric tracer transport model with prescribed global surface CO2 flux maps at a 1 × 1 km resolution. The surface fluxes used in the simulations were prepared by assembling the individual components of terrestrial, oceanic and fossil fuel CO2 fluxes. This experimental setup (i.e. a transport model running at a medium resolution, coupled to a high-resolution Lagrangian particle dispersion model together with global surface fluxes at a very high resolution), which was designed to represent high-frequency variations in atmospheric CO2 concentration, has not been reported at a global scale previously. Two sensitivity experiments were performed: (a) using the global transport model without coupling to the Lagrangian dispersion model, and (b) using the coupled model with a reduced resolution of surface fluxes, in order to evaluate the performance of Eulerian-Lagrangian coupling and the role of high-resolution fluxes in simulating high-frequency variations in atmospheric CO2 concentrations. A correlation analysis between observed and simulated atmospheric CO2 concentrations at selected locations revealed that the inclusion of both Eulerian-Lagrangian coupling and high-resolution fluxes improves the high-frequency simulations of the model. The results highlight the potential of a coupled Eulerian-Lagrangian model in simulating high-frequency atmospheric CO2 concentrations at many locations worldwide. The model performs well in representing observations of atmospheric CO2 concentrations at high spatial and temporal resolutions, especially for coastal sites and sites located close to sources of large anthropogenic emissions. While this study focused on simulations of CO2 concentrations, the model could be used for other atmospheric compounds with known estimated emissions.

A chemistry-transport model (CTM) intercomparison experiment (TransCom-CH4) has been designed to investigate the roles of surface emissions, transport and chemical loss in simulating the global methane distribution. Model simulations were conducted using twelve models and four model variants and results were archived for the period of 1990–2007. All but one model transports were driven by reanalysis products from 3 different meteorological agencies. The transport and removal of CH4 in six different emission scenarios were simulated, with net global emissions of 513 ± 9 and 514 ± 14 Tg CH4 yr−1 for the 1990s and 2000s, respectively. Additionally, sulfur hexafluoride (SF6) was simulated to check the interhemispheric transport, radon (222Rn) to check the subgrid scale transport, and methyl chloroform (CH3CCl3) to check the chemical removal by the tropospheric hydroxyl radical (OH). The results are compared to monthly or annual mean time series of CH4, SF6 and CH3CCl3 measurements from 8 selected background sites, and to satellite observations of CH4 in the upper troposphere and stratosphere. Most models adequately capture the vertical gradients in the stratosphere, the average long-term trends, seasonal cycles, interannual variations (IAVs) and interhemispheric (IH) gradients at the surface sites for SF6, CH3CCl3 and CH4. The vertical gradients of all tracers between the surface and the upper troposphere are consistent within the models, revealing vertical transport differences between models. An average IH exchange time of 1.39 ± 0.18 yr is derived from SF6 time series. Sensitivity simulations suggest that the estimated trends in exchange time, over the period of 1996–2007, are caused by a change of SF6 emissions towards the tropics. Using six sets of emission scenarios, we show that the decadal average CH4 growth rate likely reached equilibrium in the early 2000s due to the flattening of anthropogenic emission growth since the late 1990s. Up to 60% of the IAVs in the observed CH4 concentrations can be explained by accounting for the IAVs in emissions, from biomass burning and wetlands, as well as meteorology in the forward models. The modeled CH4 budget is shown to depend strongly on the troposphere-stratosphere exchange rate and thus on the model's vertical grid structure and circulation in the lower stratosphere. The 15-model median CH4 and CH3CCl3 atmospheric lifetimes are estimated to be 9.99 ± 0.08 and 4.61 ± 0.13 yr, respectively, with little IAV due to transport and temperature.

Fourteen global atmospheric transport models were evaluated by comparing the simulation of 222Rn against measurements at three continental stations in Germany: Heidelberg, Freiburg and Schauinsland. Hourly concentrations simulated by the models using a common 222Rn-flux without temporal variations were investigated for 2002 and 2003. We found that the mean simulated concentrations in Heidelberg are related to the diurnal amplitude of boundary layer height in each model. Summer mean concentrations simulated by individual models were negatively correlated with the seasonal mean of diurnal amplitude of boundary layer height, while in winter the correlation was positive. We also found that the correlations between simulated and measured concentrations at Schauinsland were higher when the simulated concentrations were interpolated to the station altitude in most models. Temporal variations of the mismatch between simulated and measured concentrations suggest that there are significant interannual variations in the 222Rn exhalation rate in this region. We found that the local inversion layer during daytime in summer in Freiburg has a significant effect on 222Rn concentrations. We recommend Freiburg concentrations for validation of models that resolve local stable layers and those at Heidelberg for models without this capability.

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