A high-resolution coastal ocean model was developed to simulate the temporal/spatial variability of the Kennebec-Androscoggin (K-A) river plume and the circulation in Casco Bay. The model results agree favorably with the moored and shipboard observations of velocity, temperature, and salinity. The surface salinity gradient was used to distinguish the plume from the ambient coastal water. The calculated plume thickness suggests that the K-A plume is surface trapped. Its horizontal scales correlate well with Q 0.25 , where Q is the volume discharge of the rivers. Directional spreading is affected by the wind with the upwelling favorable wind transporting the plume water offshore. Both the wind and the tide also enhance mixing in the plume. The inclusion of a wetting-and-drying (WAD) scheme appears to enhance the mixing and entrainment processes near the estuary. The plume becomes thicker near the mouth of the estuary, the outflow velocity of the plume is weaker, and the radius of the river plume shrinks. The flow field in the model run with the WAD is noisier, not only in shallow areas of Casco Bay but also in the plume and even on the shelf. We speculate that the WAD processes can affect much larger areas than the intertidal zones, especially via a river plume that feeds into a coastal current.