Thanks for your insightful questions. We hope that the following answers will clarify them, and we are happy to discuss them further.
1. Span semi norm $H$: As you suggest, the parameter $H$ has advantages over $t_{mix}$ as the basis for an upper bound. We included a short discussion of this in Section 5. Moreover, we have evidence suggesting that the natural analog of our algorithm should have a sample complexity upper bound of $\widetilde O (|S||A|H\epsilon^{-2})$, if we replace the uniform ergodicity assumption by one that assumes the span semi-norm of the transient value function induced by any policy is bounded by $H$. However, we didn't adopt this setup for the following reasons. First, $t_{mix}$ places a direct assumption on the primitives---namely the transition kernel---of the AMDP, yet the assumption on $H$ involves the solution to Poisson equations. Second, we are aiming at a "worst case" complexity result for the AMDP. The $H$ metric, considering the interplay between the reward and the transition kernel, has the flavor of instance-dependent complexity theory (see Khamaru et al. (2021), cited in the paper).
2. Reduction procedure: Thanks for asking about this. Before giving a direct response, we would like to point out that the idea of solving an AMDP by reducing it to a DMDP and applying dynamic programming has been considered and implemented since the 1970s, see Hordijk and Tijms (1975). In fact, the optimality of our DMDP algorithm implies that one can choose $\gamma$ of the approximating DMDP to be much smaller than that in Jin and Sidford (2021), where $(1-\gamma) = \Theta(\epsilon/t_{mix})$, and still get the optimal sample complexity. This is not the case in Jin and Sidford (2021), where enlarging $\gamma$ will worsen the sample complexity. However, note that choosing large $\gamma$ is inefficient from a computational standpoint, because solving the empirical MDP with large $\gamma$ is typically hard (consider using value iteration), and hence is not implemented in the paper. Therefore, the reduction method is mainly a computational tool. So, statistically speaking, we think it is possible to achieve the same sample complexity using a model-based approach without the reduction. However, this might lead to computational complications, e.g. slow convergence rate for non-mixing empirical AMDPs.
3. Extension to the online setting: The key ideas and bounds extend to the online setting without any issue. As our answer to your question (2) suggests, our bound is not reliant on a particular choice of the discount $\gamma$. Therefore, we do believe that our theory can be used to extend classical model-based online discounted RL algorithms to the average reward setting.