In recent times, bioconvection has numerous uses, like, biological and biotechnological problems. The present study describes the magnetic bioconvective Buongiorno's flow model with microorganisms in a stretchable area with convective heat transfer and second-order velocity slip in a non-Darcian porous regime. Here, the influence of variable viscosity, viscous dissipation, Joule heating, and heat source/sink are considered in the occurrence of higher-order chemical reactions. Employing proper similarity transformations leading equations are transformed to dimension-free form. The transformed equations are solved via MATLAB bvp4c problem solver. This study's main objective is to graphically analyze the effects of different pertinent factors on the density of motile microorganisms, velocity, concentration, temperature, number of motile microorganisms' density, skin friction, mass transport rate, and heat transport rate. The main findings drawn from this study are viscosity