Greenhouse Gas Emission Model for Tidal Flats in the Republic of Korea

: Since coastal wetlands have been severely degraded and polluted by human activities, they have increasingly become a signiﬁcant source of greenhouse gases (GHGs), so understanding the characteristics of their emissions is critical for devising future climate change mitigation strategies. This study modiﬁed a model based on carbon balance to forecast carbon stored and CO 2 , CH 4 emissions in four types of typical tidal ﬂats— Phragmites australis (PA), Spartina alterniﬂora (SA), Suaeda japonica (SJ), and Bare Tidal Flat (BTF) in Korea’s Ganghwa province from 2017 to 2047. The model was built using biomass data from salt plant species collected in different locations. The results indicate that the total annual simulated ﬂow of CH 4 increased over time in all four areas, most notably in SA, while CO 2 remained relatively stable. The mean CO 2 and CH 4 ﬂuxes in the four types of representative tidal ﬂats were in the range of 0.03 to 19.1 mg m − 2 d − 1 and 0.007 to 5.23 mg m − 2 d − 1 , respectively, across all seasons. Besides, the results indicate that the amount of carbon accumulated in the top soil increases linearly over time in nearly all areas studied, ranging from 0.01 to 0.13 (kgC m − 2 yr − 1 ). In general, the study provides a model for Korean tidal ﬂats that incorporates carbon storage and GHG emissions in the intertidal zone in order to develop potential GHG reduction scenarios.


Introduction
A tidal flat acts both as a carbon sink, by storing carbon in sediment and vegetation, and as a carbon source, by releasing major greenhouse gases such as CO 2 and CH 4 into the atmosphere [1]. CO 2 is fixed by plants and autotrophs through photosynthesis to form long-time storage organic compounds in plants and soil, a process known as carbon sequestration [2]. Recent studies have shown that the decomposition of soil organic matter is slow in the coastal wetland system because of anaerobic conditions in saturated sediments and salinity of the soil and water [3][4][5], however, decomposition still takes place through aerobic and anaerobic processes that produce CO 2 and CH 4 [6].
The median organic carbon stock sequestrated in the Korean tidal flats is 65.98 t C ha −1 in vegetated areas, and 53.91 t C ha −1 in barren areas [7], and the total amount of carbon stored in soil in the intertidal flats is approximately 6373.4 Gg [8]. Besides, statistics on GHG emissions from the tidal flats in Korea are also receiving great concern from researchers. According to a study carried out in 2016, the potential tidal flats emission corresponds to around 55.53 × 10 6 t CO 2 e [6]. Chen et al. reported CO 2 and CH 4 emissions at tidal flat types of the Yellow River Delta wetlands with mean fluxes of −20.98 to 68.12 mg m −2 h −1 and −0.12 to 0.44 mg m −2 h −1 [9]. Thus, the investigation of carbon dynamics in the tidal flat area plays an important role in management, however, Korean tidal flats are a dynamic environment in which various processes such as tides, waves and precipitation interact to control the sedimentation in a complicated manner [10]. The large variation in space and time makes it difficult to determine the amount of carbon stored in the sediment and the emission quantities from the tidal flat. Direct measurement of CO 2 and CH 4 emissions using chamber or sediment sampling at defined areas provide essential information on a location basis over a period of time and defined space. However, due to the diversity of variables, such as soil, climate, and vegetation, it is difficult to generalize this information on a larger scale [11]. Therefore, the development of a model based on the interactive reactions of the biochemical processes in the carbon cycle is necessary to calculate the amount of CO 2 and CH 4 released from tidal flats.
Currently, soil carbon dynamics are simulated using a variety of models, including RothC, Century, and ICBM (Introductory Carbon Balance Model) [12][13][14]. In the majority of treatments, these models confidently predicted long-term trends. All of these models have been validated in a variety of climate zones around the world, and the results indicate that the more complex models employ between four (RothC) and six (Century) (more or less) dynamic soil carbon pools, and have sufficient complexity to at least attempt to model short-term dynamics (on the order of weeks and months). The simplest models (such as ICBM) that are capable of simulating long-term dynamics over several decades make use of a single dynamic pool [15]. Admittedly, models for predicting carbon dynamics in tidal flat areas are limited and are attracting considerable attention from scientists.
This study modified the ICBM (Introductory Carbon Balance Model) and applied it to provide simulation results of GHG emitted from Tidal Flat over the period of 2017 and 2047. Since ICBM has been regarded as one of the excellent models to simulate carbon dynamics in soil, it has been extensively applied to estimate the amount of carbon stored in agricultural soils over a long period to establish the efficient and effective soil management plans [16][17][18]. The main objectives of this study were to improve the ICBM, to optimize the model parameters in accordance with a tidal flat area, and to evaluate the applicability of the model to predict the amount of carbon stored in the sediment and GHG emissions (CO 2 and CH 4 ) over 30 years from 2017 to 2047, using data from the vegetated tidal flat and BTF areas in Ganghwa, Korea. The main advantage of developing such a process-based model is to forecast GHG emissions quantitatively, offering suitable management practices.

