The growing interest in lignin valorization in the past decades calls for analytical techniques for lignin characterization, ranging from wet chemistry techniques to highly sophisticated chromatographic and spectroscopic methods. One of the key parameters to consider is the molecular weight profile of lignin, which is routinely determined by size-exclusion chromatography; however, this is by no means straightforward and is prone to being hampered by considerable errors. Our study expands the fundamental understanding of the bias-inducing mechanisms in gel permeation chromatography (GPC), the magnitude of error originating from using polystyrene standards for mass calibration, and an evaluation of the effects of the solvent and type of lignin on the observed bias. The developed partial least-squares (PLS) regression model for lignin-related monomers revealed that lignin is prone to association mainly via hydrogen bonding. This hypothesis was supported by functional group-based analysis of the bias as well as pulse field gradient (pfg) diffusion NMR spectroscopy of model compounds in THF-d8. Furthermore, although the lack of standards hindered drawing conclusions based on functionalities, direct infusion electrospray ionization mass spectrometry indicated that the relative bias decreases considerably for higher molecular weight species. The results from pfg-diffusion NMR spectroscopy on whole lignin samples were comparable when the same solvents were used in both experiments; in addition, the comparison between results obtained by pfg-diffusion NMR in different solvents gives some additional insights into the aggregation.