<p>Previous study has revealed the formation of the ozone valley over the Tibetan Plateau (TP) from a dynamical perspective. A realistic simulation of the Tibetan Plateau ozone valley is a prerequisite for a better understanding and prediction of the surrounding climate. Using the ERA5 reanalysis as the baseline, this study evaluates the performance of the ozone valley over the TP from CMIP6 models. Models simulate the TP ozone valley with different degrees of success and diverge in the ozone valley intensity and extent. Further, the TP ozone valley shows dual maximum variability centers across Iranian Plateau and Tibetan Plateau. High-skill models are selected to diagnose the balance of ozone transport over the TP using the Lorenz circulation decomposition. The zonal mean ozone stationary transport by the eddy (SFu3) contribute most to maintain the low ozone content center in both the high-skill models and the ERA5 reanalysis. In contrast, for models unreasonably simulating the ozone valley, ozone zonal deviation stationary transport flux by the zonal winds (SFu2) dominates, which is unfavorable for formation of the ozone low value center. In the boreal summer, the TP ozone content is strongly coupled with the giant anticyclonic circulation system, South Asian High (SAH). The intensity, area, and eastward-extending ridge point of the SAH is negatively correlated with the TP ozone content. Models better simulating the multiple facets of the SAH tend to reproduce a more realistic TP ozone valley. Additionally, climate-chemistry coupled models perform better than the remaining models in describing various parameters of this low ozone content center.</p>