skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: . Drivers of diffusive CH4 emissions from shallow subarctic lakes on daily to multi-year time scales
Lakes and reservoirs contribute to regional carbon budgets via significant emissions of climate forcing trace gases. Here, for improved modelling, we use 8 years of floating chamber measurements from three small, shallow subarctic lakes (2010–2017, n = 1306) to separate the contribution of physical and biogeochemical processes to the turbulence-driven, diffusion-limited flux of methane (CH4) on daily to multi-year timescales. Correlative data include 9 years of surface water concentration measurements (2009–2017, n = 606), total water column storage (2009–2017, n = 1593) and in situ meteorological observations. We used the latter to compute near surface turbulence based on similarity scaling and then applied the surface renewal model to compute gas transfer velocities. Chamber fluxes averaged 6.9 ± 0.3 mg CH4 m−2 d−1 and gas transfer velocities (k600) averaged 4.0 ± 0.1 cm h−1. Chamber derived gas transfer velocities tracked the power-law wind speed relation of the model. Coefficients for the model and dissipation rates depended on shear production of turbulence, atmospheric stability, and exposure to wind. Fluxes increased with wind speed until daily average values exceeded 6.5 m s−1, at which point emissions were suppressed due to rapid water column degassing reducing the water–air concentration gradient. Arrhenius-type temperature functions of the CH4 flux (Ea‘ = 0.90 ± 0.14 eV) were robust (R2 ≥ 0.93, p < 0.01) and also applied to the surface CH4 concentration (Ea‘ = 0.88 ± 0.09 eV). These results imply that emissions were strongly coupled to production and supply to the water column. Spectral analysis indicated that on timescales shorter than a month, emissions were driven by wind shear whereas on longer timescales variations in water temperature governed the flux. Long-term monitoring efforts are essential to identify distinct functional relations that govern flux variability on timescales of weather and climate change.  more » « less
Award ID(s):
1737411
PAR ID:
10211527
Author(s) / Creator(s):
Date Published:
Journal Name:
Biogeosciences
Volume:
17
ISSN:
1726-4170
Page Range / eLocation ID:
1911-1932
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Seasonally flooded forests along tropical rivers cover extensive areas, yet the processes driving air‐water exchanges of radiatively active gases are uncertain. To quantify the controls on gas transfer velocities, we combined measurements of water‐column temperature, meteorology in the forest and adjacent open water, turbulence with an acoustic Doppler velocimeter, gas concentrations, and fluxes with floating chambers. Under cooling, measured turbulence, quantified as the rate of dissipation of turbulent kinetic energy (ε), was similar to buoyancy flux computed from the surface energy budget, indicating convection dominated turbulence production. Under heating, turbulence was suppressed unless winds in the adjacent open water exceeded 1 m/s. Gas transfer velocities obtained from chamber measurements ranged from 1 to 5 cm/hr and were similar to or slightly less than predicted using a turbulence‐based surface renewal model computed with measuredεandεpredicted from wind and cooling. 
    more » « less
  2. null (Ed.)
    The hydrodynamics within small boreal lakes have rarely been studied, yet knowing whether turbulence at the air-water interface and in the water column scales with metrics developed elsewhere is essential for computing metabolism and fluxes of climate-forcing trace gases. We instrumented a humic, 4.7 ha, boreal lake with 2 meteorological stations, 3 thermistor arrays, an infra-red (IR) camera to quantify surface divergence, obtained turbulence as dissipation rate of turbulent kinetic energy (ε) using an acoustic Doppler velocimeter and a temperature-gradient microstructure profiler, and conducted chamber measurements for short periods to obtain fluxes and gas transfer velocities (k). Near-surface ε varied from 10-8 m2 s-3 to 10-6 m2 s-3 for the 0 to 4 m s-1 winds and followed predictions from Monin-Obukhov similarity theory. The coefficient of eddy diffusivity in the mixed layer was up to 10-3 m2 s-1 on the windiest afternoons, an order of magnitude less other afternoons, and near molecular at deeper depths. The upper thermocline upwelled when Lake numbers (LN) dropped below 4 facilitating vertical and horizontal exchange. k computed from a surface renewal model using ε agreed with values from chambers and surface divergence and increased linearly with wind speed. Diurnal thermoclines formed on sunny days when winds were < 3 m s-1, a condition that can lead to elevated near-surface ε and k. Results extend scaling approaches developed in the laboratory and for larger water bodies, illustrate turbulence and k are greater than expected in small wind-sheltered lakes, and provide new equations to quantify fluxes. 
    more » « less
  3. Autochamber-based CH4 fluxes and δ13C values measured with a Tunable Infrared Laser Direct Absorption Spectrometer (TILDAS, Aerodyne Research Inc.); and ancillary data, including CO2 fluxes (measured with a LGR Greenhouse Gas Analyzer), temperatures, atmospheric pressure, and photosynthetically active radiation (PAR). In addition to the data published here, data from 2011 is also available in the supplementary files to McCalley et al. (2014) under the Source data to Fig. 1 link.   METHODS: Methane fluxes were measured using a system of 8 automatic gas-sampling chambers made of transparent Lexan (n=3 each in the palsa and bog habitats, and n=2 in the fen habitat). Chambers were initially installed in the three habitat types at Stordalen Mire in 2001 (Bäckstrand et al., 2008) and the chamber lids were replaced in 2011 with the current design, similar to that described by Bubier et al 2003. Chambers cover an area of 0.2 m2 (45 cm x 45 cm), with a height ranging from 15-75 cm depending on habitat vegetation. At the Palsa and bog site the chamber base is flush with the ground and the chamber lid (15 cm in height) lifts clear of the base between closures. At the fen site the chamber base is raised 50–60 cm on lexon skirts to accommodate large stature vegetation. The chambers are instrumented with thermocouples measuring air and surface ground temperature, and water table depth and thaw depth are measured manually 3–5 times per week. The chambers are connected to the gas analysis system, located in an adjacent temperature-controlled cabin, by 3/8” Dekoron tubing through which air is circulated at approximately 2.5 L min-1. Each chamber lid is closed once every 3 hours for a period of 8 min, with a 5 min flush period before and after lid closure. We measured methane concentration using a Tunable Infrared Laser Direct Absorption Spectrometers (TILDAS, Aerodyne Research Inc.) connected to the main chamber circulation using ¼” Dekoron tubing (McCalley et al 2014). Calibrations were done every 90 min using 3 calibration gases spanning the observed concentration range (1.8–10 ppm). For each autochamber closure we calculated flux using a method consistent with that detailed by Bäckstrand et al 2008 for CO2 and total hydrocarbons, using a linear regression of changing headspace CH4 concentration over a period of 2.5 min. Eight 2.5 min regressions were calculated, staggered by 15 sec, and the most linear fit (highest r2) was then used to calculate flux. Daily average flux for each chamber was used to calculate daily flux and standard error for each cover type. References: Bäckstrand, K., Crill, P. M., Mastepanov, M., Christensen, T. R. & Bastviken, D. Total hydrocarbon flux dynamics at a subarctic mire in northern Sweden. Journal of Geophysical Research 113, (2008). Bubier, J. L., Crill, P. M., Mosedale, A., Frolking, S. & Linder, E. Peatland responses to varying interannual moisture conditions as measured by automatic CO2 chambers. Global Biogeochemical Cycles 17, (2003). McCalley, C.K., B.J. Woodcroft, S.B. Hodgkins, R.A. Wehr, E-H. Kim, R. Mondav, P.M. Crill, J.P. Chanton, V.I. Rich, G.W. Tyson, S.R. Saleska (2014), Methane dynamics regulated by microbial community response to permafrost thaw, Nature, 514:478-481, doi:10.1038/nature13798.   FILES: Files are named with the year or date range, followed by a suffix indicating data resolution: *_CH4output_clean_ckm.txt - Individual measurements of CH4 fluxes (CH4Flux), CO2 fluxes (CO2flux; for select years), and δ13C signature of emitted CH4 (Flux13CH4) for each chamber closure. CH4FluxRsq is the R2 value of the linear fit used to calculate CH4 flux, CO2Rsq is the R2 value of the linear fit used to calculate CO2 flux, and Flux13CH4_stdev is the standard deviation of the δ13C signature (standard deviation of the intercept of the Keeling plot). *_DailyCH4output_ckm.txt - Daily average CH4 fluxes (CH4Flux) and δ13C values (13CH4), grouped by site: Palsa, Bog, Fen, and Chamber 9 (bog/fen transition); along with standard deviations (stdev) and standard errors (se) of the flux or δ13C for each site type. For the Palsa, Bog, and Fen sites, these averages are calculated by chamber (n=3 for Palsa and Bog, n=2 for Fen), so each chamber's daily average is calculated, and then a daily average for that site is calculated as the average of the chambers. For Chamber 9 (bog/fen intermediate; n=1 chamber), averages are calculated by day as there are no chamber replicates. MEASUREMENT UNITS (same for both file types): CH4 flux: mg CH4 m−2 hr−1 CO2 flux: mg C m−2 h−1 δ13C: ‰ Temperature: °C Air pressure: mbar PAR: µmol photons m−2 s−1   FUNDING: This research is a contribution of the EMERGE Biology Integration Institute, funded by the National Science Foundation, Biology Integration Institutes Program, Award # 2022070.We thank the Swedish Polar Research Secretariat and SITES for the support of the work done at the Abisko Scientific Research Station. SITES is supported by the Swedish Research Council's grant 4.3-2021-00164.This study was also funded by the Genomic Science Program of the United States Department of Energy Office of Biological and Environmental Research, grant #s DE-SC0004632, DE-SC0010580, and DE-SC0016440.Autochamber measurements between 2013 and 2017 were supported by a grant from the US National Science Foundation MacroSystems program (NSF EF 1241037, PI Varner). 
    more » « less
  4. Abstract The wind‐blown flux of sand generates dunes, wind erosion, and mineral dust aerosols. Existing models predict sand flux using the wind friction velocity that characterizes near‐surface turbulent momentum fluxes. However, these models struggle to accurately predict sand fluxes. Here we analyze root causes of these model discrepancies using high‐frequency field measurements of winds and sand fluxes. We find that friction velocity is only predictive of sand fluxes on long timescales, when it correlates with horizontal wind speed. On shorter timescales, and for non‐ideal surface conditions, friction velocity is much less predictive, likely because the near‐surface wind momentum budget is dominated by other, less predictable terms. We furthermore find that variability in 30‐min averaged sand fluxes at a given friction velocity is not driven by changes in turbulence but by changes in surface conditions, raising a challenge for models. These findings can improve sand flux models and clarify their limitations. 
    more » « less
  5. Jan Vymazal (Ed.)
    Invasive species management typically aims to promote diversity and wildlife habitat, but little is known about how management techniques affect wetland carbon (C) dynamics. Since wetland C uptake is largely influenced by water levels and highly productive plants, the interplay of hydrologic extremes and invasive species is fundamental to understanding and managing these ecosystems. During a period of rapid water level rise in the Laurentian Great Lakes, we tested how mechanical treatment of invasive plant Typha × glauca shifts plant-mediated wetland C metrics. From 2015 to 2017, we implemented large-scale treatment plots (0.36-ha) of harvest (i.e., cut above water surface, removed biomass twice a season), crush (i.e., ran over biomass once mid-season with a tracked vehicle), and Typha-dominated controls. Treated Typha regrew with approximately half as much biomass as unmanipulated controls each year, and Typha production in control stands increased from 500 to 1500 g-dry mass m−2 yr−1 with rising water levels (~10 to 75 cm) across five years. Harvested stands had total in-situ methane (CH4) flux rates twice as high as in controls, and this increase was likely via transport through cut stems because crushing did not change total CH4 flux. In 2018, one year after final treatment implementation, crushed stands had greater surface water diffusive CH4 flux rates than controls (measured using dissolved gas in water), likely due to anaerobic decomposition of flattened biomass. Legacy effects of treatments were evident in 2019; floating Typha mats were present only in harvested and crushed stands, with higher frequency in deeper water and a positive correlation with surface water diffusive CH4 flux. Our study demonstrates that two mechanical treatments have differential effects on Typha structure and consequent wetland CH4 emissions, suggesting that C-based responses and multi-year monitoring in variable water conditions are necessary to accurately assess how management impacts ecological function. 
    more » « less