The main objective of the present article is to provide an analytical simulation for time dependent boundary layer flow of non-Newtonian Williamson fluid due to stretchable curved oscillatory Riga surface. Also the characteristics of heat and mass transport are studied under the influence of variable thermal conductivity and diffusivity along with convective heat and mass boundary conditions. Additionally, energy equation is also characterized with the impact of heat production. Curvilinear coordinate scheme is followed to attain the boundary layer expressions for the flow model. The governing nonlinear partial differential equations are solved analytically via homotopy analysis method (HAM). Graphs are plotted to examine comprehensively the consequences of various concerned parameters like modified magnetic parameter and radius of curvature, Williamson fluid parameter, relation of the surface's oscillating frequency to its stretching rate constant, Prandtl number, variable conductivity and heat production parameters, Schmidt number and variable diffusivity parameter on concentration, temperature, pressure and velocity profile. Also the outcomes of afore said variables on surface drag force, rate of temperature and mass transmission (Nusselt and Sherwood numbers) are shown in tabular form. The liquid velocity amplitude is enhanced with modified magnetic parameter and shows opposite behavior for Williamson fluid parameter.