Meta-analysis reveals that the effects of precipitation change on soil and litter fauna in forests depend on body size.

Anthropogenic climate change is altering precipitation regimes at a global scale. While precipitation changes have been linked to changes in the abundance and diversity of soil and litter invertebrate fauna in forests, general trends have remained elusive due to mixed results from primary studies. We used a meta-analysis based on 430 comparisons from 38 primary studies to address associated knowledge gaps, (i) quantifying impacts of precipitation change on forest soil and litter fauna abundance and diversity, (ii) exploring reasons for variation in impacts and (iii) examining biases affecting the realism and accuracy of experimental studies. Precipitation reductions led to a decrease of 39% in soil and litter fauna abundance, with a 35% increase in abundance under precipitation increases, while diversity impacts were smaller. A statistical model containing an interaction between body size and the magnitude of precipitation change showed that mesofauna (e.g. mites, collembola) responded most to changes in precipitation. Changes in taxonomic richness were related solely to the magnitude of precipitation change. Our results suggest that body size is related to the ability of a taxon to survive under drought conditions, or to benefit from high precipitation. We also found that most experiments manipulated precipitation in a way that aligns better with predicted extreme climatic events than with predicted average annual changes in precipitation and that the experimental plots used in experiments were likely too small to accurately capture changes for mobile taxa. The relationship between body size and response to precipitation found here has far-reaching implications for our ability to predict future responses of soil biodiversity to climate change and will help to produce more realistic mechanistic soil models which aim to simulate the responses of soils to global change.

Projected soil carbon loss with warming in constrained Earth system models.

The soil carbon-climate feedback is currently the least constrained component of global warming projections, and the major source of uncertainties stems from a poor understanding of soil carbon turnover processes. Here, we assemble data from long-term temperature-controlled soil incubation studies to show that the arctic and boreal region has the shortest intrinsic soil carbon turnover time while tropical forests have the longest one, and current Earth system models overestimate intrinsic turnover time by 30 percent across active, slow and passive carbon pools. Our constraint suggests that the global soils will switch from carbon sink to source, with a loss of 0.22-0.53 petagrams of carbon per year until the end of this century from strong mitigation to worst emission scenarios, suggesting that global soils will provide a strong positive carbon feedback on warming. Such a reversal of global soil carbon balance would lead to a reduction of 66% and 15% in the current estimated remaining carbon budget for limiting global warming well below 1.5 °C and 2 °C, respectively, rendering climate mitigation much more difficult.

Including land management in a European carbon model with lateral transfer to the oceans.

The use of cover crops (CCs) is a promising cropland management practice with multiple benefits, notably in reducing soil erosion and increasing soil organic carbon (SOC) storage. However, the current ability to represent these factors in land surface models remains limited to small scales or simplified and lumped approaches due to the lack of a sediment-carbon erosion displacement scheme. This precludes a thorough understanding of the consequences of introducing a CC into agricultural systems. In this work, this problem was addressed in two steps with the spatially distributed CE-DYNAM model. First, the historical effect of soil erosion, transport, and deposition on the soil carbon budget at a continental scale in Europe was characterized since the early industrial era, using reconstructed climate and land use forcings. Then, the impact of two distinct policy-oriented scenarios for the introduction of CCs were evaluated, covering the European cropping systems where surface erosion rates or nitrate susceptibility are critical. The evaluation focused on the increase in SOC storage and the export of particulate organic carbon (POC) to the oceans, compiling a continental-scale carbon budget. The results indicated that Europe exported 1.95 TgC/year of POC to the oceans in the last decade, and that CCs can contribute to reducing this amount while increasing SOC storage. Compared to the simulation without CCs, the additional rate of SOC storage induced by CCs peaked after 10 years of their adoption, followed by a decrease, and the cumulative POC export reduction stabilized after around 13 years. The findings indicate that the impacts of CCs on SOC and reduced POC export are persistent regardless of their spatial allocation adopted in the scenarios. Together, the results highlight the importance of taking the temporal aspect of CC adoption into account and indicate that CCs alone are not sufficient to meet the targets of the 4‰ initiative. Despite some known model limitations, which include the lack of feedback of erosion on the net primary productivity and the representation of carbon fluxes with an emulator, the current work constitutes the first approach to successfully couple a distributed routing scheme of eroded carbon to a land carbon model emulator at a reasonably high resolution and continental scale. SHORT ABSTRACT: A spatially distributed model coupling erosion, transport, and deposition to the carbon cycle was developed. Then, it was used to simulate the impact of cover crops on both erosion and carbon, to show that cover crops can simultaneously increase organic carbon storage and reduce particulate organic carbon export to the oceans. The results seemed persistent regardless of the spatial distribution of cover crops.

