Numerical simulations were conducted to examine the effect of episodic rainfall on nearshore groundwater dynamics in a tidally influenced unconfined coastal aquifer, with a focus on both long-term (yearly) and short-term (daily) behavior of submarine groundwater discharge (SGD) and seawater intrusion (SWI). The results showed nonlinear interactions among the processes driven by rainfall, tides, and density gradients. Rainfall-induced infiltration increased the yearly averaged fresh groundwater discharge to the ocean but reduced the extents of the saltwater wedge and upper saline plume as well as the total rate of seawater circulation through both zones. Overall, the net effect of the interactions led to an increase of the SGD. The nearshore groundwater responded to individual rainfall events in a delayed and cumulative fashion, as evident in the variations of daily averaged SGD and salt stored in the saltwater wedge (quantifying the extent of SWI). A generalized linear model (GLM) along with a Gamma distribution function was developed to describe the delayed and prolonged effect of rainfall events on short-term groundwater behavior. This model validated with results of daily averaged SGD and SWI from the simulations of groundwater and solute transport using independent rainfall data sets, performed well in predicting the behavior of the near-shore groundwater system under the combined influence of episodic rainfall, tides, and density gradients. The findings and developed GLM form a basis for evaluating and predicting SGD, SWI, and associated mass fluxes from unconfined coastal aquifers under natural conditions, including episodic rainfall
