In the research presented here, a new ground water model, FEMCoast, was developed to simulate ground water discharge to the intertidal zone of estuarine systems. This research may be the first attempt to model the ground water discharge process in a tidal estaurine system. The development of FEM Coast was undertaken as no existing ground model was capable of directly simulating the dynamic boundary conditions along the sediment water interface of the intertidal zone. Reproducing the dynamic tidal boundary conditions along the sediment water interface was determined to be essential to replicating the complex salinity gradients observed in the ground water within the intertidal zone. Field data and model results confirmed the presence of a region of ground water where an inverted salinity gradient existed. In this region the concentration of salinity decreased with depth from the ground surface. FEMCoast was also able to reproduce field data on the movement of the near shore water table and ground water discharge rates and patterns. However, the model was not able to replicate the short-term fluctuation in the concentration of salinity within the aquifer due to changes in the concentration of salinity within Cherrystone Inlet. It is believed that the inability to account for the wave action of the tides within the intertidal zone is responsible for this difficulty. The use of FEMCoast integrated with field studies provided a new method to investigate ground water discharge to tidal estuarine systems.