Impact of numerical stability in Bayesian noise wave calibration on global 21-cm experiments

Detecting the global 21-cm signal from the Cosmic Dawn and Epoch of Reionization requires calibration accuracy far below the level of astrophysical foregrounds. REACH models its receiver using the noise wave formalism, with five frequency-dependent low-noise amplifier parameters fitted jointly to multiple calibration sources. We identify a numerical instability in this Bayesian calibration pipeline: the condition number of the posterior covariance matrix reaches $\kappa(\mathbf{V}^*) \sim 10^{9}$--$10^{11}$, making solutions non-reproducible across computing environments. Singular value decomposition shows that the instability is driven by near-collinearity between the design-matrix columns associated with the excess noise source temperature, $X_\mathrm{NS}$, and the load temperature, $X_\mathrm{L}$. Using a Chebyshev basis, we develop a two-step mitigation. First, fixing $T_\mathrm{NS}$ to a scalar removes the degeneracy and reduces $\kappa(\mathbf{V}^*)$ to $\sim 60$. Second, to retain frequency dependence, we recover $T_\mathrm{NS}(\nu)$ directly from the hot-load calibration measurement. On mock data, this method preserves the stability of the reduced model while achieving comparable calibration accuracy. Masking narrow channels around cable standing-wave degeneracies further removes local artefacts in the design matrix. These steps provide a stable, reproducible, and data-driven calibration procedure. Because the $X_\mathrm{NS}$--$X_\mathrm{L}$ degeneracy is inherent to the noise wave formalism, the method is relevant to other global 21-cm experiments.