Water-sediment interactions triggered by the injection of oxidized aqueous solutions into anoxic groundwater systems usually modify both the aquifer matrix and control the final aqueous composition. The identification and quantification of these reactions in complex heterogeneous systems remains a challenge for the analysis and prediction of water quality changes. Driven by the proposed injection of large quantities of oxic water into a deep anoxic heterogeneous pyritic aquifer; this study was undertaken to quantify the reactivity of aquifer sediments with respect to oxidant consumption and to characterize the variability of the reaction rates across different lithological units. A total of 53 samples were incubated for periods of 14, 37, and 50 days, during which the gas-phase was continuously monitored and the aqueous composition analyzed. A geochemical modeling framework was developed that incorporated a mixed set of equilibrium and kinetic reactions and supported the interpretation and quantification of the geochemical controls. The good agreement between simulated and experimental results of O2 consumption, CO2 production, pH, major ions, and trace metals suggests that the framework was able to successfully quantify reaction rates of competing redox and buffering reactions for the different lithological aquifer material.