The hyperspectral image super-resolution (HSI-SR) problem aims to improve the spatial quality of a low spatial resolution hyperspectral image (LR-HSI) by fusing the LR-HSI combined with the corresponding high spatial resolution multispectral image (HR-MSI). The generated hyperspectral image with high spatial quality, i.e., the target HR-HSI, generally has some fundamental latent properties, e.g., the sparsity and the piecewise smoothness along with the three modes (i.e., width, height, and spectral mode). However, limited works consider both properties in the HSI-SR problem. In this work, a novel unidirectional total variation-based approach is been proposed. On the one hand, we consider the target HR-HSI exhibits both the sparsity and the piecewise smoothness on the three modes, and they can be depicted well by the 1-norm and total variation (TV), respectively. On the other hand, we utilize the classical Tucker decomposition to decompose the target HR-HSI (a 3-mode tensor) as a sparse core tensor multiplied by the dictionary matrices along with the three modes. Especially, we impose the 1-norm on core tensor to characterize the sparsity, and the unidirectional TV on three dictionaries to characterize the piecewise smoothness. The proximal alternating optimization (PAO) scheme and the alternating direction method of multipliers (ADMM) are used to iteratively solve the proposed model. Experiments on three common datasets illustrate the proposed approach has better performance than some current state-of-the-art HSI-SR methods.