Pulsar dispersion measures (DMs) have been used to model the electron density of the interstellar medium (ISM) in the Galactic disk as a plane-parallel medium, despite significant scatter in the DM-distance distribution and strong evidence for inhomogeneities in the ISM. We use a sample of pulsars with independent distance measurements to evaluate a model of the local ISM in the thick disk of the Galaxy that incorporates turbulent fluctuations, clumps, and voids in the electron density. The latter two components are required because $\sim 1/3$ of the lines of sight are discrepant from a strictly plane parallel model. A likelihood analysis for smooth components of the model yields a scale height $z_0=1.57^{+0.15}_{-0.14}$ kpc and a mid-plane density $n_0=0.015 \pm 0.001$ cm$^{-3}$. The scatter in the DM-distance distribution is dominated by clumps and voids but receives significant contributions from a broad spectrum of density fluctuations, such as a Kolmogorov spectrum. The model is used to identify lines of sight with outlier values of DM. Three of these pulsars, J1614$-$2230, J1623$-$0908, and J1643$-$1224, lie behind known HII regions, and the electron density model is combined with H$\alpha$ intensity data to constrain the filling factors and other substructure properties of the HII regions (Sh 2-7 and Sh 2-27). Several pulsars also exhibit enhanced DM fluctuations that are likely caused by their lines of sight intersecting the superbubble GSH 238+00+09.

Scintillation arcs have become a powerful tool for exploring scattering in the ionized interstellar medium. There is accumulating evidence that the scattering from many pulsars is extremely anisotropic resulting in highly elongated, linear brightness functions. We present a three-frequency (327~MHz, 432~MHz, 1450~MHz) Arecibo study of scintillation arcs from one nearby, bright, high-velocity pulsar, PSR~B1133+16. We show that a one-dimensional (1D), linear brightness function is in good agreement with the data at all three observing frequencies. We use two methods to explore the broadening of the 1D brightness function $B(\theta)$ as a function of frequency: 1) crosscuts of the forward arc at constant delay and 2) a 1D modeling of $B(\theta)$ using a comparison between model and observed secondary spectrum as a goodness-of-fit metric. Both methods show that the half-power width of $B(\theta)$ deviates from the expected dependence $\propto \nu^{-a}$, where $\nu$ is the observing frequency . Our estimates of $a$ have moderately large uncertainties but imply $a \lesssim1.8$, and so are inconsistent with the expected $a = 2.0$ for plasma refraction or $a = 2.2$ for Kolmogorov turbulence. In addition the shape of $B(\theta)$ cuts off more steeply than predicted for Kolmogorov turbulence. Ultimately, we conclude that the underlying physics of the broadening mechanism remains unexplained. Our results place the scattering screen at a distance that is broadly consistent with an origin at the boundary of the Local Bubble.