Comment on "Scrutinizing the carbon cycle and CO2residence time in the atmosphere" by H. Harde

Harde (2017) proposes an alternative accounting scheme for the modern carbon cycle and concludes that only 4.3% of today's atmospheric CO2 is a result of anthropogenic emissions. As we will show, this alternative scheme is too simple, is based on invalid assumptions, and does not address many of the key processes involved in the global carbon cycle that are important on the timescale of interest. Harde (2017) therefore reaches an incorrect conclusion about the role of anthropogenic CO2 emissions. Harde (2017) tries to explain changes in atmospheric CO2 concentration with a single equation, while the most simple model of the carbon cycle must at minimum contain equations of at least two reservoirs (the atmosphere and the surface ocean), which are solved simultaneously. A single equation is fundamentally at odds with basic theory and observations. In the following we will (i) clarify the difference between CO2 atmospheric residence time and adjustment time, (ii) present recently published information about anthropogenic carbon, (iii) present details about the processes that are missing in Harde (2017), (iv) briefly discuss shortcoming in Harde's generalization to paleo timescales, (v) and comment on deficiencies in some of the literature cited in Harde (2017)

Investigation of a lattice Boltzmann model with a variable speed of sound

A lattice Boltzmann model is considered in which the speed of sound can be varied independently of the other parameters. The range over which the speed of sound can be varied is investigated and good agreement is found between simulations and theory. The onset of nonlinear effects due to variations in the speed of sound is also investigated and good agreement is again found with theory. It is also shown that the fluid viscosity is not altered by changing the speed of sound

Finite difference lattice Boltzmann model with flux limiters for liquid-vapor systems

In this paper we apply a finite difference lattice Boltzmann model to study the phase separation in a two-dimensional liquid-vapor system. Spurious numerical effects in macroscopic equations are discussed and an appropriate numerical scheme involving flux limiter techniques is proposed to minimize them and guarantee a better numerical stability at very low viscosity. The phase separation kinetics is investigated and we find evidence of two different growth regimes depending on the value of the fluid viscosity as well as on the liquid-vapor ratio.Comment: 10 pages, 10 figures, to be published in Phys. Rev.

Front pinning in capillary filling of chemically coated channels

The dynamics of capillary filling in the presence of chemically coated heterogeneous boundaries is investigated, both theoretically and numerically. In particular, by mapping the equations of front motion onto the dynamics of a dissipative driven oscillator, an analytical criterion for front pinning is derived, under the condition of diluteness of the coating spots. The criterion is tested against two dimensional Lattice Boltzmann simulations, and found to provide satisfactory agreement as long as the width of the front interface remains much thinner than the typical heterogeneity scale of the chemical coating.Comment: 7 pages, 4 figures, submitted to Physical Review

Mesoscopic modeling of a two-phase flow in the presence of boundaries: the Contact Angle

We present a mesoscopic model, based on the Boltzmann Equation, for the interaction between a solid wall and a non-ideal fluid. We present an analytic derivation of the contact angle in terms of the surface tension between the liquid-gas, the liquid-solid and the gas-solid phases. We study the dependency of the contact angle on the two free parameters of the model, which determine the interaction between the fluid and the boundaries, i.e. the equivalent of the wall density and of the wall-fluid potential in Molecular Dynamics studies. We compare the analytical results obtained in the hydrodynamical limit for the density profile and for the surface tension expression with the numerical simulations. We compare also our two-phase approach with some exact results for a pure hydrodynamical incompressible fluid based on Navier-Stokes equations with boundary conditions made up of alternating slip and no-slip strips. Finally, we show how to overcome some theoretical limitations connected with a discretized Boltzmann scheme and we discuss the equivalence between the surface tension defined in terms of the mechanical equilibrium and in terms of the Maxwell construction.Comment: 29 pages, 12 figure

Ballasting by cryogenic gypsum enhances carbon export in a Phaeocystis under-ice bloom

Mineral ballasting enhances carbon export from the surface to the deep ocean; however, little is known about the role of this process in the ice-covered Arctic Ocean. Here, we propose gypsum ballasting as a new mechanism that likely facilitated enhanced vertical carbon export from an under-ice phytoplankton bloom dominated by the haptophyte Phaeocystis. In the spring 2015 abundant gypsum crystals embedded in Phaeocystis aggregates were collected throughout the water column and on the sea floor at a depth below 2 km. Model predictions supported by isotopic signatures indicate that 2.7 g m-2 gypsum crystals were formed in sea ice at temperatures below -6.5 °C and released into the water column during sea ice melting. Our finding indicates that sea ice derived (cryogenic) gypsum is stable enough to survive export to the deep ocean and serves as an effective ballast mineral. Our findings also suggest a potentially important and previously unknown role of Phaeocystis in deep carbon export due to cryogenic gypsum ballasting. The rapidly changing Arctic sea ice regime might favour this gypsum gravity chute with potential consequences for carbon export and food partitioning between pelagic and benthic ecosystems.Publisher PDFPeer reviewe

Capillary filling with wall corrugations] Capillary filling in microchannels with wall corrugations: A comparative study of the Concus-Finn criterion by continuum, kinetic and atomistic approaches

We study the impact of wall corrugations in microchannels on the process of capillary filling by means of three broadly used methods - Computational Fluid Dynamics (CFD), Lattice-Boltzmann Equations (LBE) and Molecular Dynamics (MD). The numerical results of these approaches are compared and tested against the Concus-Finn (CF) criterion, which predicts pinning of the contact line at rectangular ridges perpendicular to flow for contact angles theta > 45. While for theta = 30, theta = 40 (no flow) and theta = 60 (flow) all methods are found to produce data consistent with the CF criterion, at theta = 50 the numerical experiments provide different results. Whilst pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle heigth.Comment: 25 pages, 8 figures, Langmuir in pres