We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $\Sigma$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$\mu$m [CII] archival data toward the northern part of Orion A are also compared to the $^{13}$CO data to test whether the position and velocity offsets between the emission from these two transitions resemble those predicted by a cloud-cloud collision model. We find that the northern part of Orion A, including regions ONC-OMC-1, OMC-2, OMC-3 and OMC-4, shows qualitative agreements with the cloud-cloud collision scenario, while in one of the southern regions, NGC1999, there is no indication of such a process in causing the birth of new stars. On the other hand, another southern cluster, L1641N, shows slight tendencies of cloud-cloud collision. Overall, our results support the cloud-cloud collision process as being an important mechanism for star cluster formation in Orion A.

We identify 45 protostellar outflows in CO maps of the Orion A giant molecular cloud from the CARMA-NRO Orion survey. Our sample includes 11 newly detected outflows. We measure the mass and energetics of the outflows, including material at low-velocities by correcting for cloud contributions. The total momentum and kinetic energy injection rates of outflows is comparable to the turbulent dissipation rate of the cloud. We also compare the outflow position angles to the orientation of C$^{18}$O filaments. We find that the full sample of outflows is consistent with being randomly oriented with respect to the filaments. A subsample of the most reliable measurements shows a moderately perpendicular outflow-filament alignment which may reflect accretion of mass across filaments and onto the protostellar cores.

Context. The [CII] 158 {\mu}m far-infrared (FIR) fine-structure line is one of the most important cooling lines of the star-forming interstellar medium (ISM). High spectral resolution observations have shown complex structures in the line profiles of the [CII] emission. Aims. Our aim is to determine whether the complex profiles observed in [^{12}CII] are due to individual velocity components along the line-of-sight or to self-absorption based on a comparison of the [^{12}CII] and isotopic [^{13}CII] line profiles. Methods. Deep integrations with the SOFIA/upGREAT 7-pixel array receiver in M43, Horsehead~PDR, Monoceros~R2, and M17~SW allow for the detection of optically thin [^{13}CII] emission lines, along with the [^{12}CII] emission lines, with a high signal-to-noise ratio. We first derived the [^{12}CII] optical depth and the [CII] column density from a single component model. However, the complex line profiles observed require a double layer model with an emitting background and an absorbing foreground. A multi-component velocity fit allows us to derive the physical conditions of the [CII] gas: column density and excitation temperature. Results. We find moderate to high [^{12}CII] optical depths in all four sources and self-absorption of [^{12}CII] in Mon R2 and M17 SW. The high column density of the warm background emission corresponds to an equivalent Av of up to 41 mag. The foreground absorption requires substantial column densities of cold and dense [CII] gas, with an equivalent Av ranging up to about 13 mag. Conclusions. The column density of the warm background material requires multiple photon-dominated region (PDR) surfaces stacked along the line of sight and in velocity. The substantial column density of dense and cold foreground [CII] gas detected in absorption cannot be explained with any known scenario and we can only speculate on its origins

The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

[13CII] observations in several Galactic sources show that the fine-structure [12CII] emission is often optically thick (the optical depths around 1 to a few). The aim of our study is to test whether this also affects the [12CII] emission from nearby galaxies like the Large Magellanic Cloud (LMC). We observed three star-forming regions in the LMC with upGREAT on board SOFIA at the frequency of the [CII] line. The 4GHz band width covers all three hyperfine lines of [13CII] simultaneously. For the analysis, we combined the [13CII] F=1-0 and F=1-1 hyperfine components, as they do not overlap with the [12CII] line in velocity. Three positions in N159 and N160 show an enhancement of [13CII] compared to the abundance-ratio-scaled [12CII] profile. This is likely due to the [12CII] line being optically thick, supported by the fact that the [13CII] line profile is narrower than [12CII], the enhancement varies with velocity, and the peak velocity of [13CII] matches the [OI] 63um self-absorption. The [12CII] line profile is broader than expected from a simple optical depth broadening of the [13CII] line, supporting the scenario of several PDR components in one beam having varying [12CII] optical depths. The derived [12CII] optical depth at three positions (beam size of 14arcsec, corresponding to 3.4pc) is 1--3, which is similar to values observed in several Galactic sources shown in previous studies. If this also applies to distant galaxies, the [CII] intensity will be underestimated by a factor of approximately 2.

