Groundwater under agricultural lands is often contaminated by nitrate. In southern Japan, aquifers are covered by volcanic ash soils that can leach nitrate to the groundwater. In this study, column experiments using two volcanic ash soils (Kuroboku and Akahoya) were carried out. A mixed solution of KNO 3 and K 2 SO 4 was used in the leaching experiments. Based on the experiments, a reactive transport model was developed using the Constrained Interpolation Profile method for ion transport and used to calculate chemical equilibrium for the cation exchange reactions. Anion adsorption was modeled by retardation in the numerical model. The developed simulation model results were compared to results obtained by the reactive transport model PHAST. The developed model was shown to quite well reproduce general characteristics of the experimental results. Also, the developed model results agreed well with results from PHAST. Slight discrepancy between observed and calculated breakthrough curves was probably caused by ignoring the kinetic reaction in the model calculations.