Compositional changes can occur in the osteochondral junction (OCJ) during the early stages and progressive disease evolution of knee osteoarthritis (OA). However, conventional magnetic resonance imaging (MRI) sequences are not able to image these regions efficiently because of the OCJ region's rapid signal decay. The development of new sequences able to image and quantify OCJ region is therefore highly desirable. We developed a comprehensive ultrashort echo time (UTE) MRI protocol for quantitative assessment of OCJ region in the knee joint, including UTE variable flip angle technique for T1 mapping, UTE magnetization transfer (UTE-MT) modeling for macromolecular proton fraction (MMF) mapping, UTE adiabatic T1ρ (UTE-AdiabT1ρ) sequence for T1ρ mapping, and multi-echo UTE sequence for T2* mapping. B1 mapping based on the UTE actual flip angle technique was utilized for B1 correction in T1, MMF, and T1ρ measurements. Ten normal and one abnormal cadaveric human knee joints were scanned on a 3T clinical MRI scanner to investigate the feasibility of OCJ imaging using the proposed protocol. Volumetric T1, MMF, T1ρ, and T2* maps of the OCJ, as well as the superficial and full-thickness cartilage regions, were successfully produced using the quantitative UTE imaging protocol. Significantly lower T1, T1ρ, and T2* relaxation times were observed in the OCJ region compared with those observed in both the superficial and full-thickness cartilage regions, whereas MMF showed significantly higher values in the OCJ region. In addition, all four UTE biomarkers showed substantial differences in the OCJ region between normal and abnormal knees. These results indicate that the newly developed 3D quantitative UTE imaging techniques are feasible for T1, MMF, T1ρ, and T2* mapping of knee OCJ, representative of a promising approach for the evaluation of compositional changes in early knee OA.