In the present study, a novel theoretical model is developed for the analysis of rotating thermal-fluid flow characteristics in the presence of electrokinetic effects in the microscale gap region between two parallel disks under specified electrostatic, rotational, and thermal boundary conditions. The major flow configuration considered is a rotor-stator disk system. Axisymmetric Navier-Stokes equations with consideration of electric body force stemming from streaming potential are employed in the momentum balance. Variations of the fluid viscosity and permittivity with the local fluid temperature are considered. Between two disks, the axial distribution of the electric potential is determined by the Poisson equation with the concentration distributions of positive and negative ions obtained from Nernst-Planck equations for convection-diffusion of the ions in the flow field. Effects of disk rotation and electrostatic and thermal conditions on the electrokinetic flow and thermal characteristics are investigated. The electrohydrodynamic mechanisms are addressed with an interpretation of the coupling nature of the electric and flow fields. Finally, solutions with electric potential determined by employing nonlinear or linearized Poisson-Boltzmann equation and/or invoking assumptions of constant properties are compared with the predictions of the present model for justification of various levels of approximation in solution of the electrothermal flow behaviors in rotating microfluidic systems. (C) 2003 Elsevier Inc. All rights reserved.