When looking at futures options such as CME's Gold options or many equity index futures options, the underlying is not the index but to be precise the closest to delivery futures contract. That means, across a hypothetical surface, the underlying changes in theory.
That's not assumed in BS. I know there is the Black-76 model but I'm suddenly a bit uncertain about how to build or define a volatility surface for options which are in that way not "directly" related? I know this is being pedantic but I would like to understand the potential convexity effects involved.