Extracted spectra of the investigated sources. Left: full NIRSpec spectra with major ice features labeled. The original spectra before gas emission subtraction are plotted in gray. The spectral region of the strongest HDO absorption, the O–D stretching mode, is indicated in magenta. The spectral region plotted on the right is indicated in gray. Right: zoomed-in view of the spectral region of interest for HDO ice.