The first analysis of the photometric observation in BVR filters of a W UMa type binary system BQ Ari was performed. Light curve analysis was performed using Wilson-Devinney (W-D) code combined with a Monte Carlo (MC) simulation to determine its photometric and geometric elements and their uncertainties. These results show that BQ Ari is a contact binary system with a photometric mass ratio q=0.548pm0.019, a fillout factor f=24pm0.8 percent, and an orbital inclination of i=85.09pm0.45. We used the parallax from Gaia EDR3 for calculating the absolute parameters of the binary system. This study suggested a new linear ephemeris for BQ Ari, combining our new mid-eclipse times and the previous observations, which we analyzed using the Monte Carlo Markov Chain (MCMC) method. We present the first analysis of the systems orbital period behavior by analyzing the O-C diagram using the Genetic Algorithm (GA) and the MCMC approaches in OCFit code. We attempted to explain the analysis of the residuals of linear fit in the O-C diagram with two approaches; LiTE + Quadratic and Magnetic activity + Quadratic. Although we consider the magnetic activity to be probable, the system should be studied further in order to reveal the nature of orbital period variations.