Estimating the Air-Water Gas Transfer Velocity during Low Wind Conditions
The abundances of atmospheric carbon dioxide, CO2, and methane, CH4, are increasing. These increases affect e.g., the global carbon cycle and the climate both regionally and globally. To better understand the present and future atmospheric CO2 and CH4 concentrations and their climate impact, the gas exchange between water and the atmosphere is important. This exchange can occur in two directions. Oceans take up approximately one third of the anthropogenic CO2 release (the ocean carbon sink). At the same time coastal waters and inland waters emit large amounts of CO2 and CH4, altogether corresponding to a similar amount as the ocean sink. The interfacial gas-flux for CO2 and CH4 is controlled by the water-side. The gas-flux, F_g, is for such gases typically estimated as F_g=k_g(C_wb-ϑC_as) where k_g is the gas transfer velocity, C_wb and C_as are the gas concentrations in the water bulk and in the air at the surface, and ϑ is the dimensionless Ostwald solubility coefficient. The subject of this thesis is to describe and estimate k_g for gases that have a water-side controlled gas-flux (e.g., CO2, and CH4). Besides being important for the geophysical sciences, k_g is also used to design and optimize many applications in e.g., chemical and environmental engineering. The transfer velocity is influenced by interfacial shear stress from wind, natural convection due to surface heat flux, microscale breaking waves at moderate wind speeds, breaking waves at high wind speeds, bubbles, surfactants, and rain. This thesis focuses on the low wind condition where the forcings due to shear stress, natural convection, and surfactants are important. The relative importance of buoyancy and shear forcing is characterized via a Richardson number Ri=Bν⁄(u_*^4 ). Here B, ν, and u_* are the buoyancy flux, kinematic viscosity, and friction velocity, respectively. The thesis summarizes three papers where k_g has been studied numerically with direct numerical simulations (DNS) and one paper where field observations have been used. The results from the field measurements show close relationships for the method using flux-chambers and the parameterization using the rate of turbulent kinetic energy dissipation, and the quantities surface rms velocity and the significant wave height. A parameterization of area-integrated values of k_g from wave measurements was proposed. The DNS comprise flow conditions ranging from convection-dominated to shear-dominated cases. The results are used to: (i) evaluate different parameterizations of the air-water gas-exchange, (ii) determine, for a given buoyancy flux, the wind speed at which gas transfer becomes primarily shear driven, (iii) find an expression for the gas-transfer velocity for flows driven by both convection and shear, and (iv) investigate the influence of surfactants on gas transfer velocity. Parameterizations using either the rate of turbulent kinetic energy dissipation or the horizontal surface flow-divergence show a larger disadvantageous dependence on the type of forcing than the parameterization using the surface-normal heat-flux. Two parametrizations using the wind-speed above the surface give reasonable estimates for the transfer-velocity, depending however on the surface heat-flux. The transition from convection- to shear-dominated gas-transfer-velocity is shown to be at Ri≈0.004. This means that buoyancy fluxes in natural conditions are not important for gas exchange at wind velocities U_10 above approximately 3 ms^(-1). Below this wind speed the buoyancy fluxes should be taken into account. The transfer velocity is shown to be well represented by two different approaches: (i) Additive forcing as k_(g,sum)=A_Shear u_* (Ri⁄Ri_c +1)^(1⁄4)Sc^(-n), where Ri_c=(A_Shear⁄A_Buoy)^4 is a critical Richardson number, and (ii) either buoyancy or shear-stress forcing that gives k_g=A_Buoy (Bν)^(1⁄4)Sc^(-n) for Ri>Ri_c and k_g=A_shear u_* Sc^(-n) for Ri<Ri_c. Here A_Buoy=0.4 and A_Shear=0.1 are constants, Sc=v⁄D is the Schmidt number, D is the gas diffusivity in water, and n is an exponent that depends on the water-surface characteristics.
Parts of work
Fredriksson, S.T., Handler, R.A., Nilsson, H., Zhang, Q., and Arneborg, L. (2016) An Evaluation of Gas Transfer Velocity Parameterizations During Natural Convection using DNS. Journal of Geophysical Research. ::doi::10.1002/2015JC011112Zhang, Q., Handler, R.A., and Fredriksson, S.T. (2013) Direct numerical simulation of turbulent free convection in the presence of a surfactant. International Journal of Heat and Mass Transfer. ::doi::10.1016/j.ijheatmasstransfer.2013.01.031Fredriksson, S.T., Handler, R.A., Nilsson, H., and Arneborg, L. (2016) Surface Shear Stress Dependence of Gas Transfer Velocity Parameterizations using DNS. Submitted 2016.Gålfalk, M., Bastviken, D., Fredriksson, S.T., and Arneborg, L. (2013) Determination of the piston velocity for water-air interfaces using flux chambers, acoustic Doppler velocimetry, and IR imaging of the water surface. Journal of Geophysical Research: Biogeosciences. ::doi::10.1002/jgrg.20064
Doctor of Philosophy
Göteborgs universitet. Naturvetenskapliga fakulteten
Department of Marine Sciences ; Institutionen för marina vetenskaper
Fredagen den 29 april 2016, kl. 10.00, Stora Hörsalen, Institutionen för marina vetenskaper, Guldhedsgatan 5C, Göteborg,
Date of defence
Fredriksson, Sam T.
Air.sea gas exchange
Direct numerical simulations
Gas transfer velocity