Variation of XCO2 anomaly patterns in the Middle East from OCO-2 satellite data

ABSTRACT The anthropogenic CO2 emission is contributed to the rapid increase in CO2 concentration. In the current study the anthropogenic CO2 emission in the Middle East (ME) is investigated using 6 years column-averaged CO2 dry air mole fraction (XCO2) observation from Orbiting Carbon Observatory-2 (OCO-2) satellite. In this way, the XCO2 anomaly ( XCO2) as the detrended and deseasonalized term of OCO-2XCO2 product, was computed and compared to provide the direct space-based anthropogenic CO2 emission monitoring. As a result, the high positive and negative XCO2 values have corresponded to the major sources such as oil and gas industries, and growing seasons over ME, respectively. Consequently, the Open-source Data Inventory for Anthropogenic CO2 (ODIAC) emission and the gross primary productivity (GPP) were utilized in exploring the XCO2 relation with human and natural driving factors. The results showed the capability of XCO2 maps in detecting CO2 emission fluctuations in defined periods were detectible in daily to annual periods. The simplicity and accuracy of the method in detecting the man-made and natural driving factors including the main industrial areas, megacities, or local changes due to COVID-19 pandemic or geopolitical situations as well as the vegetation absorption and biomass burning is the key point that provides the environmental managers and policymakers with valuable and accessible information to control and ultimately reduce the CO2 emission over critical regions.


