To investigate the potential impact of emission reduction measures on ozone (O3) formation under the carbon neutrality target, we examined the changes in O₃ concentration and their sensitivity to various parameters in the urban and suburban areas of the Guangdong-Hong Kong-Macao Greater Bay Area (GBA). In this study, we used the Weather Research and Forecasting (WRF), the Sparse Matrix Operator Kernel Emissions (SMOKE) and the Community Multi-scale Air Quality Modeling system (CMAQ) air quality model to simulate O₃ formation in three key years of 2020, 2030 and 2060, based on the Ambitious-pollution-Neutral-goal scenario data from the Dynamic Projection for Emissions in China (DPEC) model. The decoupled direct method (DDM) module embedded in CMAQ was used to calculate the first-order sensitivity coefficients of O₃ to nitrogen oxides (S_{O3_NOx}) and volatile organic compounds (S_{O3_VOC}). The results show several important trends in the O₃ concentrations and sensitivity. (1) For the changes in O₃ concentrations, in terms of different seasons, the O₃ concentration in the GBA region shows an increasing trend in winter in both 2030 and 2060 compared to 2020. In terms of different cities, the O₃ concentration in Shenzhen shows a significant increasing trend compared to the other cities. (2) For changes in O₃ sensitivity, S_{O3_NOx} shows an increasing trend, with the negative area declining and the positive area increasing. In 2030, the negative absolute value of S_{O3_NOx} is reduced, indicating that the NO_x titration effect will be weakened. In 2060, S_{O3_NOx} becomes positive in most areas of the GBA region. For S_{O3_VOC}, the future scenario shows positive values throughout the study area for all years, but a decreasing trend.