Let $\Phi_1,\Phi_2 \colon S(\mathcal{H}) \to S(\mathcal{H})$ be CPTP maps on the same Hilbert space $\mathcal{H}$ which are $\varepsilon$-close in diamond norm, and let $U_1,U_2$ be respective unitary dilations on some larger space $\mathcal{H} \otimes \mathcal{K}$ (i.e. so that $\mathrm{Tr}_{\mathcal{K}}(U_i (\rho \otimes |0\rangle\langle0|) U_i^{\dagger}) = \Phi_i(\rho)$).
Now let $A$ be a unitary operating on $\mathcal{H}$, and let $$\rho_i := \mathrm{Tr}_{\mathcal{K}}(U_i^{\dagger} (A \otimes I_{\mathcal{K}}) U_i (\rho \otimes |0\rangle\langle0|) U_i^{\dagger} (A^{\dagger} \otimes I_{\mathcal{K}}) U_i).$$ That is, $\rho_i$ is obtained from $\rho$ by applying $U_i$, then $A$, then $U_i^{\dagger}$, then tracing out $\mathcal{K}$. Since all dilations are equivalent up to a local unitary, $\rho_i$ does not depend on the choice of $U_i$. Can we bound $\| \rho_1 - \rho_2 \|_1$?
I have tried to get a bound from results on the distance between Stinespring isometries, but this seems to not be enough here.