SummaryNuclear magnetic resonance (NMR) relaxation responses in porous media provide a sensitive probe of the microstructure and yet are influenced by a number of factors which are not easily detangled. Low-field T2 transverse relaxation measurements can be carried out quickly and are frequently used as pore size distributions, while adding T1 longitudinal relaxation measurements provides additional insights into surface properties and fluid content. Here we present an inverse solution workflow extracting related intrinsic physical parameters of the system by fitting experiment and numerical simulation(s). An efficient NMR forward solver for the simultaneous calculation of T1 and T2 responses is introduced, which honors existing inequality relationships between T1 and T2 parameters. We propose a Bayesian optimization approach that jointly identifies T1- and T2-related properties satisfying physical constraints by simultaneously fitting T1 and T2 experiments to simulations. This dual-task inverse solution workflow (DT-ISW) identifies the solution by minimizing the sum of the L2 norm of the fitting residuals of both T1 and T2 distributions into a single objective and jointly models the two highly correlated objectives with high efficiency using the vector-valued Gaussian process (GP) kernel for transfer learning. A multimodal search strategy is used to identify nonunique solution sets of the problem. The workflow is demonstrated on Bentheimer sandstone, identifying five intrinsic physical parameters. The performance of the joint DT-ISW (DT-ISW-J) is compared to a sequential DT-ISW (DT-ISW-S) approach as well as an independent single-task ISW (ST-ISW) of the T1 and T2 responses. Both dual-task versions converge more than two times faster than the single-task version. DT-ISW-J equally minimizes the L2 norm of T1 and T2 fitting residuals whereas DT-ISW-S only preferentially minimizes the objective assigned higher importance. A Pareto optimal solution (POS) is provided to allow operators to subjectively balance the preference of T1 and T2 data fits for the slightly conflicting objectives. The ability to extract five intrinsic physical parameters simultaneously provides new techniques for tracking wettability alteration and assessing the influence of clay amount and distribution on petrophysical property estimates.