Introduction
Compared to the preindustrial era, our planet Earth is showing a sharp increase in atmospheric carbon dioxide (CO 2 ) concentration which is contributing to global warming (IPCC 2014;Ahlström et al. 2015). Carbon dioxide has many natural sources and sinks, while the anthropogenic CO 2 emission particularly fossil fuel (including coal, oil, and natural gas) combustion is attributed as the main source of CO 2 increase (Hakkarainen, Ialongo, and Tamminen 2016;Buchwitz et al. 2021). In addition, the other sources of anthropogenic CO 2 emission such as forest fires, road traffic, oil and gas companies or electricity generations have resulted in spatial and temporal variations which are depending on the time of day or year and can be located across regions of vastly different sizes (Kumar and Nagendra 2016). These variations have prevented climate scientists from gaining a clear enough picture of emission sources to keep track of their emissions (Oda, Maksyutov, and Andres 2018). Thus, to verify the global CO 2 emission and to check the effectiveness of reduction policies (such as the Paris agreement) accurate estimates are a vital need. The primarily CO 2 emission estimation is based on bottom-up inventories using statistical data (Reuter et al. 2020), however, inaccuracies in statistical data, as well as inaccuracies in the use of these data, make error quantification a critical step (Janardanan et al. 2016).
In this way, ground-based observation is the most accurate source of data, but the sparseness of these measurements has strengthened the need for space-based observation to compensate for this limitation. Global distribution of satellite measurements with the acceptable spatial and temporal resolution are the main advantages of space-based observation which are utilized in several studies (Gavrilov et al. 2014;Nguyen et al. 2014;Kountouris et al. 2018;Mousavi, Falahatkar, and Farajzadeh 2017;Falahatkar, Mousavi, and Farajzadeh 2017;Golkar and Shirvani 2020;Mousavi and Falahatkar 2020;Mousavi et al. 2022). Column-averaged dry air mole fraction of CO 2 (XCO 2 ) is the product of satellites including the Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY, from 2002 to 2012) on board the Environmental Satellite (Envisat) with a spatial resolution of about 60 × 30 km at nadir (Houweling et al. 2005;Gloudemans et al. 2006), the Thermal And Near-infrared Sensor for Carbon Observation (TANSO) onboard the Greenhouse Gases Observing Satellite (GOSAT, 2009-present) with a footprint of approximately 10.5 km diameter at nadir (Yokota et al. 2009;Yoshida et al. 2013), and the second Orbiting Carbon Observatory (OCO-2, 2014-present) with a spatial resolution of approximately 3 km 2 Kiel et al. 2019). The above mentioned satellites are grouped in SWIR (short-wavelength infrared) category which are sensitive to near-surface CO 2 concentration (Ishizawa et al. 2016) and contain information on anthropogenic CO 2 emission. Accordingly, in recent years the XCO 2 data has been widely considered in detecting and extracting the anthropogenic CO 2 emission from space. In this way, using SCIAMACHY data, Reuter et al. (2014) analyzed the emission ratio of CO 2 -NO 2 in East Asia. Bovensmann et al. (2010) and Kort et al. (2012) estimated the fossil fuel CO 2 emissions from large point sources such as power plants and localized areas of high emissions such as large cities using GOSAT data. Janardanan et al. (2016) used GOSAT data to employ an atmospheric transport model to attribute XCO 2 to emissions in megacities. Using OCO-2 data in estimating CO 2 emission from individual power plants, and different energy-intensive industrial sources as well as detecting anthropogenic CO 2 emission for megacities such as Los Angeles has been reported in Nassar et al. (2017), Schwandner et al. (2017), Zheng et al. (2020(, Yang et al. (2020, Buchwitz et al. (2021), andWu et al. (2018) studies.
From the above studies, it was concluded that satellite observations with finer spatial resolution have improved the capability of anthropogenic CO 2 emissions monitoring as well as analyzing the localized sources and their small fraction of the background values. In this regard, several recent studies have reported the capacity of OCO-2 XCO 2 data in estimating CO 2 emissions at urban, megacity, and continental scales (Hakkarainen et al. 2019;Buchwitz et al. 2021;Zheng et al. 2020). Thus, in the current study, the OCO-2 XCO 2 dataset (with the finest footprint of 2.25 km and precision of ∼1 ppm ) has been considered in detecting the variation of CO 2 emission in the Middle East (ME). The body of this work is based on a novel and simple method developed by Tamminen (2016, 2019) to isolate the CO 2 emission signals in the study region from the satellite XCO 2 by extracting the background corrections and calculating anomalies in different time scales. The XCO 2 anomalies (DXCO 2 ) are based on computing and subtracting XCO 2 background (CO 2-bg ) values which simultaneously deseasonalized and detrended the original XCO 2 data. As a result, the △XCO 2 (including anthropogenic emissions and seasonal variability related to vegetation and biomass burning (Buchwitz et al. 2021)) with positive and negative values are often associated with critical carbon sources and sinks respectively . However, despite the subtraction of CO 2-bg concentration in DXCO 2 , the combined impact of anthropogenic emissions and environmental variables such as biosphere and biomass burning (He et al. 2020;Shi et al. 2019;Pan, Xu, and Ma 2021) have resulted in various DXCO 2 characteristics in different regions, Such that, a significant decrease in △XCO 2 is reported in the Nile Delta in summertime (Shekhar et al. 2020), while the highest values are seen in the same season in the southwest China (Wang et al. 2021). Regarding the association of DXCO 2 with a complex interaction of human activity as anthropogenic CO 2 emission and environmental variables mainly the biosphere absorption (related to photosynthesis activity), exploring the sources of different variations of DXCO 2 in different regions, requires detailed consideration of these driving factors. Consequently, by taking ME as a research area, the Open Data Inventory for Anthropogenic Carbon Dioxide (ODIAC) related to fossil fuel CO 2 emission (Oda, Maksyutov, and Andres 2018) and the gross primary productivity (GPP) data as the absorptive capacity of vegetation were collected for the study period.
The Middle East is known as one of the main air polluted regions in the world and some countries in this region are among the top 10 CO 2 producers. Despite the importance of this region, only general space-based CO 2 emissions investigations have been conducted by a few global (Hakkarainen et al. 2019;Janardanan et al. 2016) or regional studies such as Yang et al. (2020) and Wu et al. (2018) studies for Mecca and Baghdad cities or Shekhar et al. (2020) study for Nile Delta region. Thus, the purpose of this study is focused on investigating the amplitude and variation patterns of DXCO 2 and its correspond to human and natural driving factors over ME from space including anthropogenic CO 2 emission and biosphere impacts, which is conducted as follows: Section 2 introduces the study area, the datasets, and the method to derive the DXCO 2 and statistical analysis, Section 3, investigating the annual and seasonal patterns of CO 2-bg and DXCO 2 and discussing them in relation with fossil fuel anthropogenic CO 2 emission and vegetation cover, and Sections 4, the main conclusions.

Study area
The Middle East region is located in the west of Asia and includes sixteen countries as Iran, Turkey, Egypt, Saudi Arabia, Iraq, Syria, Emirates, Kuwait, Qatar, Lebanon, Bahrain, Israel, Jordan, Oman, Palestine, and Yemen ( Figure 1). The ME has a wild range of climate types but mainly arid and semi-arid with mild winter and low precipitation (Nandi, Bhattacharjee, and Chakraborty 2019). The population of ME is about 447.8 million, with an average of 1.65% population growth in 2020 (WDI 2021). Approximately more than 46% and 40% of the world's proven oil and gas reserves are located in the ME region, respectively (BP 2021), and more than 8% of the total world's anthropogenic CO 2 emissions are emitted from this area (Edgar 2020).

OCO-2XCO 2 observations
NASA's OCO-2 is the first spacecraft devoted to CO 2 monitoring, which was successfully launched in July 2014 . OCO-2 measures the CO 2 mole fraction in column-averaged dry air with the accuracy, resolution, and coverage needed to identify CO 2 sources and sinks on a regional and global scale (Crisp et al. 2004). This satellite operates in a Sun-synchronous orbit that flies at 705 km altitude with a local overpass time of 13:30 . OCO-2 have a repeat cycle of 16 days and complete global XCO 2 coverage twice per month, with a spatial resolution of approximately 3 km 2 . The validation results of OCO-2 data against high-precision measurements (ground-based or aircraft observations) indicate that precision is less than 0.5 ppm (Bi et al. 2018;Wang et al. 2021). In this research, the bias-corrected OCO-2 daily lite files (V10) from January 2015 to December 2020, with the quality flag as '0' were utilized.

ODIAC CO 2 emissions
The monthly fossil fuel CO 2 emission from the open-source Data Inventory for Anthropogenic CO 2 (ODIAC) for January 2003 to December 2019 has been utilized to be compared with DXCO 2 variations. ODIAC is a monthly global anthropogenic CO 2 emission dataset based on estimates made by the Carbon Dioxide Information and Analysis Center (CDIAC) and satellite-observed nightlights and the power plant profile dataset, with a spatial resolution of 1 km × 1 km (Oda, Maksyutov, and Andres 2018). The climatology mean (CM) of annual ODIAC CO 2 emissions from 2003 to 2019 over ME is presented in Figure 1(a).

Vegetation absorption
Gross primary production (GPP) is a key measure of carbon mass flux in carbon cycle research, which is the total carbon fixation by terrestrial ecosystems through the photosynthesis of vegetation (Wang et al. 2012). In this study, GPP from Moderate-Resolution Imaging Spectroradiometer (MODIS; MOD17A2H) with a spatial resolution of 1 km × 1 km has been utilized as a parameter related to carbon uptake by vegetation photosynthesis (Figure 1b shows the CM of annual GPP in ME). The main vegetation cover in ME is observed in the west and north of Syria, Iran, and Iraq, along the Nile River in Egypt, and most parts of Turkey, with the maximum (minimum) average index in spring (winter).

Method
For extracting the CO 2 emission information from XCO 2 , subtracting XCO 2-bg , to isolate the CO 2 source and sink signals, was considered based on the method described in Hakkarainen et al. (2019) study (Equation 1). The sensitivity of the method to the area for extracting medians, motivated the authors to calculate and compare the backgrounds in two different defined areas, the ME domain and the latitude bands ( Figure 1). The day-to-day median from daily OCO-2XCO 2 in the ME domain was calculated in three zones with 10 • latitude widths from 11°to 22°(Zone I), 22°to 33°( Zone II) and 33°to 44°(Zone III), the daily median was also considered in three latitude bands from 11°to 22°(LB I), 22°to 33°(LB II) and 33°to 44°(LB III). The XCO 2-bg values in Zone I from 2015 to 2020 are presented in Figure 2. The XCO 2-bg in defined areas (ME domain and latitude bands) are compared in Section 3.1, and the daily DXCO 2 are then extracted from the individual XCO 2 . Afterwards, the spatial daily DXCO 2 distribution maps with 0.25°× 0.25°grid cells were generated over ME using a space-time kriging prediction. The monthly, seasonal and annual DXCO 2 distribution maps to detect the spatial-temporal variations have been then extracted from daily maps.
The non-parametric Mann-Kendall trend test (MK) (Mann 1945;Kendall 1975) in analyzing time series trends and Mann-Whitney change point methods (MW) in detecting change points (CP) based on the differences in mean value (Ross 2015) are utilized for more detailed study in DXCO 2 , GPP and ODIAC CO 2 emission time series. The assumption that the observations are independent (Von Storch and Zwiers 2001) is tested by the autocorrelation function (ACF) and if the time series of data are serially correlated, uncorrelated data is utilized by applying the pre-whitened data series (Shirvani, Arpe, and Jahandideh 2020). The autocorrelation function is the coefficient of correlation between two values in a time series, which describes the autocorrolation between an observation and the same observation at a prior time step (Box, Jenkins, and Reinsel 1994). Given measurements Y 1 , Y 2 , … , Y N at time X 1 , X 2 , … , X N , the lag k autocorrelation function with confidence bands of ( + 1.96/ N √ ) is defined as (Equation 2): 3. Results and discussions

XCO 2 backgrounds
The time series of the monthly mean of CO 2-bg shows a strong seasonal cycle with an increasing trend (varying from 393 to 418 ppm) in all three zones as well as the whole ME domain and LBs (Figure 3). At the same time, the differences of CO 2-bg in ME (whole LB of 11°to 44°) from the CO 2-bg of study zones (LBs) were computed (Figures 4a and b). This comparison clearly shows the diversity in the CO 2-bg of Zone III against the other two zones, so that from April to September, the CO 2-bg shows fewer values than ME (positive values), and higher values in other months (negative values). However, Zone II shows a consistently higher CO 2-bg with the maximum values mostly in December and January and the minimum values in July and August. On the other hand, a similar comparison in LBs (Figure 4b) shows the same results, but with smoother fluctuations in LB II and LB III. In addition, the difference of the LBs CO 2-bg against the corresponding zones ( Figure 5) shows that in general the LBs CO 2-bg have higher values than the corresponding zones. However, the graphs show an annual cycle with mostly positive values (referring to higher values of LBs CO 2bg against the corresponding zones) especially in summertime (July and August) except for LB III vs. Zone III, which shows this maximum positive value almost in springtime (April and May). In addition, LB I vs. Zone I and then LB III vs. Zone III show the highest positive values particularly in August and April, respectively. However, the sharp negative values are seen in May in LB I vs. Zone I. Considering Figure 5 with the highest values in warm seasons (which is expected to have the lowest CO 2-bg based on Figure 3) and the lowest values in cold months (with the expectation of the highest CO 2-bg ) or May in LB I vs. Zone I, guided the authors to believe the CO 2-bg which are extracted from the study zones can distinguish the sources of emission more clearly than the CO 2-bg in LBs. Thus, the further results will be based on the CO 2-bg values which are calculated over the three defined study zones. Furthermore, the sensitivity of CO 2-bg choice in local and regional studies is a critical step which is mentioned in previous studies in China or South Africa Hakkarainen et al. 2019;Buchwitz et al. 2021;Fu et al. 2021).

XCO 2 anomalies
For investigating the daily DXCO 2 , the averages of daily DXCO 2 points were computed over the defined zones. The maximum and minimum in daily mean DXCO 2 shows both the most negative and positive values in Zone III for each month, while the fluctuations are the least in Zone II. The high fluctuation in Zone III can be explained by vegetation carbon sink ( Figure 1b) and precipitation at especial periods which intensify the negative DXCO 2 values in Zone III, while large population, cold winter and utilizing heating systems and fossil fuels as well as industrial activities are the sources of CO 2 discharge. In addition, the maximum daily DXCO 2 values in Zones I and II are mostly recorded in cold months, while in Zone III the maximum DXCO 2 values are seen in both July and August as well as the cold months. Regarding the continuous days with relatively the highest DXCO 2 values in June, July or August in eastern Turkey and northwestern Iran in relation to high vegetation cover strengthens the possibility of biomass burning and its effect on high DXCO 2 values in the late springtime (June) and early/mid summertime (July and August). However, wildfire reports in 2015 and 2019 in northwestern to northeastern Iran and 2018-2020 in east, north and northeast Turkey emphasize this suspicion. In this regard, the importance of fire season in Africa is considered in DXCO 2 patterns in Hakkarainen et al. (2019) study. In addition, the westerly prevailing winds in this zone seem to play an important role in the general accumulation of carbon dioxide in the eastern regions of these countries. Nevertheless, in general the minimum DXCO 2 values are often seen in spring and summer. In the following, the spatial mean in 0.25°× 0.25°grids as explained in Section 2.3 was considered for monthly, seasonal and annual periods. This leads to a spatial-temporal analysis of CO 2 sources and sinks over ME in defined periods as follows: The CM of DXCO 2 for March and April as the months with approximately the most increase in vegetation cover show mostly the negative DXCO 2 values in eastern Iraq, Turkey, southern Oman and northern and northwestern Iran as well as Egypt, while the center and northeast of Iran as well as southern and southwestern of Saudi Arabia, show relatively positive DXCO 2 values (see Figure 6a for April). This is while in cold months by decreasing the role of photosynthesis and increasing the fossil fuel use, positive values are highlighted in Turkey and northwest Iran, however, Iraq and southwest Saudi Arabia, as well as southwest Iran, are the widest areas with constant positive DXCO 2 values during the year (see Figure 6b for December). These vast areas with constant positive DXCO 2 refer to the role of oil and gas extraction facilities as well as the role of megacities and the geopolitical situation of countries such as Iraq and Syria.
The CM for seasonal DXCO 2 shows the spatial displacement of positive and negative DXCO 2 values in different seasons, particularly in Turkey, western Iran, Iraq and Egypt. In this way, comparing the DXCO 2 in winter and spring shows a clear increase in the number of grids with positive anomalies in winter (Figure 7) for the mentioned areas. However, in Figure 8(a) the CM for annual DXCO 2 shows the highest values in Iraq as well as south regions of the Persian Gulf and south and southwest of Saudi Arabia, while despite the positive DXCO 2 values in Iran, generally less DXCO 2 values in comparison to the above mentioned areas are seen. At the same time, besides the areas with denser vegetation such as Turkey as well as north and northwest of Iran (Figure 1b), in the most defined periods, the minimum DXCO 2 or even the negative values are shown in Egypt, this is while by excluding the Nile Delta, the rest of the country is mostly classified as desert. However,  more DXCO 2 values in Nile Delta in comparison to the surrounding desert go back to urbanization, dense population and irrigated farming which were mentioned in Shekhar et al. (2020) study. The latter reason refers to winter and summer as the peak of growing seasons which are receiving a high amount of carbon-rich water (Shekhar et al. 2020). This carbon-rich water will be then recycled and as a result causes CO 2 emission and thus higher DXCO 2 values in this region particularly in wintertime than in the surrounding desert (see Figure 7b). Recognizing the distinction in the Nile Delta emission by DXCO 2 in comparison to the surrounding area, as shown in Figures 6 and 7, clarifies the spatial capability and accuracy of the median method in determining the CO 2-bg and thus CO 2 emission, as resulted in Hakkarainen et al. (2019) study. However, the importance of CO 2 background choice in recognizing the local and regional emissions is mentioned in Buchwitz et al. (2021) study. In this regard, natural driving factors such as prevailing wind speed and direction as well as topography can be considered in a more detailed local investigation (Fu et al. 2021).
In the ME, deserts cover large areas including the Arabian Desert in Saudi Arabia, Dasht e Lut and Dasht e Kavir in the central part of Iran, which generally show less DXCO 2 values than the surrounding areas (Figures 6-8). However, in these countries, considering the significant role of wind (Fu et al. 2021), the transfer of CO 2 by the wind from major emission sources to surrounding areas reduces the DXCO 2 difference between deserts and the rest of the country.  On the other hand, as a general effort to detect the changes in CO 2 emissions related to the COVID-19 pandemic, the DXCO 2 in March and April 2020 as the period with widespread lockdown were compared to the same month in the previous year. The results showed an increase in the number of pixels with less DXCO 2 values in March in comparison to the same month in 2019. In addition, for megacities such as Tehran, Istanbul, Ankara, Baghdad and Abu Dhabi in April 2020 relatively less DXCO 2 values were seen, compared to the same month in 2019 (Table  1). Figure  Nevertheless, detailed information about the anthropogenic as well as natural CO 2 surface fluxes for more accurate and sophisticated analysis is suggested (Buchwitz et al. 2021).
As described in the above paragraphs the combined effect of both human activities including the oil and gas industry, urbanization, population and transportation as well as natural features such as vegetation cover, is considerable and diagnostic in the spatial-temporal distribution of DXCO 2 . So that in the following section, anthropogenic CO 2 emissions from ODIAC datasets and GPP is considered to further clarify the DXCO 2 variation in ME.

The DXCO 2 interaction with natural and anthropogenic driving factors
The CM for annual DXCO 2 (Figure 8a) comparison with the CM for annual anthropogenic CO 2 emission from the ODIAC dataset (Figure 8b) in the study period shows acceptable coordination between the regions with relatively the highest DXCO 2 values and the corresponding pixels with the highest anthropogenic emissions. On the other hand, the inter-annual DXCO 2  fluctuations over ME (Figures 6 and 7) particularly in Turkey, the Nile Delta region, the east of Iraq and the north and northwestern Iran emphasize on natural driving factors. In this way, the time series of monthly total ODIAC CO 2 emission, as well as the monthly mean of GPP in  Saudi Arabia, Iran, Iraq, Turkey, Egypt and the regions of Syria, Lebanon, Israel and Palestine (as the countries with a complex interaction of geopolitical situations, oil and gas industry or vegetation cover), were extracted to compare the natural and anthropogenic driving factors against DXCO 2 fluctuations. In general, no significant trend is detected in DXCO 2 and GPP in mentioned countries as well as the whole ME area based on the MK test, while the ODIAC CO 2 emission shows a general upward trend. At the same time, the inverse relationship between DXCO 2 and ODIAC CO 2 emission with GPP is mostly observed in Turkey, Iran, Iraq and the regions of Syria, Lebanon, Israel and Palestine, while no clear inverse relationship is considered in Saudi Arabia and Egypt (with poor vegetation cover) particularly in GPP and ODIAC CO 2 emission comparison. However, a general one to two months lag is seen in a monthly maximum of GPP values versus monthly minimum values of DXCO 2 and ODIAC CO 2 emission (see Figure 10 for Turkey). The ODIAC CO 2 emission seems to show two breaks in the general upward trend of Saudi Arabia in January 2007 and January 2016 respectively (Figure 11(. The monthly total ODIAC CO 2 emission versus the monthly mean of DXCO 2 and monthly mean of GPP for Saudi Arabia are shown in Figures 12(a and b), respectively. The short period in DXCO 2 values (available since 2015) in Figure 13(a), does not provide sufficient evidence in explaining these sharp changes. On the other hand, in Figure 12(b) despite the equal period of two datasets, the GPP time series with no significant upward nor downward trend does not explain the ODIAC CO 2 emission breaks. Accordingly, the proceedings such as pollutant filtering of oil and gas extraction industry, improving the public transportation, and changes in cooling facilities as well as planning to plant trees, probably is the reason for changes in ODIAC CO 2 emission changes since 2016. However, what is relevant to the present study is the insufficient length of DXCO 2 in justifying the changes in the anthropogenic CO 2 emission in Saudi Arabia. At the same time, CP detection in DXCO 2 , GPP and ODIAC CO 2 emission time series are considered in ME as well as the mentioned countries. The ODIAC CO 2 emission shows the only significant CP in December 2015 for Saudi Arabia and January 2016 for ME. This is while the CP in GPP was detected only in Saudi Arabia (October 2015), and no significant CP in DXCO 2 was detected in studied countries as well as the whole ME area. The almost Figure 11. The monthly total ODIAC CO 2 emission in some selected countries in ME.
simultaneous CP in ME and Saudi Arabia for ODIAC CO 2 emission can refer to the significant and critical role of some industrial regions or countries in anthropogenic CO 2 emission, however, no significant upward nor downward trend in DXCO 2 emphasizes on an overall moderating effect of some driving factors against the intensifying effect of other factors in CO 2 emission.
However, for a more detailed study the regions with relatively the highest DXCO 2 values (corresponded with the highest ODIAC CO 2 emissions (Figure 8)) were considered as individual polygons independent of political boundaries. Figure 13 shows the approximate location of selected polygons which are named Polygon I to VI. The time series of the monthly mean of CO 2 emissions from the ODIAC dataset is presented in Figure 14(a), which shows the highest level of the monthly mean of ODIAC CO 2 emission in Polygon VI (the southern margin of the Persian Gulf and Oman Sea). In this polygon, despite the small portion of the area in ME and as a result the lower total emissions, the monthly mean of CO 2 emissions show a clear higher value than other polygons. In addition, in ODIAC CO 2 emissions the only CP was detected in January 2015 in Polygon VI, while DXCO 2 shows sharp fluctuations in this polygon (Figure 14b). This confirms the capability of DXCO 2 values to detect the regions with high emissions and as a result detect the carbon cycle which is contributed to perturbation in CO 2 concentrations.

Conclusions
In our study using 6 years of OCO-2XCO 2 product, the XCO 2 anomalies over the Middle East were computed based on the daily medians in each defined zones. The daily DXCO 2 values varied in the range of −4-5.5 ppm which were smoothed in monthly, seasonal or yearly time scales. The results showed the capability of space-based observation and as a result DXCO 2 in spatial and temporal analysis of regional CO 2 sources and sinks, including mainly the anthropogenic CO 2 emission as well as biomass burning and vegetation absorption.
The DXCO 2 variation pattern over ME detected the anthropogenic CO 2 emission mainly in areas with the oil and gas industry as well as populated areas and megacities. Furthermore, biomass burning in some regions of Turkey and Iran was detectible in months with high vegetation cover along with high air temperature and arid situation. In addition, the reduction in CO 2 emission during widespread lockdown in March and April 2020 was recognized in densely populated cities such as Tehran, Istanbul, Ankara, Baghdad and Abu Dhabi based on the DXCO 2 comparison with the same months in 2019. The impact of biomass absorption and photosynthesis activity as CO 2 sinks are highlighted in seasonal patterns of DXCO 2 , particularly in Zone III with higher vegetation cover. In this regard, small countries such as Qatar and the United Arab Emirates or large countries with relatively small populations such as Saudi Arabia are recognized as the regions with constant and high DXCO 2 values, while the Nile Delta region, Figure 13. The six selected polygons in ME.
Turkey, north and northwestern Iran, and some parts of Iraq with high vegetation cover show more seasonal fluctuations in CO 2 anomalies.
The high consistency between DXCO 2 variation patterns and ODIAC as well as GPP datasets as fossil fuel CO 2 emission and biospheric contribution emphasizes on the spatial-temporal capability of the CO 2 anomaly method for detecting anthropogenic CO 2 emission patterns. However, according to the regional variation patterns of DXCO 2 and regional driving factors, additional study considering wind patterns along with topography is suggested.

Data availability statement
The data that support the findings of this study are available from the corresponding author upon reasonable request.