We documented the step-like increases of water level, flow rate, and water temperatures in a confined aquifer following multiple earthquakes. By employing tidal analysis and a coupled temperature and flow rate model, we find that post-seismic vertical permeability changes and recharge model could explain the co-seismic response. And co-seismic temperature changes are caused by mixing of different volumes of water, with the mixing ratio varying according to each earthquake.