Scale dependence of electrostatic and magnetostatic properties is investigated in the setting of spatially random linear lossless materials with statistically homogeneous and spatially ergodic random microstructures. First, from the Hill-Mandel homogenization conditions adapted to electric and magnetic fields, uniform boundary conditions are formulated for a statistical volume element (SVE). From these conditions, there follow upper and lower mesoscale bounds on the macroscale (effective) electrical permittivity and magnetic permeability. Using computational electromagnetics methods, these bounds are obtained through numerical simulations for composites of two types: (i) two-dimensional (2D) random checkerboard (two-phase) microstructures and (ii) analogous 3D random (two-phase) media. The simulation results demonstrate a scale-dependent trend of these bounds toward the properties of a representative volume element (RVE). This transition from SVE to RVE is described using a scaling function dependent on the mesoscale δ, the volume fraction v_{f}, and the property contrast k between two phases. The scaling function is calibrated through fitting the data obtained from extensive simulations (∼10000) conducted over the aforementioned parameter space. The RVE size of a given microstructure can be estimated down to within any desired accuracy using this scaling function as parametrized by the contrast and the volume fraction of two phases.