Estimation of emissions from crop residue burning in Türkiye using remotely sensed data and the Google Earth Engine platform

Abstract This research aims to find crop-based residue burning (CRB) amounts for the Southeastern Anatolia Region of Türkiye, where versatile and diverse of agricultural fields are available. We used the time series of Sentinel-2 images to identify spatial distribution of different crop types and to determine the burned agricultural fields in 2019 using Google Earth Engine platform. We implemented Intergovernmental Panel on Climate Change standards to calculate the emissions from Greenhouse gases and Particulate Matter. We also analyzed Sentinel-5P images to generate Nitrogen Dioxide maps and compared them with crop-based emission calculations. Our analysis illustrated that 14,444.307 Gg of Greenhouse Gases and 117.809 Gg of Particulate Matters were released in 2019 due to CRB practices. Our results can potentially improve the national statistics by providing spatio-temporal field-based CRB emissions and supporting country-level agricultural decision-making processes. Researchers could implement our proposed data and methodology in other regions suffering from the CRB problem.


Introduction
Agricultural production provides an essential resource to the economies of the world countries.Countries examine farming practices and create annual evaluation reports to develop effective and sustainable agricultural policies.Various communities address agricultural activity problems yearly, and goals are set to tackle these problems.One of the leading agricultural problems also reported in T€ urkiye is crop residue burning (CRB).CRB is the burning of agricultural wastes on agricultural land after the harvest season.Farmers prefer CRB since it is fast, supports the timely sowing of the next crop, and increases yield by killing emphasized the impacts of different values used for the emissions calculations resulting in different weights.They underlined the importance of identifying local values.Bahsi et al. (2019) used Sentinel-2A images for a small pilot region of 15 Â 15 km to generate crop maps and find out burned agricultural parcels to calculate total emissions of CO 2 , CO, CH 4 , NMHC, N 2 O, NH 3 , SO 2 , and NO X .McCarty (2011) reported emission factors of different crop types caused by CRB.Glushkov et al. (2021) used Sentinel-2A images obtained between 1 January and 15 May 2020 to map the burned areas in Russia accurately and generated a connection between fires and different land practices such as land abandonment, fires in natural areas, and burning of agricultural fields.They found that spring fires were mainly associated with grasslands and agricultural lands.Kalluri et al. (2020) used multiple satellite data and evaluated the effects of biomass burning on atmospheric aerosol loading.They calculated spatial and temporal atmospheric aerosols for the 2006-2017 period using Aerosol optical depth (AOD), fire pixel counts from MODIS, carbon monoxide (CO) from Measurements of Pollution in The Troposphere (MOPITT), aerosol index (AI) and tropospheric NO 2 from Ozone Monitoring Instrument (OMI), and aerosol layer depths from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO).They concluded that agricultural activity was the main factor for the high AOD values with the support of multi-sensor satellite data Lee et al. (2021) investigated the impact of crop burning in northeastern China on the PM 2.5 concentration in South Korea.Saxena et al. (2021) assessed the effects of PM 10 , PM 2.5 , NO 2 , and SO 2 , emitted during CRB activities in Haryana, on the air quality of Delhi.Shen et al. (2021) estimated the field-level NOx emissions from crop residue burning using remote sensing data in Hubei, China.Singh et al. (2021) investigated the potential of Sentinel-2A imageries in the spatio-temporal assessment of rice residue burnt areas and generated pollutants in Kurukshetra, districts of Haryana, India.Zhang et al. (2022) compared the ability of burned area products (MCD64A1, FireCCI 5.1, and the Copernicus Burnt Area) to detect crop residue burning in China.It is crucial to accurately determine the spatio-temporal distribution of CRB due to its adverse impacts on humans and the environment, such as degradation in air quality, threats to human health and increased risk of traffic accidents.Different sensors provide different advantages in spatial, spectral and temporal resolution; therefore, in this research, we aim to benefit from additional satellite data and products and mapping approaches to identify CRB-affected fields accurately.
The scientific questions addressed in this research are: (1) Can we use current remote sensing data and methods to estimate crop-based CRB emissions?(2) Can we generate accurate crop type maps based on phenological cycles of different crops?(3) What would be the contributions of new Sentinel 5-P data for the traceability of CRB? (4) How can we benefit from the Google Earth Engine (GEE) platform to generate transferable results for CRB calculations from multi-sensor data?

