We demonstrate a simple, spectrally resolved ocean color remote sensing model to estimate benthic photosynthetically active radiation (bPAR) for the waters of the Great Barrier Reef (GBR), Australia. For coastal marine environments and coral reefs, the underwater light field is critical to ecosystem health, but data on bPAR rarely exist at ecologically relevant spatio-temporal scales. The bPAR model presented here is based on Lambert-Beer's Law and uses: (i) sea surface values of the downwelling solar irradiance, Es(λ); (ii) high-resolution seafloor bathymetry data; and (iii) spectral estimates of the diffuse attenuation coefficient, Kd(λ), calculated from GBR-specific spectral inherent optical properties (IOPs). We first derive estimates of instantaneous bPAR. Assuming clear skies, these instantaneous values were then used to obtain daily integrated benthic PAR values. Matchup comparisons between concurrent satellite-derived bPAR and in situ values recorded at four optically varying test sites indicated strong agreement, small bias, and low mean absolute error. Overall, the matchup results suggest that our benthic irradiance model was robust to spatial variation in optical properties, typical of complex shallow coastal waters such as the GBR. We demonstrated the bPAR model for a small test region in the central GBR, with the results revealing strong patterns of temporal variability. The model will provide baseline datasets to assess changes in bPAR and its external drivers and may form the basis for a future GBR water-quality index. This model may also be applicable to other coastal waters for which spectral IOP and high-resolution bathymetry data exist.