We investigate the origin of the spin Seebeck effect in yttrium iron garnet (YIG) samples for film thicknesses from 20 nm to 50 μm at room temperature and 50 K. Our results reveal a characteristic increase of the longitudinal spin Seebeck effect amplitude with the thickness of the insulating ferrimagnetic YIG, which levels off at a critical thickness that increases with decreasing temperature. The observed behavior cannot be explained as an interface effect or by variations of the material parameters. Comparison to numerical simulations of thermal magnonic spin currents yields qualitative agreement for the thickness dependence resulting from the finite magnon propagation length. This allows us to trace the origin of the observed signals to genuine bulk magnonic spin currents due to the spin Seebeck effect ruling out an interface origin and allowing us to gauge the reach of thermally excited magnons in this system for different temperatures. At low temperature, even quantitative agreement with the simulations is found. The thermal excitation of a spin current by a temperature gradient is commonly called the spin Seebeck effect (SSE) which is detected by the inverse spin Hall effect (ISHE) [1,2], leading to a thermovoltage similar to the charge analogue, the Seebeck effect. Experimental evidence of the SSE, first in ferromagnetic metals [3], and later, both in semiconductors [4] and in insulators [5][6][7][8], has brought up the question about the origin of the SSE. Of particular interest for spin caloritronics is the observation of the SSE in insulators, which allows us to generate pure spin currents in insulating systems.However, the underlying mechanism, properties, and the origin of the observed signals have been highly controversial. Thermally induced magnonic spin currents have been suggested as the origin [9,10], based on the presence of the effect in magnetic insulators, which excludes charge currents as the source. Despite this explanation of the origin of the effect, direct experimental evidence has not been reported. While parasitic interface effects [11] were suggested as an alternative source of the SSE due to a polarization of the paramagnetic detector layer [12], generally, the observed effects are now primarily attributed to magnonic spin currents [13,14].Time resolved experiments trying to address the problem by probing the temporal evolution of the SSE have obtained contradictory results: For film thickness up to 61 nm, no cut-off frequency due to an intrinsic limitation by the SSE was observed [15]. In contrast, for μm thick films, a characteristic rise time was found, and a finite magnon propagation length of the order of several 100 nm was put forward as a possible explanation [16,17]. This clearly calls for study to reveal the origin of this discrepancy as it underlies the fundamental mechanism of the SSE and to determine the intrinsic length scale.To clarify the origin of the measured SSE signals, we present a detailed study of the relevant length scales of the longitudinal SSE (LSSE) coveri...