Atmospheric CH4 and its stable isotope ratios C-13/C-12 and D/H have been investigated at the Global Atmospheric Watch station Izana, Tenerife (28 degrees N, 16 degrees W, 2370 m above sea level), since late 1996. Every fortnight both spot samples and integral samples were taken, the latter collected continuously over 2-week periods. While spot samples show considerable synoptic-scale variability, the continuous samples clearly define seasonal cycles of CH4, delta(13)C, and FD, with peak-to-peak amplitudes of 30 ppb, 0.2 parts per thousand, and 3.5 parts per thousand, respectively. The measurement of the delta D seasonality is the first ever reported for atmospheric background CH4 and has been made possible by the development of a tunable diode laser based optical Methane Isotopomer Spectrometer (MISOS). The delta D is well in phase with CH4 mixing ratios, and the compact correlation between them allows to derive the average kinetic isotope effect (KIE) of the tropospheric sinks to be 1.23 +/- 0.04 (1 sigma), consistent with recent laboratory measurements of the kinetic isotope effect in the reaction of CH4 + OH. In contrast to delta D, delta(13)C is out of phase to CH4 mixing ratios, clearly indicating that delta(13)C is not only effected by the KIE of the CH4 + OH reaction, but also by seasonally varying source mixtures. Considerable short-term variations are observed in the isotopic composition and mixing ratios of the spot samples. This can be attributed to different origin of air masses arriving at Izana. Significant CH4 enhancement is observed for air masses originating from the North American continent or Europe, of 12.7 +/- 11.7 and 13.4 +/- 0.3 ppb, respectively, while air from the African continent or the North Atlantic is depleted (-11.2 +/- 2.8 and -7.0 +/- 8.2 ppb). Deviations from the mean seasonal CH4 cycle are correlated with significant deviations of delta(13)C and delta D, allowing to estimate the delta(13)C and delta D signatures of major source regions. Furthermore, deviations in CH4 mixing ratio are also clearly correlated with deviations of CO and SFL, mixing ratios. A three-dimensional inverse model is employed in order to assist with the interpretation of observational data. In general, the model shows excellent agreement with the observed mean seasonal cycles of CH4, delta(13)C and delta D, including the observed phase behavior.