<p>This study comprehensively investigates the effects of gravitational variations and throughflow of convection in anisotropic porous media under LTNE (local thermal non-equilibrium conditions). By using a normal mode analysis, the linear stability analyses are analyzed. The gravitational force is modelled in three forms- linear, parabolic, and exponential- all oriented along the spatial configuration’s vertical axis. The critical Rayleigh number for convection onset is approximated through the Galerkin approximation. The results show that the throughflow constraint consistently stabilizes the system, regardless of direction. Additionally, an increase in the thermal and mechanical conductivities leads to enhanced system stability, while greater mechanical anisotropy weakens the convection onset. Notably, the system demonstrates the highest stability under exponential gravity variations, compared to linear and parabolic variations. The initiation and consistency of convection in these systems play a crucial role in: Forecasting thermal extraction performance, reservoir temperature reduction, and crafting eco-friendly energy harvesting approaches.</p>