“…Solid line: nominal (OH flux of 10 7 cm −2 s −1 ); dash-dot and long dash -double short dash line: eddy mixing coefficient K 0 = 10 2 and 10 3 cm −2 , respectively, section 2.4; long dashed dotted line: CO fixed, no external flux. The data shown for CO is from CIRS (De Kok et al, 2007;Teanby et al, 2010), VIMS (Baines et al, 2006;Bellucci et al, 2009;Fabiano et al, 2017), Herschel SPIRE (Courtin et al, 2011), Herschel PACS (Rengel et al, 2014); for H 2 O is from ISO (Coustenis et al, 1998) (two circles at 400 km, one with small error bar is assuming constant profile, the one with large error bars is the scaled Lara profile), CIRS (Cottini et al, 2012) (lower 3 circles), INMS (Cui et al, 2009b) (upper limits are average of all the flybys, circles are individual detections from single flybys), Herschel (Moreno et al, 2012) (green line, SA profile); for CO 2 is from CIRS (Vinatier et al, 2010), INMS (Cui et al, 2009b); for H 2 CO is from CIRS (Nixon et al, 2010) Early investigations about external sources of oxygen suggested that CO 2 3658 could be produced through a chemical reaction scheme that begins with an 3659 influx of H 2 O into the upper atmosphere from micrometeorite ablation (Lara 3660 et al, 1996;Toublanc et al, 1995;Yung et al, 1984). However, these models 3661 had some difficulty in reproducing the abundance of all three observed species 3662 and a different approach for the origin of O on Titan was needed.…”