PyCO2SYS: marine carbonate system calculations in Python
PyCO2SYS: marine carbonate system calculations in Python
复制标题
PyCO2SYS:Python 中的海洋碳酸盐系统计算
DOI:
10.5281/zenodo.4133857
复制
发表时间:
2020
期刊:
影响因子:
--
通讯作者:
D. Wallace
中科院分区:
文献类型:
--
作者:
M. Humphreys;L. Gregor;D. Pierrot;S. Heuven;E. Lewis;D. Wallace
Abstract. Oceanic dissolved inorganic carbon (TC) is the largest pool of carbon that interacts considerably with the atmosphere on human timescales. Oceanic TC is increasing through uptake of anthropogenic carbon dioxide (CO2), and seawater pH is decreasing as a consequence. Both the exchange of CO2 between ocean and atmosphere and the pH response are governed by a set of parameters that interact through chemical equilibria, collectively known as the marine carbonate system. To investigate these processes, at least two of the marine carbonate system's parameters are typically measured – most commonly, two from TC, total alkalinity (AT), pH, and seawater CO2 fugacity (fCO2; or its partial pressure, pCO2, or its dry-air mole fraction, xCO2) – from which the remaining parameters can be calculated and the equilibrium state of seawater solved. Several software tools exist to carry out these calculations, but no fully functional and rigorously validated tool was previously available for Python, a popular scientific programming language. Here, we present PyCO2SYS, a Python package intended to fill this capability gap. We describe the elements of PyCO2SYS that have been inherited from the existing CO2SYS family of software and explain subsequent adjustments and improvements. For example, PyCO2SYS uses automatic differentiation to solve the marine carbonate system and calculate chemical buffer factors, ensuring that the effect of every solute and reaction is accurately included in all its results. We validate PyCO2SYS with internal consistency tests and comparisons against other software, showing that PyCO2SYS produces results that are either virtually identical or different for known reasons, with the differences negligible for all practical purposes. We discuss new insights that arose during the development process, for example that the marine carbonate system cannot be unambiguously solved from the total alkalinity and carbonate ion parameter pair. Finally, we consider potential future developments to PyCO2SYS and discuss the outlook for this and other software for solving the marine carbonate system. The code for PyCO2SYS is distributed via GitHub ( https://github.com/mvdh7/PyCO2SYS ) under the GNU General Public License v3, archived on Zenodo (Humphreys et al., 2021), and documented online ( https://PyCO2SYS.readthedocs.io ).
影响因子:
4.9
作者:
Abril, G.;Bouillon, S.;Borges, A. V.
通讯作者:
Borges, A. V.
影响因子:
3
作者:
Orr, James C.;Epitalon, Jean-Marie;Gattuso, Jean-Pierre
通讯作者:
Gattuso, Jean-Pierre