Detrended and deseasonalized annual emission anomaly and linear fitted trends for 2000–2021 (n = 22). The spread of data represents IAV of corresponding metrics. The central mark and the bottom and top edges of the box indicate the median, and the 25th and 75th percentiles, respectively. The black whiskers extend to the most extreme data points not considered outliers, which are denoted as dots. The black dots are outliers for individual years that are not 2020/2021. The region mask is defined in Supplementary Fig. 2 . GL, Global; NH, Northern Hemisphere; SH, Southern Hemisphere; NAm, North America; SAm, South America; Afr, Africa; NAs, North Asia; SAs, South Asia; SEAs, Southeast Asia.