An imaging line survey of OMC-1 to OMC-3 Averaged spectra of template regions

Context. Recently, sensitive wide-bandwidth receivers in the millimetre regime have enabled us to combine large spatial and spectral coverage for observations of molecular clouds. The resulting capability to map the distributions of lines from many molecules simul-taneously yields unbiased coverage of the various environments within star-forming regions. Aims. Our aim is to identify the dominant molecular cooling lines and characteristic emission features in the 1.3 mm window of distinct regions in the northern part of the Orion A molecular cloud. By deﬁning and analysing template regions, we also intend to help with the interpretation of observations from more distant sources which cannot be easily spatially resolved. Methods. We analyse an imaging line survey covering the area of OMC-1 to OMC-3 from 200.2 to 281.8 GHz obtained with the PI230 receiver at the APEX telescope. Masks are used to deﬁne regions with distinct properties (e.g. column density or temperature ranges) from which we obtain averaged spectra. Lines of 29 molecular species (55 isotopologues) are ﬁtted for each region to obtain the respective total intensity. Results. We ﬁnd that strong sources like Orion KL have a clear impact on the emission on larger scales. Although not spatially extended, their line emission contributes substantially to spectra averaged over large regions. Conversely, the emission signatures of dense, cold regions like OMC-2 and OMC-3 (e.g. enhanced N 2 H + emission and low HCN/HNC ratio) seem to be difﬁcult to pick up on larger scales, where they are eclipsed by signatures of stronger sources. In all regions, HCO + appears to contribute between 3 and 6% to the total intensity, the most stable value for all bright species. N 2 H + shows the strongest correlation with column density, but not with typical high-density tracers like HCN, HCO + , H 2 CO, or HNC. Common line ratios associated with UV illumination, CN/HNC and CN/HCO + , show ambiguous results on larger scales, suggesting that the identiﬁcation of UV illuminated material may be more challenging. The HCN/HNC ratio may be related to temperature over varying scales.


