Coastal aquifers may serve as major sources fur freshwater. In many coastal aquifers, intrusion of seawater has become one of the major constraints imposed on groundwater utilization. The management of groundwater in coastal acquifers means making decision as to the pumping rate and the spatial distribution of wells. Several numerical techniques for flow and solute transport simulation can provide the means to achieve this goal. As a basic study to predict the intrusion of seawater in coastal phreatic aquifers, the coupled flow and solute transport analysis was conducted by use of the 3-D finite element code, SWICHA. In order to understand how the location and the shape of freshwater-seawater transition zone were affected by the boundary conditions and hydrogeologic variables, parametric study was carried out.