Most models used to account for the hardening of nanocomposites only consider a global volume fraction of particles which is a simplified indicator that overlooks the particles size and spatial distribution. The current study aims at quantifying the effect of the real experimental particles spatial and size distribution on the strengthening of a magnesium based nanocomposites reinforced with Y2O3 particles processed by Friction Stir Processing (FSP). X-ray tomographic 3-D images allowed to identify the best FSP parameters for the optimum nanocomposite. A detailed analysis indicates that the observed hardening is mainly due to Orowan strengthening and the generation of geometrically necessary dislocations (GND) due to thermal expansion coefficients (CTE) mismatch between magnesium and Y2O3 particles. A multiscale characterization coupling 3D X-ray laboratory, synchrotron nanoholotomography and transmission electron microscopy (TEM) has been used to investigate particles size and spatial distribution over four orders of magnitude in length scales. Two dedicated micromechanical models for the two strengthening mechanisms are applied on the experimental particle fields taking into account the real particles size and spatial distribution, and compared to classical models based on average data. This required to develop a micromechanical model for CTE mismatch hardening contribution. This analysis reveals that the contribution from CTE mismatch is decreased by a factor two when taking into account the real distribution of particles instead of an average volume fraction.