Model Description
This model named as Greenhouse Emission Model for Blue Carbon (GEM_BC) is based on the carbon mass balance to simulate the carbon dynamics in tidal flat layers. The model assumes that the tidal flat layer is composed of two layers: young layer (YL), defined as the layer where organic carbons are decomposed biologically to generate CO 2 only under aerobic conditions, and old layer (OL), defined as the layer where those emit CO 2 and CH 4 under anaerobic conditions. To develop this model, it is essential to understand the variation of layer environment with time as well as organic carbon flows through the layers. Organic carbons originating from organic debris in the water layer would continuously be deposited on the YL. As already mentioned, CO 2 , one of the major GHG, is likely to be produced from the YL. Some organic carbons in the YL are resistant to biological decompositions that they may undergo humification process to accumulate carbon, leading to a gradual growth of the YL and eventually shifting to the OL. As the thickness of the YL increases, the air penetration would be limited below a certain depth, forming anaerobic conditions and thus turning the lowest part of the YL to the OL where anaerobic biological decomposition would be initiated to release CO 2 and CH 4 ( Figure 1).
The conceptual model is expressed mathematically based on the following assumptions: First-order kinetics is responsible for organic carbon dynamics: That is, the decomposition of organic matter in soil is characterized by pseudo first-order kinetics to represent carbon dynamics [19,20]. External influence coefficient is relatively constant with time, assuming a steady external input [21]. Annual amount of carbon sequestration is accounted for by that in the total biomass of the one-year plants. Lateral C flux: Currently, there are no specific research results on the lateral C flux in the tidal flat areas of Korea. Besides, the west coast of South Korea receives large amounts of sediment from rivers, it is considered a carbon sink [1]. This amount of carbon is mainly deposited in the tidal flat areas [8]. Therefore, the model assumed that the amount of carbon transported by the lateral C flux in the study areas is negligible.
YL and eventually shifting to the OL. As the thickness of the YL increases, the air penetration would be limited below a certain depth, forming anaerobic conditions and thus turning the lowest part of the YL to the OL where anaerobic biological decomposition would be initiated to release CO2 and CH4 (Figure 1). The conceptual model is expressed mathematically based on the following assumptions:  First-order kinetics is responsible for organic carbon dynamics: That is, the decomposition of organic matter in soil is characterized by pseudo first-order kinetics to represent carbon dynamics [19,20].  External influence coefficient is relatively constant with time, assuming a steady external input [21].  Annual amount of carbon sequestration is accounted for by that in the total biomass of the one-year plants. . This amount of carbon is mainly deposited in the tidal flat areas [8]. Therefore, the model assumed that the amount of carbon transported by the lateral C flux in the study areas is negligible.

Kinetic Equations
Kinetic equations are denoted for organic carbon dynamics in YL and OL as follows: where C represents the areal amount of organic carbon (kgC m −2 ), i denotes the input rate of organic carbon to the YL (kgC m −2 yr −1 ), k21 denotes the decomposition rate constant in YL (yr −1 ), is the external influence coefficient, h is humification coefficient, k22 is decomposition rate constant in OL (yr −1 ). Subscripts 1 and 2 represent YL and OL, respectively. Integration of Equations (1) and (2) allow the following expression for organic carbon content with time in YL and OL: . Where, C 1ss = total organic carbon (kgC m −2 ) of YL in a steady state, C 2ss = total amount of organic carbon in OL in the steady-state (kgC m −2 ).

