Opacity broadening and interpretation of suprathermal CO linewidths: Macroscopic Turbulence and Tangled Molecular Clouds

(Abridged) Many of the observed CO line profiles exhibit broad linewidths that greatly exceed the thermal broadening expected within molecular clouds. These suprathermal CO linewidths are assumed to be originated from the presence of unresolved supersonic motions inside clouds. Typically overlooked in the literature, in this paper we aim to quantify the impact of the opacity broadening effects on the current interpretation of the CO suprathermal line profiles. Without any additional contributions to the gas velocity field, a large fraction of the apparently supersonic (${\cal M}\sim$2-3) linewidths measured in both $^{12}$CO and $^{13}$CO (J=1-0) lines can be explained by the saturation of their corresponding sonic-like, optically-thin C$^{18}$O counterparts assuming standard isotopic fractionation. Combined with the presence of multiple components detected in our C$^{18}$O spectra, these opacity effects seem to be also responsible of the highly supersonic linewidths (${\cal M}>$8-10) detected in the broadest $^{12}$CO and $^{13}$CO spectra in Taurus. Our results demonstrate that most of the suprathermal $^{12}$CO and $^{13}$CO linewidths could be primarily created by a combination of opacity broadening effects and multiple gas velocity components blended in these saturated emission lines. Once corrected by their corresponding optical depth, each of these gas components present transonic intrinsic linewidths consistently traced by the three CO isotopologues within a factor of 2. Highly correlated and velocity-coherent at large scales, the largest and highly supersonic velocity differences inside clouds are generated by the relative motions between individual gas components. This highly discretized structure of the molecular gas traced in CO suggest that the gas dynamics inside molecular clouds could be better described by the properties of a fully-resolved macroscopic turbulence.