In this study, magnetohydrodynamic cross fluid model is used to formulate the 2D boundary layer equations of fluid moving on the parabolic surface. The surface is assumed to be in vertical shape due to the convection. The underlying effects of viscous dissipation and chemical reaction are modeled to observe the transportation of heat and mass rate. To understand fluid behavior, the thermophysical properties are considered as variable because they have huge influence in food processing, viscometers, lubricants and various industrial works. The assumed geometry of fluid flows is similar in shape of bullet, submarine, aircraft and car’s bonnet, namely paraboloid surface. The modeled equations of cross fluid with mentioned effects are obtained in form of PDEs and then converted these equations in form of ODEs by assuming set of scaling transformations. For the sake of numerical and graphical outcomes, the resulting equations are solved numerically on Matlab software via BVP4c method. Results achieved across velocity indicate that the decay happens by numerous values of Weissenberg number, power law index, viscosity parameter and Hartmann number. The temperature and concentration fields attained increasing influence by thermal conductivity coefficient, Eckert number and thermal diffusivity parameter.