Study area
Our study area lies in the South-eastern Anatolia Region of T€ urkiye, consisting of Kilis, Gaziantep, Adiyaman, Sanliurfa, Mardin, Diyarbakır, Batman, Sirnak, and Siirt provinces (Figure 1).It covers the Euphrates-Tigris Basin and upper Mesopotamia plains, including the region of the Southeast Anatolia Project (with the Turkish acronym of GAP).The GAP region involves 22 dams, 19 hydroelectric power plants (HEPPs), and a total irrigation area of 1.8 million ha.The GAP project is vital with its extensive contributions to the economy, agricultural productivity, and new employment opportunities (URL 1, Sevinc ¸et al. 2019).Approximately 4.5 million ha of land has been cultivated in the GAP region, which is 11.4% of the total cultivated area of T€ urkiye.Moreover, 55.12% of the total cotton and 32.48% of the whole wheat have been obtained from the GAP region (URL 2).The main crops grown in the region are wheat, wheat-barley, cotton, maize and others from national sources (Table 1, URL 3).
Although the land cover mainly comprises agricultural lands, it also includes mountains, small urban areas, and wetlands.There are two main reasons for the selection of this region.The first is the increase in agricultural activities in the region with the GAP project implemented in 1989.Unfortunately, the increase in agricultural production triggers the application frequency of CRB.The second is that the agricultural diversity of the area is abundant.We believe that the emission results obtained in the study would provide an advantage in evaluating different crops of some national statistics (TSI 2020) and has the potential to support decision-making processes.Crop Residue Burning (CRB) is one of the most common methods used worldwide, although it was forbidden officially.This practice has been employed not only in T€ urkiye (Figure 2) but also in other countries such as the United States (McCarty et al. 2007), China (Gao et al. 2017) and India (Singh et al. 2021).Farmers prefer CRB since it is fast and labour-free, and they also think it is an effective method to control pests and weeds and increase fertility.However, there are varying adverse effects of CRB locally and globally.Un-harvested land and healthy vegetation areas within CRB areas have been negatively affected, and traffic accidents occurred due to the smoke during the fire.Moreover, CRB results in the loss of organic matter on the topsoil, which is one of the main factors affecting the yield of agriculture, and it is essential to preserve a high amount of organic matter in the soil (C ¸ıtak et al. 2006).

