Reactive transport simulation on unstructured meshes can provide fundamental insight into the effect that geometric complexity of geologic structures has on fluid flow and development of reaction fronts. When applied to conditions ranging from ambient to hydrothermal and combined with compressible flow, accounting for geometric complexity provides an advantage for applications such as enhanced geothermal systems, carbon dioxide sequestration, hydrothermal ore formation, and radioactive waste disposal. We introduce CSMP-GEMS, a thermo-hydro and chemical multicomponent reactive transport code based on coupling of the Complex System Modeling Platform (CSMP) transport modeling framework with the GEMS3K chemical speciation solver. GEMS3K features a comprehensive suite of non-ideal activity and equation-of-state models of solution phases (aqueous electrolyte, gas and fluid mixtures, solid solutions). Current features include transient, compressible, single-phase advective and/or dispersive fluid flow, mass transport, heat transport in saturated porous media, and geochemical reactions in subsurface hydrothermal systems. We present two one-dimensional numerical experiments to compare CSMP-GEMS with the reactive transport codes OpenGeoSys-GEM and TOUGHREACT. Each experiment simulates calcite dissolution and dolomite precipitation during advection and hydrodynamic dispersion. One experiment corresponds to an existing isothermal benchmark; the second explores the applicability of the codes to non-isothermal problems. We also present a two-dimensional example that illustrates the application of CSMP-GEMS on unstructured meshes that can represent complex geologic relations. The results suggest that all three codes are well suited to predicting fluid circulation, heat transport, and mineral stability within hydrothermal systems relevant to enhanced geothermal systems and carbon dioxide sequestration in deep aquifers. Self-consistent accounting for kinetic processes is a major advantage of TOUGHREACT, but published applications are restricted to orthogonal meshes, potentially limiting the applicability of TOUGHREACT to geometrically less complex natural systems. OpenGeoSys-GEM can operate on unstructured meshes that may include multiple element types, facilitating the examination of non-orthogonal domains. However, due to its reliance on the groundwater equations, OpenGeoSys-GEM may be best suited for application to systems in which flow includes dispersion/diffusion and is not compressible. CSMP-GEMS does not currently calculate reaction kinetics, but may be useful for application to geometrically complex systems.