Why wasn cubic convolution interpolation used instead of nearest-neighbor?
During the projection from geographic to UTM coordinates, the mid-high latitudes have much overlap between scenes, and, therefore, much need for interpolation. When a lot of missing values (-32768) are present during a cubic convolution resampling in these areas the result would be less useful than a nearest neighbor resampling. Thus a nearest-neighbor resampling was adopted.