In this paper, we systematically study spacetimes of gravitational plane waves in Einstein-aether theory. Due to the presence of the timelike aether vector field, now the problem in general becomes overdetermined. In particular, for the linearly polarized plane waves, there are five independent vacuum Einstein-aether field equations for three unknown functions. Therefore, solutions exist only for particular choices of the four free parameters $c_{i}$'s of the theory. We find that there exist eight cases, in two of which any form of gravitational plane waves can exist, similar to that in general relativity, while in the other six cases, gravitational plane waves exist only in particular forms. Beyond these eight cases, solutions either do not exist or are trivial (simply representing a Minkowski spacetime with a constant or dynamical aether field.).

In this paper, we carry out a systematic analysis of the theoretical and observational constraints on the dimensionless coupling constants $c_i$ ($i=1,2,3,4$) of the Einstein-aether theory, taking into account the events GW170817 and GRB 170817A. The combination of these events restricts the deviation of the speed $c_T$ of the spin-2 graviton to the range, $- 3\times 10^{-15} < c_T -1 < 7\times 10^{-16}$, which for the Einstein-aether theory implies $\left|c_{13}\right| \le 10^{-15}$ with $c_{ij} \equiv c_{i} + c_{j}$. The rest of the constraints are divided into two groups: those on the ($c_1, c_{14}$)-plane and those on the ($c_2, c_{14}$)-plane, except the strong-field constraints. The latter depend on the sensitivities $\sigma_{\ae}$ of neutron stars, which are not known at present in the new ranges of the parameters found in this paper.