ABSTRACT Hyperspectral images are easily contaminated by multi-modal noise in the process of data acquisition. Generally, hyperspectral signatures and noise appear unexpectedly scattered in the spatial-spectral domain. Conversely, in the transform domain, the principle components of an image are gathered in the low-frequency bands, while the details and noise often occur in high-frequency parts. The existing denoising methods usually divide the whole frequency domain into two parts, that is, high frequency and low frequency, which can not make the most of the knowledge contained at different frequencies or scales, especially in the multiscale transform domain. Furthermore, aiming to restore a clean image from the noise-corrupted data, hyperspectral image denoising often amounts to an ill-posed problem, in which reshaping cubic data into matrix form results in the loss of structure information. To address the aforementioned issues, this paper incorporates the reweighted smoothing regularization in the transform domain into the noise-perturbed degradation model to reformulate a novel optimization problem. First, the multiscale curvelet is leveraged to transform the spatial matrix into frequency domain, where the coefficients are carefully designed to weight the components at different scales and orientations. Second, l 1 -norm is imposed on the abundance term to promote sparsity in the transform domain for performance enhancement. Finally, using proximal alternating optimization, an efficient algorithm is well developed to obtain the closed-form solutions to the proposed denoising problem. The experimental results on several synthetic and real datasets, illustrate that the proposed method improves PSNR by at least 1.5 dB and reduces ERGAS error by more than 35% than the state-of-the-art approaches under the multi-modal noise mixed environments, which verifies the validity of regularization terms in hyperspectral images denoising.