We present three-dimensional astrochemical simulations and synthetic observations of magnetised, turbulent, self-gravitating molecular clouds. We explore various galactic interstellar medium environments, including cosmic-ray ionization rates in the range of 𝜁 CR = 10 −17 -10 −14 s −1 , far-UV intensities in the range of 𝐺 0 = 1-10 3 and metallicities in the range of 𝑍 = 0.1-2 Z . The simulations also probe a range of densities and levels of turbulence, including cases where the gas has undergone recent compression due to cloud-cloud collisions. We examine: i) the column densities of carbon species across the cycle of C , C and CO, along with O , in relation to the H -to-H 2 transition; ii) the velocity-integrated emission of [C ] 158𝜇m, [ 13 C ] 158𝜇m, [C ] 609𝜇m and 370𝜇m, [O ] 63𝜇m and 146𝜇m, and of the first ten 12 CO rotational transitions; iii) the corresponding Spectral Line Energy Distributions; iv) the usage of [C ] and [O ] 63𝜇m to describe the dynamical state of the clouds; v) the behavior of the most commonly used ratios between transitions of CO and [C ]; and vi) the conversion factors for using CO and C as H 2 -gas tracers. We find that enhanced cosmic-ray energy densities enhance all aforementioned line intensities. At low metallicities, the emission of [C ] is well connected with the H 2 column, making it a promising new H 2 tracer in metal-poor environments. The conversion factors of 𝑋 CO and 𝑋 CI depend on metallicity and the cosmic-ray ionization rate, but not on FUV intensity. In the era of ALMA, SOFIA and the forthcoming CCAT-prime telescope, our results can be used to understand better the behaviour of systems in a wide range of galactic and extragalactic environments.