  The surfaces of neutron stars are sources of strongly polarized soft X rays
due to the presence of strong magnetic fields. Radiative transfer mediated by
electron scattering and free-free absorption is central to defining local
surface anisotropy and polarization signatures. Scattering transport is
strongly influenced by the complicated interplay between linear and circular
polarizations. This complexity has been captured in a sophisticated magnetic
Thomson scattering simulation we recently developed to model the outer layers
of fully-ionized atmospheres in such compact objects, heretofore focusing on
case studies of localized surface regions. Yet, the interpretation of observed
intensity pulse profiles and their efficacy in constraining key neutron star
geometry parameters is critically dependent upon adding up emission from
extended surface regions. In this paper, intensity, anisotropy and polarization
characteristics from such extended atmospheres, spanning considerable ranges of
magnetic colatitudes, are determined using our transport simulation. These
constitute a convolution of varied properties of Stokes parameter information
at disparate surface locales with different magnetic field strengths and
directions relative to the local zenith. Our analysis includes full general
relativistic propagation of light from the surface to an observer at infinity.
The array of pulse profiles for intensity and polarization presented highlights
how powerful probes of stellar geometry are possible. Significant
phase-resolved polarization degrees in the range of 10-60% are realized when
summing over a variety of surface field directions. These results provide an
important background for observations to be acquired by NASA's new IXPE X-ray
polarimetry mission.
