To achieve optimization operation between the energy hub operator and the load aggregators, this paper proposes a collaborative model for the energy hub and load aggregator considering the carbon intensity-driven and uncertainty. Firstly, a time-coupled carbon emission flow model for energy storage is proposed to cope with the change in carbon intensity at the injection node. Meanwhile, a carbon intensity-based incentive price is developed for interruptible loads with diverse carbon intensity. Secondly, the Wasserstein metric-based distributionally robust optimization method is leveraged to height against the energy hub uncertainty, in which the infinite-dimensional distributionally robust optimization model is reformulated into a tractable linear programming. Thirdly, considering the multi-entity characteristics of energy hub and load aggregators, a collaborative optimization method based on the alternating direction method of multipliers algorithm is constructed. Finally, to meet the convergence requirements of the alternating direction method of multipliers algorithm, the iterative convexification algorithm is applied to transform the nonconvex bi-level energy hub-users model into an iterative convex model by avoiding binary variables and big-M parameters. Consequently, the energy hub and load aggregator can operate cooperatively and protect privacy in uncertain conditions. Simulation results verify the effectiveness of the proposed model and method. Simulation results demonstrate the effectiveness of the proposed model and method. Specifically, the iterative convexification algorithm achieved a profit of ¥1653.82 for the load aggregator, which is 24.4 % higher than the profit obtained by CONOPT (¥1329.73) and 20.3 % higher than that obtained by IPOPT (¥1374.69), thus providing substantial advantages over traditional nonlinear solvers.