Using multiple observational arguments, recent work has shown that cosmological GRBs are typically viewed at angles within, or close to the cores of their relativistic jets. One of those arguments relied on the lack of tens-of-days-long periods of very shallow evolution that would be seen in the afterglow light-curves of GRBs viewed at large angles. Motivated by these results, we consider that GRBs efficiently produce $\gamma$-rays only within a narrow region around the core. We show that, on these near-core lines-of-sight, structured jets naturally produce shallow phases in the X-ray afterglow of GRBs. These plateaus would be seen by a large fraction of observers and would last between $10^2-10^5$ s. They naturally reproduce the observed distributions of time-scales and luminosities as well as the inter-correlations between plateau duration, plateau luminosity and prompt $\gamma$-ray energy. An advantage of this interpretation is that it involves no late time energy injection which would be both challenging from the point of view of the central engine and, as we show here, less natural given the observed correlations between plateau and prompt properties.

The up-to-now only binary neutron star merger gravitational wave event with detected electromagnetic counterparts, GRB170817A, occurred in a rarefied medium with a density smaller than $10^{-3}-10^{-2}~{\rm cm}^{-3}$. As neutron star binaries are imparted kicks upon formation, and due to their long delay-times before merger, such low-density circum-merger media are generally expected. However, there is some indirect evidence for a class of faster-merging binaries, which would coalesce in denser environments. Nonetheless, astronomical data is largely inconclusive on the possibility of these high-density mergers. We describe a method to directly probe this hypothetical population of high-density mergers through multi-messenger observations of binary neutron star merger afterglows, exploiting the sharp sensitivity of these to the circum-merger medium density. This method is based on a yet-to-be-collected sample of merger afterglows. Its constraining power is large even with a small sample of events. We discuss the method's limitations and applicability. In the upcoming era of 3$^{\rm rd}$-generation gravitational wave detectors, this method's potential will be fully realized as it will allow to probe mergers having occurred soon after the peak of cosmic star formation.

Following the historical observations of GW170817 and its multi-wavelength afterglow, more radio afterglows from neutron star mergers are expected in the future as counterparts to gravitational wave inspiral signals. We wish to describe these events using our current knowledge of the population of neutron star mergers coming from gamma-ray burst science, and taking into account the sensitivities of current and future gravitational wave and radio detectors. We combine analytical models for the merger gravitational wave and radio afterglow signals to a population model prescribing the energetics, circum-merger density and other relevant parameters of the mergers. We report the expected distributions of observables (distance, orientation, afterglow peak time/flux, etc.) from future events and study how these can be used to further probe the population of binary neutron stars, their mergers and related outflows during future observing campaigns. In the case of the O3 run of the LIGO-Virgo Collaboration, the radio afterglow of one third of gravitational-wave-detected mergers should be detectable (and detected if the source is localized thanks to the kilonova counterpart) by the Very Large Array, and these events should have viewing angles similar to that of GW170817. These findings confirm the radio afterglow as a powerful insight on these events, though some key afterglow-related techniques, such as very long baseline interferometry imaging of the merger remnant, may no longer be feasible as the gravitational wave horizon increases.

Following the historical observation of GW170817 and its electromagnetic follow-up, new neutron star merger afterglows are expected to be observed as counterparts to gravitational wave signals during the next science runs of the gravitational interferometer network. The diversity of the observed population of afterglows of these future events is subject to various factors, which are (i) intrinsic, such as the energy of the ejecta, (ii) environmental, such as the ambient medium density or (iii) observational, such as the viewing angle and distance of the source. Through prescribing a population of mergers and modelling their afterglows, we study the diversity of those events to be observed jointly in gravitational waves and electromagnetic bands. In the future, observables of detected events such as viewing angle, distance, afterglow peak flux or proper motion will form distributions which together with predictions from our study will provide insight on neutron star mergers and their environments.