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