Although significant progress has been made in recent years in understanding the composition of the proton's spin from its quark and gluon constituents, a complete picture has yet to emerge. Such information is encoded in spin-dependent parton distribution functions (PDFs) that, as a consequence of being inherently nonperturbative, must be extracted through global QCD analyses of polarized lepton-nucleon and proton-proton collisions. Experiments that measure a final state hadron from these reactions are particularly useful for separating the individual quark and anti-quark polarizations, but require knowledge of parton-to-hadron fragmentation functions (FFs) to describe theoretically.

In this thesis, we present a new approach to global QCD analyses, that were performed recently by the Jefferson Lab Angular Momentum (JAM) Collaboration to determine the spin PDFs and FFs from deep inelastic scattering (DIS), semi-inclusive DIS, and single inclusive electron-positron annihilation observables. While previous global QCD studies typically used a single chi2 minimization procedure, the JAM Collaboration applies a robust Monte Carlo fitting methodology to extract the central values and uncertainties of the relevant distributions. The results from these JAM global QCD analyses, which include a first ever simultaneous fit of the spin PDFs and FFs, resolve a long-standing puzzle regarding the strange quark polarization and provide new information about the proton spin structure.