We present high sensitivity 12CO and 13CO (1-0) molecular line maps covering the full extent of the parsec scale Haro~6-10 Herbig-Haro (HH) flow. We report the discovery of a molecular CO outflow along the axis of parsec-scale HH flow. Previous molecular studies missed the identification of the outflow probably due to their smaller mapping area and the confusing spectral features present towards the object. Our detailed molecular line study of the full 1.6 pc extent of the optical flow shows evidence for both blueshifted and redshifted gas set in motion by Haro~6-10 activity. The molecular outflow is centered at Haro~6-10, with redshifted gas being clumpy and directed towards the northeast, while blueshifted gas is in the southwest direction. The molecular gas terminates well within the cloud, short of the most distant HH objects of the optical flow. Contamination from an unrelated cloud along the same line of sight prevents a thorough study of the blueshifted outflow lobe and the mass distribution at the lowest velocities in both lobes. The cloud core in which Haro~6-10 is embedded is filamentary and flattened in the east-west direction. The total cloud mass is calculated from 13CO(1-0) to be ~200Msun. The lower limit of the mass associated with the outflow is ~0.25Msun.

We present high sensitivity 12/13CO(1-0) molecular line maps covering the full extent of the parsec scale L1551 molecular outflow, including the redshifted east-west (EW) flow. We also present 12CO(3-2) data that extends over a good fraction of the area mapped in the 1-0 transition. We compare the molecular data to widefield, narrow-band optical emission in H$\alpha$. While there are multiple outflows in the L1551 cloud, the main outflow is oriented at 50\arcdeg position angle and appears to be driven by embedded source(s) in the central IRS 5 region. The 3-2 data indicate that there may be molecular emission associated with the L1551 NE jet, within the redshifted lobe of main outflow. We have also better defined the previously known EW flow and believe we have identified its blueshifted counterpart. We further speculate that the origin of the EW outflow lies near HH 102. We use velocity dependent opacity correction to estimate the mass and the energy of the outflow. The resulting mass spectral indices from our analysis, are systematically lower (less steep) than the power law indices obtained towards other outflows in several recent studies that use a similar opacity correction method. We show that systematic errors and biases in the analysis procedures for deriving mass spectra could result in errors in the determination of the power-law indices. The mass spectral indices, the morphological appearance of the position-velocity plots and integrated intensity emission maps of the molecular data, compared with the optical, suggest that jet-driven bow-shock entrainment is the best explanation for the driving mechanism of outflows in L1551. The kinetic energy of the outflows is found to be comparable to the binding energy of the cloud and sufficient to maintain the turbulence in the L1551 cloud.