Introduction
The northern part of the Orion A molecular cloud is one of the most prominent regions of current low-to intermediate-mass star formation, whose close proximity of just ∼400 pc (Menten et al. 2007;Kounkel 2017) enables us to spatially resolve its physically and chemically different regions. Continuum maps show substructures which divide Orion A into morphologically different regions: the bright OMC-1 in the south, with a group of filaments radiating away from its central region (O'Dell et al. 2008 and references therein, e.g. Wiseman & Ho 1998), and the less prominent OMC-2/3 in the north.
OMC-1 hosts star formation and is heavily influenced by intense UV radiation from the young massive Trapezium stars. In addition to the Trapezium, Orion BN/KL (hereafter KL) and Orion South are sites of recent of star formation (O'Dell et al. 2008). Their positions are also indicated in the later discussed Fig. 2. OMC-1 also includes an archetypical photon-dominated region (PDR), the Orion Bar, outflows, and the eponymous hot core (e.g. Masson et al. 1984) containing a rich chemistry of complex molecules, but also more quiescent gas around it. OMC-2 and OMC-3, located northwards of OMC-1, appear to be a continuation of the gas in OMC-1, apart from a shift in velocity (e.g. O'Dell et al. 2008;. A large number of pre-stellar Class 0 and Class I objects following the filamentary structure have been discovered through submillimetre observations (e.g. Chini et al. 1997;Johnstone & Bally 1999;Lis et al. 1998), in addition to a number of brown dwarfs . In contrast to OMC-1, there are no massive O/B stars, resulting in very different conditions with outflows driven by young embedded stars  and references therein, e.g. Yu et al. 1997), but without strong UV radiation.
Investigating the characteristics of these distinct regions helps us to understand how low-to high-mass star formation influences the surrounding material and vice versa. Previous line surveys could be biased in that they have often been focused on a single or few positions (e.g. Sutton et al. 1985;Blake et al. 1987;Tercero et al. 2010Tercero et al. , 2011Esplugues et al. 2013;Johnstone et al. 2003), or spatially extended regions were mapped in selected molecular lines only (e.g. Ungerechts et al. 1997), missing out on information a complete frequency coverage offers.
New opportunities arise with sensitive, wide-bandwidth receivers which combine an extensive frequency coverage with a mapping speed sufficient for large fields of view. This enables us now to obtain a comprehensive and unbiased picture of a molecular cloud. It facilitates the deduction of a variety of physical and chemical conditions using several molecular tracers. These possibilities start being utilised in lower frequency regimes (Pety et al. 2017 and companion papers, e.g. Gratier et al. 2017;Bron et al. 2018 focusing on Orion B or the LEGO project started with Kauffmann et al. 2017). Our data set described in Sect. 2 deals with higher energies, and is thus sensitive to higher critical densities and constitutes a useful complement to these observations. This will help us to better understand the conditions (e.g. column density, temperature, strength of UV illumination) under which different molecules are excited, and in particular identify those molecules predominantly present in very specific environments (e.g. N 2 H + in dense cores; see also Pety et al. 2017). This information is crucial for understanding emission of more distant molecular clouds, including extragalactic sources.
Although CO emission is dominant in the 1.3 mm window, its share of the total intensity (and thus cooling) changes for distinct regions (see also Goldsmith 2001). Comparing the influence of different molecular species to the cooling in Orion A under various conditions will help us to develop templates. These can be used to "reproduce" spectra of other clouds when regions with different physical and chemical conditions are not resolved.

Observations
All observations were carried out with the Atacama Pathfinder Experiment (APEX) 12 m submillimetre telescope (Güsten et al. 2006) using the PI230 receiver operating in the 1.3 mm atmospheric window. We covered the area of OMC-1 to OMC-3 from 200.2 to 281.8 GHz with data collected over several observing periods from October 2015 to November 2018.
The PI230 receiver offers 32 GHz bandwidth per tuning, apportioned into two sidebands and two polarisations, with the two polarisation mixers co-aligned on sky. Eight Fourier Transform Spectrometer (FFTS4G) backends each provide 6 5536 channels for 4 GHz. An overlap of 0.2 GHz between two backends results in 7.8 GHz coverage for each sideband and polarisation per tuning. The sideband rejection is ≈20 dB. The width of the telescope beam changes from around 22 -31 over the 1.3 mm window.
The overall spatial coverage of around 400 × 1900 , corresponding to roughly 0.8 × 3.8 pc at the distance of the Orion Nebula Cluster, was achieved with on-the-fly maps scanned in the x-and y-directions in steps of 8 with a dump time of 0.3 s. The reference position was α 2000 = 5 h 31 m 14.5 s , δ 2000 = −5 • 52 29.0 . The frequency coverage was obtained with 12 overlapping frequency set-ups as shown in Fig. 1. The mean and median system temperatures in main-beam brightness temperature scale (T mb ) determined over all observations are T mean sys ≈ 266 K and T median sys ≈ 239 K, respectively. Mean and median receiver temperatures were T mean rec ≈ 81 K and T median rec ≈ 78 K, respectively.

Calibration and data reduction
All data reduction was done with the GILDAS 1 software package. This included analysis of spectra and maps of the lines' intensity distribution. In the first processing step the data were calibrated and correction factors for the beam efficiencies 1 http://www.iram.fr/IRAMFR/GILDAS Fig. 1. Twelve overlapping frequency set-ups were used to cover the 1.3 mm window, each providing 2 × 7.8 GHz bandwidth. applied separately for each observing period. These correction factors were obtained from Jupiter continuum pointings at different frequencies. The observed peak temperature of the planet was compared to the expected value calculated by the GILDAS ASTRO module, yielding the main beam efficiency. These main beam efficiencies were then fitted by the Ruze formula, facilitating a frequency dependent correction for every observing period, summarised in Table 1. The exact value of the surface accuracy of the dish has only a marginal impact on the data observed in the 1.3 mm window, however. The main beam efficiency varies by 5 % over the observed frequency range. The spectra are also converted from corrected antenna temperature, T * A , to T mb units in this step, which is thus the temperature scale used in this paper.
A substantial amount of data was collected during the commissioning phase of the PI230 receiver, entailing larger uncertainties. However, parts of the maps have been observed repeatedly in different observing periods such that efficiency corrected spectra could be compared and the earlier data validated, suggesting our calibration uncertainty to be around 30% overall.
In the next processing step we removed spikes, artefacts, and spurious lines that only appeared during the early commissioning of the receiver/backend system and affect parts of the data. Furthermore, spectra with T sys ≥ 600 K were rejected. While single channel spikes are readily identified by simple scripts, the artefacts needed to be manually identified and blanked as they showed continuous behaviour towards their edges, making it more challenging to discriminate between them and unaffected channels. A larger number of spectra were checked for each backend and each frequency setting for the different observing days.
To identify and remove spurious lines, we did a first assignment of all observed lines to molecular transitions from the JPL 2 (Pickett et al. 1998) and CDMS 3 (Endres et al. 2016) molecular spectroscopy databases using the WEEDS (Maret et al. 2011) extension for GILDAS. If a line could not be identified, we checked if it was present in all relevant backends and frequency settings. The respective channels in the originating backends were blanked if the line only appeared in some instances.

Final data set
In the following processing step, a first-order baseline was subtracted from each spectrum. The boundaries of the line windows were based on the Orion KL spectrum, where we expect the broadest lines and strongest overlap between them. To simplify the analysis and speed up both imaging and the extraction of spectra, we assembled all data for each backend and each frequency set-up and gridded it using the same template cube with a smoothed beam size of 32 . This results in 96 data cubes (8 backends ×12 frequency set-ups), each containing around 3000 spectra and each of which with 65 536 channels (corresponding to a resolution of 0.07−0.09 kms −1 or 61 kHz). These spectra are intensity corrected, (mostly) artefact-free, baseline subtracted, and arranged on the same grid, and form the basis for all further analyses. The 1.3 mm window gives us access to a variety of molecular species and transitions. In this paper we want to concentrate on 29 species (55 isotopologues), which are listed in Table 2. We can work with our data set both in terms of imaging and spectral analysis. An example of the differing spatial extensions of typical molecular tracers is given in the integrated intensity maps of Fig. 2. The maps are resampled to a resolution of 0.4 km s −1 , resulting in a typical rms noise of around 0.5 K. However, not all species contained in our data set can be mapped as for an individual image pixel their intensity may be too low. Those molecules (like CF + ) are detectable when we average over larger areas. These averaged spectra will be further described in Sects. 3.2 and 3.3.

Analysis
The final data set offers a starting point for varied analyses. In this paper, we want to concentrate on the cooling by different molecular species depending on their environment, focusing on an observational point of view.

Ancillary data
Attributing column density, temperature, or the strength of UV irradiation to different parts of Orion A is done with the help of ancillary data. Dust column density 4 and temperature based on Herschel data are from Guzmán et al. (2015), while C65α emission data (here used to define the dense PDR region) are from Wyrowski et al. (1997). Dust and gas temperature are expected to be coupled for densities ∼3 × 10 4 cm −3 (Galli et al. 2002), and Guzmán et al. (2015) found that for the examined MALT90 clumps (Jackson et al. 2013;Foster et al. 2013Foster et al. , 2011 under 22 K ammonia and dust temperatures agree within ±3 K, while the uncertainties become larger with increasing temperature. We will mainly use the temperature map to differentiate between 2 http://spec.jpl.nasa.gov/ 3 https://cdms.astro.uni-koeln.de/ 4 Converted from g cm −2 to cm −2 under the simplifying assumption that all mass comes from H 2 , 1 g ≈ 2.134 × 10 23 molecules. Notes. Summary of the species and isotopologues used for our analysis. A List of all transitions is found in Appendix A.1 colder (<25 K) and warmer (≥25 K) regions. Maps of gas kinetic temperatures based on NH 3 (Orion A) and H 2 CO (OMC-1 only) observations can be found in Friesen et al. (2017) and Tang et al. (2018), respectively. As the utilised molecular lines require higher column densities to be excited, these maps do not cover all of the ambient material within our maps. Hence we decided on the use of dust-derived temperatures. Reproductions of the column density and dust temperature maps (reprojected using our template data cube) are shown in Fig. 3.

Regions and masks
Masks were used to define regions meeting certain conditions (e.g. column density or temperature ranges). These masks were mostly based on dust column density and temperature maps, but in some cases also selected around chosen coordinates. To this end, the maps were reprojected using the same template cube as for the gridding of our observational data, such that combinations of masks (e.g. high column density concurrent with low temperature) were possible. See Table 3 for a summary of the selected regions and their basic properties. For a better visual idea of the spatial extension of these regions, images of the masks are included in Fig. B.1. Our choice of masks aims to assess the effects of temperature and column density, for example. For the reasons described in Sect. 3.3, KL and Orion South are excluded in some masks, A39, page 3 of 36 A&A 636, A39 (2020) Fig. 2. Selected integrated intensity [5,15] km s −1 maps. Black contours run from 10 to 90% of the maximum value in steps of 10%, while the white contour shows the 3σ level. The maps highlight some of the molecules typically used as tracers: CO for the bulk of molecular gas, HCO + as an indicator of high density, N 2 H + for cold dense gas, CH 3 OH for shocked material, and C 2 H associated with UV irradiation. The velocity range used for the maps does not encompass the whole line width (which differs strongly, especially comparing OMC-1 to OMC-3), but covers all central velocities.
despite matching the column density or temperature criteria. We have masks with only one condition ("high column density" vs. "low column density", "high temperature without KL and Orion South" vs. "low temperature"), but also more restrictive masks, which also define other parameters ("high column density, low temperature" vs. "low column density, low temperature" vs. "low column density, high temperature" vs. "high column density, high temperature without KL and Orion South"). The masks "high column density" and "high column density without KL and Orion South" intend to show how the spatially quite compact emission of KL and Orion South influences the larger scale emission. Additionally, UV irradiation is considered in the masks "H II" and "dense PDR". Choosing this "dense PDR" traced by radio recombination lines, instead of a more "diffuse PDR" traced by C II (Pabst et al. 2019), has the advantage of the edge-on view in Orion. This means we will only probe the actual PDR and are not affected by projection effects. Earlier papers (Turner & Thaddeus 1977) pointed out that the "radical region" may be a chemically distinct environment within OMC-1, although subsequent observations were not conclusive (Greaves & White 1992). We decided to include this region in our analysis, hoping that our unbiased data set might advance its characterisation. Column density and temperature values listed in Table 3 for the "radical region" are thus not its defining features. Lastly, using no mask at all ("all averaged") and employing the results from the other regions, we can evaluate which features or regions dominate on the largest considered scale and which might disappear, indicating that those latter regions may be hard to identify in other, spatially unresolved sources.
The exact values we chose for the column density and temperature thresholds (N = 3.6 × 10 22 cm −2 , T = 25 K) are somewhat arbitrary. We selected them such that they divide Orion A broadly into the denser parts of the filament (N ≥ 3.6 × 10 22 cm −2 ) and the ambient material (N < 3.6 × 10 22 cm −2 ), and distinguish between the historically established regions OMC-1 (T > 25 K) and OMC-2/3 (T < 25 K). Based on these basic distinctions, combinations of masks are added and complemented by the "H II", "dense PDR" and "radical region" masks.
Our masks are thus not devised to divide Orion A into strictly disjointed regions, but to explore the emission of template regions meeting physical or chemical conditions. A listing of all overlaps between the selected regions is given in Table A . Dust derived column density and temperature map, modified (reprojected and cropped, the column density unit converted) from Guzmán et al. (2015). The area around KL and Orion South is saturated and blanked.

Obtaining total intensities
The total intensity of different molecular species for the different regions was determined in three steps. Firstly, the spectrum of the region was extracted from our data cubes using the masks described in Sect. 3.2. All spectra within a mask were combined into a single spectrum, the "regional spectrum" or average. The idea is that local effects should average out in this step when the area is large, while common characteristic features will add up.
Smoothing along the frequency axis to a resolution of around 305 kHz was done to improve the signal-to-noise ratio (S/N). Three of these averaged spectra are shown as examples in Fig. 4.
In the second step, the spectral lines listed in Table A.1 were considered in Gaussian fits. Line-of-sight velocities do not vary significantly within the mapped area, from ∼9 km s −1 for most of OMC-1 (especially when excluding KL and Orion South with slightly lower velocities) and increasing to ∼12 km s −1 towards OMC-3. Different velocity components thus overlap in a spectrum. The list of lines is based on the spectrum obtained with the high column density mask (which includes Orion KL and Orion South), where we expect the greatest number of lines. Omitted are complex molecules with the exception of CH 3 OH, CH 3 CN, and CH 3 CCH. The lines of CH 3 OCH 3 or CH 3 OCHO, for example, are numerous but weak and their contribution to the overall emission as typical hot core molecules should be negligible on larger scales. The main goal is to capture most of the line emission, which in some cases means that not each individual transition is fitted. The automated fitting is more robust when strongly overlapping lines are not fitted separately, but covered in a single Gaussian fit. This favouring of robustness was abandoned when the overlapping lines belonged to different species. In that case, the lines were fitted separately. Overlap between lines becomes an issue whenever KL and Orion South are involved, hence the decision to only include them in two larger regions, as automated fitting is not feasible otherwise.
In the third and final step all lines above 5σ with reasonable widths and velocity were considered to be real detections and the area of their fit added to the total intensity of their respective species.
Noise levels vary within a single spectrum depending on frequency, and between different spectra due to the varied spatial extension of the regions. This influences the possibility to detect very weak lines and will be addressed further in Fig

Results
The averaged total intensity for each species and region is listed in Table 4, together with the detection limit for each region. Owing to space constraints, the formal fit errors are listed separately in Table A.2, but they usually amount to 2−6% for all species and regions. Due to the overall low noise levels and the high number of free parameters, the fit errors for a single line are very small. The actual uncertainties are thus dominated by the calibration uncertainties discussed in Sect. 2.2.
The high column density region has the highest averaged total intensity with ∼1360 K km s −1 . Removing just the region around KL and Orion South has a significant impact on the overall emission and reduces the averaged total intensity by 48%, although the removed area represents only 8% of the high column density region. The regions "high column density, high temperature, without KL and Orion South" and "H II" have comparable averaged total intensities (∼954 K km s −1 and ∼925 K km s −1 ), which are a factor of ∼3 above average (∼319 K km s −1 ). The averaged total CO intensity is also similar for these two regions, their differences lying in the species with less emission. The "dense PDR" and "radical" regions both emit a factor of ∼2 above average, but their share of CO emission differs by around 20%, indicating that their emission profiles (meaning the breakdown of the total emission into contributions from different species) are distinct from each other. The "high temperature without KL and Orion South" region is most similar to the "averaged" region not only in terms of total intensity, but also regarding the most prominent species. The remaining regions, including "high column density, low temperature", emit below average.
For easier comparison, hereafter the data is presented in three complementing ways, each highlighting different aspects. Figure 5 shows the visualisation for the different species with the example of our average data (where no mask was used). The normalised intensity of all species in descending order is shown in Fig. 5a with a logarithmic plot, down to 0.1% of the CO intensity. This order of species is used in the corresponding total intensity plots for all other regions (see Fig. B.2). For the comparison of regions, this representation of the data helps to make shifts in the influence of species more apparent. The pie charts (Figs. 5b and B.4) visualise the percentage of the total intensity originating A39, page 5 of 36 A&A 636, A39 (2020)  from different species, with shares under 2% summed under "other". Plots in this form are more conducive from an observational point of view, as they highlight the most dominant and accessible species. The third visualisation (Fig. B.3) concentrates more on a comparison between species. The plot again shows the percentaged share of each molecule, while the colour bar helps to highlight in which region each species has its largest or smallest share and which species are overall stronger or weaker emitters.
Of the 55 species considered in our analysis, 15 are seen to contribute over 2% each of the total intensity for at least one region. Members of these 15 species account for around 88-94% of line emission in the 1.3 mm window in all cases (or 5.6-53% when excluding CO and its isotopologues). In addition to CO and 13 CO, only HCO + and H 2 CO are prominent in every region, each contributing between 3-6% in all cases.
In terms of averaged total intensity, there is a factor of ∼8 between the lowest ("low column density, low temperature") and highest ("high column density") regions. The high column density region has the lowest share of CO and the highest diversity of species noticeably involved in cooling. SO 2 , with a share of 11.4% being the most important coolant after CO in the high column density region, is no longer relevant when KL and Orion South are removed (its share drops to around 1%; see also Figs. B.4a,B.4b,B.2a and B.2b). SO and CH 3 OH show similar trends. In the high column density region, their total intensities exceed that of the typically strong emitters HCN, HCO + , and H 2 CO, while they lose importance without KL and Orion South.
The regions "high column density, low temperature" and "high column density, high temperature, without KL and Orion South" are parts of "high column density without KL and Orion South" (see also Fig. B.1), coinciding with OMC-2/OMC-3 and OMC-1, respectively. Both in terms of averaged total intensity and as suggested by the pie charts and detailed normalised intensities (Figs. B.4b, B.4c, B.4d, B.2b, B.2c and B.2d), OMC-1 dominates the emission. Higher N 2 H + emission, indicative of OMC-2/OMC-3, is not seen in the "high column density without KL and Orion South" region, which is instead very similar to OMC-1 alone.
The regions "low column density", "low temperature", "low column density, low temperature" and "low column density, high temperature" are similar in the sense of high CO shares (60-67% for CO, 15-21% for 13 CO) and only two to four species over 2% (Figs 4. Spectra obtained from different masks as described in Sects. 3.2 and 3.3. The restriction of the y-axis to 5 K for the middle panel is done to make the weaker lines visible (the CO line has a peak temperature of roughly 57 K), while the zoomed-in boxes illustrate some of the molecular variety.
The regions "high temperature, without KL and Orion South" and "dense PDR" are similar to each other andwith the exception of SO emission and slightly different CO shares -to the "averaged" region (Figs. B.4i, B.4l and 5b). The averaged total intensity is noticeably higher in the "dense PDR" region than in the two others. The differences between regions lie again in the fainter species (Figs. B.2i, B.2l and 5a).
The "H II" region is the only other region, in addition to "high column density", with relevant SO 2 emission and also has higher shares of CH 3 OH and SO, which is the dominant coolant after CO and 13 CO. The "radical region" has the highest shares A39, page 7 of 36 A&A 636, A39 (2020)    Notes. The detection limit is based on the median line width for each region and assumes a 5σ median feature. Real detections below this limit may occur, for example when σ local < σ median or when the line is narrow. Normalised total intensity and pie chart for our averaged data. The order of species in the normalised plot is also used for every other region (see Fig. B.2). The pie chart gives the percentage of the total intensity originating from different species. Shares under 2% are summed under "other".
of C 2 H and CN and is the only other region besides "high column density, low temperature" (dense part of OMC-2/OMC-3) with notable shares of N 2 H + .

Approximating the emission of KL and Orion South
While automated fitting procedures are problematic for high column density and high temperature environments with heavily overlapping lines, we can approximate the emission around KL and Orion South from our existing regions. With the known pixel sizes of "high column density" and "high column density without KL and Orion South" we can gauge the emission of the region and compute its average. The result is not as robust as the others, as more lines are expected for some species which are not accounted for in our routine (especially for SO 2 , SO, and CH 3 OH), and omitted complex organic molecules like CH 3 OCH 3 presumably have nonnegligible impact here. However, the overall results, as presented in Table 5 and Fig. B.5, are in general agreement with Schilke et al. (1997) for the 325-360 GHz frequency regime. The dominant species are in both cases SO 2 , followed by CO, SO, CH 3 OH, and HCN.

Correlations and line ratios on large scales
In the following analysis of correlations and line ratios, the integrated intensity for a given species always refers to the sum of all of its considered transitions in the 1.3 mm window listed in Table A.1. For some species (CO, 13 CO, C 18 O, HCO + , H 13 CO + , HCN, HNC, N 2 H + , and CS) this involves only one transition, but includes several for others (C 2 H, CN, CH 3 OH, H 2 CO, and SO). For the discussion of line ratios, we use simplified quantum numbers for species with one transition, and mark those species with multiple transitions with the letter Σ. Not all energy levels, for example of CH 3 OH, will be populated in all regions, but its integrated intensity still carries meaning as a measure of the cooling in the examined frequency range. Furthermore, neighbouring transitions from complex species like CH 3 OH will strongly overlap in the broad lines from extragalactic sources, thus also inherently limiting the transferability of results concerning single or a few transitions only. Correlations between species do not only point to similarities in their physical and chemical behaviour, but also reduce the number of necessary transitions for gauging the conditions in a molecular cloud. If the emission of two species correlates strongly it is possible, for example, to limit the frequency coverage to one species and save observation time. We limit our analysis to 14 species that are typically used as tracers.
All emission correlates to the first order with column density (see also Fig. 6), which will thus influence the correlation between species. To reduce this effect in our analysis of the 13 regional spectra, we divided the total intensity for each species and region by the respective median column density for the region (see Table 3). This is not as good as a pixel-by-pixel normalisation, but it is more feasible, especially for SO with its numerous partially overlapping lines. Furthermore, the number of data points is small (13 regions), meaning we do not have enough statistical data for a truly reliable correlation coefficient. However, we mainly hope to distinguish between species with strong correlation and those without correlation. Linear correlations were measured with the Pearson correlation coefficient 5 and are shown in Fig. 6, while selected plots are presented in Figs. B.6 and B.7. Using our averaged (and therefore unresolved) regional data adds uncertainty for the interpretation compared to a pixel-by-pixel analysis (e.g. the discussed column density value, which is only an approximation for the whole region). On the other hand, our averaged data enables us to include species like SO and CH 3 OH even in low column density and low-temperature environments.

Correlations and tracers
We expect chemical effects to influence correlation, and also different optical depths of the species. We do not correct for the latter, as our aim is to stay as close to the data an observer might receive from an unresolved source, where the assumptions needed to correct for optical depth may add additional uncertainty.
Strong correlations ( 0.90) are found between the typical high-density tracers HCN, HCO + , H 2 CO, HNC, CS, but also CN. There is generally little spread between the data points in the correlation plots, for example between H 2 CO and HCO + (Fig. B.6a) or H 2 CO and HCN (Fig. B.6b), but the "H II" region with its high integrated intensities (after normalisation) constitutes a more isolated data point.  To a lesser degree, correlations are also found between H 13 CO + and the other high-density tracers (Fig. B.6c). However, the integrated intensity of the former is often an order of magnitude smaller than for the latter, making its observation more challenging.
Not only is CN correlated with C 2 H, which is often associated with UV irradiation (Nagy et al. 2015 and references therein), but it is also correlated with high-density tracers (e.g. Fig. B.6d, e). Another strong correlation is found between the shock tracers (e.g. Bachiller & Pérez Gutiérrez 1997;Sakai et al. 2012;Nagy et al. 2015 and references therein, e.g. Wakelam et al. 2004) SO and CH 3 OH (Fig. B.6f).
For the evaluation of correlations between species and column density in this case we have to factor in that the intensities are not normalised (with the median column density). Hence, it is possible to see two species correlating strongly with column density, but not with each other after normalisation (see C 18 O and N 2 H + ). Of all the considered species, N 2 H + shows the strongest correlation with the median column density (Fig. B.7a), the outlier being the non-detection (<5σ) in the "dense PDR" region, where N 2 H + is most likely expected to be destroyed. We also see a strong correlation between C 18 O and the median column density (Fig. B.7b), but the data points show a wider spread. When considering slope (and intercept) of the linear fit between species and median column density, we see that N 2 H + reacts more sensitively to changes in column density. From the more luminous high-density tracers, HNC may correlate the most with the median column density, in agreement with Pety et al. (2017).
CO and its isotopologues do not show clear correlations with other species (but they do with each other). This may be explained to some degree by optical depth effects, meaning that CO mainly traces the surface of the cloud, while other species probe deeper layers. Especially for C 18 O, which is unlikely to be optically thick, depletion may be an important factor. Highdensity tracers generally profit from higher (volume) densities, while C 18 O may freeze out in these environments if the temperatures are low (as the chance to collide with and stick to dust grains is increased). N 2 H + only correlates with the median column density.

Line ratios
A related important diagnostic, especially for extragalactic observations, but also for molecular clouds in the Milky Way, are line or integrated intensity ratios. They are more reliable tracers of physical or chemical conditions than lines from a single species as calibration errors can cancel out for ratios, depending on the observation technique. These ratios can also be examined with our data.
We will examine some integrated intensity ratios which have also been discussed in Pety et al. (2017), Gratier et al. (2017), or Bron et al. (2018, among others, in the 3 mm window. The considered transitions are thus different, but general trends (e.g. HCN/HNC value increasing with temperature) may still be seen with our data. Selected integrated intensity ratios are compiled in Fig. 7.
Our sample supports the notion from Pety et al. (2017) that the CN/HCN ratio is not a reliable tracer of UV illuminated gas. As the photodissociation of HCN produces CN, one might expect the highest ratio in the "dense PDR" region of our data (CN(Σ)/HCN(3−2) ≈ 0.9). Instead we find it in the "low A39, page 11 of 36 column density" and "radical region" environments (both ≈ 1.3). The lowest ratio is found for the "high column density" region (≈0.5). The CN(1−0)/HCO + (1−0) ratio, discussed in Bron et al. (2018), is suggested to help distinguish UV-illuminated gas from shielded gas, with higher values associated with higher illumination. We find the lowest value (CN(2−1)/HCO + (3−2) ≈ 0.4) for our "low column density, low temperature" region, the highest (≈1.2) in the "radical region", not with the expected "dense PDR" (≈0.7). So at least on our examined large scales and without correction for optical depth, using this ratio in the 1.3 mm window to trace UV illumination seems difficult too.
The ratio HCN/HNC should increase with temperature (Pety et al. 2017 and references therein, e.g. Graninger et al. 2014) due to HNC reacting with H to form HCN at temperatures 30 K. Indeed we find the three lowest values (HCN(3−2)/HNC(3−2) = 1.7−1.8) for our regions associated with the lowest median temperatures, while the highest two values of 5.7−6.7 coincide with the highest temperatures, found in our case in the "H II" and "dense PDR" regions. Larger deviations from this trend are found for the "radical region", for example (see also Fig. B.7c).
In their PCA of Orion B data in the 3 mm window, Gratier et al. (2017) found that a higher N 2 H + (J = 1−0, F1 = 2−1, F = 3 − 2, )/CH 3 OH(J = 2−1, K = 0−0, A + ) ratio possibly highlights the chemistry of the densest cores. This cannot be meaningfully examined for the 1.3 mm window with our spatially unresolved data, and we do not see a correlation with column density. While the emission of both species is higher in regions of enhanced column density in our data, their ratio is not. The highest ratios are not found in the regions encompassing the dense filament, but instead in the "low temperature" and "low column density, low temperature" regions (N 2 H + (3−2)/CH 3 OH(Σ) ≈ 3.8−4.1), while the spectra from the high column density regions have notably lower values (≈0.1 for the "high column density", ≈2.4 for the "high column density, low temperature", and ≈0.4 for the "high column density, high temperature, without KL and Orion South" region). However, we see a potential correlation of N 2 H + (3−2)/CH 3 OH(Σ) with temperature (Fig. B.7d).

Line luminosities
For the interpretation of emission from unresolved and/or extragalactic sources, information on line luminosities [K km s −1 pc 2 ] is important. Regions with low overall emission may still contribute considerably to the emission of some species if these regions are spatially extended. Conversely, strong but compact emission may be diluted on larger scales. Figure 8 lists the total luminosity (sum over all considered species in the 1.3 mm window) and the values for selected species for all our examined regions, including the approximations for KL and Orion South (see also Sect. 4.1).
Despite its overall lower emission, the extended "low column density" region is the source of ∼75% of the CO, 13 CO, and C 18 O emission, while ∼25% can be allotted to the "high column density" region. For the high-density tracers HCO + and HNC, but also N 2 H + , the allocation between low and high column density is around 50%/50%. For HCN and CS it is less of an even split between the two column density regimes and ∼65% originate from the high column density region.
A lot of the emission from the high column density region can be attributed to the environment of KL and Orion South: with a size of <1% of the total examined area, it emits ∼20% of the HCO + and HNC, and ∼30% of HCN and CS. When considering all emission in the 1.3 mm window, KL and Orion South are the source of roughly 25%.

Discussion
By averaging over the area, we highlight the individual emission profile of the regional spectra, including lower intensity species, while ignoring their actual spatial extent (which here mainly influences the respective noise level). This approach might help to characterise sources for which the spatial resolution is not good enough to differentiate between distinct regions. It is also a different approach compared to the analysis of the Orion B data set (Pety et al. 2017;Gratier et al. 2017;Bron et al. 2018) in the 3 mm window, where the focus is more on selected, typically more luminous transitions, not necessarily on a complete inventory of species. We discuss our findings further in this section and compare them with results from other authors. We find that temperature has a significant impact on the total intensities of our selected regions; considering column density alone is not sufficient. This is illustrated most obviously in the cases of the "high column density, low temperature" and "H II" regions. The former has higher column densities but low overall emission, while the latter has lower median column density but a factor of ≈2.6 higher total intensity (total intensity summed over all considered species). We mostly consider temperatures and column densities, but it also seems instructive to keep the main feedback processes in different regions in mind. While the column densities in OMC-2/3 are similar to values found in parts of OMC-1, the energy input in the former is mainly driven by outflows from low-mass protostars, while in the latter it is shared between outflows and feedback from the H II region, as discussed in Berné et al. (2014) and inferred from CO and 13 CO (2−1) emission. The kinetic energy in KL was found to be dominated by feedback of massive protostars (outflows, jets, explosive motion). In contrast, the Bar exhibits very little outflow activity; its kinetic energy is mainly caused by the expanding H II region.
We find very similar total CO intensities for the "high column density" and "dense PDR" regions (∼393 and ∼392 K km s −1 ), but their total intensities summed over all species vary by more than 50% (∼1360 and ∼603 K km s −1 ). While the different column densities and temperatures in the two regions still result in similar CO intensities, other species react very differently and are much more enhanced in the "high column density" region. This might again be partly due to opacity effects and CO only tracing the outer layers.
The median column densities for the "dense PDR" and "H II" regions are about a factor of 3 lower than for "high column density" and "high column density, high temperature without KL and Orion South", but their median temperatures are the highest for all considered regions. This might explain their high emission despite overall lower column densities.

Correlations
For the further examined most luminous species, strong correlations are typically found between high-density tracers, excluding N 2 H + . This may be related to the optical depth effects mentioned before, but also to the influence of temperature, as shown in the case of the warm "dense PDR" region, where N 2 H + does not emit over 5σ despite enhanced column densities. N 2 H + shows no correlation with any other prevalent species in the 1.3 mm window, but instead with column density. If N 2 H + data is not available, HNC might overall be a better tracer of column density than HCN or HCO + , for example. HNC shows similar intensities to N 2 H + , however, so it should not be a question of observation time. That HCN does not exclusively trace dense gas was also shown in Kauffmann et al. (2017) for the 1−0 transition in Orion A, where HCN was found to trace lower densities ∼10 3 cm −3 in cold sections of the cloud. Additionally, they found that the cold dense gas emits too little HCN to explain the luminosities observed in extragalactic sources. This is consistent with the below average emission we see for the OMC-2/3 region in our data ("high column density, low temperature").
Different papers discuss additional mechanisms which may excite "dense gas" tracers. As elaborated in Goldsmith & Kauffmann (2017), electron excitation may be important for high-dipole moment molecules in regions where the fraction of ionised carbon is significant. The low-J transitions of HCN (but also HCO + , CN, and CS) could thus be observed in lower density environments and may not qualify as indicators of high density. Another mechanism could be radiative trapping (Shirley 2015;Pety et al. 2017). It is argued that fundamental lines of A39, page 13 of 36 A&A 636, A39 (2020) HCN, HNC, and HCO + could be excited in regions well below their critical density as the latter is computed assuming optically thin emission only. Both mechanisms could explain our observations. Particularly in OMC-1, where C + is abundant (Pabst et al. 2019), the emission of HCN, for example, may originate in part from electron excitation. Additionally, some transitions may have a high optical depth, increasing contributions from photontrapping. An alternative approach to the detection of low-lying rotational lines is discussed in Liszt & Pety (2016). They discuss the observability of transitions for HCO + , HNC, and CS in regions where the density is far below the critical density. The authors find, in the limit of weak collisional excitation, that there is a column density (not a volume density) that will produce a given output.
We find that the typical UV tracers CN and C 2 H correlate not only with each other, but also strongly with high-density tracers. This means that tracing UV illumination on large scales could be challenging, as enhanced CN or C 2 H emission might not necessarily be indicative of higher UV illumination.

Emission on larger scales
The stark difference in the emission profile of KL and Orion South (see Sect. 4.1) compared to the other regions was also found for some species in Ungerechts et al. (1997), where the integrated intensity maps of SO or HC 3 N showed distributions strongly peaked towards KL. While we find that line emission from KL and Orion South has a noticeable impact on spectra averaged over larger scales (compare Figs. B.4a and b, but also appreciate that SiO and SO 2 are detectable in the averaged spectrum despite not being detected for the majority of regions), this influence of high column density regions on the averaged spectrum was not seen by Watanabe et al. (2017) for the spectrum of W51 in the 3 mm window. Their spatial coverage of 39 pc × 39 pc is significantly larger than for our data set, so emission from high-density regions is expected to be smeared out more, such that they find a dominance of the quiescent material over the averaged spectrum.
In their analysis of a line survey of the central parts of the starburst galaxy M82 in the 1.3 and 2 mm window, Aladro et al. (2011) found that the physical processes are dominated by PDRs. Arguing that feedback from young OB stars leaves an imprint on molecular composition in the form of an overabundance of CO + , HCO, c-C 3 H 2 , and CH 3 CCH, for example, they find M82 to match these criteria. While we can confirm the high c-C 3 H 2 emission in our dense PDR region, CH 3 CCH is below the 5σ limit there. This is also interesting because their spatial resolution varies between 158 pc and 333 pc, which is 2 to 3 orders of magnitude higher than our averaged spectra, but we observe CH 3 CCH emission as quite compact and thus strongly diluted on larger scales. The noise levels are comparable in both studies, ∼3−8 mK for M82, around 8−20 mK median noise for the different regions in our Orion data.

Line ratios on larger scales
Line ratios seem to be more ambiguous on large scales, where both CN(Σ)/HCN(3−2) and CN(Σ)/HCO + (3−2) do not clearly highlight regions with enhanced UV irradiation. We can confirm a correlation between HCN(3−2)/HNC(3−2) and temperature. On even larger scales, however, as examined by Meier & Turner (2005), among others, for the nuclear region of IC 342 in the 3 mm window (∼50 pc resolution), the ratio seems to be fairly constant and not related to (kinetic) temperature. Their value of HCN(1−0)/HNC(1−0) ≈ 1−2 is found in our data set for those regions associated with low temperature (see Fig. B.7c). The considered transitions are different to ours, but both species still have comparable upper energies, such that a comparison seems meaningful.
The correlation of N 2 H + (3−2)/CH 3 OH(Σ) with temperature we found on larger scales might be explained both with the temperature sensitivity of N 2 H + and the association of CH 3 OH with shocks. If temperatures rise, more CO can enter the gas phase and subsequently destroy N 2 H + , while more CH 3 OH evaporates.

Summary
We have conducted an imaging line survey of OMC-1 to OMC-3 from 200.2 to 281.8 GHz and examined the emission of distinct regions. These were selected to represent regimes of low or high column density and differing temperature, but also to cover the influence of UV irradiation. By contrasting their emission with each other, we aim to provide templates for the interpretation of other more distant, spatially unresolved sources. Comparing spectra from these unresolved sources with our different templates might reveal similarities that can help to further characterise the distant object. The transitions of the 29 species (55 isotopologues) listed in Table A.1 were considered for the analysis. The integrated intensity of a given species is considered here to be the sum over all of their respective transitions. From our analysis we would like to highlight the following results: 1. Line emission from an Orion KL-like source can contribute significantly to spectra averaged over larger regions, both in terms of averaged total intensity and chemical diversity. In terms of line luminosities, KL and Orion South contribute around 25% of all emission in the 1.3 mm window in the area of OMC-1 to OMC-3.
2. Regions like OMC-2/3, with pre-stellar Class 0 and Class I objects and enhanced column density but low temperature, have a much lower total intensity. Their signatures (e.g. high N 2 H + emission coincident with a low HCN(3−2)/HNC(3−2) ratio) would be difficult to pick up in a non-resolved source. 3. While the contribution of CO to the share of the total intensity can vary more strongly in the examined 1.3 mm window, HCO + (3−2) contributes ∼3−6% in all cases (average ≈4%). This seems to hold true even in the case of the emission around KL and Orion South. 4. The emissions of the high-density tracers HCN, HCO + , H 2 CO, HNC, CS, but also CN are strongly correlated with each other, but not with N 2 H + . Of all the examined species, N 2 H + shows the strongest correlation with column density. If N 2 H + observations are not available, HNC seems to trace the column density most reliably. 5. Around 50% of the line luminosity of HCO + and HNC in the 1.3 mm window comes from lower column density material, for CS and HCN ∼35%. 6. The ratios HCN(3−2)/HNC(3−2) and N 2 H + (3−2)/ CH 3 OH(Σ) may be related to temperature. 7. Identifying UV illuminated material on large scales seems to be challenging, as both CN(Σ)/HCN(3−2) and CN(Σ)/HCO + (3−2) show ambiguous results. A&A 636, A39 (2020) Appendix A: Additional tables     Notes. To be read from left to right, e.g. ∼39% of the "high column density" region overlaps with the "low temperature" region. Conversely, ∼9% of the "low temperature" region overlaps with the "high column density" region. The majority of these overlaps are discussed in this paper as they correspond to regions themselves (e.g. the overlap between the "high column density" and "low temperature" regions is identical to the region "high column density, low temperature"