Data
We used two different Sentinel satellites in this study, namely Sentinel-2 and Sentinel-5P belonging to Earth observation missions from the Copernicus Programme.Sentinel-2 satellites include two identical satellites, which are Sentinel-2A (launched on 23 June 2015) and Sentinel-2B (launched on 23 June 2015).These satellites provide 10 to 60 m spatial resolution data for 13 spectral bands within 5-day periods (ESA 2020).Most of the agricultural areas in T€ urkiye are fragmented and small; therefore, high spatial resolution is necessary to determine the boundaries of the agricultural parcels accurately.Another factor in the selection of the satellite is the temporal resolution since the burning effect of a field exposed to CRB disappears in different time periods.In our study area, for some of the agricultural parcels, CRB effect disappeared within 15 days; however, there are several cases in different parts of the world in which shorter time windows occurred even up to one day (Vadrevu et al. 2022).Therefore, frequent satellite images obtained from multi-sensor system are very important for CRB research.Afterwards, the related agricultural parcel seems like a traditional harvested area; therefore, satellites having a temporal resolution of one week or shorter would be desirable to find the CRB-affected regions precisely.We examined the spectral change of some agricultural parcels exposed to CRB during and after the practice using Sentinel-2A and 2B data (Figure 3).We observed that up to ten days after the residue burning, there are still indications of burning in the field represented by low spectral reflectance values.The image captured after the fifteenth day of CRB exhibits similar reflectance with an unburned field, making CRB detection impossible.Moreover, the cultivation and harvesting periods of crops vary for some regions and farmers, which affects the CRB application date.Therefore, frequent monitoring of agricultural areas is essential.
The high spectral resolution of Sentinel-2 allows the opportunity to observe the crops' phenological stages and use various remote sensing indices such as Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Normalized Burned Ration Index (NBR), and difference of NBR (dNBR) (Table 2) (Bahsi et al. 2019;Garc ıa-Llamas et al. 2019).We used Level 2 A products to minimize the geometric and atmospheric distortions (ESA 2020).Copernicus program has provided Level 2 A images globally since December 2018, in which atmospheric and radiometric corrections were implemented to eliminate or minimize the scattering, absorption, cloud, and haze effects.Level-2A products are also geometrically corrected and presented in a projection system (ESA 2020; Aksoy et al. 2022).We also used the Sentinel-5 Precursor (P) satellite data in this research.Sentinel-5P is the first Copernicus mission, developed to monitor the atmosphere with its TROPOspheric Monitoring Instrument (TROPOMI) instrument.TROPOMI measures the ultraviolet and visible (270-495 nm), near-infrared (675-775 nm) and shortwave infrared (2305-2385 nm) light with A temporal resolution of less than one day (URL 6).It was launched on 13 October 2017, and it became fully operational on 5 March 2019 after the commissioning phase ended (URL 7).The native spatial resolution of Sentinel-5P is 5.5 Ã 3.5 km.However, we used NO 2 productS of Sentinel-5P available in the GEE platform, and this product is resampled to approximately 1 Ã 1 km (URL 8).Sentinel-5P Level 2 products are streamed as the near-real-time (NRTI), the Non-Time Critical or Offline (OFFL) and the Reprocessing (RPRO).NRTI data becomes available within 3 h after data acquisition and is widely used for rapid operational processing (URL 9).We used Near Real-Time NO 2 Level 3 (NRTI/L3_NO 2 ) data available in the GEE platform, including gridded high-resolution imagery of NO 2 concentrations (URL 10).In this study, we generated Sentinel-5P NO 2 density maps, showing the total vertical column of NO 2 (ratio of the slant column density of NO 2 and the total air mass factor) (URL 11).
We visually analysed the active fire data from NASA's Fire Information for Resource Management System (FIRMS), consisting of data from the NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the Aqua and Terra satellites and the Visible Infrared Imaging Radiometer Suite (VIIRS) aboard S-NPP and NOAA 20 (URL 14; Schroeder et al. 2014;Giglio et al. 2016).We evaluated the compliancy of FIRMS data with Sentinel-5P and burned agricultural fields that we identified from Sentinel-2 data.

Methods
We conducted four steps in this study in the GEE platform: pre-possessing, classification, emission calculation, and satellite-based emission (Figure 4).First, we generated a complete mosaic image of the study area for August 2019.Then, we identified water and nonagricultural areas using various remote sensing indices and masked these regions from the mosaic image.After the masking, we used multi-temporal Sentinel-2 images selected based on the phenologic cycles of different crops and implemented the Support Vector Machines (SVM) classification (Vapnik et al. 1997).As a result, we created a crop type map consisting of maize, cotton, wheat-barley, and other classes.Besides, we used an extra set of Sentinel-2 satellite images, selected based on the harvest dates of different crop types to identify burned areas.Afterwards, we determined the burned agricultural lands in the region for 2019 using the dNBR index.Integrated usage of crop type maps and burned area maps let us conduct crop-specific emission calculations for the agricultural parcels exposed to residue burning.We employed the emission calculations following the Intergovernmental Panel on Climate Change (IPCC) standards.We estimated the emissions per crop and per province as spatial extent.Finally, to better analyse the emission calculation, we created maps of the NO 2 product of the Sentinel-5P satellite in the GEE platform and compared these maps with calculated values.3.1.Pre-processing In the pre-processing, we first created the study area's ortho-mosaic image covering 31 Sentinel-2A/B image frames for August 2019 (Figure 3).We then masked the non-agricultural areas using NDWI and NDVI indices to deal with only agricultural areas for the later stages and minimize uncertainties.First, we applied the NDVI index to all images for the April-November period and performed the visual analysis with different threshold values.As a result, we masked the non-agricultural areas based on a threshold value of 0.65, which was determined by evaluating the average NDVI value for well-known agricultural fields.Spectral reflectance values vary for different land categories; data acquisition conditions and diversity of geographical regions also impact the NDVI threshold.We employed the second masking process for water areas.We applied NDWI to images in a similar period and used the accepted value of 0.30 in the literature as the threshold value (Topalo glu et al. 2021).Thus, we masked natural water bodies and waterways in the study area.We generated the crop type map after masking the water and non-agricultural areas.