Although the 2P3/2-2P1/2 transition of [CII] at 158um is known to be an excellent tracer of active star formation, we still do not have a complete understanding of where within star formation regions the emission originates. Here, we use SOFIA upGREAT observations of [CII] emission toward the HII region complex Sh2-235 (S235) to better understand in detail the origin of [CII] emission. We complement these data with a fully-sampled Green Bank Telescope radio recombination line map tracing the ionized hydrogen gas. About half of the total [CII] emission associated with S235 is spatially coincident with ionized hydrogen gas, although spectroscopic analysis shows little evidence that this emission is coming from the ionized hydrogen volume. Velocity-integrated [CII] intensity is strongly correlated with WISE 12um intensity across the entire complex, indicating that both trace ultra-violet radiation fields. The 22um and radio continuum intensities are only correlated with [CII] intensity in the ionized hydrogen portion of the S235 region and the correlations between the [CII] and molecular gas tracers are poor across the region. We find similar results for emission averaged over a sample of external galaxies, although the strength of the correlations is weaker. Therefore, although many tracers are correlated with the strength of [CII] emission, only WISE 12um emission is correlated on small-scales of the individual HII region S235 and also has a decent correlation at the scale of entire galaxies. Future studies of a larger sample of Galactic HII regions would help to determine whether these results are truly representative.

Understanding the matter cycle in the interstellar medium of galaxies from the assembly of clouds to star formation and stellar feedback remains an important and exciting field in comtemporary astrophysics. Many open questions regarding cloud and structure formation, the role of turbulence, and the relative importance of the various feedback processes can only be addressed with observations of spectrally resolved lines. We here stress the importance of two specific sets of lines: the finestructure lines of atomic carbon as a tracer of the dark molecular gas and mid-J CO lines as tracers of the warm, active molecular gas in regions of turbulence dissipation and feedback. The observations must cover a wide range of environments (i.e., physical conditions), which will be achieved by large scale surveys of Galactic molecular clouds, the Galactic Center, the Magellanic clouds, and nearby galaxies. To date, such surveys are completely missing and thus constitute an important science opportunity for the next decade and beyond. For the successful interpretation of the observations, it will be essential to combine them with results from (chemical) modelling and simulations of the interstellar medium.

We investigate the relationship between turbulence and feedback in the Orion A molecular cloud using maps of $^{12}$CO(1-0), $^{13}$CO(1-0) and C$^{18}$O(1-0) from the CARMA-NRO Orion survey. We compare gas statistics with the impact of feedback in different parts of the cloud to test whether feedback changes the structure and kinematics of molecular gas. We use principal component analysis, the spectral correlation function, and the spatial power spectrum to characterize the cloud. We quantify the impact of feedback with momentum injection rates of protostellar outflows and wind-blown shells as well as the surface density of young stars. We find no correlation between shells or outflows and any of the gas statistics. However, we find a significant anti-correlation between young star surface density and the slope of the $^{12}$CO spectral correlation function, suggesting that feedback may influence this statistic. While calculating the principal components, we find peaks in the covariance matrix of our molecular line maps offset by 1-3 km s$^{-1}$ toward several regions of the cloud which may be produced by feedback. We compare these results to predictions from molecular cloud simulations.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

The ubiquitous presence of filamentary structures in the interstellar medium asks for an unbiased characterization of their properties including a stability analysis. We propose a novel technique to measure the spectrum of filaments in any two-dimensional data set. Using anisotropic wavelets we can quantify and distinguish local and global anisotropies and measure the size distribution of filaments. The wavelet analysis does not need any assumptions on the alignment or shape of filaments in the maps, but directly measures their typical spatial dimensions. In a rigorous test program, we calibrate the scale-dependence of the method and test the angular and spatial sensitivity. We apply the method to molecular line maps from magneto-hydrodynamic (MHD) simulations and observed column density maps from Herschel observations. When applying the anisotropic wavelet analysis to the MHD data, we find that the observed filament sizes depend on the combination of magnetic-field dominated density-velocity correlations with radiative transfer effects. This can be exploited by observing tracers with different optical depth to measure the transition from a globally ordered large-scale structure to small-scale filaments with entangled field lines. The unbiased view to Herschel column density maps does not confirm a universal characteristic filament width. The map of the Polaris Flare shows an almost scale-free filamentary spectrum up to the size of the dominating filament of about 0.4pc. For the Aquila molecular cloud the range of filament widths is limited to 0.05-0.2pc. The filaments in Polaris show no preferential direction in contrast to the global alignment that we trace in Aquila. By comparing the power in isotropic and anisotropic structures we can measure the relative importance of spherical and cylindrical collapse modes and their spatial distribution.

The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{{\rm CII}}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{{\rm CII}}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{{\rm CII}}$ to change in environments with a lower or higher radiation field than simulated here.

In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

