  We perform a next to leading order QCD global analysis of nuclear deep
inelastic scattering and Drell-Yan data using the convolution approach to
parameterize nuclear parton densities. We find both a significant improvement
in the agreement with data compared to previous extractions, and substantial
differences in the scale dependence of nuclear effects compared to leading
order analyses.