Classification
Pixel-based (supervised and unsupervised) classification methods have been widely used to extract meaningful information from remote sensing data (such as crop type maps, yield maps, land cover and land use maps) (Garc ıa-Llamas et al. 2019).Recently, machine learning methods have been used more frequently in the literature due to their ability to produce more accurate results than traditional methods (Demidova et al. 2019;Aksoy et al. 2022).We used the SVM method for the classification, which provides more accurate results than other pixel -based classification methods and requires very few samples (Mountrakis et al. 2011;Alganci et al. 2013).As with most machine learning methods, the SVM method uses unique training sets to adjust its parameters to solve the uncertainties in the classification.Thus, it learns how to classify unknown samples from the training data set.SVM is widely used in agricultural areas that are likely to show similar reflections since it makes classes separable (Mucherino et al. 2009;Alganci et al. 2013).
Besides, SVM algorithms use mathematical parameters defined as kernels.The kernel takes the data as input and converts it into the required form.It uses different kernel types for different SVM algorithms.Linear, nonlinear, polynomial, sigmoid, and radial basis functions (RBF) are examples different kernel types.We used RBF in this study with gamma and a cost value of 128, and 64, respectively.After applying the grid search function in the GEE platform, we found these values, a hyperparameter optimization technique (Aksoy et al. 2022).Cross-validation methods are such as Jack Knife Test, Boot-strapping, Monte Carlo Test, k-fold etc (Aksoy et al. 2022).We implemented the 5-fold cross-validation method to analyse the training and validation accuracy which divides the data set into two separate groups, training and validation, to quantify the accuracy of the classification (Saud et al. 2020).

Crop type mapping
The main crops grown in the region are maize, cotton, wheat, and barley (Bahsi et al. 2019).We carried out the classification considering the phenological stages of these crops.We specifically analyzed different phenological stages of maize and cotton while constructing our methodological steps and selecting the data acquisition times for satellite images.The first crucial phenological phase of the maize, tasselling stage, occurs before the flowering stage and maize exhibits higher spectral reflectance at this phase than other crops.As a result we used the satellite images of the tasselling stage for the classification.However, maize was planted in the region in two periods: the main and the second periods.Therefore, it is crucial to capture highly vegetated maize conditions for both periods, which also vary in different provinces of the region from north to south.We found appropriate times for each province in the literature for these two periods (Alganci et al. 2013;URL 12).In this context, we used July images to determine the main maize crop and mid-September images for the second maize crop.Besides, the tasselling phase is generally parallel to the pesticide phase that enables cotton to come out of the cones.Thus, cotton fields with pesticides can be easily distinguished from other crops due to the their whitish appearance on satellite images.Wheat-barley phenology is the opposite of the phenology of maize and cotton.Wheat and barley classes appear as empty fields when maize and cotton are cultivated.Finally, we classified the remaining pixels as the other class, those which do not exhibit any of the spectral characteristics mentioned above.
Considering the phenological characteristics, for each crop, we selected 25 samples from provinces with small areas and 50 from regions with larger provinces.Then, we produced a crop type map using the SVM method with a total of 1775 training data sets.We also conducted an independent accuracy assessment using test polygons for each class.

