Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced for the Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. We report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

We conduct a comprehensive study of the effects of incorporating galaxy morphology information in photometric redshift estimation. Using machine learning methods, we assess the changes in the scatter and catastrophic outlier fraction of photometric redshifts when galaxy size, ellipticity, S\'{e}rsic index and surface brightness are included in training on galaxy samples from the SDSS and the CFHT Stripe-82 Survey (CS82). We show that by adding galaxy morphological parameters to full $ugriz$ photometry, only mild improvements are obtained, while the gains are substantial in cases where fewer passbands are available. For instance, the combination of $grz$ photometry and morphological parameters almost fully recovers the metrics of $5$-band photometric redshifts. We demonstrate that with morphology it is possible to determine useful redshift distribution $N(z)$ of galaxy samples without any colour information. We also find that the inclusion of quasar redshifts and associated object sizes in training improves the quality of photometric redshift catalogues, compensating for the lack of a good star-galaxy separator. We further show that morphological information can mitigate biases and scatter due to bad photometry. As an application, we derive both point estimates and posterior distributions of redshifts for the official CS82 catalogue, training on morphology and SDSS Stripe-82 $ugriz$ bands when available. Our redshifts yield a 68th percentile error of $0.058(1+z)$, and a catastrophic outlier fraction of $5.2$ per cent. We further include a deep extension trained on morphology and single $i$-band CS82 photometry.