Response of a mesoscopic superconducting disk to perpendicular magnetic fields is studied by using the multiple-small-tunnel-junction method, in which transport properties of several small tunnel junctions attached to the disk are measured simultaneously. This allows us for the first experimental distinction between the giant vortex states and multivortex states. Moreover, we experimentally find magnetic field induced rearrangement and combination of vortices. The experimental results are well reproduced in numerical results based on the nonlinear Ginzburg-Landau theory.The appearance of vortices in various quantum systems, such as superconductors, superfluids and BoseEinstein condensates, is an intriguing phenomenon in nature. A conventional quantum vortex is singly quantized, having a core where the value of the order parameter decreases to zero, while its phase changes by 2π when encircling the core. Recently, an important breakthrough was established by the observation of doubly quantized vortex lines in superfluid 3 He-A[1]. For superconductors expectations are even more spectacular. In macroscopic type-II superconductors a triangular lattice of single flux quanta is formed, whereas two kinds of fundamentally new vortex states have theoretically been predicted in mesoscopic superconductors where the sample size approaches the size of Cooper-pairs [2,3,4,5]; (i) multivortex states (MVSs) with a unique spatial arrangement of singly quantized vortices, and (ii) multiply quantized or giant vortex states (GVSs) with a single core in the center [6,7].Although several experimental techniques have been developed for observing these novel states [7,8,9,10,11,12,13], none of them has been able to make a clear distinction between MVSs and GVSs. In this Letter, we present the first experimental evidence for the existence of GVSs and MVSs in a circular disk, and demonstrate magnetic-field induced MVS-GVS and MVS-MVS transitions. Our results are in good agreement with the theoretical prediction based on the nonlinear GinzburgLandau (G-L) theory.Here we used the multiple-small-tunnel-junction (MSTJ) method, in which several small tunnel junctions with high tunnel resistance are attached to a mesoscopic superconductor to simultaneously detect small changes in the local density of states (LDOS) under the junctions [14,15]. Since the LDOS depends on the local supercurrent density, the MSTJ method gives us information on the distribution of the supercurrent, which reflects the detailed vortex structure inside the disk. Figure 1 shows a schematic drawing and a scanning electron microscopy (SEM) image of the sample. Four normal-metal (Cu) leads are connected to the periphery of the superconducting Al disk through highly resistive small tunnel junctions, A, B, C, and D. The sample is designed to be symmetrical with respect to the central axis SS ′ . The angles AOD and BOC are 120 and 32 degrees, respectively. Although junctions A and D and junctions B and C ideally have the same area and tunnel resistance, small differences actu...