diff --git a/specreduce/calibration_data.py b/specreduce/calibration_data.py
index 20052ed..f90a4a4 100644
--- a/specreduce/calibration_data.py
+++ b/specreduce/calibration_data.py
@@ -526,7 +526,7 @@ class AtmosphericTransmission(AtmosphericExtinction):
         two columns, wavelength and transmission (unscaled dimensionless). If
         this isn't provided, a model is built from a pre-calculated table of values
         from 0.9 to 5.6 microns. The values were generated by the ATRAN model,
-        https://atran.arc.nasa.gov/cgi-bin/atran/atran.cgi (Lord, S. D., 1992, NASA
+        https://ntrs.nasa.gov/citations/19930010877 (Lord, S. D., 1992, NASA
         Technical Memorandum 103957). The extinction is given as a linear transmission
         fraction at an airmass of 1 and 1 mm of precipitable water.
 
diff --git a/specreduce/extract.py b/specreduce/extract.py
index f3dab4b..cabe856 100644
--- a/specreduce/extract.py
+++ b/specreduce/extract.py
@@ -7,6 +7,7 @@
 from astropy import units as u
 from astropy.modeling import Model, models, fitting
 from astropy.nddata import NDData, VarianceUncertainty
+from scipy.integrate import trapezoid
 from scipy.interpolate import RectBivariateSpline
 from specutils import Spectrum1D
 
@@ -794,7 +795,7 @@ def __call__(self, image=None, trace_object=None,
             else:  # interpolated_profile
                 fitted_col = interp_spatial_prof(col_pix, xd_pixels)
                 kernel_vals.append(fitted_col)
-                norms.append(np.trapz(fitted_col, dx=1)[0])
+                norms.append(trapezoid(fitted_col, dx=1)[0])
 
         # transform fit-specific information
         kernel_vals = np.vstack(kernel_vals).T