Mapping risks associated with soil copper contamination using availability and bio-availability proxies at the European scale.

Soil contamination by trace elements like copper (Cu) can affect soil functioning. Environmental policies with guidelines and soil survey measurements still refer to the total content of Cu in soils. However, Cu content in soil solution or free Cu content have been shown to be better proxies of risks of Cu mobility or (bio-)availability for soil organisms. Several empirical equations have been defined at the local scale to predict the amount of Cu in soil solution based on both total soil Cu content and main soil parameters involved in the soil/solution partitioning. Nevertheless, despite the relevance for risk assessment, these equations are not applied at a large spatial scale due to difficulties to perform changes from local to regional. To progress in this challenge, we collected several empirical equations from literature and selected those allowing estimation of the amount of Cu in solution, used as a proxy of available Cu, from the knowledge of both total soil Cu content and soil parameters. We did the same for the estimation of free Cu in solution, used as a proxy of bio-available Cu. These equations were used to provide European maps of (bio-)available Cu based on the one of total soil Cu over Europe. Results allowed comparing the maps of available and bio-available Cu at the European scale. This was done with respective median values of each form of Cu to identify specific areas of risks linked to these two proxies. Higher discrepancies were highlighted between the map of bio-available Cu and the map of soil total Cu compared to the Cu available map. Such results can be used to assess environmental-related issues for land use planning.

Lowering water table reduces carbon sink strength and carbon stocks in northern peatlands.

Peatlands at high latitudes have accumulated >400 Pg carbon (C) because saturated soil and cold temperatures suppress C decomposition. This substantial amount of C in Arctic and Boreal peatlands is potentially subject to increased decomposition if the water table (WT) decreases due to climate change, including permafrost thaw-related drying. Here, we optimize a version of the Organizing Carbon and Hydrology In Dynamic Ecosystems model (ORCHIDEE-PCH4) using site-specific observations to investigate changes in CO2 and CH4 fluxes as well as C stock responses to an experimentally manipulated decrease of WT at six northern peatlands. The unmanipulated control peatlands, with the WT <20 cm on average (seasonal max up to 45 cm) below the surface, currently act as C sinks in most years (58 ± 34 g C m-2  year-1 ; including 6 ± 7 g C-CH4 m-2  year-1 emission). We found, however, that lowering the WT by 10 cm reduced the CO2 sink by 13 ± 15 g C m-2  year-1 and decreased CH4 emission by 4 ± 4 g CH4 m-2  year-1 , thus accumulating less C over 100 years (0.2 ± 0.2 kg C m-2 ). Yet, the reduced emission of CH4 , which has a larger greenhouse warming potential, resulted in a net decrease in greenhouse gas balance by 310 ± 360 g CO2-eq  m-2  year-1 . Peatlands with the initial WT close to the soil surface were more vulnerable to C loss: Non-permafrost peatlands lost >2 kg C m-2 over 100 years when WT is lowered by 50 cm, while permafrost peatlands temporally switched from C sinks to sources. These results highlight that reductions in C storage capacity in response to drying of northern peatlands are offset in part by reduced CH4 emissions, thus slightly reducing the positive carbon climate feedbacks of peatlands under a warmer and drier future climate scenario.

Responses of soil nitrification activities to copper after a moisture stress.

Some steps of the soil nitrogen (N) cycle are sensitive to environmental pressures like soil moisture or contamination, which are expected to evolve during the next decades. Individual stresses have been well studied, but their combination is not yet documented. In this work, we aimed at assessing the importance of the soil moisture on the impact of copper (Cu) contaminations on the N cycling soil function using the potential nitrification activities (PNA) as bioindicator. A two-step experiment was performed. First, a loamy soil was incubated 5 weeks in either 30, 60, or 90% of its water holding capacity (WHC) or alternating drought and rewetting periods. Thereafter, soil samples were exposed to a gradient of Cu concentrations through a bioassay involving nitrification. The dose-response curves of PNA in function of added Cu were modeled to calculate the effective Cu concentrations, namely ECx with x being the percentage of PNA inhibition. These values were then compared between experimental conditions to highlight differences in threshold values. The preincubation moisture treatments significantly affected the PNA responses to the secondary Cu stress with, for instance, hormetic responses in all cases except for the dry-rewetting treatment. Small PNA inhibitions were estimated for high Cu doses in the soils with low water contents (30% WHC) or submitted to dry-rewetting cycles, contrarily to the patterns observed for the soils with high water contents (90% WHC) or submitted to a single period of drought. Overall, significant differences were found in estimated ECx values between moisture treatments.

