Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assumptions, such as the masses of the Kepler-62 planets and the possibility of outer bodies. Significant obliquity variation occurs when the rotational precession frequency overlaps with one or more of the secular orbital frequencies, but most variations are limited to $\lesssim$10$^\circ$. Moderate variations ($\sim$10$^\circ - 20^\circ$) can occur over a broader range of initial obliquities when the relative nodal longitude ($����$) overlaps with the frequency and phase of a given secular mode. However, we find that adding outer gas giants on long period orbits ($\gtrsim$ 1000 days) can produce large ($\sim$60$^\circ$) variations in obliquity if Kepler-62f has a very rapid (4 hr) rotation period. The possibility of giant planets on long period orbits impacts the climate and habitability of Kepler-62f through variations in the latitudinal surface flux, where the timescale for large variation can occur on million year timescales.

15 pages, 17 figures, 5 tables; Revised relative to Referee Reports. Submitted to Astrobiology