Tidally-induced head fluctuation is a natural phenomenon in coastal regions. The discharge of groundwater through sediments will occur anywhere that the aquifer is hydraulically connected to a surface water body and the time averaged tidally-influenced water level in the aquifer is higher than sea level, and almost all coastal regions are subject to such flow. With the development of coastal areas, the discharge of groundwater contaminants into tidally affected coastal water bodies has become a significant problem. Biota that live in the benthic region are known to be sensitive to the concentration of discharging anthropogenic chemical compounds. Thus the contaminant concentration entering the benthos is of very significant practical importance and its study is the focus of this dissertation. An investigation of the effect of tides on the concentration of groundwater contaminants discharging to a surface water body is studied using a one-dimensional homogeneous sand column. Results of the experiment are confirmed using a three-dimensional heterogeneous groundwater tank model. A constant water level is imposed upgradient, and the downgradient water level is controlled by a wave generator that controls the hydraulic head to mimic a 12 hour tidal fluctuation. The experimental results demonstrate that the tidal fluctuations in the downgradient reservoir result in a decrease in average contaminant concentration at the point of groundwater discharge to the surface-water body. The further upstream the well is located, the smaller the amplitude of the concentration oscillation. In addition, upstream migration of concentration oscillations is observed in spite of a net downgradient flow. Fourier analysis suggests that the dominant frequency of the peaks of pressure and chemical data at different locations along the length of the column is identically two cycles per day and that the amplitudes of the concentration oscillations increase with time at measurement locations at the upstream responding probes. As the classical groundwater flow and transport model cannot reproduce the phenomena we observed, an innovative multi-mobility model, is proposed with one highly mobile liquid phase, one less mobile liquid phase and a solid phase. Averaging theory is applied to develop the mass conservation equation from the microscale to the macroscale and facilitate the reduction of dimensionality to obtain one-dimensional governing equations with closure relations. A new finite volume method is utilized to solve the resulting equations. The simulation confirms the existence of the enhanced tidally-induced mixing process.

The prediction of coastal groundwater movement is necessary in coastal management. However, the study in this field is still a great challenge due to the involvement of tidal-groundwater interactions and the phenomena of hydrodynamic dispersion between salt-fresh water in the coastal region. To date, numerous theories for groundwater dynamic have been made available in analytical, numerical and also experimental forms. Nevertheless, most of them are based on the zeroth-order shallow flow, i.e. Boussinesq approximation. Two main components for coastal unconfined aquifer have been completed in this Thesis: the vertical beach model and the sloping beach model. Both solutions are solved in closed-form up to higher order with shallow water parameter ([epsilon]) and tidal amplitude parameter ([alpha]). The vertical beach solution contributes to the higher-order tidal fluctuations while the sloping beach model overcomes the shortcomings in the existing solutions. From this study, higher-order components are found to be significant especially for larger value of [alpha] and [epsilon]. Other parameters such as hydraulic conductivity (K) and the thickness of aquifer (D) also affect the water table fluctuations. The new sloping solution demonstrated the significant influence of beach slope ([beta]) on the water table fluctuations. A comprehensive comparison between previous solution and the present sloping solution have been performed mathematically and numerically and the present solution has been demonstrated to provide a better prediction