Global analysis of carbon disulfide (CS2) using the 3-D chemistry transport model STOCHEM

Carbon disulfide (CS2), a precursor to the long-lived carbonyl sulphide (OCS) is one of the main contributors to the atmospheric sulfate layer. The annual fluxes from its sources and sinks are investigated using a 3-D chemistry transport model, STOCHEM-CRI. In terms of the flux analysis, the oxidation of CS2 by OH is found to be the main removal process (76–88% of the total loss) and the dry deposition loss contributes 11–24% to the total loss of CS2. The global burden of CS2 was calculated, varying between 6.1 to 19.2 Tg and the lifetime of CS2 was determined to be within the range of 2.8–3.4 days. The global distribution of CS2 found the Northern Hemisphere (NH) continental landmasses to be the areas of concentration maxima with peak concentrations reaching up to 20 ppt during June-July-August (J-J-A) season and 40 ppt during December-January-February (D-J-F) season in anthropogenic source regions. Oceanic regions returned low CS2 levels of less than 2 ppt. The vertical profile of CS2 shows higher levels up to 3 ppt at 30°N–45°N during J-J-A and up to 4 ppt at 30°N–55°N during D-J-F. The oxidation of CS2 by OH can produce a substantial amount (0.58 Tg/yr) of atmospheric OCS and the annual average surface distribution of this flux shows up to 12 Tg/yr OCS formed in the regions with highest anthropogenic pollution (e.g., South east Asia). In general, the model-measurement comparison reveals an underprediction of model CS2 compared with measured CS2 for most of the regions. It is likely that the emissions of CS2 are being underestimated and there are likely much larger emission sources of atmospheric CS2 than previously suggested.


