Abstract
Because the Japanese archipelago is surrounded by sea and important social-economic activities concentrate in areas below sea-level, it is prerequisite to predict impacts of sea-level rise on groundwater systems in coastal aquifers. For applying numerical models to areas with inclined land surface which allows storage of surface water, we developed a simulation code, SIFEC, based on coupling of finite elements for density-dependent flow and method of characteristics linked with finite elements for mass transport. To improve accuracy for interpolating concentrations in solving transport equation, single-step reverse particle tracking is used with biquadratic interpolation scheme. Accuracy of the technique was compared with results of experiments and other simulation techniques, and excellent results were obtained.