Burned area mapping
Determining the burn severity, defined as the magnitude of the ecological change caused by the fire (Petropoulos and Islam 2017), enables a more sensitive analysis of the burned area.Most satellite-based burn severity studies use methods based on remote sensing indices because they are easy to calculate and simple to implement.Spectral indices using the NIR and Short-Wave Infrared (SWIR) bands, such as Normalized Burn Ratio (NBR) and its bi-temporal version, the difference in Normalized Burn Ratio (dNBR), were used to detect the burned areas (Keeley 2009;Garc ıa-Llamas et al. 2019).
We generated NBR images for the selected periods and the formed dNBR images by taking the difference of May and November NBR images.We then spotted the burnt agricultural fields by taking the minimum NBR index value from the different days in May and created a complete pre-NBR mosaic image of the region.Then, we took the maximum NBR index value for each pixel from the different days in November and created a complete post-NBR mosaic image.Finally, we generated the dNBR index image by taking the difference between pre-NBR and post-NBR images to find the CRB exposed agricultural parcels and implemented threshold values widely accepted in the literature (Fang and Yang 2014).
We performed an area-based accuracy assessment to evaluate the accuracy of the generated burned area map.First, we divided the burned area map into grids of 1 hectare (100 m Â 100 m).Then, we randomly selected 200 grids, purely representing the related classes.We used the post-NBR image as a reference for visual comparison.

Emission calculation
We overlapped the burned area map and the crop type map to determine the crop type of burned agricultural fields.Then, we calculated crop-based emissions from CRB using Equations 1 and 2 proposed by the IPCC (IPCC 2020).
Burned area (A) (decare) and Yield (Y) (kg/decare) information of each crop for 2019 were taken from the burned area map, and the national statistical report of the Turkish Statistical Institute (TSI 2020), respectively.The emission factors (EF) for NO 2 , SO 2 , CO, CO 2 , CH 4 gases, and PM 2.5 , PM 10 particles are taken from literature (Table 3).The EF of the wheat-barley class was calculated by taking the average of the values determined separately for both.Also, the combustion factor (Cf), dry matter ratio (d), and the ratio between the mass of crop residues and the crop yield (s) values were obtained from the literature.These values were found for maize, cotton, and wheat-barley classes.However, when calculating the emission from the wheatbarley class, the values of wheat were used since the amounts determined for the two crops were similar.On the other hand, in the other class, these values could not be determined from past studies, and calculations were made by taking the average of the other three crop's values (Table 4).Finally, since the proportion of burning residues (pb) value in the formula is about 1 (EEA 2016), it is not considered.Greenhouse gases and PMs emissions from CRB based on crop type were determined by calculation.Estimations were made for each province and values were collected to calculate the total emission in the study area.

Satellite based emission estimation
We created NO 2 density maps for June, July, and August when the harvesting is intense for different crop types.Besides, we produced a NO 2 density map for November following the local news on CRB applications in November.We generated monthly average values for NO 2 concentrations averaged over 30-31 days.After analysing minimum and maximum values of NO 2 concentrations for different months and operational NO maps provided by the Copernicus program (URL 13), we used the 0-150 micromole/m 2 range for our maps.

Crop type mapping
We created a crop type map consisting of maize, cotton, wheat-barley, and other classes using the SVM method.Our results showed a total maize area of 2233.42 km 2 , cotton of 10,378.85km 2 , 19,644.47 km 2 (wheat and barley), and 3090.82km 2 other in the study area.
Maize and cotton crops are mostly grown in the Sanliurfa province, and the wheat-barley class is grown in Diyarbakir (Figure 5).
We obtained the highest cross-validation result for the province of Sanliurfa (0.95) and the lowest values for Diyarbakır (0.81) (Table 5).We also conducted independent accuracy assessment using 25 test polygons for each class for the whole region.These polygons  were different from the train and validation dataset used in cross-validation.We obtained high class-wise accuracy for different crops (Table 6).

