Mass stratification in the globular cluster system revealing the assembly history of the nearest S0 galaxy NGC 3115

Galaxy formation and evolution is hierarchical. The most massive galaxies are thought to form their central regions early through violent dissipational processes, then grow inside-out by accreting smaller satellites. While widely supported, direct observational confirmation of this process in individual galaxies remains lacking, except for the Milky Way. We present a detailed analysis of globular cluster (GC) candidates within a $70^\prime$ ($\sim190$ kpc) radius around the nearest S0 galaxy, NGC 3115, using images in \textit{g,r,z} bands from the DESI Legacy Imaging Surveys and data from Gaia. We report the discovery of mass stratification in the GC system (GCS), evident in two ways: first, the effective radius of the GCS increases monotonically from the bright to faint end, up to the detection limit near the turnover magnitude of the GC luminosity function (GCLF); second, the GCLF shows fainter turnover magnitudes and smaller standard deviations at larger galactocentric radii. This stratification cannot be readily explained by radial migration or tidal dissolution, but most likely reflects the hierarchical assembly of NGC 3115's stellar halo, with later-accreted satellites deposited across broader galactocentric distances. This interpretation is supported by cosmological simulations of subhalos with comparable mass and bulge-to-total mass ratios and is consistent with the negative color gradients observed in the GCS. Additionally, we identify several substructures within the GCS, indicating ongoing assembly of NGC 3115. This work highlights the power of GCS as tracers of galaxy assembly and sets the stage for upcoming space-based wide-field imaging surveys to constrain the assembly of massive galaxies.