WL measurements have well-known shear estimation biases, which can be partially corrected for with the use of image simulations. We present an analysis of simulated images that mimic HST/ACS observations of high-redshift galaxy clusters, including cluster specific issues such as non-weak shear and increased blending. Our synthetic galaxies have been generated to match the observed properties of the background-selected samples in the real images. First, we used simulations with galaxies on a grid to determine a revised signal-to-noise-dependent correction for multiplicative shear measurement bias, and to quantify the sensitivity of our bias calibration to mismatches of galaxy or PSF properties between the real data and the simulations. We studied the impact of increased blending and light contamination from cluster and foreground galaxies, finding it negligible for $z>0.7$ clusters, whereas there is an effect at the $\sim 1\%$ level for lower redshift clusters. Finally, we studied the impact of fainter neighbours and selection bias mimicking the positions and magnitudes of galaxies in CANDELS data. The initial SExtractor object detection causes a selection bias of $-0.028 \pm 0.002$, reduced to $-0.016 \pm 0.002$ by further cuts. We compared our CANDELS-based estimate to a grid-based analysis, with added clustered galaxies reaching even fainter magnitudes, yielding a refined estimate of $\sim -0.013$. Our pipeline is calibrated to an accuracy of $\sim 0.015$, which is fully sufficient for current and near-future weak lensing studies of high-redshift clusters. As an application, we used it for a refined analysis of three highly relaxed clusters from the SPT-SZ survey, including measurements down to $r>200$ kpc. Compared to previously employed scales, this tightens the cluster mass constraints by a factor 1.38 on average.

In modern weak-lensing surveys, the common approach to correct for residual systematic biases in the shear is to calibrate shape measurement algorithms using simulations. These simulations must fully capture the complexity of the observations to avoid introducing any additional bias. In this paper we study the importance of faint galaxies below the observational detection limit of a survey. We simulate simplified Euclid VIS images including and excluding this faint population, and measure the shift in the multiplicative shear bias between the two sets of simulations. We measure the shear with three different algorithms: a moment-based approach, model fitting, and machine learning. We find that for all methods, a spatially uniform random distribution of faint galaxies introduces a shear multiplicative bias of the order of a few times $10^{-3}$. This value increases to the order of $10^{-2}$ when including the clustering of the faint galaxies, as measured in the Hubble Space Telescope Ultra-Deep Field. The magnification of the faint background galaxies due to the brighter galaxies along the line of sight is found to have a negligible impact on the multiplicative bias. We conclude that the undetected galaxies must be included in the calibration simulations with proper clustering properties down to magnitude 28 in order to reach a residual uncertainty on the multiplicative shear bias calibration of a few times $10^{-4}$, in line with the $2\times10^{-3}$ total accuracy budget required by the scientific objectives of the Euclid survey. We propose two complementary methods for including faint galaxy clustering in the calibration simulations.