Spatial association of clumps from different tracers turns out to be a valuable tool to determine the physical properties of molecular clouds. It provides a reliable estimate for the $X$-factors, serves to trace the density of clumps seen in column densities only and allows to measure the velocity dispersion of clumps identified in dust emission. We study the spatial association between clump populations, extracted by use of the GAUSSCLUMPS technique from $^{12}$CO (1-0), $^{13}$CO (1-0) line maps and Herschel dust-emission maps of the star-forming region Rosette, and analyse their physical properties. All CO clumps that overlap with another CO or dust counterpart are found to be gravitationally bound and located in the massive star-forming filaments of the molecular cloud. They obey a single mass-size relation $M_{\rm cl}\propto R_{\rm cl}^\gamma$ with $\gamma\simeq3$ (implying constant mean density) and display virtually no velocity-size relation. We interpret their population as low-density structures formed through compression by converging flows and still not evolved under the influence of self-gravity. The high-mass parts of their clump mass functions are fitted by a power law ${\rm d}N_{\rm cl}/{\rm d}\,\log M_{\rm cl}\propto M_{\rm cl}^{\Gamma}$ and display a nearly Salpeter slope $\Gamma\sim-1.3$. On the other hand, clumps extracted from the dust-emission map exhibit a shallower mass-size relation with $\gamma=2.5$ and mass functions with very steep slopes $\Gamma\sim-2.3$ even if associated with CO clumps. They trace density peaks of the associated CO clumps at scales of a few tenths of pc where no single density scaling law should be expected.

Photon Dominated Regions (PDRs) are interfaces between the mainly ionized and mainly molecular material around young massive stars. Analysis of the physical and chemical structure of such regions traces the impact of far-ultraviolet radiation of young massive stars on their environment. We present results on the physical and chemical structure of the prototypical high UV-illumination edge-on Orion Bar PDR from an unbiased spectral line survey with a wide spectral coverage. A spectral scan from 480-1250 GHz and 1410-1910 GHz at 1.1 MHz resolution was obtained by the HIFI instrument onboard the Herschel Space Observatory. For molecules with multiple transitions we used rotational diagrams to obtain excitation temperatures and column densities. For species with a single detected transition we used an optically thin LTE approximation. In case of species with available collisional rates, we also performed a non-LTE analysis to obtain kinetic temperatures, H2 volume densities, and column densities. About 120 lines corresponding to 29 molecules (including isotopologues) have been detected in the Herschel/HIFI line survey, including 11 transitions of CO, 7 transitions of 13CO, 6 transitions of C18O, 10 transitions of H2CO, and 6 transitions of H2O. Most species trace kinetic temperatures in the range between 100 and 150 K and H2 volume densities in the range between 10^5 and 10^6 cm^-3. The species with temperatures and / or densities outside of this range include the H2CO transitions tracing a very high temperature (315 K) and density (1.4x10^6 cm^-3) component and SO corresponding to the lowest temperature (56 K) measured as a part of this line survey. The observed lines/species reveal a range of physical conditions (gas density /temperature) involving structures at high density / high pressure, obsoleting the traditional 'clump / interclump' picture of the Orion Bar.

Context. Models of photon-dominated regions (PDRs) still fail to fully reproduce some of the observed properties, in particular the combination of the intensities of different PDR cooling lines together with the chemical stratification, as observed e.g. for the Orion Bar PDR. Aims. We aim to construct a numerical PDR model, KOSMA-\tau 3D, to simulate full spectral cubes of line emission from arbitrary PDRs in three dimensions (3D). The model is to reproduce the intensity of the main cooling lines from the Orion Bar PDR and the observed layered structure of the different transitions. Methods. We build up a 3D compound, made of voxels ("3D pixels") that contain a discrete mass distribution of spherical "clumpy" structures, approximating the fractal ISM. To analyse each individual clump the new code is combined with the KOSMA-\tau PDR model. Probabilistic algorithms are used to calculate the local FUV flux for each voxel as well as the voxel-averaged line emissivities and optical depths, based on the properties of the individual clumps. Finally, the computation of the radiative transfer through the compound provides full spectral cubes. To test the new model we try to simulate the structure of the Orion Bar PDR and compare the results to observations from HIFI/Herschel and from the Caltech Submillimetre Observatory (CSO). In this context new Herschel data from the HEXOS guaranteed-time key program is presented. Results. Our model is able to reproduce the line integrated intensities within a factor 2.5 and the observed stratification pattern within 0.016 pc for the [Cii] 158 \mu m and different 12/13 CO and HCO+ transitions, based on the representation of the Orion Bar PDR by a clumpy edge-on cavity wall. In the cavity wall, a large fraction of the total mass needs to be contained in clumps. The mass of the interclump medium is constrained by the FUV penetration. Furthermore, ...