Burned area mapping
We used Sentinel-2A/B satellite images obtained before and after the harvest from May and November for the classification step.We used dNBR image to find out the spatial extent of burned agricultural parcels.dNBR values are generally suitable for the study area; but, we spotted some dark-coloured fields in the moderate low severity class.After carefully analysing each burnt severity class, we decided that Moderate-high severity (0.44-0.66) and High Severity (0.66-) classes successfully represented the burned areas (Figure 6).During our visual analyses, we used some agricultural parcels with burned smoke, easily seen in the Sentinel-2 images, and also benefitted from FIRMS active fire products.As a result of the study, we found 5063.52 km 2 of Moderate-high severity and 80.53 km 2 of high-severity burned area, and 31% of these were in Diyarbakır province, which was the highest.For Sanliurfa, Mardin, Sirnak, Siirt, Adiyaman, Batman, Gaziantep, and Kilis; the burnt area percentages were 19%, 15%, 10%, 5%, 4%, 3% and 1%, respectively.We found that burned areas spatially well matched with the agriculturally intensive areas in the crop-type map (Figure 6).We found some agricultural parcels for which simultaneous Sentinel-2 images acquired on the day of CRB practice were available.We concluded that high-severity burned areas represent these occasions based on these images.On the other hand, moderate-high severity areas were mostly for those agricultural parcels whose image was taken about a week after being burned.Since the possibility of having a simultaneous image with the CRB date is low, the area of the high severity class was less than the other class (Figure 6).The total emission caused by the maize residue burning was higher than that of other crops.The same result was also obtained by Dennis et al. (2002).According to the results, it was determined that the burned area map had 87% overall accuracy and 0.73 kappa value.We downloaded FIRMS active fire data of the study region for 1 March À 30 November 2019 period and visually compared this map with the determined CRBexposed areas (Figure 7).These two maps, in general, spatially matched and we observed that Sanliurfa, Diyarbakir and Mardin provinces had higher similarities in both maps.NASA's Fire information for resource management systems (FIRMS) consists of satellite fire products from MODIS (Giglio et al. 2016) and Visible Infrared Imaging Radiometer Suite (VIIRS) (Schroeder et al. 2014) active fire detection products.The spatial resolutions of MODIS and VIIRS are 1000 m and 375 m, respectively.

Emission estimation
In the study area, we found that 14,444.307Gg greenhouse gas and 117.809Gg PMs were released due to CRB in 2019.Our results showed that CO 2 and CO values were higher than other gases and PMs (Table 7).For the whole region, PM emissions were primarily due to maize and cotton burning.Most of the CH 4 emissions were caused by cotton burning.In addition, our crop-based total emission analysis illustrated that residue burning of maize caused the highest emission (5863.58Gg), while wheat-barley caused the lowest emission (2129.40Gg).The total emission from the cotton burning was 4437.98 Gg.