Introduction
Carbon disulfide (CS 2 ) exerts a significant influence on the global atmospheric sulfur budget [1,2]. The principal source of CS 2 is supposed to be oceanic [3][4][5][6][7][8], but spatial analysis of CS 2 indicates large concentrations over continental masses, suggesting the main sources of CS 2 are of anthropogenic origin [9][10][11]. Blake et al. [11] reported four major industrial sources of CS 2 as black carbon production, rayon manufacture, CS 2 chemical production and use, and sulfur recovery. Ren [12] found that carbonyl sulphide (OCS) converted to CS 2 in the roots and shoots of barley and chickpeas, suggesting it may be a precursor to CS 2 in vegetation and soil. Soil has no consistent flux direction for CS 2 [13]. Soil and wetland fluxes of CS 2 are sparse [10] because of the large variety of soil types (e.g., oxic soils, anoxic soils) and wetland types (e.g., freshwater mangrove swamp, freshwater grassy marsh, and paddy fields) which make them extremely hard to estimate.
The main removal process of CS 2 is the photochemical oxidation with OH radicals producing both sulfur dioxide (SO 2 ) and OCS as the major products [14,15].
The oxidation of CS 2 and the further oxidation of OCS form sulphate aerosols which can influence the radiative properties of the Earth's atmosphere, potentially contributing to climate change and the stratospheric ozone concentration [9,[16][17][18][19][20]. CS 2 can also undergo removal by O( 3 P) radicals although this is negligible due to low O( 3 P) radical abundances in the troposphere [9]. The photolysis of CS 2 is possible only in the upper parts of the troposphere by shortwave radiation, but very small at altitudes below 6 km [9].
The dry deposition through vegetation can act as a significant loss process of CS 2 [21], but this is not consistent throughout the studies. Taylor et al. [22] reported CS 2 deposition velocities in the range of 0-1.5 mm/s. However, Xu et al. [21] reported CS 2 deposition velocity of 5.6 mm s −1 with large associated uncertainties at a spruce forest in the Solling Mountains, Germany. Wet deposition loss is considered to be minimal due to CS 2 being relatively insoluble [23].
The climatic significance of the atmospheric sulfur cycle is difficult to assess without accurate estimation of the global atmospheric reduced sulfur compounds (e.g., OCS, DMS, and CS 2 ). The budget for DMS is relatively secure and balanced [10,24,25]; the budget of OCS was previously considered balanced [10,26] but recent studies have shown large uncertainties in its global budget with large missing sources [27][28][29][30][31][32]. There are very few modeling studies to estimate the global budget and surface distribution of tropospheric CS 2 [3,9,[33][34][35], many uncertainties exist in the assessments of its global budget [10]. In this study, we employ STOCHEM-CRI, a 3-D global chemistry transport model, to evaluate the global budget and distribution of CS 2 in the troposphere. We show a comparison between model CS 2 concentrations and a judicious selection of observed CS 2 concentrations from individual field measurements and some flight data set.

Methodology
The global 3-D chemistry transport model, STOCHEM used in this study is from the UK Meteorological Office. This model has been extensively evaluated and assessed in a number of model inter-comparison exercises [36,37], which have shown it to be ideally suited to modeling atmospheric chemistry at a level comparable with Eulerian models. The model uses a Lagrangian approach, where the troposphere is divided into 50,000 air parcels which are advected every three hours by winds from the UK Meteorological Office Hadley Centre general circulation model (GCM). The Lagrangian cells are based on a grid of 1.25° longitude, 0.8333° latitude and 12 unevenly spaced (with respect to altitude) vertical levels with an upper boundary of 100 hPa [38,39]. Each air parcel contains the complete 231 trace gases in the code, which undergo chemical processes. The chemical processes that occur within the parcel, together with emission, dry and wet deposition, mixing and removal processes, are generally uncoupled from advection [39]. During the advection process, the Lagrangian cells are considered to be isolated parcels of air and a treatment of mixing is achieved by imposing a fixed Eulerian grid over the model domain [40]. At each timestep, trace species in Lagrangian air parcels within a Eulerian grid cell are relaxed to the average of the trace species from that group of parcels. Thus, the Eulerian mixing processes are artificially imposed on the Lagrangian parcels to simulate the effects of diffusion in the atmosphere. Turbulent mixing in the boundary layer is achieved by randomly re-assigning the vertical co-ordinates of air parcels over the depth of the planetary boundary layer height. Small-scale convective processes are treated by randomly mixing a fraction of the air parcels between the surface and the cloud top, depending on the convective cloud top, convective cloud cover, and convective precipitation rate. Boundary layer and convection parameterisations were fine-tuned by using 222 Rn observations [27]. Other meteorological variables, such as cloud distributions, are kept fixed over the 3-hour chemistry model time-step [41]. More detailed description of the meteorological parameterisations in STOCHEM can be found in Collins et al. [38] with updates described by Derwent et al. [42].
The chemical mechanism used in STOCHEM, is the common representative intermediates mechanism version 2 and reduction 5 (CRI v2-R5). The CRI v2-R5 was developed on a compound-by-compound basis using 5-day box model simulations with the performance of the chemistry for each compound being compared and optimised with the master chemical mechanism (MCM v3.1) with ozone production being the primary criterion. The detail of the CRI v2-R5 mechanism is given by Jenkin et al. [43], Watson et al. [44], and Utembe et al. [45] with updates highlighted in Utembe et al. [39] and Utembe et al. [46]. In the study, the mechanism consists of 231 chemical species competing in 630 reactions and gives excellent agreement with the MCM v3.1 over a full range of NO x levels [43,44].
Emissions are treated as an additional term to the source fluxes of each species during each integration time step, rather than a step change in species concentration [38,47]. Emission data for oceans and soil is mapped onto a monthly, 5° longitude  5° latitude resolution, two-dimensional source map [48]. The anthropogenic surface distribution of CS 2 has been used as like SO 2 surface distribution which were based on EDGAR v3.2 data for 1995 [49]. Total emission values for CO, NO x , and NMVOCs are obtained from the Precursor of Ozone and their Effects in the Troposphere (POET) inventory [50] for the year 1998 with added emissions of CS 2 from anthropogenic (0.34 Tg/yr), soil + volcano (0.14 Tg/yr) and oceans (0.18 Tg/yr) taken from Watts [10]. Dry deposition is one of the loss processes for CS 2 which is accounted in the model via gravitational movement process. Deposition velocity depends on the location of the air parcel and in the STOCHEM model, there is a distinction made between air parcels over the sea and the land. The dry deposition velocities over land and sea used in the model as 1.5 mm/s [22] and 0.1 mm/s (calculated from the gas transfer velocity based on global average wind speed of ~8 m/s), respectively. The uptake of CS 2 by vegetation with the deposition velocity of 5.4 mm/s [21] was used as dry deposition over land with unchanging deposition over sea (0.1 mm/s) in the model to investigate how the increment of dry deposition over land can affect the global budget of tropospheric CS 2 . The photolysis rate of CS 2 in STOCHEM is calculated explicitly for each air parcel at a time resolution of one hour. The photolysis rate of CS 2 is calculated in STOCHEM using the following integral: represents the photolysis rate constant of CS 2 ; represents the spherically integrated actinic flux at a given wavelength; represents the cross section of CS 2 , at a given wavelength; represents the quantum yield for dissociation of CS 2 at a given wavelength.
The cross-section and quantum yields data were taken from Atkinson et al. [51]. More details about the photolysis reactions implemented in STOCHEM can be found in Khan et al. [52].
The concentrations produced from the model simulation is mapped onto a Eulerian grid resolution 5°  5° with 9 vertical levels. Summing the 50,000 air parcels produces a global burden for CS 2 , which is broken down into the respective source and sink fluxes. The flux outputs are calculated within each grid square by dividing the averaged emissions per air parcel (molecules s −1 ) by its volume, which gives volume-averaged fluxes with units of molecules cm −3 s −1 . Five simulations, STOCHEM-base case detailed in Utembe et al. [39], STOCHEM-DEPO with increasing deposition loss of CS 2 over land, STOCHEM-PHOT with including photolysis loss of CS 2 , STOCHEM-OCSL with increasing global oceanic and soil emissions of CS 2 and STOCHEM-ANTH with increasing global anthropogenic emissions of CS 2 (Table 1) were conducted with meteorology from 1998 for a period of 24 months with the first 12 allowing the model to spin up. Analysis is performed on the subsequent 12 months of data.
The amount of measurement data from both surface stations and flight campaigns is very limited and thus a judicious selection of data is used to evaluate the model simulations. The simulated data for model-measurement comparison was extracted from the closest grid square to the region in question. Some experimental data was taken for a large area and thus encompasses many grid boxes and in such circumstances, data were extracted from multiple grid squares and the average was taken. The modeled data are monthly averages for a given region adjusted to represent the months studied.  Table 2 details the production and loss processes of CS 2 . In this study, only the direct emissions of CS 2 from anthropogenic sources, oceans, soils had been integrated into the model. Three sinks of CS 2 were identified in the STOCHEM-base case: the reaction with OH, the reaction with O( 3 P) and the dry deposition. Oxidation by OH radicals is found to be the main sink of CS 2 (0.58 Tg/yr, 88% of the total sinks), which is concordant with the data estimates (0.57 Tg/yr) from Chin and Davis [9] and Watts [10] who estimated that oxidation by OH accounts for 56-100% of CS 2 removal. The simulated flux of the CS 2 conversion to OCS contributes 43% of the total production flux of OCS (1.31 Tg/yr) estimated by Watts [10]. The concentrations of O( 3 P) are very low across the troposphere and hence the contribution of the loss process, CS 2 + O( 3 P) is found to be very small (<1%) in the study. The dry deposition is another sink of CS 2 (0.07 Tg/yr, 11% of the total loss) which has been found to be significant (0.15 Tg/yr, 23% to the total loss) after a 5-fold increment of deposition velocity over land in the simulation, STOCHEM-DEPO. The emission fluxes of CS 2 from the previous four modeling studies [3,[33][34][35] with a range of global emissions (0.4-2.4 Tg/yr) of CS 2 used in their models (Table 2) gave a range of global burdens (7.8-80 Gg) and lifetimes (4-13 days). However, in our model study, the global burden of CS 2 is found to be small compared with the studies of Pham et al. [33], Weisenstein et al. [34] and Khalil and Rasmussen [3], the discrepancy in emission data is the likely reason for the disagreement between the results. The simulated lifetime of CS 2 for the STOCHEM-base is found to be 3.4 days, which is at the lower end of the range of literature estimates shown in Table 2. The photolysis of CS 2 has been ignored in all previous modeling studies [35], but we included the photolysis of CS 2 in STOCHEM-PHOT and found this loss process very negligible (0.01 Gg/yr) compared with other loss processes. The inclusion of emission classes of Khalil and Rasmussen [3] in STOCHEM-OCSL (increased global soil and oceanic emissions in the model) and Lee and Brimblecombe [54] in STOCHEM-ANTH (increased global anthropogenic emissions in the model) increases the global burden of CS 2 to 16.2 Gg and 19.2 Gg, respectively which are 212 and 271% increments relative to the STOCHEM-PHOT. Thus, when emissions are the only atmospheric source of a compound, a range of emission classes makes it difficult to ensure an accurate and valid output from a modeling study. The oxidation of CS 2 by OH in STOCHEM-OCSL and STOCHEM-ANTH can provide a substantial amount of atmospheric OCS, 1.6 Tg/yr and 1.9 Tg/yr, respectively. n/a n/a n/a n/a Notes: *n/a represents no available data.

Global surface distribution of CS 2
The surface distributions of CS 2 in STOCHEM-base ( Figure 1) show a difference between the continental regions in the northern hemisphere (NH) and southern hemisphere (SH) with increased CS 2 concentrations of up to 20 pptv during June-July-August (J-J-A) and up to 40 pptv during December-January-February (D-J-F) over anthropogenic emission regions (e.g., North America, central Europe and South East Asia) in the NH. The surface distribution pattern of CS 2 in our model is comparable with the studies of Pham et al. [33] and Kjellström et al. [35] who also found highest concentrations of CS 2 over areas of industrial activity. However, the magnitude of surface CS 2 levels in our study is found to be more variable than their studies which are most likely due to the different emission classes used in the models as well as because of the different resolutions of the models. For example, the emissions total used in STOCHEM (this study) is slightly higher than that used in IMAGES by Pham et al. [33], but 2-fold lower than that used in ECHAM3 by Kjellström et al. [35] (see Table 2), but the vertical resolution of STOCHEM (9 vertical levels) is coarser than that of IMAGES (25 vertical levels) and ECHAM3 (19 vertical levels). The surface layer is taken between the surface and 1 km in STOCHEM compared with the models IMAGES and ECHAM3 which were taken from 30-35 m and therefore much greater dilution by mixing and removal that can take place in STOCHEM. Thus, the decreased vertical resolution of STOCHEM may give rise to the lower mixing ratios of CS 2 from areas of increased anthropogenic emissions.
The east coast of the US and South East Asia reach elevated concentrations of up to 20 pptv during J-J-A ( Figure 1A) which reflects the increased industrial activities of these areas. Europe has background CS 2 concentrations of between 3-7 ppt with several hot spots of high CS 2 concentrations reaching up to 15 ppt over Eurasia (e.g., Turkey) and the Middle East (e.g., Iraq, Saudi Arabia) ( Figure 1A) most likely attributed to the oil, gas, and petrochemical industries as well as increased anthropogenic pollution. There are a few hot spots of CS 2 (e.g., 10 ppt) in the southern continental regions such as South Africa ( Figure 1A). The elevated CS 2 concentrations around South Africa are expected because of its significant petrochemical, plastic, and synthetic fibre industries [55]. In South America, there appear to be two areas (e.g., near Rio de Janeiro and Northern Brazil) of increased CS 2 concentrations (2-3 ppt) are most likely to be due to increased anthropogenic emissions partly due to the high population densities in these places.

Figure 1. The average surface level distribution of CS 2 simulated by the STOCHEM-base for (A) J-J-A and (B) D-J-F seasons.
During D-J-F, there is a significant increment in concentrations over the continental landmasses (e.g., up to 40 ppt in south East Asia, Figure 1B). This is to be expected as removal by OH accounts for 88% of  Figure 2B), respectively. These increased amounts of CS 2 can enhance the oxidation products (e.g., OCS, SO 2 ) which can alter the formation of the global atmospheric sulfate layer.

Figure 2. The percentage change of surface CS 2 level (A) from STOCHEM-PHOT to STOCHEM-OCSL and (B) from STOCHEM-PHOT to STOCHEM-ANTH.
CS 2 is the most important precursor gas to OCS, for which a large missing source is needed to balance sources and sinks [28][29][30][31][32]. The flux analysis using the STOCHEM-base model shows that up to 12 Tg OCS /yr per grid box can be formed from the reaction CS 2 + OH in the regions ( Figure 3A) with very high anthropogenic pollution (e.g., South east Asia). Incorporating the global anthropogenic emissions of CS 2 [54] in the STOCHEM-ANTH model produces up to 40 Tg OCS /yr per grid box in South East Asia ( Figure 3B). This revised production of OCS from CS 2 can play an especially important role over South East Asia regions, where comparisons between measurements and anthropogenic emission inventories show a significant missing source of OCS [11].

Global vertical distribution of CS 2
The vertical profile data in STOCHEM-base (Figure 4) shows the decay of CS 2 concentrations as altitude increases. The J-J-A season vertical distribution shows higher levels of CS 2 up to 3 ppt at 30°N-45°N ( Figure 4A) due to the larger surface concentrations over the latitude due to areas of high anthropogenic emissions from Europe and North America. A very small area (0-15°S) of higher CS 2 (1.5 ppt) is visible in the tropics ( Figure 4A) because of the increased oceanic emissions. Heading south from 30°S latitude to 90°S, CS 2 is found in small amounts throughout the troposphere. In general, the model OH concentration decreases from the equator for each altitude bin, and so the loss of CS 2 via reaction with OH decreases in a similar way ( Figure 4A). During D-J-F, the CS 2 concentration is also found to be higher (up to 4 ppt) in the region of 30-55 o N ( Figure 4B) due to the larger anthropogenic activity, but the vertical distribution of CS 2 for D-J-F mirrors that of J-J-A; lower OH concentration in the NH during D-J-F leads to lower removal of CS 2 and therefore it can reach higher altitudes in the region of 60-90 o N. During D-J-F, the peak in the tropics disappears due to the decreased emissions from ocean but another peak at the region of 50-60 o S appears due to the increased emissions from SH oceans.

. The average zonal distribution of CS 2 simulated by the STOCHEM-base for (A) J-J-A and (B) D-J-F seasons.
The vertical CS 2 distribution shows similarities to the vertical profile found by Kjellström [35]. Both profiles display peaks around 60 o N and 60 o S, but the concentrations seen in this study are found to be smaller than that in Kjellström [35]. Ultimately, as the emission data and the resolution of the two studies are different, it is unsurprising that the absolute concentrations are different, yet what is encouraging is that despite these factors, the annual vertical distributions are similar.

Model-measurement comparison
In order to develop an understanding of the model accuracy, data points were extracted from the model and compared with 12 field measurements and 9 flight data sets ( Table 3). The model uses meteorological data from 1998, but the field and flight data sets were taken from different studies for 1979-2005 which were combined in the same plot panels in Figures 5 and 6, this could involve potential artefactual patterns from using different measurement methods/sampling times. The simulated CS 2 data shows no clearly defined seasonal cycle ( Figure 5). The seasonal cycle of CS 2 is controlled by seasonality in their emissions, losses (by OH, deposition) and atmospheric transport. The loss process due to higher OH concentrations during summer would produce the seasonal cycle of summer minimum with winter maximum. However, the oceanic emissions would be highest in summer months and lowest in winter months. Little is known about the production mechanisms of CS 2 in the ocean (one of the major sources of CS 2 ), but evidence suggests that the production can be through photochemical and biological processes [6,7]. The competition between oceanic emissions and loss by OH resulted in a complex seasonal cycle of CS 2 which is reflected in Figure 5.  The area best modeled is the marine boundary layer likely due to this area being subject to insignificant amounts of anthropogenic emissions, but an underprediction with a mean bias of −2.7 ppt for the modeled data compared with the measured data is found when looking at regions within marine clean environments (e.g., North Atlantic, Southern Ocean, Pacific, eastern Indian and Southern Ocean). The global oceanic emissions of CS 2 are highly variable ranging from 3.3 to 2346 Gg/yr for different marine environments [8], but we used the global oceanic emissions of 180 Gg/yr in the STOCHEM-base case which can explain the discrepancies between model and measured CS 2 in the marine environment. Increasing the global soil and oceanic emissions [3] in the STOCHEM-OCSL brings the model into closer agreement (mean bias reduces to −0.6 ppt) with measurements for these marine campaigns. The oceanic emissions of CS 2 are unclear because of lacking process understanding, a further investigation into how the discrepancies can be improved would be highly desirable.
The correlation between modeled and measured CS 2 is found to be very good for some clean areas (e.g., Jeju Island in Korea, Fremantle in Australia), but the CS 2 concentrations in approximately background air (relatively unpolluted air) arriving at the east coast of the USA, Philadelphia and Wallops Island and in air arriving at rural areas Harwell, UK and Tai O, Hong Kong are reported as 36 pptv in November, 40 pptv in December, 15 pptv in January and 47 pptv in March-April-May, respectively [56,59,60] which equates to a large underprediction (mean bias of −24.7 ppt) of the simulated values by STOCHEM-base ( Figure 6). Harwell and Tai O are located in the rural environment, but there are nearby industrialized areas which may contain large anthropogenic sources of CS 2 which can be advected into the rural areas increasing the ambient rural conditions.
Increasing the global anthropogenic emissions of CS 2 in the STOCHEM-ANTH model improves the agreement with measurements (the mean bias reduces to −1.0 ppt) for these locations. Relatively low concentrations of simulated CS 2 are found in Fiji which are likely to be due to the absence of significant anthropogenic sources surrounding the area. The measured data collected near Fiji recorded concentrations as high as 12 pptv, because there are significant biomass burning events in the area [62], which may have skewed the results and therefore may be unrepresentative of Fiji's natural ambient condition. Due to limited data within the remote SH and hence clean atmospheric conditions, this hypothesis couldn't be tested any further. Increased experimental data within the remote SH would help to validate global atmospheric models and perhaps filling this void should be a focus of future expeditions.
Air samples collected from industrial areas (e.g., Beijing, China and Ansan, South Korea) found very high concentrations of CS 2 due to the anthropogenic emissions, giving very large underprediction of the STOCHEM-ANTH modeled values with a mean bias of −490 ppt compared with the measured values. Due to the large difference in magnitude it is likely that there may be additional sources of CS 2 currently being under-estimated or not being accounted for, suggesting primarily that anthropogenic emissions are likely to be much larger than current estimates. In addition, the coarse resolution of STOCHEM model (5° latitude  5° longitude) makes it difficult to accurately model CS 2 in highly polluted areas which make up only a small fraction of the grid. The meteorological data incorporated in STOCHEM was from the year 1998 which was one of the strongest El Niño years in the past 50 years, this could lead to differences in the model data when compared with measured data from other years. Consideration of inter-annual variation give better representation of the model data which could improve the disagreement between model and measured data. Figure 6 shows the vertical distributions of CS 2 over the North Atlantic, South tropical Atlantic, Western Pacific, Central/Eastern Pacific, Wallops Island-Virginia to Natal-Brazil, and Japan sea coast with the simulated data overlaid. The modeled CS 2 at different altitude levels for North Atlantic and South tropical Atlantic fit within one standard deviation of the measured mean except one outlier at 3.5 km. The STOCHEM-base case simulation underpredicted CS 2 by on average 1.5 ppt in North Atlantic and South tropical Atlantic and increasing the global ocean and soil emissions in the model, STOCHEM-OCSL reduces the underprediction by 47%. In the North Atlantic, the concentration at the surface is almost four times higher than that found in the South Atlantic because of the possible higher anthropogenic emissions of CS 2 in NH than that in SH. The global distribution of CS 2 ( Figure 1) indicates that Europe and the east coast of the US have much higher concentrations compared with South America and Africa and thus through mixing and advection might have influenced CS 2 concentrations above the ocean [64]. In addition, the area of the South Atlantic that was sampled in this study was closer to the equator, thus higher OH concentrations yielding increased CS 2 removal and lower atmospheric concentrations.
In extremely polluted environments e.g., the Japan sea-coast (part of the NASA Transport and Chemical Evolution over the Pacific (TRACE-P) campaign to investigate outflow from polluted Asian regions and clean air from the Pacific Ocean), STOCHEM-base underpredicted the CS 2 concentrations in the Japan sea-coast with a mean bias of −64 ppt suggesting that the emissions of CS 2 integrated in the model are likely to be underestimated. The long-range transport of CS 2 from the source regions on the Asian continental outflow mixing with marine air over the Japan Sea can be responsible for the higher CS 2 over Japan sea coast [63]. The coarse model resolution (5°  5°) of STOCHEM is responsible for model underprediction over the Japan sea-coast because it cannot represent typical Asian continental pollution with sufficient spatial accuracy. Increasing the global anthropogenic emissions of CS 2 in the STOCHEM-ANTH model reduces the model-measurement mean bias to −55 ppt. There are some other marine flight campaigns as part of TRACE-P e.g., Western Pacific and Central/Eastern Pacific, where the mean bias for STOCHEM-base is found to be −3.2 ppt. STOCHEM-OCSL and STOCHEM-ANTH simulations reduces the mean biases to 22% and 53%, respectively. In the study, a more comprehensive regional comparison with the literature was intended but there is a lack of up-to-date measurement data, more field and flight measurement data would be invaluable to validate our model accuracy.

Conclusion
The global analysis for CS 2 was performed using a 3-D chemistry and transport model, STOCHEM-CRI which was concordant with several literature sources suggesting that the model results show a tendency towards underestimation. A suit of simulation results give the global burden and lifetime of CS 2 as 6.1 to 19.2 Tg and 2.8-3.4 days, respectively. The global surface distribution of CS 2 displayed the expected characteristics; with highest concentrations over areas of industrial activities. Comparing the global surface and vertical distribution of CS 2 simulated by the STOCHEM-CRI with other model studies suggested that the trends are found to be in good agreement with the concentrations simulated by IMAGES and ECHAM3. The oxidation of CS 2 produces a large amount (0.5-1.9 Tg/yr) of OCS which can reduce the imbalance between the sources and the sinks of OCS. Increasing oceanic emissions in the model reduces the model-measurement disagreement significantly for the marine campaigns. The major underprediction by the model in polluted areas is likely due to the level of anthropogenic emissions being much greater than previously thought but also due to the coarse resolution of STOCHEM. Currently many of the campaigns have looked at very similar areas, for example the Global Tropospheric Experiment/Chemical Instrumentation Test and Evaluation (GTE/CITE 3) and Global Atmospheric Measurements Experiment of Tropospheric Aerosols and Gases (GAMETAG) expeditions have carried out extensive sampling around the East coast of USA, the Atlantic Ocean and down to South America. The TRACE-P study looked in detail at China, South East Asia, Korea and Japan, yet there has been very little sampling from Europe, the Middle East and the SH in general. Such experimental data will be the key to validate our model results.