AbstractWhile benchmark dose (BMD) methodology is well‐established for settings with a single exposure, these methods cannot easily handle multidimensional exposures with nonlinear effects. We propose a framework for BMD analysis to characterize the joint effect of a two‐dimensional exposure on a continuous outcome using a generalized additive model while adjusting for potential confounders via propensity scores. This leads to a dose–response surface which can be summarized in two dimensions by a contour plot in which combinations of exposures leading to the same expected effect are identified. In our motivating study of prenatal alcohol exposure, cognitive deficits in children are found to be associated with both the frequency of drinking as well as the amount of alcohol consumed on each drinking day during pregnancy. The general methodological framework is useful for a broad range of settings, including combinations of environmental stressors, such as chemical mixtures, and in explorations of the impact of dose rate rather than simply cumulative exposure on adverse outcomes.