Seventeen accreting neutron star pulsars, mostly high mass X-ray binaries with half of them Be-type transients, are known to exhibit Cyclotron Resonance Scattering Features (CRSFs) in their X-ray spectra, with characteristic line energies from 10 to
60 keV. To date about two thirds of them, plus a few similar systems without known CRSFs, have been observed with Suzaku. We present an overview of results from these observations, including the discovery of a CRSF in the transient 1A 1118-61 and pulse phase resolved spectroscopy of GX 301-2. These observations allow for the determination of cyclotron line parameters to an unprecedented degree of accuracy within a moderate amount of observing time. This is important since these parameters vary - e.g., with orbital phase, pulse phase, or luminosity - depending on the geometry of the magnetic field of the pulsar and the properties of the accretion column at the magnetic poles. We briefly introduce a spectral model for CRSFs that is currently being developed and that for the first time is based on these physical properties. In addition to cyclotron line measurements, selected highlights from the Suzaku analyses include dip and flare studies, e.g., of 4U 1907+09 and Vela X-1, which show clumpy wind effects (like partial absorption and/or a decrease in the mass accretion rate supplied by the wind) and may also display magnetospheric gating effects.
We present results from two observations of the wind-accreting X-ray pulsar 4U 1907+09 using the Suzaku observatory. The broadband time-averaged spectrum allows us to examine the continuum emission of the source and the cyclotron resonance scattering
feature at ~19 keV. Additionally, using the narrow CCD response of Suzaku near 6 keV allows us to study in detail the Fe K bandpass and to quantify the Fe K beta line for this source for the first time. The source is absorbed by fully-covering material along the line of sight with a column density of NH ~2e22 /cm^2, consistent with a wind accreting geometry, and a high Fe abundance (~3-4 x solar). Time and phase-resolved analyses allow us to study variations in the source spectrum. In particular, dips found in the 2006 observation which are consistent with earlier observations occur in the hard X-ray bandpass, implying a variation of the whole continuum rather than occultation by intervening material, while a dip near the end of the 2007 observation occurs mainly in the lower energies implying an increase in NH along the line of sight, perhaps indicating clumpiness in the stellar wind.
