Freeman J. Cook

A complete model using analytical solutions for one‐dimensional oxygen transport from the atmosphere into soil with microbial and root sinks that builds on work over 30 years is developed. This new model uses concepts from a previously published model for one distributed sink and two sinks with a distributed (microbial) and line (root) sink. It removes the problem, in previous publications, of matching the flux at the joining point between the two sink solution where the root sink ceases and the single sink at a finite depth. Analytical solutions are developed for integer values of p = Z _r / Z _m , where Z _m is the scaling depth for microbial respiration and Z _r is the scaling depth for root length density. The solutions allow two critical diffusivity ( D ) values to be defined ( D _c ) and ( D _{c 2} ). When D _{c 2} ≤ D < D _c , a procedure is presented to calculate the depth, z ₁ , where C ^′ = 0 and this is the depth where root uptake of oxygen ceases and is shown to be related to D / D _c . When D < D _{c 2} , a procedure is presented to estimate the depth, Z ₀ , at which the oxygen concentration = 0 and is shown to be related to D/D _{c 2} . These results have useful applications in determining soil aeration, soil biogeochemical reactions, soil surface flux of oxygen and carbon dioxide, and the effect of climate change on these processes through the temperature dependence of the solution. These results suggest the oxygen diffusion rate (ODR) is likely to be the best estimator of soil aeration but there will not be a universal value for all plants. The surface flux density of oxygen into the soil for both the microbial sink ( S _m ) and total sink ( f ₀ ) are presented and the ratio is shown to be related to ( is the diffusivity in air). The possible range in S _m / f ₀ is shown to be compatible with measured values from the literature. The solutions have many applications in environmental science.

This manuscript is mainly well written. My only issue with it is the explanation of the time component and hence the mass of amendments applied. The manuscript states that amendments were applied at 0, 60, 120, and 180 days, but does not stipulate wh...

Key Points: Techniques for estimating evaporation were compared to measured values Penman‐Monteith model performed well with local site data and wind function Model predictions are sensitive to local variation in meteorological data