Nodal bubble tower solutions to slightly subcritical elliptic problems with Hardy terms

We study the possible blow-up behavior of solutions to the slightly subcritical elliptic problem with Hardy term $$\{\begin{array}{rlrl}{\displaystyle -\mathrm{\Delta }u-\mu \frac{u}{\mathrm{\mid }x{\mathrm{\mid }}^2}}& {\displaystyle =\mathrm{\mid }u{\mathrm{\mid }}^{2^{\ast }-2-\epsilon }u}& {\displaystyle }& {\displaystyle \text{in }\mathrm{\Omega },}\\ {\displaystyle u}& {\displaystyle =0}& {\displaystyle }& {\displaystyle \text{on }\mathrm{\partial }\mathrm{\Omega },}\end{array}$$ in a bounded domain $\mathrm{\Omega }\subset {\mathbb{R}}^N(N\ge 7)$ with $0\in \mathrm{\Omega }$, as $\mu ,\epsilon \to 0^{+}$. In \cite{BarGuo-ANS}, we obtained the existence of nodal solutions that blow up positively at the origin and negatively at a different point as $\mu =O({\epsilon }^{\alpha })$ with $\alpha >\frac{N-4}{N-2}$, $\epsilon \to 0^{+}$. Here we prove the existence of nodal bubble tower solutions, i.e.\ superpositions of bubbles of different signs, all blowing up at the origin but with different blow-up order, as $\mu =O(\epsilon )$, $\epsilon \to 0^{+}$.