<p>An efficient scheme is applied to generate a nonisospectral Botie-Pempinelli-Tu (BPT) integrable hierarchy under the case where <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(\lambda_{t} = \sum_{j=0}^n k_{j}(t){\lambda}^{-j}\)</EquationSource> <EquationSource Format="MATHML"><math display="block"> <msub> <mi>λ</mi> <mrow> <mi>t</mi> </mrow> </msub> <mo>=</mo> <munderover> <mo>∑</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>k</mi> <mrow> <mi>j</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>t</mi> <mo stretchy="false">)</mo> <msup> <mrow> <mi>λ</mi> </mrow> <mrow> <mo>−</mo> <mi>j</mi> </mrow> </msup> </math></EquationSource> </InlineEquation>. Based on an expanding higher-dimensional Lie algebra, we obtain a nonisospectral BPT integrable coupling hieararchy. It follws that some new nonisospectral nonlinear systems are obtained by reducing these two non-isospectral BPT hierarchies. Actually, these nonisospectral integrable models that we obtained can enrich the existing integrable models and possibly describe new nonlinear phenomena.</p>