Empirical estimates of regional carbon budgets imply reduced global soil heterotrophic respiration.

Resolving regional carbon budgets is critical for informing land-based mitigation policy. For nine regions covering nearly the whole globe, we collected inventory estimates of carbon-stock changes complemented by satellite estimates of biomass changes where inventory data are missing. The net land-atmospheric carbon exchange (NEE) was calculated by taking the sum of the carbon-stock change and lateral carbon fluxes from crop and wood trade, and riverine-carbon export to the ocean. Summing up NEE from all regions, we obtained a global 'bottom-up' NEE for net land anthropogenic CO2 uptake of -2.2 ± 0.6 PgC yr-1 consistent with the independent top-down NEE from the global atmospheric carbon budget during 2000-2009. This estimate is so far the most comprehensive global bottom-up carbon budget accounting, which set up an important milestone for global carbon-cycle studies. By decomposing NEE into component fluxes, we found that global soil heterotrophic respiration amounts to a source of CO2 of 39 PgC yr-1 with an interquartile of 33-46 PgC yr-1-a much smaller portion of net primary productivity than previously reported.

Large historical carbon emissions from cultivated northern peatlands.

When a peatland is drained and cultivated, it behaves as a notable source of CO2 However, we lack temporally and spatially explicit estimates of carbon losses from cultivated peatlands. Using a process-based land surface model that explicitly includes representation of peatland processes, we estimate that northern peatlands converted to croplands emitted 72 Pg C over 850-2010, with 45% of this source having occurred before 1750. This source surpassed the carbon accumulation by high-latitude undisturbed peatlands (36 to 47 Pg C). Carbon losses from the cultivation of northern peatlands are omitted in previous land-use emission assessments. Adding this ignored historical land-use emission implies an 18% larger terrestrial carbon storage since 1750 to close the historical global carbon budget. We also show that carbon emission per unit area decrease with time since drainage, suggesting that time since drainage should be accounted for in inventories to refine land-use emissions from cultivated peatlands.

A mixed-effect model approach for assessing land-based mitigation in integrated assessment models: A regional perspective.

Given the prospects of low short-term emissions reduction, carbon removals (CDRs) are expected to play an important role in achieving ambitious mitigation targets in future scenarios of integrated assessment models (IAMs), particularly Bioenergy with Carbon Capture and Storage (BECCS). In this paper, we explore the IAMC 1.5℃ database to depict the characteristics of the two main CDR options present in mitigation scenarios: BECCS and afforestation/reforestation. We apply a linear mixed-effect model to capture the specific regional and cross-IAM effects. Results reveal that the distribution of BECCS and afforestation deployment differs across IAMs and regions and, to a second extent, time. BECCS is preferred in the scenarios not for its ability to expand energy use but actually because it appears as an alternative to afforestation, which is associated with a decrease in energy use. However, the regional distribution of CDR deployment does not show a common pattern across scenarios and IAMs. Therefore, a more comprehensive investigation is needed before it can support policy proposals.

Feasibility of the 4 per 1000 aspirational target for soil carbon: A case study for France.

Increasing soil organic carbon (SOC) stocks is a promising way to mitigate the increase in atmospheric CO2 concentration. Based on a simple ratio between CO2 anthropogenic emissions and SOC stocks worldwide, it has been suggested that a 0.4% (4 per 1000) yearly increase in SOC stocks could compensate for current anthropogenic CO2 emissions. Here, we used a reverse RothC modelling approach to estimate the amount of C inputs to soils required to sustain current SOC stocks and to increase them by 4‰ per year over a period of 30 years. We assessed the feasibility of this aspirational target first by comparing the required C input with net primary productivity (NPP) flowing to the soil, and second by considering the SOC saturation concept. Calculations were performed for mainland France, at a 1 km grid cell resolution. Results showed that a 30%-40% increase in C inputs to soil would be needed to obtain a 4‰ increase per year over a 30-year period. 88.4% of cropland areas were considered unsaturated in terms of mineral-associated SOC, but characterized by a below target C balance, that is, less NPP available than required to reach the 4‰ aspirational target. Conversely, 90.4% of unimproved grasslands were characterized by an above target C balance, that is, enough NPP to reach the 4‰ objective, but 59.1% were also saturated. The situation of improved grasslands and forests was more evenly distributed among the four categories (saturated vs. unsaturated and above vs below target C balance). Future data from soil monitoring networks should enable to validate these results. Overall, our results suggest that, for mainland France, priorities should be (1) to increase NPP returns in cropland soils that are unsaturated and have a below target carbon balance and (2) to preserve SOC stocks in other land uses.

