There are few methods to calculate PAR underwater from PAR subsurface values. There are all based on the exponential decrease of light within th water column, which is controled by the vertical attenuation coefficient KdPAR. You can get such KdPAR from papers of Morel and also the one by Lee
Morel, A. (1988). Optical modeling of the upper ocean in relation to its biogenous matter content (case 1 water), Journal of Geophysical Research, 93, 10,749-10,768.
Lee Zhongping (2009) KPAR: An optical property associated with ambiguous values. Journal of Lake Sciences.