In the process of urban agglomeration, water-food security can be deemed as a key to support urban development and human living, but which can be challenged by expanded population growth, accelerated industrialization, unbalance regional economic development and diversity of weather (due to climate changes). In this study, a water resources allocation and food production (WF) optimization is developed for regional sustainability under multiple uncertainties. A hybrid two-stage fuzzy programming with Laplace criterion (TSFL) is proposed into a WF optimization to handle hybrid indeterminacies, which can increase the robustness of decision-making. The WF optimization with proposed TSFL method can be applied to a practical case of Beijing-Tianjin-Hebei (BTH) region. The obtained results associated with water deficits, optimal water allocations, inadequate capacities of food production, rational irrigation schedules, sound livestock scales, optimized agricultural possessing layouts and system benefits under various population-economy regulation scenarios can be obtained. The results can reflect the tradeoff between economic development and water-food safety; meanwhile, they can display risk violation of WF plan under various credibility levels and Laplace criterions (based on TSFL method). All above results can facilitate to produce an optimized water-food plan to support the synergetic development of BTH region in a robust manner.