Quantifying spatial variability in intertidal benthic productivity is necessary to guide management of estuaries and to understand estuarine ecological processes, including the amount of benthic organic carbon available for grazing, burial and transport to the pelagic zone. We developed a model to assess microphytobenthic (MPB) primary production using (1) remotely sensed information on MPB biomass and remotely sensed information on sediment mud content, (2) surface irradiance and ambient temperature (both from local meteorological observations), (3) directly-measured photosynthetic parameters and (4) a tidal model. MPB biomass was estimated using the normalized-difference vegetation index (NDVI) and mud content was predicted using surface reflectance in the blue and near-infrared, both from Landsat 8 satellite imagery. The photosynthetic capacity (maximum photosynthesis rate normalized to MPB chl-a) was estimated from ambient temperature, while photosynthetic efficiency and the light saturation parameter were derived from in situ fluorometry-based production measurements (PAM). The influence of tides (submergence by turbid water) on MPB production was accounted for in the model. The method was validated on several locations in two temperate tidal basins in the Netherlands (Oosterschelde and Westerschelde). Model based production rates (mg C m−2 h−1) matched well with an independent set of in situ (PAM) measurement based production rates (Oosterschelde: RMSE = 9.7, mean error = 1.5, χ = 0.57; Westerschelde: RMSE = 46.7, mean error = −17.6, χ = 0.9). The relationship between photosynthetic capacity and temperature shows considerable variation and may be improved by using sediment surface temperature instead of ambient temperature. A sensitivity analysis revealed that emersion duration and mud content determine most of the variability in MPB production. Our results demonstrate that it is possible to derive a satellite remote sensing-based overview of average hourly and daily MPB primary production rates at the macro scale. As the proposed model is generic, the model can be applied to other tidal systems to assess spatial variability in MPB primary production at the macro scale after calibration at the site of interest. Model calibration, results and possible applications for regular monitoring of MPB production are discussed below.