A two-stage stochastic integer programming model is developed to address the joint capacity planning and distribution network optimization of multiechelon coal supply chains (CSCs) under uncertainty. The proposed model not only introduces the uses of compound real options in sequential capacity planning, but also considers uncertainty induced by both risks and ambiguities. Both strategic decisions (i.e., facility locations and initial investment, service assignment across the entire CSC, and option holding status) and scenario-based operational decisions (i.e., facility operations and capacity expansions, outsourcing policy, and transportation and inventory strategies) can be simultaneously determined using the model. By exploiting the nested decomposable structure of the model, we develop a new distributed parallel optimization algorithm based on nonconvex generalized Bender decomposition and Lagrangean relaxation to mitigate the computation resource limitation. One of the main CSCs in China is studied to demonstrate the applicability of the proposed model and the performance of the algorithm. (c) 2017 American Institute of Chemical Engineers AIChE J, 64: 1246-1261, 2018