Climate warming from managed grasslands cancels the cooling effect of carbon sinks in sparsely grazed and natural grasslands.

Grasslands absorb and release carbon dioxide (CO2), emit methane (CH4) from grazing livestock, and emit nitrous oxide (N2O) from soils. Little is known about how the fluxes of these three greenhouse gases, from managed and natural grasslands worldwide, have contributed to past climate change, or the roles of managed pastures versus natural grasslands. Here, global trends and regional patterns of the full greenhouse gas balance of grasslands are estimated for the period 1750 to 2012. A new spatially explicit land surface model is applied, to separate the direct effects of human activities from land management and the indirect effects from climate change, increasing CO2 and regional changes in nitrogen deposition. Direct human management activities are simulated to have caused grasslands to switch from a sink to a source of greenhouse gas, because of increased livestock numbers and accelerated conversion of natural lands to pasture. However, climate change drivers contributed a net carbon sink in soil organic matter, mainly from the increased productivity of grasslands due to increased CO2 and nitrogen deposition. The net radiative forcing of all grasslands is currently close to neutral, but has been increasing since the 1960s. Here, we show that the net global climate warming caused by managed grassland cancels the net climate cooling from carbon sinks in sparsely grazed and natural grasslands. In the face of future climate change and increased demand for livestock products, these findings highlight the need to use sustainable management to preserve and enhance soil carbon storage in grasslands and to reduce greenhouse gas emissions from managed grasslands.

Can N2 O emissions offset the benefits from soil organic carbon storage?

To respect the Paris agreement targeting a limitation of global warming below 2°C by 2100, and possibly below 1.5°C, drastic reductions of greenhouse gas emissions are mandatory but not sufficient. Large-scale deployment of other climate mitigation strategies is also necessary. Among these, increasing soil organic carbon (SOC) stocks is an important lever because carbon in soils can be stored for long periods and land management options to achieve this already exist and have been widely tested. However, agricultural soils are also an important source of nitrous oxide (N2 O), a powerful greenhouse gas, and increasing SOC may influence N2 O emissions, likely causing an increase in many cases, thus tending to offset the climate change benefit from increased SOC storage. Here we review the main agricultural management options for increasing SOC stocks. We evaluate the amount of SOC that can be stored as well as resulting changes in N2 O emissions to better estimate the climate benefits of these management options. Based on quantitative data obtained from published meta-analyses and from our current level of understanding, we conclude that the climate mitigation induced by increased SOC storage is generally overestimated if associated N2 O emissions are not considered but, with the exception of reduced tillage, is never fully offset. Some options (e.g. biochar or non-pyrogenic C amendment application) may even decrease N2 O emissions.

Ensemble modelling, uncertainty and robust predictions of organic carbon in long-term bare-fallow soils.

Simulation models represent soil organic carbon (SOC) dynamics in global carbon (C) cycle scenarios to support climate-change studies. It is imperative to increase confidence in long-term predictions of SOC dynamics by reducing the uncertainty in model estimates. We evaluated SOC simulated from an ensemble of 26 process-based C models by comparing simulations to experimental data from seven long-term bare-fallow (vegetation-free) plots at six sites: Denmark (two sites), France, Russia, Sweden and the United Kingdom. The decay of SOC in these plots has been monitored for decades since the last inputs of plant material, providing the opportunity to test decomposition without the continuous input of new organic material. The models were run independently over multi-year simulation periods (from 28 to 80 years) in a blind test with no calibration (Bln) and with the following three calibration scenarios, each providing different levels of information and/or allowing different levels of model fitting: (a) calibrating decomposition parameters separately at each experimental site (Spe); (b) using a generic, knowledge-based, parameterization applicable in the Central European region (Gen); and (c) using a combination of both (a) and (b) strategies (Mix). We addressed uncertainties from different modelling approaches with or without spin-up initialization of SOC. Changes in the multi-model median (MMM) of SOC were used as descriptors of the ensemble performance. On average across sites, Gen proved adequate in describing changes in SOC, with MMM equal to average SOC (and standard deviation) of 39.2 (±15.5) Mg C/ha compared to the observed mean of 36.0 (±19.7) Mg C/ha (last observed year), indicating sufficiently reliable SOC estimates. Moving to Mix (37.5 ± 16.7 Mg C/ha) and Spe (36.8 ± 19.8 Mg C/ha) provided only marginal gains in accuracy, but modellers would need to apply more knowledge and a greater calibration effort than in Gen, thereby limiting the wider applicability of models.

