This paper proposes a low-carbon operation model for an energy hub (EH) that combines the distributionally robust optimization (DRO) method with the Stackelberg game. Firstly, a bilevel single-leader-multi-follower Stackelberg game model is presented where the EH is the leader while users and electric vehicles (EVs) are regarded as two followers. Then, the Kullback-Leibler (KL) divergence-based DRO model is developed to deal with the uncertainty of renewable generation (RG) in the EH. Besides, Karush–Kuhn–Tucker (KKT) conditions, strong duality theory, and big- M approach are combined to transform the bilevel model into a single-level model. The reformulated single-level operation model is incorporated into the KL-based DRO approach. Furthermore, since the crafted column and constraint generation (C&CG) algorithm can prevent possible numerical problems caused by the exponential function and accelerate the solution speed, the crafted C&CG algorithm with linearization for the upper-level slave problem is proposed to iteratively solve the KL-based DRO integrated with Stackelberg game. Finally, numerical case studies are conducted with all simulation results confirming the effectiveness of the proposed model and method.