We use the Separate Universe technique to calibrate the dependence of linear and quadratic halo bias $b_1$ and $b_2$ on the local cosmic web environment of dark matter haloes. We do this by measuring the response of halo abundances at fixed mass and cosmic web tidal anisotropy $\alpha$ to an infinite wavelength initial perturbation. We augment our measurements with an analytical framework developed in earlier work which exploits the near-Lognormal shape of the distribution of $\alpha$ and results in very high precision calibrations. We present convenient fitting functions for the dependence of $b_1$ and $b_2$ on $\alpha$ over a wide range of halo mass for redshifts $0\leq z\leq1$. Our calibration of $b_2(\alpha)$ is the first demonstration to date of the dependence of non-linear bias on the local web environment. Motivated by previous results which showed that $\alpha$ is the primary indicator of halo assembly bias for a number of halo properties beyond halo mass, we then extend our analytical framework to accommodate the dependence of $b_1$ and $b_2$ on any such secondary property which has, or can be monotonically transformed to have, a Gaussian distribution. We demonstrate this technique for the specific case of halo concentration, finding good agreement with previous results. Our calibrations will be useful for a variety of halo model analyses focusing on galaxy assembly bias, as well as analytical forecasts of the potential for using $\alpha$ as a segregating variable in multi-tracer analyses.

The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $\alpha$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represent two distinct statistical links: one between the internal property and $\alpha$, and the other between $\alpha$ and large-scale ( $>30h^{-1}$Mpc) halo bias $b_1$. We focus on scalar internal properties of haloes related to formation time (concentration $c_{\rm vir}$), shape (mass ellipsoid asphericity $c/a$), velocity dispersion structure (velocity ellipsoid asphericity $c_v/a_v$ and velocity anisotropy $\beta$) and angular momentum (dimensionless spin $\lambda$) in the mass range $8\times10^{11}< M_{\rm vir}/(h^{-1}M_\odot)<5\times10^{14}$. Using conditional correlation coefficients and other detailed tests, we show that the joint distribution of $\alpha$, $b_1$ and any of the internal properties $c\in\{\beta,c_v/a_v,c/a,c_{\rm vir},\lambda\}$ is consistent with $p(\alpha,b_1,c)\simeq p(\alpha)p(b_1|\alpha)p(c|\alpha)$, at all but the largest masses. $\textit{Thus, the assembly bias trends $c-b_1$ reflect the two fundamental correlations $c-\alpha$ and $b_1-\alpha$.}$ Our results are unaffected by the exclusion of haloes with recent major merger events or splashback objects, although the latter are distinguished by the fact that $\alpha$ does not explain their assembly bias trends. The overarching importance of $\alpha$ provides a new perspective on the nature of assembly bias of distinct haloes, with potential ramifications for incorporating realistic assembly bias effects into mock catalogs of future large-scale structure surveys and for detecting galaxy assembly bias.