Kinetic Equations
Kinetic equations are denoted for organic carbon dynamics in YL and OL as follows: where C represents the areal amount of organic carbon (kgC m −2 ), i denotes the input rate of organic carbon to the YL (kgC m −2 yr −1 ), k 21 denotes the decomposition rate constant in YL (yr −1 ), is the external influence coefficient, h is humification coefficient, k 22 is decomposition rate constant in OL (yr −1 ). Subscripts 1 and 2 represent YL and OL, respectively. Integration of Equations (1) and (2) allow the following expression for organic carbon content with time in YL and OL: where C 10 and C 20 are initial values for C 1 and C 2 , respectively. In YL, soil carbon decomposes aerobically through microbial respiration, releasing CO 2 into the atmosphere. The humification coefficient (h) controls the fraction of YL that enters OL and (1−h) then represents the fraction of the outflow from YL that directly becomes CO 2 -C ( Figure 1). The CO 2 emissions from YL is calculated using Equation (5): where Q YLCO 2 is CO 2 emission from YL. In OL, organic carbons are anaerobically decomposed to release CO 2 and CH 4 . The study assumed a ratio of carbon sequestration (C s ) to carbon in CH 4 (CH 4 _C) for PA, SA, SJ, and BTF areas in OL (see Table 1), based on the correlation between annual C s and CH 4 emissions investigated in tidal flat areas [7,[22][23][24]. The CO 2 and CH 4 emissions from OL are calculated using Equations (6) and (7): where Q OLCO 2 is CO 2 emission from OL. Q OLCH 4 is CH 4 emission from OL. 1/a is the ratio of CH 4 _C: C s from OL.

Determination of Model Parameters
Five parameters ( , k 21 , k 22 , i, and h) are needed to calculate the kinetic equation described in Section 2.1.1.
, denoting the eternal influences on soil organic matter decompositions rates of the tidal flat, depends on three main factors: (i) soil temperature factor ( t ), (ii) soil moisture factor ( m ), and (iii) soil salinity factor ( s ) of the tidal flat. The daily values for were calculated using the following equation: = t · m · s . In this equation, t and m were calculated from daily meteorological data and soil data following the approach in Andrén et al. [21] and Fortin et al. [25]. The s factor depends on the osmotic potential at the actual soil water content and electrical conductivity. The exponential relationship between microbial activity and salinity is determined by the formula given in Yuan et al. [26].
Where EC 1:5 is the electrical conductivity (dS.m −1 ) of 1:5 soil/water extract. The EC 1:5 was converted from EC e using the equation proposed by Shaw et al. [27].
This study assumed the value of k 21 and k 22 as 0.8 yr −1 and 0.006 yr −1 following Andrén et al. [16].
The fourth parameter i in vegetated area is calculated following the formula , where, M i is the biomass of vegetation (kg m −2 ), CF i is carbon conversion factor for grasses, and A t is the area occupied by vegetation community (m −2 ). In BTF area, the carbon source is mainly provided from the deposition of terrestrial sediments and suspended particles from seawater [8]. Thus, the annual carbon input was determined based on the mean C sequestration rate per year (kgC m −2 yr −1 ).
Lastly, h, denoting the proportion of young soil carbon that becomes sediment carbon (humus), is estimated following the formula: h = (h s · i s + h r · i r )/i [16]. Where h s and h r are humification coefficient of stem and roots, while i s and i r are carbon amount of stem and roots, respectively. In BTF area, the carbon can be considered as sewage sludge, so humification coefficient was set 0.41 [28].

Model Validation
The model's simulation of topsoil carbon mass, CO 2 and CH 4 emissions were validated against observed data available from a study conducted by the Korea Marine Environment Management Corporation as part of the "Development of Blue Carbon Information System and its Assessment for management" project from Ganghwa-do (37.7 • N, 126.5 • E) [29], which encompasses the southern portion of the tidal flat. The study area is sandwiched between the Sukmo and Yeomha channels to the west and east, respectively, and thus receives a significant amount of carbon sediment annually. A layer of mud and salt plants covers most of its surface area. Numerous salt plants are found in the area, with three dominant species: Phragmites australis (PA), Spartina alterniflora (SA), and Suaeda japonica (SJ). The sampling areas, Yeochari, Dongmakri, and Donggeomdo which are homes to PA, SA, and SJ species, respectively, are located south of Ganghwado along the coastline. At Dongmakri, bare tidal flat soil carbon data were also collected ( Figure 2). roots, respectively. In BTF area, the carbon can be considered as sewage sludge, so humification coefficient was set 0.41 [28].

Model Validation
The model's simulation of topsoil carbon mass, CO2 and CH4 emissions were validated against observed data available from a study conducted by the Korea Marine Environment Management Corporation as part of the "Development of Blue Carbon Information System and its Assessment for management" project from Ganghwa-do (37.7° N, 126.5° E) [29], which encompasses the southern portion of the tidal flat. The study area is sandwiched between the Sukmo and Yeomha channels to the west and east, respectively, and thus receives a significant amount of carbon sediment annually. A layer of mud and salt plants covers most of its surface area. Numerous salt plants are found in the area, with three dominant species: Phragmites australis (PA), Spartina alterniflora (SA), and Suaeda japonica (SJ). The sampling areas, Yeochari, Dongmakri, and Donggeomdo which are homes to PA, SA, and SJ species, respectively, are located south of Ganghwado along the coastline. At Dongmakri, bare tidal flat soil carbon data were also collected ( Figure 2).

Experimental Data
During the growing season (May to November) from 2017 to 2019, this study collected PA, SA, and SJ biomass data. The above-ground biomass estimation method makes use of a quadrat to represent sampling areas within plots following Howard et al. [1]. In this experiment, quadrats measured 50 cm × 50 cm in size. Stems were collected in each quadrat to determine the wet and dry weights. In the laboratory, wet weight stem biomass was oven-dried for approximately 72 h at 60 degrees Celsius to obtain a constant dry weight stem biomass.

Experimental Data
During the growing season (May to November) from 2017 to 2019, this study collected PA, SA, and SJ biomass data. The above-ground biomass estimation method makes use of a quadrat to represent sampling areas within plots following Howard et al. [1]. In this experiment, quadrats measured 50 cm × 50 cm in size. Stems were collected in each quadrat to determine the wet and dry weights. In the laboratory, wet weight stem biomass was oven-dried for approximately 72 h at 60 degrees Celsius to obtain a constant dry weight stem biomass.
This study used direct sampling to determine the biomass of the roots (below-ground biomass) and the samples were taken by extracting a core as described in Howard et al. [1]. Further, following methodology in Saunders et al. (2006), 10 cm diameter cores for sampling 2.5 cm wide segments cut at various depths of 2.5, 5, 7.5, 12.5, 15, 22.5, and 25 cm, and another sampling 5 cm wide segments cut at depths of 35, 45, 55, and 65 cm were used [30]. The segments were then visually separated into root, rhizome, and dead litter components after washing them over a 1 mm screen. The living and dead things were distinguished based on color [30]. Following that, living root and rhizome materials were oven-dried to a constant dry weight at 60 degrees Celsius. The below-ground biomass was then determined in the manner described by Howard et al.
[1]. Table 1 shows the data input for the model. Total biomass was calculated in vegetated areas using stem, roots, and rhizodeposition (the amount of biomass allocated to the soil during the growing season as root change and exudates). This study used Pausch et al. estimate of 68.75% root turnover for rhizodeposition, which is calculated by multiplying the estimated root biomass C by 0.6875 and adding that as C input from root turnover [31]. The halophytes that inhabit in Ganghwa are highly dependent on the weather conditions during the winter months. Despite minor differences, most juveniles emerge in early spring and grow until autumn. In the winter, they turn brown and cease to grow. Therefore, this study assumed that the entire C amount in biomass would be converted into annual carbon input to soil. SJ and PA are two common native plants in the Ganghwa area. Meanwhile, SA was first discovered near Dongmak Beach in 2012 [32], however, it has expanded, and the density of individuals is high (334 no. stem/m 2 ), so the amount of carbon stored in this plant is greater than that of the other two species-SJ and PA. Optimization of the parameter h was done in accordance with the relative contribution of root versus shoot material to the formation of soil organic matter. The relative contribution of root versus stem material to soil was estimated using Kätterer et al. results [28]. According to the findings of Kätterer et al., root-derived material humifies 2.3 times faster than above-ground plant material [29]. The estimation factors for h s and h r are 0.15 and 0.35, respectively.

Input Data and Parameter Values in Ganghwa-do, Korea
Among the most important factors influencing decomposition rates in soil are soil temperature and water content. Also, various studies show that increased salinity slows microbial degradation [26,33,34]. After having gathered the raw input data (see Table 1), this study analyzed the raw data to obtain the external factors parameter . For example, to calculate the effect of m , evaporation (ET o ) was calculated using equation 6 of Allen et al. and has an annual mean value of 2380 (mm) [35]; Pedotransfer functions were used to estimate the soil's field capacity and wilting point [36], which were determined as 35% and 22%, respectively. Furthermore, the EC e data was gathered from the literature to determine the impact of soil salinity on the rate of decomposition [37]. Due to the study areas proximity and similar climatic conditions, the difference in mean value of is very close, and almost equal. In general, the value of represents the peak during the summer months, which is attributable to soil temperature and soil water content, and is usually low during the cold and dry winter months (see Figure A1).  Figure 3 illustrates the model's prediction results for the carbon content of the topsoil (0-25 cm depth) corresponding to halophyte types and BTF in the Ganghwa area from 2017 to 2047. As shown, the amount of carbon accumulated in the soil increases linearly over time in almost all of the areas studied. However, PA and SA predicted a significant increase in the amount of carbon accumulated in the soil. Specifically, the amount of carbon in the topsoil increased by an annual average of 0.1, 0.13, 0.01, and 0.02 (kgC m −2 yr −1 ), in the PA, SA, SJ, and BTF study areas, respectively. This depends largely on the amount of carbon that enters the ecosystem each year. According to the Blue Carbon Manual published by the IUCN (International Union for Conservation of Nature), the composition and density of vegetation in coastal wetlands have an impact on the total carbon stock [39]. Studies regarding the vegetated tidal flat area indicate that high sedimentary organic carbon content in salt marshes is primarily due to increased carbon input via plant litter and roots, as well as their location in saturated and, presumably, more anaerobic environments as opposed to other types of marshes [40]. For the BTF region, the increase in carbon in topsoil sediments was primarily due to organic carbon deposition from marine and terrestrial sources. The Ganghwa tidal flat, which is influenced by terrestrial inputs such as freshwater discharge, municipal sewage, and discharge, contributes to the burial of organic carbon in sediments [22]. All the relevant data related to Figure 3 is provided in Appendix A, Table A1.

Simulation of Topsoil Carbon Mass
For the purpose of validating the simulation results of the model, Figure 4 shows a comparison of the simulated and measured carbon stored in sediments at a depth of 25 cm in Ganghwa from 2017 to 2020, which was conducted at four different locations throughout the year. As illustrated in Figure 4, the annual carbon stored in the simulated sediment is generally consistent with the data collected from field sites at various points in time. Also, at some points the simulated magnitude is slightly different from the observed data, with observed carbon in sediments being higher or lower. This could be because tidal energy in a variety of topographical and morphological settings could be a critical factor in increasing not only the supply of waterborne sediment particles, but also organic carbon burial in intertidal sediments. Lee et al. demonstrated that the amount of sedimentary carbon varies over time and space due to a variety of factors such as wave, tidal, and erosion and coastal environments [22]. Thus, the unique parameter affecting the carbon mass in soil prediction results is i. Besides that, ϵ is a parameter that has a direct effect on the decomposition process. It is worth noting that this study did not account for carbon deposits from terrestrial inputs caused by riparian discharge into salt marshes, which could have a significant impact. Additionally, because of the scarcity of data, the model's test results do not account for wave, tidal, or erosion effects, all of which contribute to annual carbon uncertainty.   Thus, the unique parameter affecting the carbon mass in soil prediction results is i. Besides that, ϵ is a parameter that has a direct effect on the decomposition process. It is worth noting that this study did not account for carbon deposits from terrestrial inputs caused by riparian discharge into salt marshes, which could have a significant impact. Additionally, because of the scarcity of data, the model's test results do not account for wave, tidal, or erosion effects, all of which contribute to annual carbon uncertainty.  Thus, the unique parameter affecting the carbon mass in soil prediction results is i. Besides that, is a parameter that has a direct effect on the decomposition process. It is worth noting that this study did not account for carbon deposits from terrestrial inputs caused by riparian discharge into salt marshes, which could have a significant impact. Additionally, because of the scarcity of data, the model's test results do not account for wave, tidal, or erosion effects, all of which contribute to annual carbon uncertainty.

Simulation of GHG Emissions
The simulation results for CO 2 and CH 4 emissions in all four test sites from 2017 to 2047 were as shown in Figure 5. Litter decomposition is the primary source of CO 2 , which is influenced by temperature and humidity. The results demonstrate that the annual carbon input and the external coefficient have a strong influence on the simulation results of CO 2 emitted. The comparison of monthly simulated and observed soil CO 2 emissions for the years 2019 at all four sites were as shown in Figure 6. Both prediction and measurement of CO 2 emissions in the fields followed a similar pattern, with the highest levels of emissions being recorded during the May to August growing season. This could be because during the summer temperatures are high, resulting in increased respiration rates for roots and soil microorganisms. Although the data are only for seasonally specific months, the CO 2 emissions predictions are comparable to the measured. However, the observed and simulated daily fluxes differ in terms of peak timing and magnitude. The model's CO 2 emission peak occurs in May for all regions as follows; 7.61 mg m −2 day −1 for PA, 7.92 mg m −2 day −1 for SA, 6.03 mg m −2 day −1 for SJ, and 0.83 mg m −2 day −1 for BTF. This is because the applicable sites are in close proximity to one another, implying that the natural conditions are similar. Meanwhile, the areas PA and BTF have emission peaks of 12 mgCO 2 -C m −2 day −1 , and 11 mgCO 2 -C m −2 day −1 in August, while the SA and BTF areas have emission peaks of 11 mgCO 2 -C m −2 day −1 in May.

Conclusions
This study developed the GHG Emissions Model for tidal flat areas to make estimates of CO2 and CH4 emissions based on the carbon balance of the carbon sink. The model of Tidal Flat GHG emissions has been verified using annual carbon accumulation in sediment from 2017 to 2020 and daily CO2 and CH4 emissions in 2019 from four different locations in Ganghwa, Korea. The results indicate that, in comparison to 2017, CH4 emissions in the PA and SA areas have increased over time due to increased carbon accumulation in the sediment. In comparison, the results of the simulations in the SJ and BTF regions remained relatively stable over time. CO2 emissions from sites also remained relatively stable over time. The study demonstrated that the model is capable of capturing significant trends in soil CO2 and CH4 emissions under natural conditions in Ganghwa, Korea. According to the results of the simulated and observed comparisons, CO2 and CH4 emissions occur primarily during the growing season, while are near zero in most areas CH 4 is produced during the degradation of organic matter in anaerobic (OL) conditions. The major factors affecting the amount of CH 4 released are dissolved organic carbon (DOC) and soil moisture. Carbon input and methanogenic activity are correlated positively, i.e., increase in carbon input from biomass and exudates results to increase in methane production [41]. Thus, the more DOC sediment present, the more CH 4 is released. According to Figure 5, the amount of CH 4 released in the PA and SA areas increases over time, owing to the increase in carbon accumulated in the sediment. In comparison, CH 4 emissions decreased in the BTF and SJ regions, owing to the region's low annual carbon input, which is insufficient to compensate for the amount of sediment degraded by microbial decomposition activities. Figure 6 illustrates the comparison of simulated and observed values in various vegetative areas in 2019. Although the DOC has a significant positive correlation with CH 4 emissions, Kim et al. assert that it has no direct effect on CH 4 emissions [41]. The dependence on DOC and soil moisture results in a proportional increase in bacteria that reduce sulfate, thereby affecting CH 4 emissions [41]. Sulfate-reducing bacteria competitively inhibits methanogenic activities which are the primary source of CH 4 [42]. The composition and quality of carbon sources have a significant effect on the structure of the methanogen-generating community [43], and thus the amount of CH 4 released varies based on salt plant regions. Recent studies report that the potential for methane production is greater in SA regions than in PA, SJ, and BTF regions, owing to the influence of microbiological content, water table, and depth [23,40,41]. In comparison to the observed values, the simulated CH 4 concentration in winter months was highly consistent, whereas the simulated results for the growing season in the PA, SJ, and BTF areas were generally higher.
In general, the simulation results of the model show a significant increase in CH 4 emissions in the future, while CO 2 emissions will remain relatively stable. The main reason is due to increased sequestration of carbon in the sediments. The difference between simulation and observed results is due to cross factors that contribute to the uncertainty of the model such as tides, changes in salinity, content of microbial enzymes, etc. as well as factors affecting the change in carbon content in sediments.

Conclusions
This study developed the GHG Emissions Model for tidal flat areas to make estimates of CO 2 and CH 4 emissions based on the carbon balance of the carbon sink. The model of Tidal Flat GHG emissions has been verified using annual carbon accumulation in sediment from 2017 to 2020 and daily CO 2 and CH 4 emissions in 2019 from four different locations in Ganghwa, Korea. The results indicate that, in comparison to 2017, CH 4 emissions in the PA and SA areas have increased over time due to increased carbon accumulation in the sediment. In comparison, the results of the simulations in the SJ and BTF regions remained relatively stable over time. CO 2 emissions from sites also remained relatively stable over time. The study demonstrated that the model is capable of capturing significant trends in soil CO 2 and CH 4 emissions under natural conditions in Ganghwa, Korea. According to the results of the simulated and observed comparisons, CO 2 and CH 4 emissions occur primarily during the growing season, while are near zero in most areas during the winter. Following an examination of the relationship between carbon sequestration and GHG emissions, the tidal flat ecosystem is critical for carbon sequestration.
The simulation results were generally consistent with the corresponding experimental data under the model application conditions. However, when peak time and magnitude differences between observed and daily flux simulations were compared, the modeling simulation did not completely match the measured values. This is because the model's experimental results did not account for the amount of carbon deposited by terrestrial currents that accumulate in salt marshes areas, the effects of tides, erosion on soil carbon content, and salinity variation (limitations of study). This causes uncertainty in the annual carbon input and also affect the model's simulation results. Therefore, a detailed study of lateral C flow and carbon deposition in the experimental areas is required to determine the amount of carbon added or eroded in the topsoil layer, from which the model's simulations can be run more precisely. To obtain more accurate predictions of CO 2 and CH 4 emissions, the model must be refined further by incorporating wave energy, water tide, and the correlation between gas volume emission and the enzyme content of the microorganisms involved. Also, salinity is a dynamic property of the soil that changes in response to environmental factors, however, this study assumed it as a constant. Due to the complexity of the processes involved in the development of salinity and the scarcity of data, it is impossible to calculate the uncertainty associated with electrical conductivity (EC) measurements.
In conclusion, the model forecasted carbon dynamics until 2047 for Korea's tidal flats, including carbon storage and GHG emissions in the intertidal zone, for future management scenarios. Further research is recommended to consider additional methods (such as DMI, GLUE etc.) for reducing the model's uncertainty, and it can also be extended by applying it to other areas of the Korean coastal ecosystem (salt marshes and BTF). To continue developing the model, it could be extended in the future to fully establish the processes involved in the carbon cycle in the blue carbon field. Data Availability Statement: All the raw data used in this study was obtained from Korea Marine Environment Management Corporation (https://www.koem.or.kr/site/eng/main.do, accessed on 6 Janurary 2021). All the remaining data is provided in the tables and figures of the main article and appendix. Further inquiries can be directed to the corresponding author.

Acknowledgments:
We would like to express our gratitude to Korea Marine Environment Management Corporation for providing the data used to validate the model as part of the Development of Blue Carbon Information System and its Assessment for Management program. The authors would also like to thank the reviewers and the journal's editor for their valuable comments, which have greatly improved the study.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. Table A1. Model projections (2017-2047) of soil carbon mass (kgC m −2 yr −1 ) at 25 cm depth at the four experimental sites in Ganghwa province of Korea. Table A2. Simulation of CO 2 emissions (kg ha −1 yr −1 ) over years at the four experimental sites in Ganghwa province of Korea. Table A3. Simulation of CH 4 emissions (kg ha −1 yr −1 ) over years at the four experimental sites in Ganghwa province of Korea. Figure A1. Mean within 4 years (2017-2020) dynamics of in Ganghwa, Korea.