Weak lensing catalogue and image simulations

UNIONS-3500 Weak Lensing, Rencontres de Moriond 2026, La Thuile


Sacha Guerrini

on behalf of the UNIONS collaboration

March 27, 2026

The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS, Gwyn et al., 2025)

\(r\), \(u\) bands
\(i\) band
\(g\), \(z\) bands

Euclid photometry in the North

Scaramella+2025

An excellent photometry

Depth of the photometric bands (Gwyn et al., 2025)

FWHM in the \(r\)-band: \(\sim 0.7"\) (DES: \(\sim 0.98"\))

Outline of the talk

  1. Selection and weak lensing sample properties
  2. Systematic tests and mitigation strategy
  3. Image simulations and shear bias calibration
  4. Weak-lensing science with UNIONS

1. Selection and weak lensing sample properties

Wright+2025

Technical specifications of UNIONS:

  • Target area: \(\sim 6,250\) deg\(^2\)
  • Current area: \(\sim 3,500\) deg\(^2\)
  • Depth: \(24.5\) (\(r\)-band)
  • Seeing: \(\sim 0.7"\) (\(r\)-band)
  • Shape measurement only in the \(r\)-band
  • Image reduction with MegaPipe (Gwyn+2008)
  • Shape processing with ShapePipe (Farrens+2022)

1. Source detection

Credit: F. Hervas-Peters

1. Masking

3 steps of masking:

  • Pre-processing (MegaCam)
    • chip defects
    • cosmic rays
  • ShapePipe masks
    • Blended objects
    • Bright stars
    • Messier or NGC objects
    • ngmix quality flags
    • others…
  • Post-processing (manual, maximask Paillassa+2020)
    • remaining bright star haloes
    • satellite trails
    • number of pointings/exposures

Only keep \(\sim 25 \%\) of the detected objects. Reduces the \(3,500 \textrm{deg}^2\) imaging area to an effective footprint of \(2,894 \textrm{deg}^2\).

1. Final galaxy selection

  • \(15 \leq r_{\rm mag} \leq 30\) (SExtractor magnitude)
  • \(10 \leq \mathrm{SNR} \leq 500\) (ngmix ratio of flux to flux error)
  • \(0.707 \leq r_{\rm h} / r_{\rm h, PSF} \leq 3\) (\(r_{\rm h, .}\) galaxy or PSF half-light radius)

Last two selections applied in Metacalibration (Huff & Mandelbaum, 2017) to determine the selection bias.

1. Weak lensing sample properties

  • Number of galaxies: \(N_{\rm gal} = 61,378,891\)
  • Effective area: \(A_{\rm eff} = 2,894 \textrm{deg}^2\)
  • Effective number density: \(n_{\rm eff} = 5.09\) galaxies per arcmin\(^2\)
  • Shape noise per component: \(\sigma_{\rm e} = 0.27\)

2. Systematic tests and mitigation strategy

galaxy clustering

Mandelbaum 2018

2. Calibrating biases

Point spread function (PSF) error model: \(\delta \boldsymbol{e}^\mathrm{sys}_\mathrm{PSF} = \alpha \underbrace{\boldsymbol{e}_\mathrm{PSF}}_{\text{Leakage}} + \beta \underbrace{(\boldsymbol{e}_* - \boldsymbol{e}_\mathrm{PSF})}_{\text{Ellipticity error}} + \eta \underbrace{\boldsymbol{e}_\mathrm{PSF} \left(\frac{T_* - T_\mathrm{PSF}}{T_*} \right)}_{\text{Size error}}\).

\(\alpha\), \(\beta\) and \(\eta\) free parameters.

Estimate leakage in bins of size ratio and SNR doing linear regressions.

Subtract the contribution in each objects ellipticity from PSF leakage.

2. Galaxy-PSF cross-correlations (Guerrini+2025)

Galaxy-Star correlations: \(\tau_0(\vartheta) = \langle e^\mathrm{obs} e^\mathrm{PSF} \rangle(\vartheta)\), \(\tau_2(\vartheta) = \langle e^\mathrm{obs} \delta e^\mathrm{PSF} \rangle(\vartheta)\), \(\tau_5(\vartheta) = \langle e^\mathrm{obs} \delta T^\mathrm{PSF} \rangle(\vartheta)\)

2. Galaxy-PSF cross-correlations (Guerrini+2025)

Galaxy-PSF correlation helps to:

  • Select our weak lensing sample
  • Choose scale cuts for the inference

\[ \xi_\mathrm{sys}(\vartheta) = \alpha \tau_0(\vartheta) + \beta \tau_2(\vartheta) + \eta \tau_5(\vartheta) \]

2. Other systematics check

No significant correlation of shapes around CCD and tile centers.

2. Other systematics check

No significant correlation of shapes around Gaia stars.

3. Image simulations and shear bias calibration

Multiplicative bias origin:

  • Model bias
  • Noise bias
  • Selection effect
  • Blending (not accounted for by Metacalibration)

Image simulations on a grid (left) and with realistic positions (right). Credit: F. Hervas-Peters

3. Image simulations

3. \(m\)-bias estimations

  • Estimation of the multiplicative bias with respect to nearest neighbour distance in the input catalogue.
  • We can see the improved seeing of UNIONS compared to DES.
  • FWHM:
    • UNIONS: \(\sim 0.7"\)
    • DES: \(\sim 0.98"\)

3. \(m\)-bias estimations

  • Evolution of the \(m\)-bias with selection effects, offset bug in the pipeline, and realistic positionning.

4. Weak-lensing science with UNIONS

  • UNIONS has a unique overlap with spectroscopic surveys.
  • Opens exciting opportunities for cosmic shear and galaxy-galaxy lensing analyses with DESI.

4. Intrinsic Alignment of galaxies

  • Strong detection of intrinsic alignment with SDSS CMASS and LRGs
  • Direct measurement used in the cosmological analysis (see Cail’s talk)
  • More intrinsic alignment studies in preparation (see Calum’s talk)

Hervas-Peters+2025

4. Cluster lensing

  • Cluster selection using DLIS + DESI DR1 spectroscopic redshifts.
  • Lensing signal is obtained from UNIONS weak lensing catalog.
  • Selection of evolved and evolving clusters.

Ahad+2025

4. Void lensing

  • Void lensing using BOSS foreground lenses and UNIONS background sources.
  • \(6.2 \, \sigma\) detection of void lensing.
  • Galaxy bias relative to underdense regions, \(b_{\rm g} = 2.45 \pm 0.36\)

Martin+2026

Conclusion

  • UNIONS is a unique photometric survey with promising cross-correlations with spectroscopic surveys in the northern hemisphere
  • Current weak lensing sample of \(5.09\) galaxies per arcmin\(^2\) over \(2,894 \textrm{deg}^2\)
  • A target area of \(6,250 \textrm{deg}^2\) will be reached in the next years.
  • First cosmic shear cosmological release of the collaboration: UNIONS-3500 weak lensing
    • Paper I: Weak lensing catalogues (Hervas-Peters et al., in prep.)
    • Paper II: Image simulations for cosmic shear (Hervas-Peters et al., in prep.)
    • Paper III: \(B\)-mode validation and comparison for cosmic shear (Daley et al., in prep.)
    • Paper IV: 2D cosmological constraints in configuration space (Goh et al., in prep.)
    • Paper V: 2D cosmological constraints in harmonic space (Guerrini et al., in prep.)

Thank you for your attention. Special acknowledgements to F. Hervas-Peters for his contribution to this work and this presentation.
Image: JC Cuillandre

Backup slides: star selection

  • Size-magnitude distribution of SExtractor detections on a tile. Violet points are selected stars along the stellar locus.

Backup slides: star selection

  • Mean ellipticities of the selected galaxies in the MegaCam focal plane. No significant structure on the focal plane.

Backup slides: systematic maps

Backup slides: systematic deprojection

Backup slides: noise bias in shape measurement