Energy densities of six dominant benthic groups (Actinopterygii, Amphipoda, Bivalvia, Cumacea, Isopoda, and Polychaeta) and total prey energy were modeled for the nearshore western gray whale feeding area, Sakhalin Island, Russia, as part of a multi-disciplinary research program in the summer of 2015. Energy was modeled using generalized additive mixed models (GAMM) with accommodations for zero-inflation (logistic regression and hurdle models) and regression predictions combined with kriging to interpolate energy densities across the nearshore feeding area. Amphipoda energy density was the highest nearshore and in the south whereas Bivalvia energy density was the highest offshore and in the northern portion of the study area. Total energy was the highest in mid-range distances from shore and in the north. Amphipoda energy density was higher than minimum energy estimates defining gray whale feeding habitats (312–442 kJ/m2) in 13% of the nearshore feeding area whereas total prey energy density was higher than the minimum energy requirement in 49% of the habitat. Inverse distance-weighted interpolations of Amphipoda energy provided a broader scale representation of the data whereas kriging estimates were spatially limited but more representative of higher density in the southern portion of the study area. Both methods represented the general trend of higher Amphipoda energy density nearshore but with significant differences that highlight the value of using multiple methods to model patterns in highly complex environments.