On convex lower-level black-box constraints in bilevel optimization with an application to gas market models with chance constraints