Box modeling is a versatile tool to explore earth systems processes, ranging from transient changes in the marine carbonate system to the long-term evolution of biogeochemical cycles. The Earth Science Box Modeling Toolkit is a python based toolkit that allows for the rapid creation and deployment of box models. It abstracts typical modeling tasks, e.g., air-sea gas exchange, weathering, seafloor carbonate precipitation/dissolution, kinetic isotope fractionation, etc., to python classes. Class instances can then simply be combined to build a model. While there is no graphical interface, this approach significantly reduces coding complexity and model development time. Crucially, the model structure is independent of the numerical implementation. Instead the model is parsed to dynamically create the necessary equation systems that can be passed to ode solver libraries like ODEPACK. Separating model description from numerical implementation results in well-documented model code, and combines the computational efficiency of state-of-the-art numerical libraries with the ease of use of python. The efficiency of this process is demonstrated by a 12-box model with air-sea gas exchange, tracers for carbon isotopes, and water column carbonate chemistry that requires about 1 CPU second to calculate the model evolution over 30 million years.