The mechanical response of multiphase metallic materials is governed by the strain and stress partitioning behavior among their phases, crystals, and subgrains. Despite knowledge about the existence of these complex and long-ranging interactions, the experimental characterization of such materials is often limited to surface observations of microstructure evolution and strain partitioning, i.e. ignoring the influence of the underlying features. Hence, for the interpretation of the observed surface behavior it is imperative to understand how it might be influenced by the subsurface microstructure. In the present study, we therefore systematically change the subsurface microstructure of synthetic dual-phase polycrystals and investigate the altered response of a 2D region of interest. The series of high-resolution crystal plasticity simulations are conducted with a fast and efficient spectral-based iterative scheme for calculating the mechanical response of complex crystalline materials. To overcome the slow convergence of the conventional spectral-based solver when dealing with heterogeneous materials of large contrast in stiffness (or strength), direct and mixed variational conditions for mechanical equilibrium and strain compatibility have been formulated such that they can be combined with a general class of non-linear solution methods. The different solution techniques have been implemented into DAMASK, the Dusseldorf Advanced Material Simulation Kit, and the ones showing the best performance are used in this study. The results show that the subsurface microstructure has a dominant influence on the observed stress and strain partitioning. Additionally, it can be seen that the zone of influence increases with increasing heterogeneity of the microstructure.