New software enables global ocean biogeochemical modeling in Python

Science Highlights

5 September 2025

Have you ever wondered what life would be like if you could write and run complex biogeochemical models easily and conveniently in Python? Wonder no more. In a paper published in J. Adv. Model. Earth Syst., Samar Khatiwala (2025; see reference below) describes tmm4py, a new software to enable efficient, global scale biogeochemical modelling in Python.

tmm4py is based on the Transport Matrix Method (TMM), an efficient numerical scheme for “offline” simulation of tracers driven by circulations from state-of-the-art physical models and state estimates. tmm4py exposes this functionality in Python, providing the tools needed to implement complex models in pure Python using standard modules such as NumPy, and run them interactively on hardware ranging from laptops to supercomputers. No knowledge of parallel computing required! tmm4py even extends the interactivity to models written in Fortran, allowing the many existing models coupled to the TMM, e.g., MITgcm, to be used from the familiar comfort of Python. Whether you’re a seasoned modeler, just want to try out an idea, or illustrate a concept in your teaching, tmm4py is designed to make biogeochemical modeling more widely accessible.

Download the code from: https://github.com/samarkhatiwala/tmm

Click here to open the image in a new window for a larger view.
Figure: Schematic illustrating the structure of tmm4py and its relationship with the various libraries and components it is built on or interacts with. Outlined boxes represent user‐supplied code (such as the “Hello World” example of the ideal age tracer shown on the left). Other low-level libraries on which tmm4py depends, for example, BLAS and LAPACK for linear algebra, MPI for parallel communication, and CUDA for GPUs, are not shown.

Reference:

Khatiwala, S. (2025). tmm4py: Global Ocean Biogeochemical Modeling in Python With the Transport Matrix Method. Journal of Advances in Modeling Earth Systems, 17. doi: 10.1029/2025ms005028

Joint Science Highlight with US Ocean Carbon and Biogeochemistry Programme (OCB).

Latest highlights

Continuous record of the Antarctic Circumpolar Current latitude over the last glacial-interglacial cycles

The meridional positions of the oceanic fronts separating subtropical and Antarctic waters are key to constraining the mechanisms that drive the degassing of deeply-stored CO2 at the end of the glacial periods…

23 Million Years of iron sources reconstructed in the Northwest Pacific Ocean using iron isotopes in a ferromanganese crust

Using iron isotopes and elemental records, Chu and co-workers analysed an iron-manganese crust recovered from the Hongyan Seamount, and reconstructed the sources of dissolved Fe to the Northwest Pacific Ocean…

Variability of the past deep-water masses in the Arabian Sea over the last 41ka

Shukla and colleagues propose a high-resolution authigenic ƐNd record measured from a sediment core located in the eastern Arabian Sea…

Arctic outflows to Davis Strait and the Labrador Sea traced by radionuclide distributions

This study demonstrates that Arctic outflows through the Canadian Arctic Archipelago play a much more significant role than previously recognised.

Rechercher