Province based emission estimation results
The province-based emissions analysis showed that the highest amount was released in Diyarbakır, whereas the lowest was in Kilis.We examined the emission amounts of all gases and PMs for each provinces based on crop types (Table 8).In Adiyaman, we observed that the highest emission of 263.42 Gg was released from CO2 gas due to cotton residue burning.In contrast, CO 2 emissions of maize, wheat-barley and other crops were 64.81 Gg, 7.14 Gg, and 36.52 Gg, respectively (Table 8).
We determined that the total amount of emission due to CRB was 559.28 Gg in Gaziantep; the highest contribution was for CO 2 (531.25 Gg) primarily caused by cotton burning, similar to the Adiyaman case.In Gaziantep, most of the emissions were caused  by cotton burning, followed by maize and other crops burning.Emission amounts caused by CO and CO 2 were higher compared to other gases and PMs (Table 9).
In Sanliurfa province, the primary cause of emissions was maize burning which released a total of 1863.74Gg, with the highest proportion from CO2 (1804.43Gg) and CO (42.59 Gg) (Table 10).Although the emissions caused by cotton burning were lower than the maize in the Sanliurfa region, total emissions caused by cotton burning were still higher compared to neighbourhood provinces such as Adiyaman, Gaziantep, Mardin, and Siirt.
In Mardin province, the primary source of emissions were maize and cotton residue burnings, respectively; whereas the wheat barley residue burning resulted in less CO 2 and CO gases produced higher emissions than other gases and PMs (Table 11).
When we analyzed the total emissions caused by different crop residue burnings in selected provinces, we concluded that Mardin and Sanliurfa were exposed to higher emissions primarily by maize and cotton residue burnings.Kilis found as the cleanest province in terms of CRB emissions and total emissions of CH 4 , NO 2 , and SO 2 were too low to be ignored.Total emissions caused by maize residue burning were lower for Siirt, Sirnak and Kilis, whereas it reached the maximum for Mardin (2860.46Gg) followed by Sanliurfa (1863.74Gg), and Diyarbakir (800.16Gg) (Table 12).
We found that wheat-barley residue burning resulted in the highest CO 2 emission in Diyarbakir compared to all other provinces.
In general, CO 2 emissions were high for all crops, followed by CO emissions, and in most cases, CH4-related emissions were higher than SO 2 and NO 2 emissions.
We found that maize and cotton are the primary crops causing the highest CO emissions.The NO 2 density was determined to be the lowest value in the wheat-barley and the highest value in the other class.The emission from cotton was very dominant in the distribution of SO 2 .Finally, although the emissions of PMs from cotton, maize, and wheat-barley were similar, the PM 25 density resulting from the burning of the other was lower than PM 10 (Table 12).

Satellite based emission estimation
As mentioned in the previous sections, we calculated the monthly average NO 2 density maps for June, July, August and November using Sentinel-59 data by applying 0-150 micromole/m 2 range (Figure 8).We observed that the high-density areas in the NO 2 maps were consistent with where the CRB was applied (Figure 9).
Finally, we generated a Sentinel-5P NO 2 map for November to check the claims by the local people presented in the media regarding to the CRB cases in November.NO 2 concentration was found to be high in November (Figure 10), however it is not possible to say that it is precisely due to the CRB because of the cooling of the weather in November and the emission increase due to warming (Figure 10).

Discussion
The human impact on climate change has been clearly stated in the IPCC 6th Report (Marques et al. 2022, Nakaegawa andMurazaki 2022).T€ urkiye is located in the Mediterranean Basin, where vulnerability to climate change is very high.Therefore, human-induced fires such as residue burning and wildfires should be closely monitored since the extreme temperatures due to climate change trigger the occurrence and spread of fires.Within the scope of climate change adaptation and environmental policies, reducing emissions caused by greenhouse gases and particulate matter is significantly essential.Therefore, accurate identification of burned agricultural areas and calculation of CRBrelated emission amounts are vital to cope with climate change and correctly determine air quality effects.
Moreover, it is also essential to obtain crop-based emission factors and phenological information to generate national-level reliable CRB-related emission estimations.We used   the EF values published by Santiago-De La Rosa et al. (2018) and McCarty et al. (2012).Still, it is important to determine the region and crop-specific values for T€ urkiye in future studies.We used Sentinel-2 images in this research to find the spatial distribution of agricultural areas, crop types and burned areas.We selected Sentinel-2 due to its 10 m spatial resolution compliance to determine small and fragmented agricultural parcels in our study region.In case of availability of large agricultural parcels, Landsat-8 and Landsat-9 images could also be integrated into implementation.Most of the studies in the literature used MODIS data for CRB determination due to its daily temporal resolution (Chen et al. 2017;McCarty et al. 2017;Lee et al. 2021;Kalluri et al. 2020).Although MODIS data could produce good results for large agricultural parcels, its spatial resolution is not sufficient to obtain high quality details, specifically for small and fragmented agricultural areas.Sentinel-2 satellites, with their high spectral and spatial resolution, provide an essential data set for CRB estimation; however, the number of CRB research with Sentinel-2 is quite limited (Ding et al. 2020).Deshpande et al. (2022) successfully used Sentinel-2 images for Central India to find burned agricultural areas; however, they did not conduct specific crop-based emission calculations.Instead, they estimated emissions assuming that all agricultural burning was due to wheat residue since wheat was one of the major winter crops for their area.Singh et al. (2021) used Sentinel-2 data to estimate CRB for the rice fields in the northern part of India.Our study region is large, 76458,9 km2, and includes different crop types and we estimated emissions due to different crop residue burning.The usage of the GEE platform provided an advantage to this work in terms of data processing time, dynamic access to the massive amount of satellite images, availability of several image processing functions and no requirement to download data.
The NO 2 density product of the Sentinel-5P satellite is generally helpful in finding hot spot areas exposed to high NO 2 releases.Although the satellite provides daily data, the loss of some data sets (CO, CH 4 , and SO 2 ) in the study area is a disadvantage for detecting emissions of all greenhouse gases.

