Skip to content

Stellar Reddening VAC

Overview

Y1 spectra and catalog data for creating the DESI dust reddening map of Zhou et al. (2024).

While the paper uses the first two years of DESI data, only the Y1 portion of the data is released here. Also note that while both the MWS VAC and this VAC use RVSpecfit results and are based on the same DESI data, there are some key differences: 1) the RVSpecfit model spectra in this VAC are the zero-extinction spectra (thus suitable for estimating the “intrinsic” stellar colors), whereas the MWS VAC model spectra include an additional multiplicative term to fit the smooth-varying component of the observed spectra; 2) the DR1 MWS VAC uses a different (newer) version of rvspecfit than the version used here. For general purposes (other than for studying dust), the MWS VAC should be used instead of this VAC.

Code for creating the dust map: https://github.com/rongpu/desidust

Data Access

Data URL: https://data.desi.lbl.gov/public/dr1/vac/dr1/stellar-reddening/

NERSC access for DESI collaborators:

/global/cfs/cdirs/desi/public/dr1/vac/dr1/stellar-reddening/

Spectra

The directory spectra/ contains the observed spectra (input for RVSpecfit) from DESI Y1 observations and the corresponding RVSpecfit model spectra. It includes standard stars (STD_FAINT and STD_BRIGHT) from the dark and bright programs, MWS stars (MWS_MAIN_BLUE, MWS_MAIN_RED and MWS_BROAD) from the bright program, and standard and MWS stars (GAIA_STD_FAINT, GAIA_STD_BRIGHT, BACKUP_BRIGHT, BACKUP_FAINT) from the backup program. This selection is purely based on targeting, and it includes the spectra of stars that are ultimately dropped by the quality cuts. They follow the same healpix-based directory structure and file naming convention as the DESI spectra. Example file: /stellar-reddening/v1.0/spectra/main/bright/0/0/rvspecfit-main-bright-0.fits

The first 15 HDUs (FIBERMAP, _WAVELENGTH, _FLUX, _IVAR, _MASK, SCORES, REDSHIFTS) are copied from the DESI pipeline output, specifically the iron reduction of the coadd files and the redrock files.

The RVTAB HDU contains the RVSpecfit parameters The last two HDUs, SPECMOD and WAVELENGTH, containing the RVSpecfit model spectra and its wavelength grid, are used to obtain the zero-extinction model photometry of each star. The model spectrum SPECMOD, in \(f_{\lambda}\) with arbitrary absolute normalization, is obtained by evaluating the interpolation of PHOENIX templates on the best-fit stellar atmospheric parameters.

HDUs

HDU# EXTNAME Type Contents
HDU1 FIBERMAP BINTABLE fibermap table
HDU2 B_WAVELENGTH IMAGE Wavelength array of b-channel spectra
HDU3 B_FLUX IMAGE Flux of b-channel spectra
HDU4 B_IVAR IMAGE Inverse variance of b-channel spectra
HDU5 B_MASK IMAGE Mask of b-channel spectra
HDU6 R_WAVELENGTH IMAGE Wavelength array of r-channel spectra
HDU7 R_FLUX IMAGE Flux of r-channel spectra
HDU8 R_IVAR IMAGE Inverse variance of r-channel spectra
HDU9 R_MASK IMAGE Mask of r-channel spectra
HDU10 Z_WAVELENGTH IMAGE Wavelength array of z-channel spectra
HDU11 Z_FLUX IMAGE Flux of z-channel spectra
HDU12 Z_IVAR IMAGE Inverse variance of z-channel spectra
HDU13 Z_MASK IMAGE Mask of z-channel spectra
HDU14 SCORES BINTABLE Scores table
HDU15 REDSHIFTS BINTABLE Table with redshifts and spectral classifications
HDU16 RVTAB BINTABLE RVSpecfit results
HDU17 WAVELENGTH IMAGE Wavelength array of SPECMOD spectra, in angstrom
HDU18 SPECMOD IMAGE RVSpecfit model spectra of the stars obtained by evaluating the interpolation of PHOENIX templates on the best-fit stellar atmospheric parameters. The spectra are in \(f_{\lambda}\), with arbitrary absolute normalization