Simulating Erosion-Induced Soil and Carbon Delivery From Uplands to Rivers in a Global Land Surface Model.

Global water erosion strongly affects the terrestrial carbon balance. However, this process is currently ignored by most global land surface models (LSMs) that are used to project the responses of terrestrial carbon storage to climate and land use changes. One of the main obstacles to implement erosion processes in LSMs is the high spatial resolution needed to accurately represent the effect of topography on soil erosion and sediment delivery to rivers. In this study, we present an upscaling scheme for including erosion-induced lateral soil organic carbon (SOC) movements into the ORCHIDEE LSM. This upscaling scheme integrates information from high-resolution (3″) topographic and soil erodibility data into a LSM forcing file at 0.5° spatial resolution. Evaluation of our model for the Rhine catchment indicates that it reproduces well the observed spatial and temporal (both seasonal and interannual) variations in river runoff and the sediment delivery from uplands to the river network. Although the average annual lateral SOC flux from uplands to the Rhine River network only amounts to 0.5% of the annual net primary production and 0.01% of the total SOC stock in the whole catchment, SOC loss caused by soil erosion over a long period (e.g., thousands of years) has the potential to cause a 12% reduction in the simulated equilibrium SOC stocks. Overall, this study presents a promising approach for including the erosion-induced lateral carbon flux from the land to aquatic systems into LSMs and highlights the important role of erosion processes in the terrestrial carbon balance.

Leaching of dissolved organic carbon from mineral soils plays a significant role in the terrestrial carbon balance.

The leaching of dissolved organic carbon (DOC) from soils to the river network is an overlooked component of the terrestrial soil C budget. Measurements of DOC concentrations in soil, runoff and drainage are scarce and their spatial distribution highly skewed towards industrialized countries. The contribution of terrestrial DOC leaching to the global-scale C balance of terrestrial ecosystems thus remains poorly constrained. Here, using a process based, integrative, modelling approach to upscale from existing observations, we estimate a global terrestrial DOC leaching flux of 0.28 ± 0.07 Gt C year-1 which is conservative, as it only includes the contribution of mineral soils. Our results suggest that globally about 15% of the terrestrial Net Ecosystem Productivity (NEP, calculated as the difference between Net Primary Production and soil respiration) is exported to aquatic systems as leached DOC. In the tropical rainforest, the leached fraction of terrestrial NEP even reaches 22%. Furthermore, we simulated spatial-temporal trends in DOC leaching from soil to the river networks from 1860 to 2010. We estimated a global increase in terrestrial DOC inputs to river network of 35 Tg C year-1 (14%) from 1860 to 2010. Despite their low global contribution to the DOC leaching flux, boreal regions have the highest relative increase (28%) while tropics have the lowest relative increase (9%) over the historical period (1860s compared to 2000s). The results from our observationally constrained model approach demonstrate that DOC leaching is a significant flux in the terrestrial C budget at regional and global scales.

Microbial dynamics and soil physicochemical properties explain large-scale variations in soil organic carbon.

