In this study, an efficient hybrid continuum-atomistic method is proposed to study electrokinetic transport of aqueous solutions in nanofluidics. The aqueous phase is considered as a continuous phase containing immersed ion particles. The behavior of the system is then simulated through utilization of an improved hybrid continuum-atomistic four-way coupled approach, including the MultiPhase Particle-In-Cell method for the short-ranged interaction between the ion particles, the Brownian force for the collision between the aqueous phase molecules and the ion particles, and a wall force accounting for the short-ranged interaction of ions and walls. The validation of the proposed model with the results of Molecular Dynamics simulations suggests that this model can be a promising approach for studying the electrokinetic phenomena in more complicated geometries where the Molecular Dynamics approach is computationally prohibitive. Finally, the effects of electrokinetic parameters, such as the height of the channel, the external electric field, and bulk ionic concentration, on the electroosmotic flow in a nanochannel are investigated and discussed.