Low-salinity waterflooding is an emerging enhanced-oil-recovery (EOR) technique in which the salinity of the injected water is substantially reduced to improve oil recovery over conventional higher-salinity waterflooding. Although there are many low-salinity experimental results reported in the literature, publications on modeling this process are rare. Although there remains some debate regarding the mechanisms of low salinity waterflooding process (LoSal EOR V R ) * , the geochemical reactions that control the wetting of crude oil on the rock are likely to be central to a detailed description of the process. Because no comprehensive geochemical-based modeling has been applied in this area, it was decided to couple a state-of-the-art geochemical package, IPhreeqc (Charlton and Parkhurst 2011), developed by the US Geological Survey, with UTCOMP (Chang 1990), the compositional reservoir simulator developed by The University of Texas at Austin.A step-by-step algorithm is presented for integrating IPhreeqc with UTCOMP. Through this coupling, we are able to simulate homogeneous and heterogeneous (mineral dissolution/precipitation), irreversible, and ion-exchange reactions under nonisothermal, nonisobaric, and both local-equilibrium (away from the wellbore) and kinetic (near wellbore) conditions. Consistent with the literature, there are significant effects of water-soluble hydrocarbon components-e.g., carbon dioxide (CO 2 ), methane (CH 4 ), and acidic/basic components of the crude-on buffering the aqueous pH value and more generally, on the crude oil, brine, and rock reactions. Thermodynamic constraints are used to explicitly include the effect of these water-soluble hydrocarbon components. Hence, this combines the geochemical power of IPhreeqc with the important aspects of hydrocarbon flow and compositional effects to produce a robust, flexible, and accurate integrated tool capable of including the reactions needed to mechanistically model low-salinity waterflooding.Different geochemical-based approaches to modeling wettability change in sandstones (e.g., interpolation on the basis of total ionic strength and multicomponent ion exchange through surface complexation of the organometallic components) were implemented in UTCOMP-IPhreeqc, and the integrated tool is then used to match and interpret a low-salinity experiment published by Kozaki (2012) and the field trial performed by BP at the Endicott field.