We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or M92A, has a pulse frequency of 316.5~Hz (3.16~ms) and a dispersion measure of 35.45 pc cm$^{-3}$. The pulsar is a member of a binary system with an orbital period of 0.20~days around a low-mass companion which has a median mass of $\sim$0.18~\Ms. From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ~5000~s in the orbital phase range 0.1--0.5. The other lasted for ~500~s and occurred between 1000--2000~s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a ``red-back'' system with a low-mass main sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing.

The Green Bank North Celestial Cap (GBNCC) pulsar survey will cover the entire northern sky ($\delta > -40^\circ$) at 350 MHz, and is one of the most uniform and sensitive all-sky pulsar surveys to date. We have created a pipeline to re-analyze GBNCC survey data to take a 350MHz census of all pulsars detected by the survey, regardless of their discovery survey. Of the 1413 pulsars in the survey region, we were able to recover 661. For these we present measured signal-to-noise ratios (S/N), flux densities, pulse widths, profiles, and where appropriate, refined dispersion measurements (647 out of 661) and new or improved spectral indices (276 out of 661 total, 15 new, 261 improved). Detection scans for several hundred sources were reanalyzed in order to inspect pulsars' single pulse behavior and 223 were confirmed to exhibit evidence of nulling. With a detailed analysis of measured and expected S/N values and the evolving radio frequency interference environment at 350MHz, we assess the GBNCC survey's sensitivity as a function of spin period, dispersion measure, and sky position. We find the sky-averaged limiting flux density of the survey to be 0.74mJy. Combining this analysis with PsrPopPy pulsar population simulations, we predict 60/5 non-recycled/millisecond pulsar discoveries in the survey's remaining 21,000 pointings, and we begin to place constraints on population model parameters.

We present a detailed analysis of the International Pulsar Timing Array (IPTA) Second Mock Data Challenge. We tested our analysis methods using the open datasets, and then analyzed the closed datasets. In both the open and the closed datasets, we were able to detect some, but not all, of the injected gravitational wave signals. This work presents two search cases that are not well explored in the pulsar timing array (PTA) literature: a simultaneous search for a stochastic GW background and an individual loud super-massive black hole binary (SMBHB) and a simultaneous search for two SMBHB sources. While we have constructed a cohesive framework for performing these GW searches, our analyses required fine-tuning of the sampling method used in order to appropriately converge. Given the nature of real PTA data in which multiple sources will be present in data, improved techniques will be required in the future to accurately detect and characterize these GW signals.