Massive multiple-input multiple-output (MIMO) technology and millimeter wave (mmWave) communication are key technologies for the fifth generation (SG) wireless communications. The combination of mmWave and massive MIMO has the potential to dramatically improve wireless access and throughput performance. Such systems benefit from large available signal bandwidths and small antenna form factor. In the literature, most of the massive MIMO channel measurements are carried out at sub-6 GHz frequency bands, while the effects caused by large antenna arrays at mmWave bands have not been studied yet. In this paper, we conduct channel measurements at 11, 16, 28, and 38 GHz frequency bands combined with large antenna arrays in an indoor office environment. The space-alternating generalized expectation-maximization (SAGE) algorithm is applied to obtain the multipath component (MPC) parameters. New propagation characteristics like spherical wavefront, cluster birth-death, and non-stationarity over antenna array axis are validated for the four mmWave bands by investigating the temporal-spatial channel characteristics like power delay profile (PDP), power azimuth profile (PAP), power elevation profile (PEP), root mean square (RMS) delay spread (DS), and azimuth and elevation angular spread (AS). The results indicate that massive MIMO effects should be fully considered for mmWave channel models under systems with large antenna arrays.