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.).

We numerically study the gravitational collapse of a massless scalar field with spherical symmetry in Einstein-aether theory, and show that apparent, spin-0 and dynamical universal horizons (dUHs) can be all formed. The spacetime and the aether field are well-behaved and regular, including regions nearby these horizons (but away from the center of spherical symmetry). The spacetime outside the apparent and spin-0 horizons settles down to a static configuration, and some of such resulting static black holes were already found numerically in the literature. On the other hand, the proper distance of the outermost dUH from the apparent (or spin-0) horizon keeps increasing on $\mbox{\ae}$ther-orthogonal time slices. This indicates that the outermost dUH is evolving into the causal boundary, even for excitations with large speeds of propagation.

In this paper, we study the non-projectable 2d Ho\v{r}ava gravity coupled with a non-relativistic scalar field, where the coupling is in general non-minimal and of the form $f(\phi)R$, where $f(\phi)$ is an arbitrary function of the scalar field $\phi$, and $R$ denotes the 2d Ricci scalar. In particular, we first investigate the Hamiltonian structure, and show that there are two-first and two-second class constraints, similar to the pure gravity case, but now the local degree of freedom is one, due to the presence of the scalar field. Then, we present various exact stationary solutions of this coupled system, and find that some of them represent black holes but now with universal horizons as their boundaries. At these horizons, the Hawking radiations are thermal with temperatures proportional to their surface gravities, which normally depend on the non-linear dispersion relations of the particles radiated, similar to the (3+1)-dimensional case.