The collective mechanical behavior of multilayer colloidal arrays of hollow silica nanoparticles (HSNP) is explored under spherical nanoindentation through a combination of experimental, numerical, and theoretical approaches. The effective indentation modulus E(ind) is found to decrease with an increasing number of layers in a nonlinear manner. The indentation force versus penetration depth behavior for multilayer hollow particle arrays is predicted by an approximate analytical model based on the spring stiffness of the individual particles and the multipoint, multiparticle interactions as well as force transmission between the layers. The model is in good agreement with experiments and with detailed finite element simulations. The ability to tune the effective indentation modulus, E(ind), of the multilayer arrays by manipulating particle geometry and layering is revealed through the model, where E(ind) = (0.725m(-3/2) + 0.275)E(mon) and E(mon) is the monolayer modulus and m is number of layers. E(ind) is seen to plateau with increasing m to E(ind_plateau) = 0.275E(mon) and E(mon) scales with (t/R)(2), t being the particle shell thickness and R being the particle radius. The scaling law governing the nonlinear decrease in indentation modulus with an increase in layer number (E(ind) scaling with m(-3/2)) is found to be similar to that governing the indentation modulus of thin solid films E(ind_solid) on a stiff substrate (where E(ind_solid) scales with h(-1.4) and also decreases until reaching a plateau value) which also decreases with an increase in film thickness h. However, the mechanisms underlying this trend for the colloidal array are clearly different, where discrete particle-to-particle interactions govern the colloidal array behavior in contrast to the substrate constraint on deformation, which governs the thickness dependence of the continuous thin film indentation modulus.