We study 4 dimensional (4d) gravitational waves (GWs) with compact wavefronts, generalizing Robinson-Trautman (RT) solutions in Einstein gravity with an arbitrary cosmological constant. We construct the most general solution of the GWs in the presence of a causal, timelike, or null boundary when the usual tensor modes are turned off. Our solution space besides the shape and topology of the wavefront which is a generic compact, smooth, and orientable 2d surface Σ, is specified by a vector over Σ satisfying the conformal Killing equation and two scalars that are arbitrary functions over the causal boundary, the boundary modes (soft hair). We work out the symplectic form over the solution space using covariant phase space formalism and analyze the boundary symmetries and charges. The algebra of surface charges is a Heisenberg algebra. Only the overall size of the compact wavefront and not the details of its shape appears in the boundary symplectic form and is canonical conjugate to the overall mass of the GW. Hence, the information about the shape of the wavefront can’t be probed by the boundary observer. We construct a boundary energy-momentum tensor and a boundary current, whose conservation yields the RT equation for both asymptotically AdS and flat spacetimes. The latter provides a hydrodynamic description for our RT solutions.