First-order organic matter decomposition models are used within most Earth System Models (ESMs) to project future global carbon cycling; these models have been criticized for not accurately representing mechanisms of soil organic carbon (SOC) stabilization and SOC response to climate change. New soil biogeochemical models have been developed, but their evaluation is limited to observations from laboratory incubations or few field experiments. Given the global scope of ESMs, a comprehensive evaluation of such models is essential using in situ observations of a wide range of SOC stocks over large spatial scales before their introduction to ESMs. In this study, we collected a set of in situ observations of SOC, litterfall and soil properties from 206 sites covering different forest and soil types in Europe and China. These data were used to calibrate the model MIMICS (The MIcrobial-MIneral Carbon Stabilization model), which we compared to the widely used first-order model CENTURY. We show that, compared to CENTURY, MIMICS more accurately estimates forest SOC concentrations and the sensitivities of SOC to variation in soil temperature, clay content and litter input. The ratios of microbial biomass to total SOC predicted by MIMICS agree well with independent observations from globally distributed forest sites. By testing different hypotheses regarding (using alternative process representations) the physicochemical constraints on SOC deprotection and microbial turnover in MIMICS, the errors of simulated SOC concentrations across sites were further decreased. We show that MIMICS can resolve the dominant mechanisms of SOC decomposition and stabilization and that it can be a reliable tool for predictions of terrestrial SOC dynamics under future climate change. It also allows us to evaluate at large scale the rapidly evolving understanding of SOC formation and stabilization based on laboratory and limited filed observation.

Fine-Root Turnover, Litterfall, and Soil Microbial Community of Three Mixed Coniferous-Deciduous Forests Dominated by Korean Pine (Pinus koraiensis) Along a Latitudinal Gradient.

Carbon dynamics in forests and in particular in soils are of primary importance in the context of climate change. A better understanding of the drivers controlling carbon storage is needed to improve climate mitigation strategies. Carbon storage is the result of a balance between inputs and outputs. Carbon inputs in the soil come from plant residues and root exudates, which are further transformed by microorganisms and stored in the long term. Here, we measured litter and fine-root production in three mixed forests dominated by Pinus koraiensis along a latitudinal gradient and performed a litterbag experiment to better understand the driving factors of decomposition. We found that over the three sites litter production was controlled by climatic factors, soil properties, and forest stand characteristics, whereas decay rates were mainly controlled by microbial community structure and soil stoichiometry. For fine roots, production differed among sites, and higher production was consistently observed in the top soil layers compared to deep soil, although the root distribution along the soil profile differed among sites. Fine-root decay rates were mainly controlled by fine-root stoichiometric characteristics. This article emphasizes the complexity of fine-root dynamics even for a single species. Environmental drivers impact on both production and decay, and we suggest performing manipulative field experiments to better identify the mechanisms involved in soil carbon cycling.

Decadal biomass increment in early secondary succession woody ecosystems is increased by CO2 enrichment.

Increasing atmospheric CO2 stimulates photosynthesis which can increase net primary production (NPP), but at longer timescales may not necessarily increase plant biomass. Here we analyse the four decade-long CO2-enrichment experiments in woody ecosystems that measured total NPP and biomass. CO2 enrichment increased biomass increment by 1.05 ± 0.26 kg C m-2 over a full decade, a 29.1 ± 11.7% stimulation of biomass gain in these early-secondary-succession temperate ecosystems. This response is predictable by combining the CO2 response of NPP (0.16 ± 0.03 kg C m-2 y-1) and the CO2-independent, linear slope between biomass increment and cumulative NPP (0.55 ± 0.17). An ensemble of terrestrial ecosystem models fail to predict both terms correctly. Allocation to wood was a driver of across-site, and across-model, response variability and together with CO2-independence of biomass retention highlights the value of understanding drivers of wood allocation under ambient conditions to correctly interpret and predict CO2 responses.

Global regulation of methane emission from natural lakes.

Methane (CH4) emissions from lakes are the largest of the emissions from freshwater ecosystems. We compile open water CH4 emission estimates from individual lakes from all over the world and consider the three main emission pathways: diffusive; ebullitive; and storage. The relationships between emissions, environmental variables, lake characteristics and methodological approaches are investigated for the measurements from 297 lakes. We show that environmental factors, such as temperature and precipitation, act as important driving factors for CH4 emissions, with higher emissions occurring where air temperature and precipitation are high. The diffusive flux of CH4 was found to be positively related to dissolved organic carbon concentration. Diffusive flux is the most frequently estimated component of the total flux, while the other emission pathways are often neglected. Based on the cases where all three components of the total flux were measured (30 lakes), we estimate that measuring the diffusive emission only, and then assuming that the value obtained is a good surrogate for the total emission, would have led to a 277% underestimation of the real total flux. In addition we show that the estimation of fluxes is method-dependent with substantial differences revealed between the flux estimates obtained from different measurement techniques. Some of this uncertainty is due to technical constraints which should not be neglected, and lake CH4 flux measurement techniques require thorough re-evaluation.