Galaxy archaeology for wet mergers: Globular cluster age distributions in the Milky Way and nearby galaxies

Context. Identifying past wet merger activity in galaxies has been a longstanding issue in extragalactic formation history studies. Gaia's 6D kinematic measurements in our Milky Way (MW) have vastly extended the possibilities for Galactic archaeology, leading to the discovery of a multitude of early mergers in the MW's past. As recent work has established a link between younger globular clusters (GCs; less than about 10-11 Gyr old) and wet galaxy merger events, the MW provides an ideal laboratory for testing which GC properties can be used to trace extragalactic galaxy formation histories. Aims. To test the hypothesis that GCs trace wet mergers, we relate the measured GC age distributions of the MW and three nearby galaxies, M 31, NGC 1407, and NGC 3115, to their merger histories and interpret the connection with wet mergers through an empirical model for GC formation. Methods. The GC ages of observed galaxies are taken from a variety of studies to analyze their age distributions side-by-side with the model. For the MW, we additionally cross-match the GCs with their associated progenitor host galaxies to disentangle the connection to the GC age distribution. For the modeled GCs, we take galaxies with similar GC age distributions as observed to compare their accretion histories with those inferred through observations. Results. We find that the MW GC age distribution is bimodal, mainly caused by younger GCs (10-11 Gyr old associated with Gaia-Sausage/Enceladus (GSE) and in part by unassociated high-energy GCs. The GSE GC age distribution also appears to be bimodal. We propose that the older GSE GCs (12-13 Gyr old) were accreted together with GSE, while the younger ones formed as a result of the merger. For the nearby galaxies, we find that clear peaks in the GC age distributions coincide with active early gas-rich merger phases. Even small signatures in the GC age distributions agree well with the expected wet formation histories of the galaxies inferred through other observed tracers. From the models, we predict that the involved cold gas mass can be estimated from the number of GCs found in the formation burst. Conclusions. Multimodal GC age distributions can trace massive wet mergers as a result of GCs being formed through them. From the laboratory of our own MW and nearby galaxies we conclude that the ages of younger GC populations of galaxies can be used to infer the wet merger history of a galaxy.