We present results from a comprehensive submillimeter spectral survey toward the source Orion South, based on data obtained with the HIFI instrument aboard the \textit{Herschel Space Observatory}, covering the frequency range 480 to 1900 GHz. We detect 685 spectral lines with S/N $>$ 3$\sigma$, originating from 52 different molecular and atomic species. We model each of the detected species assuming conditions of Local Thermodynamic Equilibrium. This analysis provides an estimate of the physical conditions of Orion South (column density, temperature, source size, \& V$_{LSR}$). We find evidence for three different cloud components: a cool (T$_{ex} \sim 20-40$ K), spatially extended ($> 60"$), and quiescent ($\Delta V_{FWHM} \sim 4$ km s $^{-1}$) component; a warmer (T$_{ex} \sim 80-100$ K), less spatially extended ($\sim 30"$), and dynamic ($\Delta V_{FWHM} \sim 8$ km s $^{-1}$) component, which is likely affected by embedded outflows; and a kinematically distinct region (T$_{ex}$ $>$ 100 K; V$_{LSR}$ $\sim$ 8 km s $^{-1}$), dominated by emission from species which trace ultraviolet irradiation, likely at the surface of the cloud. We find little evidence for the existence of a chemically distinct "hot core" component, likely due to the small filling factor of the hot core or hot cores within the \textit{Herschel} beam. We find that the chemical composition of the gas in the cooler, quiescent component of Orion South more closely resembles that of the quiescent ridge in Orion-KL. The gas in the warmer, dynamic component, however, more closely resembles that of the Compact Ridge and Plateau regions of Orion-KL, suggesting that higher temperatures and shocks also have an influence on the overall chemistry of Orion South.

The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 mm, H2 S(3) rotational line at 9.3 mm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9x10^(-10) and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.

The CH$^+$ ion is a key species in the initial steps of interstellar carbon chemistry. Its formation in diverse environments where it is observed is not well understood, however, because the main production pathway is so endothermic (4280 K) that it is unlikely to proceed at the typical temperatures of molecular clouds. We investigation CH$^+$ formation with the first velocity-resolved spectral mapping of the CH$^+$ $J=1-0, 2-1$ rotational transitions, three sets of CH $\Lambda$-doubled triplet lines, $^{12}$C$^+$ and $^{13}$C$^+$, and CH$_3$OH 835~GHz E-symmetry Q branch transitions, obtained with Herschel/HIFI over $\approx$12 arcmin$^2$ centered on the Orion BN/KL source. We present the spatial morphologies and kinematics, cloud boundary conditions, excitation temperatures, column densities, and $^{12}$C$^+$ optical depths. Emission from C$^+$, CH$^+$, and CH is indicated to arise in the diluted gas, outside of the explosive, dense BN/KL outflow. Our models show that UV-irradiation provides favorable conditions for steady-state production of CH$^+$ in this environment. Surprisingly, no spatial or kinematic correspondences of these species are found with H$_2$ S(1) emission tracing shocked gas in the outflow. We propose that C$^+$ is being consumed by rapid production of CO to explain the lack of C$^+$ and CH$^+$ in the outflow, and that fluorescence provides the reservoir of H$_2$ excited to higher ro-vibrational and rotational levels. Hence, in star-forming environments containing sources of shocks and strong UV radiation, a description of CH$^+$ formation and excitation conditions is incomplete without including the important --- possibly dominant --- role of UV irradiation.

Probability distribution functions (PDFs) of column densities are an established tool to characterize the evolutionary state of interstellar clouds. Using simulations, we show to what degree their determination is affected by noise, line-of-sight contamination, field selection, and the incomplete sampling in interferometric measurements. We solve the integrals that describe the convolution of a cloud PDF with contaminating sources and study the impact of missing information on the measured column density PDF. The effect of observational noise can be easily estimated and corrected for if the root mean square (rms) of the noise is known. For $\sigma_{noise}$ values below 40% of the typical cloud column density, $N_{peak}$, this involves almost no degradation of the accuracy of the PDF parameters. For higher noise levels and narrow cloud PDFs the width of the PDF becomes increasingly uncertain. A contamination by turbulent foreground or background clouds can be removed as a constant shield if the PDF of the contamination peaks at a lower column or is narrower than that of the observed cloud. Uncertainties in the definition of the cloud boundaries mainly affect the low-column density part of the PDF and the mean density. As long as more than 50% of a cloud are covered, the impact on the PDF parameters is negligible. In contrast, the incomplete sampling of the uv plane in interferometric observations leads to uncorrectable distortions of the PDF of the produced maps. An extension of ALMA's capabilities would allow us to recover the high-column density tail of the PDF but we found no way to measure the intermediate and low column density part of the underlying cloud PDF in interferometric observations.