TORAX incorporates modular physics-based and machine-learning models, leveraging JAX for fast runtimes with JIT compilation and automatic differentiation, enabling gradient-based optimisation workflows, Jacobian-based PDE solvers, and coupling to machine-learning surrogate models of physics.