Conclusion
In this research, we determined the spatio-temporal distribution of crop types and burned areas due to CRB using Sentinel-2 images for a large region with diverse, fragmented agricultural parcels.We calculated crop-specific emissions caused by CRB using the EF values in the literature for maize, cotton, wheat and barley.The fact that the EF value, which varies depending on the country, region, and climate, had to be taken from the literature made it difficult to evaluate this study on a national scale; therefore, region-based determination of EF via field measurements should be conducted in future studies to have more precise emission estimations.We also used Sentinel-5P NO 2 density maps, MODIS and VIIRS active fire data to check the spatial matching of our results with these satellitebased products.MODIS and VIIRS spatial resolutions are insufficient to find small agricultural parcels exposed to burning.
Moreover, in some locations, fires did not occur, although there were indications in MODIS and VIIRS active fire data.Knowledge of the crop's phenologic stages and periods and harvest time is also essential to accurately generate crop type maps and identify burned areas.Both high spatial and temporal resolution sensors are needed for CRB detection since fires might occur on small agricultural parcels and start quickly and last a short time.Therefore, integrating multi-sensor data to have high spatio-temporal coverage might be necessary to identify CRB-affected fields.

Figure 1 .
Figure 1.Location and mosaic image of the study area.

Figure 2 .
Figure 2. Some examples of CRB in the study area (Source: URL 4, URL 5).

Figure 3 .
Figure 3.The temporal and spectral change of a field in Diyarbakır where CRB was applied.

Figure 5 .
Figure 5. Crop-type map of the study area.

Figure 6 .
Figure 6.Burned area mapping of the study.

Figure 7 .
Figure 7. FIRMS active fire data of the study region for 1 March-30 November 2019 period (a) Modis based active fire data (b) VIRS based active fire data (URL 14).

Figure 10 .
Figure 10.The N02 concentration of the region in November.

Table 1 .
The main crops grown in the study area.

Table 2 .
Sentinel-2 bands used in various application parts of the study.

Table 3 .
Emission factors of greenhouse gases and PM from literature.

Table 4 .
The values used for the parameters in the emission calculation.De LaRosa et al. (2018).bMcCarty et al. (2012).c Average of other crops.

Table 6 .
Accuracy assessment of crop-type map.

Table 7 .
Crop based emission amount of different gases and PMs in GAP region (Gg).

Table 8 .
Crop based emission amount of different gases and PMs (Gg) in Adiyaman.

Table 9 .
Crop based emission amount of different gases and PMs (Gg) in Gaziantep.

Table 10 .
Crop based emission amount of different gases and PMs (Gg) in S¸anlıurfa.

Table 11 .
Crop based emission amount of different gases and PMs (Gg) in Mardin.

Table 12 .
Crop based total emisions (Gg) for each province.
Figure 8.The NO 2 maps for June, July, and August.