Scientific progress in the context of seismic precursors reveals a systematic mechanism, namely lithosphere–atmosphere–ionosphere coupling (LAIC), to elaborate the underlying physical processes related to earthquake preparation phases. In this study, a comprehensive analysis was conducted for two earthquakes that occurred on the sea coast through tidal force fluctuation to investigate ocean–lithosphere–atmosphere–ionosphere coupling (OLAIC), based on oceanic parameters (i.e., sea potential temperature and seawater salinity), air temperature and electron density profiles. The interrupted enhancement and diffusion process of thermal anomalies indicate that the intensity of seismic anomalies in the atmosphere is affected by the extent of land near the epicenter. By observing the evolution of the ocean interior, we found that the deep water was lifted and formed upwelling, which then diffused along the direction of plate boundaries with an “intensification-peak-weakening” trend under the action of the accelerated subduction of tectonic plates. Furthermore, the analysis shows that the seismic anomalies have two propagation paths: (i) along active faults, with the surface temperature rising as the initial performance, then the air pressure gradient being generated, and finally the ionosphere being disturbed; (ii) along plate boundaries, upwelling, which is the initial manifestation, leading to changes in the parameters of the upper ocean. The results presented in this study can contribute to understanding the intrinsic characteristics of OLAIC.