We report interferometric observations tuned to the redshifted neutral hydrogen (HI) 21cm emission line in three strongly lensed galaxies at $z \sim 0.4$ with the Giant Metrewave Radio Telescope (GMRT). One galaxy spectrum (J1106+5228 at z=0.407) shows evidence of a marginal detection with an integrated signal-to-noise ratio of 3.8, which, if confirmed by follow-up observations, would represent the first strongly lensed and most distant individual galaxy detected in HI emission. Two steps are performed to transcribe the lensed integrated flux measurements into HI mass measurements for all three target galaxies. First, we calculate the HI magnification factor $\mu$ by applying general relativistic ray-tracing to a physical model of the source-lens system. The HI magnification generally differs from the optical magnification and depends largely on the intrinsic HI mass $M_{\rm HI}$ due to the HI mass-size relation. Second, we employ a Bayesian formalism to convert the integrated flux, amplified by the $M_{\rm HI}$-dependent magnification factor $\mu$, into a probability density for $M_{\rm HI}$, accounting for the asymmetric uncertainty due to the declining HI mass function (Eddington bias). In this way, we determine a value of $\log_{\rm 10} (M_{\rm HI}/M_\odot) = 10.2^{+0.3}_{-0.7}$ for J1106+5228, consistent with the estimate of $9.4\pm0.3$ from the optical properties of this galaxy. The HI mass of the other two sources are consistent with zero within a 95 per cent confidence interval however we still provide upper limits for both sources and a $1\sigma$ lower limit for J1250-0135 using the same formalism.

The Event Horizon Telescope (EHT) aims to spatially resolve the silhouette (or shadow) of the supermassive black holes in the Galactic Centre (Sgr A$^\star$) and M87. The primary scientific objectives are to test general relativity in the strong-field regime and to probe accretion and jet-launch physics at event-horizon scales. This is made possible by the technique of Very Long Baseline Interferometry (VLBI) at (sub)millimetre wavelengths, which can achieve angular resolutions of order $\sim10~\mu$-arcsec. However, this approach suffers from unique observational challenges, including scattering in the troposphere and interstellar medium; rapidly time-variable source structure in both polarized and total intensity; as well as non-negligible antenna pointing errors. In this, the first paper in a series, we present the MeqSilhouette software package which is specifically designed to accurately simulate EHT observations. It includes realistic descriptions of a number of signal corruptions that can limit the ability to measure the key science parameters. This can be used to quantify calibration requirements, test parameter estimation and imaging strategies, and investigate systematic uncertainties that may be present. In doing so, a stronger link can be made between observational capabilities and theoretical predictions, with the goal of maximising scientific output from the upcoming order of magnitude increase in EHT sensitivity.