Catalog

The catalog file desi_dust_stellar_catalog_y1.fits contains stars used in the dust map (after the quality cuts are applied). The columns (listed below) include those copied from the DESI output, RVSpecfit output, Gaia DR3, Legacy Surveys DR9, and model spectrophotometry and reddening measurements.

Data model

EXTNAME = STELLAR_REDDENING

Name Type Units Description
Columns PROGRAM to DELTACHI2 are copied from the DESI redrock output
PROGRAM str DESI program type - BRIGHT, DARK, BACKUP, OTHER
TARGETID_DESI int64 Unique DESI target ID assigned by DESI
COADD_FIBERSTATUS int32 bitwise-AND of input FIBERSTATUS
TARGET_RA float64 deg Barycentric right ascension in ICRS
TARGET_DEC float64 deg Barycentric declination in ICRS
FA_TARGET int64 Targeting bit internally used by fiberassign (linked with FA_TYPE)
FA_TYPE uint8 Fiberassign internal target type (science, standard, sky, safe, suppsky)
OBJTYPE str Object type: TGT, SKY, NON, BAD
SUBPRIORITY float64 Random subpriority [0-1) to break assignment ties
OBSCONDITIONS int32 Bitmask of allowed observing conditions
PHOTSYS_DESI str Photometric system for DESI targeting: ‘N’ for the MzLS/BASS photometric system, ‘S’ for DECaLS
PRIORITY_INIT int64 Target initial priority from target selection bitmasks and OBSCONDITIONS
NUMOBS_INIT int64 Initial number of observations for target calculated across target selection bitmasks and OBSCONDITIONS
DESI_TARGET int64 DESI (dark time program) target selection bitmask
BGS_TARGET int64 BGS (Bright Galaxy Survey) target selection bitmask
MWS_TARGET int64 Milky Way Survey targeting bits
SCND_TARGET int64 Target selection bitmask for secondary programs
PLATE_RA float64 deg Barycentric Right Ascension in ICRS to be used by PlateMaker
PLATE_DEC float64 deg Barycentric Declination in ICRS to be used by PlateMaker
COADD_NUMEXP int16 Number of exposures in coadd
COADD_EXPTIME float32 s Summed exposure time for coadd
COADD_NUMNIGHT int16 Number of nights in coadd
COADD_NUMTILE int16 Number of tiles in coadd
MEAN_DELTA_X float32 mm Mean (over exposures) fiber difference requested - actual CS5 X location on focal plane
RMS_DELTA_X float32 mm RMS (over exposures) of the fiber difference between measured and requested CS5 X location on focal plane
MEAN_DELTA_Y float32 mm Mean (over exposures) fiber difference requested - actual CS5 Y location on focal plane
RMS_DELTA_Y float32 mm RMS (over exposures) of the fiber difference between measured and requested CS5 Y location on focal plane
MEAN_FIBER_RA float64 deg Mean (over exposures) RA of actual fiber position
STD_FIBER_RA float32 arcsec Standard deviation (over exposures) of RA of actual fiber position
MEAN_FIBER_DEC float64 deg Mean (over exposures) DEC of actual fiber position
STD_FIBER_DEC float32 arcsec Standard deviation (over exposures) of DEC of actual fiber position
MEAN_PSF_TO_FIBER_SPECFLUX float32 Mean of input exposures fraction of light from point-like source captured by 1.5 arcsec diameter fiber given atmospheric seeing
TSNR2_GPBDARK float32 template (S/N)^2 for dark targets in guider pass band
TSNR2_ELG float32 ELG template (S/N)^2 summed over B,R,Z
TSNR2_GPBBRIGHT float32 template (S/N)^2 for bright targets in guider pass band
TSNR2_LYA float32 LYA template (S/N)^2 summed over B,R,Z
TSNR2_BGS float32 BGS template (S/N)^2 summed over B,R,Z
TSNR2_GPBBACKUP float32 template (S/N)^2 for backup targets in guider pass band
TSNR2_QSO float32 QSO template (S/N)^2 summed over B,R,Z
TSNR2_LRG float32 LRG template (S/N)^2 summed over B,R,Z
Z float64 Redshift measured by Redrock
ZERR float64 Redshift error from redrock
ZWARN int64 Redshift warning bitmask from Redrock
CHI2 float64 Best fit chi squared
NPIXELS int64 Number of unmasked pixels contributing to the Redrock fit
SPECTYPE str Spectral type of Redrock best fit template (e.g. GALAXY, QSO, STAR)
SUBTYPE str Spectral subtype
DELTACHI2 float64 chi2 difference between first- and second-best redrock template fits
Columns VRAD to SUCCESS are copied from the RVSpecfit output
VRAD float64 km s^-1 Radial velocity
VRAD_ERR float64 km s^-1 Radial velocity error
VRAD_SKEW float64 Radial velocity posterior skewness
VRAD_KURT float64 Radial velocity posterior kurtosis
LOGG float64 Log of surface gravity
TEFF float64 K Effective temperature
ALPHAFE float64 [alpha/Fe] from template fitting
FEH float64 [Fe/H] from template fitting
LOGG_ERR float64 Log of surface gravity uncertainty
TEFF_ERR float64 K Effective temperature uncertainty
ALPHAFE_ERR float64 [alpha/Fe] uncertainty from template fitting
FEH_ERR float64 [Fe/H] uncertainty from template fitting
VSINI float64 km s^-1 Stellar rotation velocity
CHISQ_TOT float64 Total chi-square for all arms
CHISQ_C_TOT float64 Total chi-square for all arms for polynomial only fit
CHISQ_B float64 Chi-square in the B arm
CHISQ_C_B float64 Chi-square in the B arm after fitting continuum only
CHISQ_R float64 Chi-square in the R arm
CHISQ_C_R float64 Chi-square in the R arm after fitting continuum only
CHISQ_Z float64 Chi-square in the Z arm
CHISQ_C_Z float64 Chi-square in the Z arm after fitting continuum only
RVS_WARN int64 RVSpecfit warning flag
SN_B float32 Median S/N in the B arm
SN_R float32 Median S/N in the R arm
SN_Z float32 Median S/N in the Z arm
SUCCESS bool Did we succeed or fail
Columns GAIA_SOURCE_ID to GAIA_PMDEC_ERROR are copied from Gaia DR3
GAIA_SOURCE_ID int64 Gaia unique source designation
GAIA_REF_EPOCH float32 yr Gaia reference epoch (Julian Years)
GAIA_RA float64 deg Gaia ICRS right ascension
GAIA_RA_ERROR float64 mas Gaia standard error of right ascension
GAIA_DEC float64 deg Gaia ICRS declination
GAIA_DEC_ERROR float64 mas Gaia standard error of declination
GAIA_PHOT_G_MEAN_MAG float32 mag Gaia G band magnitude
GAIA_PHOT_G_MEAN_FLUX_OVER_ERROR float32 Gaia G band signal-to-noise
GAIA_PHOT_BP_MEAN_MAG float32 mag Gaia BP band magnitude
GAIA_PHOT_BP_MEAN_FLUX_OVER_ERROR float32 Gaia BP band signal-to-noise
GAIA_PHOT_RP_MEAN_MAG float32 mag Gaia RP band magnitude
GAIA_PHOT_RP_MEAN_FLUX_OVER_ERROR float32 Gaia RP band signal-to-noise
GAIA_PHOT_BP_RP_EXCESS_FACTOR float32 Gaia BP/RP excess factor
GAIA_PHOT_G_N_OBS int32 Gaia G band number of observations
GAIA_ASTROMETRIC_EXCESS_NOISE float32 Gaia astrometric excess noise
GAIA_ASTROMETRIC_EXCESS_NOISE_SIG float32 Gaia significance of excess noise
GAIA_DUPLICATED_SOURCE bool Gaia duplicated source flag
GAIA_ASTROMETRIC_SIGMA5D_MAX float32 mas Gaia longest semi-major axis of the 5-d error ellipsoid
GAIA_ASTROMETRIC_PARAMS_SOLVED bool which astrometric parameters were estimated for a Gaia source
GAIA_RUWE float32 Gaia renormalised unit weight error
GAIA_IPD_GOF_HARMONIC_AMPLITUDE float32 Gaia amplitude of the IPD GoF versus position angle of scan
GAIA_IPD_FRAC_MULTI_PEAK bool Gaia percent of successful-IPD windows with more than one peak
GAIA_PARALLAX float32 mas Gaia parallax
GAIA_PARALLAX_ERROR float32 mas Gaia standard error of parallax
GAIA_PMRA float32 mas yr^-1 Gaia proper motion in right ascension direction
GAIA_PMRA_ERROR float32 mas yr^-1 Gaia standard error of proper motion in right ascension direction
GAIA_PMDEC float32 mas yr^-1 Gaia proper motion in declination direction
GAIA_PMDEC_ERROR float32 mas yr^-1 Gaia standard error of proper motion in declination direction
Columns RELEASE to PHOTSYS_PHOTOM are copied from Legacy Surveys DR9
RELEASE int16 Imaging surveys release ID
BRICKID int32 Brick ID from tractor input
BRICKNAME str Brick name from tractor input
OBJID int32 Imaging Surveys OBJID on that brick (renamed BRICK_OBJID)
TYPE str Morphological Model type from Tractor; renamed MORPHTYPE in most files
DCHISQ float32 Difference in chi-squared between Tractor model fits
EBV_SFD float32 mag Galactic extinction E(B-V) reddening from SFD98
FLUX_G float32 nanomaggy Flux in the Legacy Survey g-band (AB)
FLUX_R float32 nanomaggy Flux in the Legacy Survey r-band (AB)
FLUX_Z float32 nanomaggy Flux in the Legacy Survey z-band (AB)
FLUX_W1 float32 nanomaggy WISE flux in W1 (AB)
FLUX_W2 float32 nanomaggy WISE flux in W2 (AB)
FLUX_W3 float32 nanomaggy WISE flux in W3 (AB)
FLUX_W4 float32 nanomaggy WISE flux in W4 (AB)
FLUX_IVAR_G float32 nanomaggy^-2 Inverse variance of FLUX_G (AB)
FLUX_IVAR_R float32 nanomaggy^-2 Inverse variance of FLUX_R (AB)
FLUX_IVAR_Z float32 nanomaggy^-2 Inverse variance of FLUX_Z (AB)
FLUX_IVAR_W1 float32 nanomaggy^-2 Inverse variance of FLUX_W1 (AB)
FLUX_IVAR_W2 float32 nanomaggy^-2 Inverse variance of FLUX_W2 (AB)
FLUX_IVAR_W3 float32 nanomaggy^-2 Inverse variance of FLUX_W3 (AB)
FLUX_IVAR_W4 float32 nanomaggy^-2 Inverse variance of FLUX_W4 (AB)
MW_TRANSMISSION_G float32 Milky Way dust transmission in LS g-band assuming SFD98
MW_TRANSMISSION_R float32 Milky Way dust transmission in LS r-band assuming SFD98
MW_TRANSMISSION_Z float32 Milky Way dust transmission in LS z-band assuming SFD98
MW_TRANSMISSION_W1 float32 Milky Way dust transmission in WISE W1 assuming SFD98
MW_TRANSMISSION_W2 float32 Milky Way dust transmission in WISE W2 assuming SFD98
MW_TRANSMISSION_W3 float32 Milky Way dust transmission in WISE W3 assuming SFD98
MW_TRANSMISSION_W4 float32 Milky Way dust transmission in WISE W4 assuming SFD98
NOBS_G int16 Number of images for central pixel in g-band
NOBS_R int16 Number of images for central pixel in r-band
NOBS_Z int16 Number of images for central pixel in z-band
NOBS_W1 int16 Number of images that contribute to the central pixel in W1: filter for this object (not profile-weighted)
NOBS_W2 int16 Number of images that contribute to the central pixel in W2: filter for this object (not profile-weighted)
NOBS_W3 int16 Number of images that contribute to the central pixel in W3: filter for this object (not profile-weighted)
NOBS_W4 int16 Number of images that contribute to the central pixel in W4: filter for this object (not profile-weighted)
RCHISQ_G float32 Profile-weighted
RCHISQ_R float32 Profile-weighted
RCHISQ_Z float32 Profile-weighted
RCHISQ_W1 float32 Profile-weighted
RCHISQ_W2 float32 Profile-weighted
RCHISQ_W3 float32 Profile-weighted
RCHISQ_W4 float32 Profile-weighted
FRACFLUX_G float32 Profile-weighted fraction of the flux from other sources divided by the total flux in g (typically [0,1])
FRACFLUX_R float32 Profile-weighted fraction of the flux from other sources divided by the total flux in r (typically [0,1])
FRACFLUX_Z float32 Profile-weighted fraction of the flux from other sources divided by the total flux in z (typically [0,1])
FRACFLUX_W1 float32 Profile-weighted fraction of the flux from other sources divided by the total flux in W1 (typically [0,1])
FRACFLUX_W2 float32 Profile-weighted fraction of the flux from other sources divided by the total flux in W2 (typically [0,1])
FRACFLUX_W3 float32 Profile-weighted fraction of the flux from other sources divided by the total flux in W3 (typically [0,1])
FRACFLUX_W4 float32 Profile-weighted fraction of the flux from other sources divided by the total flux in W4 (typically [0,1])
FRACMASKED_G float32 Profile-weighted fraction of pixels masked from all observations of this object in g, strictly between [0,1]
FRACMASKED_R float32 Profile-weighted fraction of pixels masked from all observations of this object in r, strictly between [0,1]
FRACMASKED_Z float32 Profile-weighted fraction of pixels masked from all observations of this object in z, strictly between [0,1]
FRACIN_G float32 Fraction of a source’s flux within the blob in g, near unity for real sources
FRACIN_R float32 Fraction of a source’s flux within the blob in r, near unity for real sources
FRACIN_Z float32 Fraction of a source’s flux within the blob in z, near unity for real sources
ANYMASK_G int16 Bitwise mask set if the central pixel from any image satisfies each condition in g
ANYMASK_R int16 Bitwise mask set if the central pixel from any image satisfies each condition in r
ANYMASK_Z int16 Bitwise mask set if the central pixel from any image satisfies each condition in z
ALLMASK_G int16 Bitwise mask set if the central pixel from all images satisfy each condition in g
ALLMASK_R int16 Bitwise mask set if the central pixel from all images satisfy each condition in r
ALLMASK_Z int16 Bitwise mask set if the central pixel from all images satisfy each condition in z
WISEMASK_W1 uint8 Bitwise mask for WISE W1 data
WISEMASK_W2 uint8 Bitwise mask for WISE W2 data
PSFSIZE_G float32 arcsec Median PSF size evaluated at the BRICK_PRIMARY objects in this brick in g-band
PSFSIZE_R float32 arcsec Median PSF size evaluated at the BRICK_PRIMARY objects in this brick in r-band
PSFSIZE_Z float32 arcsec Median PSF size evaluated at the BRICK_PRIMARY objects in this brick in z-band
PSFDEPTH_G float32 nanomaggy^-2 PSF-based depth in g-band
PSFDEPTH_R float32 nanomaggy^-2 PSF-based depth in r-band
PSFDEPTH_Z float32 nanomaggy^-2 PSF-based depth in z-band
GALDEPTH_G float32 nanomaggy^-2 Galaxy model-based depth in LS g-band
GALDEPTH_R float32 nanomaggy^-2 Galaxy model-based depth in LS r-band
GALDEPTH_Z float32 nanomaggy^-2 Galaxy model-based depth in LS z-band
PSFDEPTH_W1 float32 nanomaggy^-2 PSF-based depth in WISE W1
PSFDEPTH_W2 float32 nanomaggy^-2 PSF-based depth in WISE W2
WISE_COADD_ID str unWISE coadd file name for the center of each object
SHAPE_R float32 arcsec Half-light radius of galaxy model (>0)
SHAPE_R_IVAR float32 arcsec^-2 Inverse variance of SHAPE_R
SHAPE_E1 float32 Ellipticity component 1 of galaxy model for galaxy type MORPHTYPE
SHAPE_E1_IVAR float32 Inverse variance of SHAPE_E1
SHAPE_E2 float32 Ellipticity component 2 of galaxy model for galaxy type MORPHTYPE
SHAPE_E2_IVAR float32 Inverse variance of SHAPE_E2
FIBERFLUX_G float32 nanomaggy Predicted g-band flux within a fiber of diameter 1.5 arcsec from this object in 1 arcsec Gaussian seeing
FIBERFLUX_R float32 nanomaggy Predicted r-band flux within a fiber of diameter 1.5 arcsec from this object in 1 arcsec Gaussian seeing
FIBERFLUX_Z float32 nanomaggy Predicted z-band flux within a fiber of diameter 1.5 arcsec from this object in 1 arcsec Gaussian seeing
FIBERTOTFLUX_G float32 nanomaggy Predicted g-band flux within a fiber of diameter 1.5 arcsec from all sources at this location in 1 arcsec Gaussian seeing
FIBERTOTFLUX_R float32 nanomaggy Predicted r-band flux within a fiber of diameter 1.5 arcsec from all sources at this location in 1 arcsec Gaussian seeing
FIBERTOTFLUX_Z float32 nanomaggy Predicted z-band flux within a fiber of diameter 1.5 arcsec from all sources at this location in 1 arcsec Gaussian seeing
MASKBITS int16 Bitwise mask from the imaging indicating potential issue or blending
FITBITS int16 Bitwise mask detailing properties of how a source was fit
SERSIC float32 Power-law index for the Sersic profile model (MORPHTYPE=’SER’)
SERSIC_IVAR float32 Inverse variance of SERSIC
TARGETID_PHOTOM int64 Unique ID based on the imaging data
PHOTSYS_PHOTOM str Photometric system of the imaging data: ‘N’ for the MzLS/BASS photometric system, ‘S’ for DECaLS
Columns MODEL_GMAG_N_REDDENED to ERZ_ERR are derived from the RVSpecfit model spectra
MODEL_GMAG_N_REDDENED float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in BASS g-band filter with extinction assuming SFD98
MODEL_GMAG_N float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in BASS g-band filter with no extinction
MODEL_GMAG_S_REDDENED float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam g-band filter with extinction assuming SFD98
MODEL_GMAG_S float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam g-band filter with no extinction
MODEL_RMAG_N_REDDENED float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in BASS r-band filter with extinction assuming SFD98
MODEL_RMAG_N float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in BASS r-band filter with no extinction
MODEL_RMAG_S_REDDENED float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam r-band filter with extinction assuming SFD98
MODEL_RMAG_S float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam r-band filter with no extinction
MODEL_IMAG_S_REDDENED float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam i-band filter with extinction assuming SFD98
MODEL_IMAG_S float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam i-band filter with no extinction
MODEL_ZMAG_N_REDDENED float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in MzLS z-band filter with extinction assuming SFD98
MODEL_ZMAG_N float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in MzLS z-band filter with no extinction
MODEL_ZMAG_S_REDDENED float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam z-band filter with extinction assuming SFD98
MODEL_ZMAG_S float64 mag RVSpecfit model magnitude (with arbitrary absolute zero point) derived from the SPECMOD spectrum in DECam z-band filter with no extinction
EGR float64 mag Reddening in DECam g-r color
EGR_ERR float64 mag Reddening error in DECam g-r color
ERZ float64 mag Reddening in DECam r-z color
ERZ_ERR float64 mag Reddening error in DECam r-z color

Contact

Contact Rongpu Zhou for questions about this catalog.