Modeling coupled biogeochemical cycles