Data acquisition. We launched a detailed review of the literature published in peer-reviewed journals through the year 1978-2020 (cut-off date on October 20, 2020). We extracted original experimental data from publications on aquatic CH4 and N2O fluxes as well as related parameters from six non-marine waters including rivers, reservoirs, lakes, ponds streams and estuaries. A combination of search terms [“CH4” OR “methane” AND “N2O” OR “nitrous oxide” AND “flux” OR “emission” OR “release” OR “evasion” AND “river” OR “reservoir” OR “lake” OR “pond” OR “stream” OR “estuary”]. All selected data were extracted from the Web of Science and Google Scholar, and also the publication sources by gathering and reevaluating the older literature cited in the earlier reviews. The experimental locations were mostly clustered in Asia, Europe, and North America, with only several studies scattered in South and Northern Hemispheres with high latitudes (Supplementary Fig. S1). Overall, we established a solid dataset consisting of 2997 in-situ flux or concentration measurements of CH4 and N2O sourced from 277 peer-reviewed publications, of which 52 studies with 196 simultaneous flux measurements of CH4 and N2O were included (Supplementary Fig. S1).
We only included in-situ measurements of CH4 and N2O fluxes or concentrations from non-marine waters. With the focus on natural aquatic systems, the gas flux data from artificial constructed ponds used for sewage treatment or agricultural aquaculture ponds were excluded from our dataset due to intensive human impacts. We only incorporated studies in which aquatic systems were clearly defined. For studies in which the type of riverine systems was not clearly defined, we grouped them into rivers or streams according to the specific Strahler stream order44. We ensured that the data on gas fluxes and geographical information were reported or can be made available from authors. In case of absence of data on climate (e.g., temperature, precipitation) information, we obtained relevant data from the World Meteorological Organization (http://www.worldweather.cn/zh/home.html). When the same site was reported in multiple studies, we used the study that included the largest number of sampling dates, either across seasons or years. If the data were collected across multiple years, we calculated the average fluxes or concentrations over the whole measurement period. Generally, diffusive and ebullitive CH4 fluxes spanned over four orders of magnitude, ranging from 0.00 to 56.00 mg m-2 h-1 and 0.00 to 60.42 mg m-2 h-1 in inland waters, and varied from -0.15 to 17.78 mg m-2 h-1 and 0.01 to 0.18 mg m-2 h-1 in estuaries (Supplementary Table S4). Diffusive N2O fluxes and concentrations had a range from -79.00 to 1151.77 μg m-2 h-1 and 0.50 to 1500.00 nmol L-1 in inland waters, and from -11.90 to 322.67 μg m-2 h-1 and 4.35 to 210.30 nmol L-1 in estuarine open waters (Supplementary Table S5).
Upscaling and uncertainties. We used a bottom-up approach to upscale CH4 and N2O emissions from non-marine waters by multiplying averaged emission rates by the estimated global areal extents of water bodies12,18. In terms of CH4 emissions, only studies with simultaneous measurement data for both diffusive and ebullitive fluxes were included in our estimates to reduce bias. For N2O emissions, considering the ebullition was not the representative pathway of N2O release from water bodies, we chose to only estimate diffusive N2O emissions from non-marine waters, finally leading to exclusion of several sporadic ebullitive N2O fluxes from our dataset. Moreover, in order to reduce uncertainties, we only estimated diffusive CH4 emissions from estuaries due to insufficient data on ebullitive CH4 fluxes. To generate annual mean gas fluxes, we assumed the average seasonal fluxes were representative of the entire year in tropical and frigid regions. In other regions with typical seasonal differences, the seasonal flux data (collected in summer or winter) were rectified by the annual mean temperature using a temperature-dependent empirical model reviewed by Yan et al.45. However, we did not extrapolate CH4 and N2O emissions from ponds in this study due to large uncertainties in current available areal extents determination. All the data on global surface areas of aquatic ecosystems were cited from recently published literature with solid updated estimates18,46,47.
Compared to previous global estimates based on limited and localized CH4 and N2O flux rates, we have dedicated to exploring the fractions of global total freshwater CH4 emissions resulting from two different emission pathways (diffusion and ebullition) by incorporating simultaneous flux measurement data. Meanwhile, we provided a full understanding of the magnitude and drivers of diffusive N2O emissions from diverse non-marine waters, relative to previously mostly limited in a single aquatic system5,18. However, uncertainties remained existed for our estimates. First, while a wide range in diffusive fluxes has been reported for estuarine open waters, measurements of ebullition remain notably scare, especially for the simultaneous measurement data with diffusion (Supplementary Table S4). Second, except in tropical and frigid regions, flux data showed considerable variations with seasons, with general higher flux rates occurring in summer than in other seasons, although we have attempted to account for this in our analysis by rectifying the flux data using an earlier established temperature-dependent empirical model45. Third, we did not estimate the indirect EF5 of N2O based on the IPCC methodology to create a comparison in this study due to lack of detailed information on N inputs in most studies. Thus, given that future changes in climate and anthropogenic N loading are expected to increase N2O emissions from non-marine waters, more extensive direct measurements of N2O fluxes coupled with aquatic N loading rates are highly needed to make the IPCC methodology applicable to bridge the gap between global bottom-up and top-down inventories.
Calculation of indirect N2O emission factors (EF5). Indirect N2O emission factors for riverine systems (EF5r) and estuaries (EF5e) were estimated in this study to create a comparison with the recently updated IPCC default value of 0.26% (ref.34). The indirect EF5 of N2O represents N2O emissions from a given water body to the atmosphere as a fraction of N loading into the system6. IPCC defined the indirect EF5 of N2O as a ratio of N2O-N emitted from leached N and N in runoff divided by the fraction of all N added to, or mineralized within managed soils that is lost through leaching and runoff 48. Due to incomplete acquisition of the specific information (e.g., data on N leaching and runoff) required to determine the indirect EF5 based on the IPCC methodology for all aquatic systems, we therefore alternatively used the concentration method, i.e., the N2O-N/NO3--N mass ratio derived from the concentration data of N2O and nitrate (NO3-) reviewed from water bodies to estimate the indirect EF5 of N2O (ref.6) using the following equation:
Where EF5 is the indirect emission factor determined by the N2O-N/NO3--N mass ratio method, CN2O-N and NO3--N are concentrations measured at the water-air interface and dissolved in surface water of aquatic ecosystems, respectively35,44.
Estimation of CO2-equivalent emissions. Total CO2-equivalent emissions or emission intensity of CH4 and N2O from aquatic ecosystems were estimated using the conversion factors (mass basis) of 28 for CH4 and 265 for N2O over the time horizon of 100 years49.
Statistical analyses. One-way analysis of variance (ANOVA) was performed to test the difference in CH4 and N2O fluxes between two CH4 emission pathways (diffusion and ebullition), two flux-derived methods (chamber-based and model-based), and among different aquatic ecosystems. Linear or nonlinear regressions were used to examine the dependence of CH4 and N2O fluxes on potential driving factors. Linear stepwise regression models with the personality of Ordinary Least Squares (OSL) were conducted to identify the appropriate subset of environmental parameters that can best predict CH4 and N2O fluxes from inland waters or individual aquatic systems. All statistical analyses were carried out using JMP version 7.0 and R, and statistical significance was determined at the 0.05 probability level.