Photocatalytic chlorine atom production on mineral dust–sea spray aerosols over the North Atlantic

Significance Using a combination of field data and global modeling, we demonstrate a mechanism in which a mix of Sahara dust and sea spray aerosol activated by sunlight produces large amounts of active chlorine. This mechanism resolves a number of unexplained observations and significantly revises our understanding of atmospheric chlorine, reducing uncertainties in the source budget. The chlorine formed by this mechanism impacts two important greenhouse gasses, methane and tropospheric ozone, with an estimated catalytic efficiency of removing ca. 45 methane molecules per iron atom per day. The inclusion of Cl. made from the photocatalytic oxidation of ocean chloride in models will reduce critical uncertainties in estimates of methane emissions and improve our ability to predict future climate change.

atmosphere because its composition is also changed by the OH and Cl sinks. Due to the KIE, oxidation of CH 4 by OH and Cl enriches the remaining CH 4 in 13 C, leading to an observed δ 13 C-CH 4 value of, for example, −47.1‰, thus making the observed δ 13 C-CH 4 seem to include less biological methane sources than what was actually emitted. Here, we present a significant chlorine source, and if it is not included in isotope-based emission budgets, these models will underestimate the biological methane source. The Cl-KIE makes top-down models especially sensitive to the Cl + CH 4 reaction; each 1% increase in CH 4 loss via Cl increases mean δ 13 C-CH 4 by ~0.5‰ (18). To put this in perspective, the recent rapid increase of atmospheric methane (3) is associated with a change in the observed δ 13 C-CH 4 from −47.1‰ in 2007 to −47.5‰ in 2022 (21). A better understanding of the Cl sink of CH 4 will significantly improve our understanding of the sources responsible for the recent CH 4 increase. If current models underestimate the source of Cl, they would compensate by underestimating the relative contribution from 13 C-depleted biological sources and overestimating the contribution of 13 C-rich fossil sources (19,20).
Anomalous results from a number of field studies could be explained using a so-far-unknown, episodic source of Cl. Read et al. measured the intra-annual cycles of nonmethane volatile organic compounds and other species at Cape Verde and found that CO and ethane, which have similar lifetimes with respect to OH degradation, vary intra-annually with the sinusoidal variability expected due to their reaction with OH but with CO showing a smaller cycle amplitude (22). This is unusual because the OH sink will not discriminate between CO and ethane, so the ratio of CO to ethane should be stable. While Read et al. suggest that the anomaly in the CO:ethane ratio is due to substantial chlorine chemistry, there was no known chlorine source that would operate along the entire 50-d lifetime from source to detection, and which has significant seasonal variation. In separate work, Lawler et al. measured HOCl in marine air at Cape Verde and found unusual HOCl levels that could be explained only by adding an unknown Cl source that is at 100 to 1,000 times stronger than conventional acid displacement (23). Finally, Saharan dust has been observed to affect O 3 in the free troposphere at the Monte Cimone Climate Observatory, with implications for surface air quality in cities of northern Italy and southern Spain (24).
Here, we present field and modeling evidence of the mechanism of the production of atomic Cl via the photocatalytic oxidation of chloride in aerosols containing Sahara mineral dust.
By this mechanism, Cl 2 and Cl are generated when lofted iron-bearing mineral dust aerosol from North Africa descends into the marine boundary layer (MBL) over the Atlantic and mixes with sea spray aerosol to form Mineral Dust-Sea spray Aerosols (MDSA). We combine data from field with global atmospheric modeling and predict extremely low δ 13 C-CO values that match those seen in CO in air samples from Barbados (25); these results remained unexplained for 20 y. Finally, we discuss the global significance of this mechanism that is not yet included in global models.

Results
Barbados Observations. Mak et al. (25) saw the largest 13 C depletions in CO ever recorded, in air samples taken on Barbados (13.2 °N, 59.5 °W) from 1996 to 1999 (Fig. 1). The Barbados samples were depleted in 13 C by up to −3‰ relative to the other northern hemisphere locations in episodes occurring from July to December, reducing the seasonal minimum δ 13 C-CO at Barbados to −34‰. For other locations, for example, Tenerife (28.3 °N, 16.5 °W), the seasonal minimum of δ 13 C-CO remained above −31‰ (Fig. 1). As seen in the fig., the seasonal maximum at Barbados was similar to that seen at the other locations, around −26‰. Thus, the depletion events are seasonal, episodic, and location specific. Mak et al. (25) suggest the observed changes in 13 CO and C 18 O are evidence of the reaction of Cl with CH 4 . CO is the first stable product in atmospheric methane oxidation, and δ 13 C-CO is sensitive to even small CH 4 removal by Cl because atmospheric (CO) is more than an order of magnitude lower than [CH 4 ] (8) and because of the large KIE of the Cl reaction. They estimated that the data could be explained if 3 ppbv of CH 4 were destroyed by reaction with [Cl] = 5 × 10 4 atoms cm −3 over 3 d. However, the known Cl sources could not explain the seasonal depletion in 13 CO (25).
We incorporated the MDSA mechanism into the global 3-D chemistry-climate model Community Atmosphere Model with Chemistry (CAM-Chem) (Methods) (v4) to investigate whether Cl generated by MDSA photocatalysis could be responsible for the CO isotope signal observed by Mak et al. (25). CAM-Chem is included in the well-documented CESM (Community Earth System Model) framework (26,27) and includes an updated scheme for halogen chemistry (Cl, Br, and I) (13 observations from the North Atlantic region near Barbados, which makes it useful for describing the CO isotope anomaly at Barbados, although interpretation in other regions may be less reliable. Output from our MDSA model shows good agreement with the timing, episodic variability, and magnitude of δ 13 C-CO observations of Mak (Fig. 1). We will discuss this in detail after we first introduce the mechanism.

The Mechanism of Chlorine Production by Sahara Dust
Photochemistry. MDSA is formed when iron-bearing mineral dust and sea spray aerosols mix (28). Iron is released from the minerals by weathering, and according to the Pourbaix diagram (29), iron oxides and hydroxides, abundant in the crust and in mineral dust, are converted into Fe(III) ions under the oxidizing and acidic conditions found in well-mixed marine and dust aerosols found in the North Atlantic (30). Acidity favors formation of Fe(III) chlorides (FeCl n 3-n with n = 1..4) by keeping hydroxide concentrations low (31). Fe(III) chlorides have an intense ligand to metal electron transfer absorption activated by sunlight, yielding Fe(II) and oxidized chlorine Cl (14). The Cl atoms can escape directly from the aerosol to the gas phase but most likely combine to escape as Cl 2 , which is then quickly photolyzed by sunlight (31). Fe(II) is reoxidized to Fe(III) by O 3 , H 2 O 2 , and other species, making the process catalytic in iron by reforming the iron chloride species (29). When Cl reacts with methane, it forms HCl, which can be reabsorbed by the MDSA, making the process catalytic in Cl as well (see SI Appendix, Fig. S1 for a scheme of the mechanism). While this mechanism resembles the release of iodine from dust reported by Koening et al. (32), in this case, the halogen source is the large reservoir of chloride from sea spray already present in the aerosol phase in contrast to deposition of iodine-containing atmospheric gases on the mineral surface. This critical difference means that MDSA-chlorine is produced close to the surface in the MBL from abundant sea spray aerosol. Using a smog chamber, Wittmer et al. (31) demonstrated that formation of Cl atoms by photooxidation of chloride occurs in aerosols that form when iron oxide is mixed with sea salt under acidic conditions. Until now, this mechanism has not been observed in situ and its global impact was unknown, although the necessary elements are found in the troposphere: Photoactive iron is known to exist above the North Atlantic (33)(34)(35) and mineral particles in the MBL accumulate oxidation products and salts (including sea salt) on their surface within a day (28). Zhu et al. (34) collected mineral aerosol at Barbados and measured the daily variability of Fe(II). They saw a clear diurnal pattern in the concentration of soluble Fe(II), with daytime values (mean 3.7 ng m −3 ) about twice nighttime values (1.5 ng m −3 ), showing that photochemical processing of Fe takes place. The same researchers found that a photochemical steady state is established in which Fe(III) is photoreduced to Fe(II) as quickly as Fe(II) is oxidized to Fe(III) (33), and they showed that there are at most 11.4 Fe(II)-Fe(III) cycles h −1 at maximum solar intensity. This rate is at the low end of the range reported by Wittmer of 6 to 78 h −1 (31); atmospheric conditions could result in lower values, more similar to those assumed here, due to the suppression of this reaction in the presence of environmental sulfate (36). Several mechanisms for the photochemical reduction of Fe(III) are known including the oxidation of organics yielding H 2 O 2 (33), but we assume in our model parameters that the dominant photoreduction reaction in MDSA yields OH and/or Cl that both lead to Cl 2 production (31). Considering that these dust particles will be mixed with sea salt and assuming the irradiated combination gives rise to oxidized chlorine via the MDSA mechanism, each cycle would produce one Cl atom. The agenda is to i) determine whether Fe photochemistry in Sahara dust is able to produce Cl atoms, ii) model the MDSA mechanism over the North Atlantic, iii) compare the model prediction with the observed CO isotope anomaly in Barbados, and, if successful, iv) estimate the global impact of MDSA photochemistry.
Model Results. The MDSA mechanism is implemented in the CESM assuming: i) The cycling rate observed near Barbados between Fe(II) and Fe(III) is fully equated with the rate of Cl production, and ii) the influence of sulfate is reflected in these observations because FeSO 4 (which would reduce the Fe available for Cl production) does not take part in the observed photochemistry cycles (36). We ran the model with and without MDSA chemistry and used a tracer molecule to derive the additional CO generated by the Cl + CH 4 oxidation due to the addition of MDSA. We use this to calculate the change in δ 13 C-CO and add the result to the same inverse model as was used by Mak et al. (25) (see Methods for details). We then evaluate how this affects the predicted seasonal minimum and maximum of the inverse model.
The model results show that MDSA photocatalysis produces 110 ± 10 Cl atoms (active Fe atom) −1 d −1 on average at Barbados, which is 70 g Cl 2 (g active Fe) −1 d −1 . The modeled average cycling rate between Fe(II) and Fe(III) of 11 h −1 at noon, is similar to that of the observations, consistent with the implemented reaction rates observed by Zhu (33). According to the model results, the average concentration of photoactive iron in MDSA at Barbados is 2.5 ng m −3 with peak concentrations up to 34 ng m −3 [similar to observations reported by Zhu (34)], resulting in a Cl 2 production rate of up to 6 × 10 5 cm −3 s −1 (average 2 × 10 4 cm −3 s −1 ) and peak [Cl] of >10 5 cm −3 .
The model output shows good agreement with observation (25) for the location, timing, episodic variability, and magnitude of δ 13 C-CO (Fig. 1). The MDSA model reduced the seasonal minimum of weekly average δ 13 C-CO at Barbados by −2‰ in 1996 and by −4‰ in 1997. Because of the episodic nature of the dust intrusions, the springtime dust does not significantly change the seasonal maximum, changing it by only −0.1‰ and −0.3‰ in 1996 and 1997, respectively. The large episodic variations in δ 13 C-CO at Barbados are also seen in the observations. In our model, at Spitzbergen, both the seasonal minimum and maximum did not change by more than −0.3‰, confirming that the MDSA impact is local and specific to the western North Atlantic. Cl formation in the model depends on specific conditions, including dust concentration, wind speed/sea spray aerosol loading, atmospheric circulation, wind direction, and actinic flux. Mak et al. (25) suggested that chlorine is the cause of the 13 C anomaly in CO, and our CESM simulation is consistent with our hypothesis that the Cl is associated with Sahara dust episodes.
Our model results show that seasonal Cl 2 production rates up to 6 × 10 5 cm −3 s −1 are needed to explain the CO isotope anomaly at Barbados. Other known Cl production mechanisms are not able to produce a seasonality of this magnitude. For example, acid displacement in sea salt forms gaseous HCl which can then react with OH becoming a large source of Cl above the oceans, Cl (14,37). We estimate that the maximum Cl production from this source is 3 × 10 4 cm −3 s −1 (Methods), which is more than an order of magnitude too small to explain the observations. In addition, the global seasonality of acid displacement as a Cl source comes from solar-intensity-driven OH and DMS seasonality (6), which makes acid displacement seasonality more apparent away from the equator, and therefore does not match with the uniqueness of the isotope anomaly at Barbados. Finally, nonsea salt sulfate levels (and associated HCl levels) at Barbados reach their seasonal maximum in June/July (38), which is opposite to the observed CO isotope effect. Other substantial Cl sources such as ClO and HOCl have at most a similar order of magnitude and are part of the Cl radical chain with their seasonality depending on an initial input of Cl from sources like OH + HCl. We therefore conclude that MDSA photocatalysis is the necessary mechanism for producing Cl and with the right concentration and spatial and temporal structure to explain the observations. As was already observed by Mak et al. (25), the observations in Barbados are depleted compared to the observations in Tenerife (Fig. 1). In agreement with this observation, the model predicts much smaller MDSA-related dips in δ 13 C-CO at Tenerife (Fig. 1B) and less frequently. Braunlich (39) discussed the Tenerife 13 CO record in detail and noted that the summer minimum in Tenerife was lower in 1997 compared to 1996, 1998, and 1999. This was also observed in stations at Spitzbergen and Sonnblick, but no explanation could be found. This period aligns with the prediction by our model of significant Cl production by the MDSA mechanism in Tenerife for August 1997.
The modeled annual average shift in δ 13 C-CO across the North Atlantic region shows a shift of up to −4‰ at 30 °W over the North Atlantic ( Fig. 2; monthly in SI Appendix, Fig. S2). The zone of maximum shift changes through the year depending on winds in the Inter-Tropical Convergence Zone that carry Sahara dust to Barbados from June through December. Even though more dust is transported across the Atlantic in summer, the model predicted that the strongest δ 13 C-CO excursions at Barbados would occur in the fall as is observed. This is because in summer, most dust is transported at high altitude, above the MBL, which limits mixing with sea spray aerosol (SI Appendix, Fig. S3). This agrees with satellite observations (40) and in situ measurements (41) that show that the Saharan Air Layer transports dust at a higher altitude and velocity in summer than fall. Depletion in 13 C-CO is also seen in the spring (Fig. 1). While Sahara dust is farther south of Barbados this time of year, its impact is still strong enough to influence δ 13 C-CO at Barbados. However, because of the episodic nature of the dust intrusions, the springtime dust photocatalysis does not significantly change the seasonal maximum.
MDSA photochemistry increases the monthly average [Cl] to >10 4 cm −3 over the equatorial Eastern Atlantic (Fig. 3A) greatly increasing the fraction of CH 4 oxidized by Cl. Fig. 3B shows the reduction in O 3 , which is a precursor for OH and leads to a reduction in OH shown in Fig. 3C. The overall effect on CH 4 loss is the sum of the additional loss via Cl and the reduced loss via OH. For most of the North Atlantic, we find that CH 4 removal increases with the addition of the MDSA mechanism ( Fig. 3D  and SI Appendix, Fig. S7); the local CH 4 loss rate increases by 20%, even though OH is reduced by up to 10% (Fig. 3C and  SI Appendix, Fig. S6).
For the North Atlantic region (0 °N to 35° N; 90° W to 0° W), an annual average burden of 0.22 Gg photochemically active Fe in MDSA in the MBL above the ocean surface resulted in a Cl 2 production of 2.9 Tg y −1 [36 g Cl 2 (g active iron) −1 d −1 ], resulting in a net CH 4 loss of 0.6 Tg y −1 [7.5 g CH 4 (g active iron) −1 d −1 ], from an increased Cl-mediated CH 4 loss of 1.1 Tg y −1 and a reduced OH-mediated loss of 0.5 Tg y −1 (SI Appendix, Table S1). The total burden of tropospheric O 3 was reduced by 6% above the North Atlantic (SI Appendix, Table S1).
These model results suggest that MDSA produce a large amount of Cl 2 that is not considered in existing Cl budgets (11): 3.8 Tg Cl 2 y −1 across the entire North Atlantic (continental and marine), of which 2.9 Tg y −1 are over marine areas. Globally, the MDSA Cl source increased the tropospheric inorganic Cl source in our CESM simulations by ~41%, adding 13 Tg y −1 of chlorine to 23 Tg y −1 from sea-salt cycling, 8 Tg y −1 from anthropogenic emission, and 1 Tg y −1 from photochemical oxidation of organic chlorine. Because the MDSA source is very strong over the North Atlantic, this increases the Cl contribution to CH 4 oxidation by at least fourfold, making the MDSA mechanism the dominant source of atomic chlorine in this region. Outside of the North Atlantic region, our model predicted Cl 2 production by MDSA in dust from the east coast of Argentina, as well as in the Gulf of Aden, around India and around Australia.

Discussion
Li et al. (13) found that the effect of chlorine in the atmosphere today is to reduce OH production and increase its loss particularly in the free troposphere and over the open ocean, leading to a net increase in the CH 4 burden. In this work, using the same model, we obtain the same initial result for the global average behavior. In addition, we see that atmospheric behavior changes during intense MDSA episodes: Chlorine production continues to remove methane after tropospheric ozone levels (the precursor for OH) have been reduced. Thus, after an initial increase, there is a significant decrease in methane lifetime in regions with high Cl  Table S3. Many atmospheric species are sensitive to changes in chlorine. Future research may be able to show that the MDSA mechanism can explain a variety of phenomena in particular at the regional scale, including ozone suppression, VOC oxidation, VOC oxidation-driven secondary organic aerosol formation, HCl formation, etc. Additional Cl from MDSA would be able to explain the variation in the CO:ethane ratio observed by Read et al. (22), and an additional chlorine source is consistent with Lawler et al.'s observation of elevated Cl 2 and HOCl in air masses arriving at Cape Verde (23). Lawler observed HOCl up to 100 ppt at Cape Verde and estimated that the unknown Cl source responsible for the increase in HOCl is at least one order of magnitude greater than acid displacement of HCl, which is in agreement with the magnitude that we found in this study. Over the North Atlantic, the monthly average HOCl levels in our model output were increased from 0.3 to 2.0 pptv (without MDSA) to 0.9 to 52.1 pptv (with MDSA), during our simulation period. Cl from MDSA reduces the annual average tropospheric ozone burden in the North Atlantic MBL by 6%. The substantial local MDSA-mediated ozone loss helps to reconcile observations and models by providing a mechanism for dust-induced tropospheric ozone depletion in air masses rich in mineral dust from the North African desert (43). Our model includes the stratosphere but found no impact on stratospheric ozone.
Our model showed that Cl from mineral dust/salt aerosols increases CH 4 loss via Cl by 4.8 Tg y −1 globally (0.9% of total CH 4 loss). This is important given the aforementioned uncertainty in estimates of the fraction of CH 4 removed by Cl (9)(10)(11)(12)(13)(14)37). The additional 0.9% Cl oxidation of CH 4 shifts global δ 13 C-CH 4 by ~0.5‰ (18). Adding this into top-down methane emission models shifts 12 Tg y −1 of estimated CH 4 emissions (2% of the global total emissions) toward 13 C-depleted biological sources such as agriculture, potentially reducing the estimated fraction of fossil-fuel methane emissions by 7% and increasing the ruminant contribution by a corresponding amount (19).
A doubling of this large change in methane isotopes from the MDSA mechanism is likely to have occurred over the 20th century since available in situ and paleoclimate data suggest an almost-doubling of the transport of North African dust into the North Atlantic, especially since 2014 (43,44). The large changes in atmospheric dust loading seen over geological time scales show that it is highly variable. Another study predicts a 60% reduction in dust from the Holocene to today (45). These changes could be expected to affect chlorine production via the MDSA mechanism. Depending on the conditions where the chlorine is added, this could modestly increase or decrease methane lifetime. Perhaps more importantly, enhancing the Cl sink of methane shifts our estimates of methane sources. It is thus extremely important to accurately quantify this source as part of efforts to understand and control atmospheric methane. Using a box model, we calculated the change in the atmospheric 13 CH 4 isotope (21) in response to recent changes in North African dust emissions, by scaling the Cl production of MDSA linearly with the reported dust emission (44). Fig. 4 shows that these fluctuations in dust emission give 13 CH 4 isotope changes that are significant compared to typical yearly changes (16,17) and could in that way mask changes to emission sources. In conclusion, the MDSA chlorine source should be verified in future laboratory and field studies and should be included in methane emission models to improve their accuracy, and further studies should be done to evaluate the importance of the MDSA mechanism in other regions.

Materials and Methods
CESM Model Setup. The model employed in this work is a global chemistry-climate model CAM-Chem, which is included in the well-documented CESM framework (CESM) (26,27), and incorporates an updated halogen chemistry scheme for chlorine, bromine, and iodine (14). In this work, CAM-Chem was configured with a horizontal resolution of 0.9° latitude by 1.25° longitude and 56 levels from the surface to the stratosphere.

MDSA Chemistry Implementation in CESM.
For this implementation of the MDSA, we use a simplified mechanism to describe the overall process based on results from field studies. At present, there is no model that combines iron aerosols and chlorine chemistry to describe the elementary processes in detail. In our parameterization, the amount of MDSA was calculated by estimating the amount of photoactive iron. In studies of Fe(II) and photoreducible Fe(III) in North Atlantic dust aerosols, it was observed that photoreducible Fe(III) is found mainly in fine dust particles of sizes below 2.5 µm and across a wide range of dust loads (46)(47)(48). These observations imply that a substantial fraction of these dust/salt aerosols are in the acceptable pH range of the Pourbaix diagram (29) for the reoxidation of Fe(II) to Fe(III). We used these observations as input for our model. We used the size bin of 1 to 2.5 µm from the model (referred to as χ DST02 below) because at least 90% of Sahara dust <2.5 µm is in the 1 to 2.5 µm range, and χ DST02 therefore provides us with a good reference value to compare with the observations. Additionally, the smallest size fraction is often associated with high sulfate concentrations (30), which could suppress the MDSA (49). In situ observations suggest that most fine particles downwind of the Saharan plume will be acidic enough for the MDSA mechanism (30, 50).
The following equation was used to calculate the amount Fe in MDSA.
where α = 0.063% is the mass fraction of photoactive Fe in the dust, c is a pressure-dependent factor converting χ from kg kg −1 in the CESM model to kg m −3 , and γ is a parameter set to 1 or 0 depending on the availability of Cl as described below.
This mechanism should be active only in downwind regions from dust plumes and in the MBL since this will assure a high-enough Cl:Fe ratio, above 13 for MDSA concentrations up to 90 ng m −3 (31). In order to ensure these characteristics, we limit the Cl production using the constraints ρSSLT02 > 10 −9 and altitude below 900 hPa, by setting γ = 1, and otherwise setting γ = 0. Mineral particles in the MBL accumulate oxidation products and salts including sea salt on their surface (28), with 80 to 90% of the dust particles becoming mixed with sea salt within a day.
The factor α is calculated by multiplying the fraction of Fe in the dust (3.5%) with the fraction of the iron that is soluble and photoactive (1.8%). We used 3.5% as the amount of iron in the Saharan dust, based on the findings of Trapp et al. 2010 (48), as a typical value for dust-sea salt mixtures across the North Atlantic.
Observational evidence suggests an inverse relationship between dust amount and solubility (51). The two explanations for this behavior are 1) differences in particle size (smaller combustion-sourced particles having higher solubility than dust-sourced particles) and 2) downwind atmospheric processing of particles increasing the solubility and the higher gravitational settling rate of larger particles (52). However, the fraction of the total iron that is photoactive does not appear to follow this inverse relationship at locations downstream from the dust sources (48) because observations across the North Atlantic show a similar fraction of total Fe to be photoactive, for both high and low dust loads (46)(47)(48). Therefore, we assume a constant value for the fraction of the iron that is photoactive for the long range-transported dust of 1.8% based on cruises far from North Africa (35,(46)(47)(48) and average values for the model-data syntheses for dust in the North Atlantic . This observed cycling of Fe(II) and Fe(III) during the day and night demonstrated that photochemical processing of Fe does take place at Barbados, and we postulate here that it is entirely due to the Cl mechanism. Next, the MDSA number w MDSA is converted to d MDSA in mol Fe m −3 , using M w (Fe) = 55.845 g mol −1 .  (33). The MDSA aerosol chemistry is described by Wittmer (36,49). Fe(III) chlorides FeCl n 3-n (n = 1..4) have intense UV charge-transfer absorptions. On absorption of light, an electron is transferred from chloride to Fe(III), resulting in photoreduction of Fe(III) to Fe(II) and yielding chlorine atoms (Cl) in the liquid phase. Chlorine atom recombination is catalyzed by Cl − via Cl 2 − leading to Cl 2 which degasses. In summary, Cl production occurs at the rate at which iron is cycled between Fe(II) and Fe(III) within d MDSA . Assuming steady state, the fraction of d MDSA that is in the Fe(II) state is given by [3] Fe(II) fraction = k 2 In which k 1 is the reaction rate constant for oxidation of Fe(II) to Fe(III), and k 2 is the reaction rate constant for photoreduction of Fe(III) to Fe(II) at solar intensity I 0 = 600 W m −2 , and I is the temporally varying solar intensity.
At steady state, the production of Cl can be calculated from either Fe(II) oxidation or Fe(III) photoreduction. Using Fe(II) oxidation, the Cl production rate (amount per unit time; r Cl ) is Combining Eqs. 1 to 4 and using two Cl atoms to create one Cl 2 molecule gives the following rate of Cl 2 production:  limiting factor, we used a photoreduction reaction rate of 10 times as fast: k 2 = 1.9 min −1 . The resulting Cl 2 production rate is in min −1 as well. These values of k 1 and k 2 will result in a maximum of 10.2 photochemistry cycles per hour at maximum sunlight, which is well within the range of cycles measured by Wittmer (31) of 6 to 78 Cl h −1 .
Analysis. To monitor the reaction of Cl + CH 4 , we added a tracer to the model with a lifetime of 60 d, representative of the lifetime of CO in the North Atlantic region, i.e., 0° N to 35° N; 90° W to 0° W), based on Pfister et al. (53). We then ran the model with and without MDSA chemistry, using a stabilization time of 12 mo. Using a CO yield of 90%, the difference in the tracer molecule represented the extra CO generated by the Cl + CH 4 oxidation due to MDSA implementation.
To calculate the isotope shift in the CO produced via CH 4 + Cl, we used 100‰ depletion in 13 C relative to the ambient CO without the MDSA mechanism. The MDSA implementation also changed the amount of CO produced by CH 4 + OH slightly; the isotope shift from this effect was less than 0.05‰ for a 1 ppb change in CO, calculated using the Rayleigh fractionation equation, and was therefore further neglected in the analysis.
Sampling at Barbados and Inverse Modeling of CO Isotopes. Samples were collected weekly at Ragged Point in Barbados in 1996 to 1999, at the top of a 15-m tower, which sits atop a 15-m bluff overlooking the Atlantic. The vast majority of samples were collected during an east-south-east flow in the morning (8:00 to 12:00) with wind speeds in excess of 10 m s −1 . CO isotope analysis was done with an isotope ratio mass spectrometer. More details on the sample collection protocol and analysis of CO isotopes are provided in Mak et al. (25). Mak et al. (25) data were shown as 3-wk averages, but in our Fig. 1, we show the raw data that have not been published before. Fig. 1 shows the same inverse model simulations as used by Mak et al. (25), which is an inverse model simulation using TM2 from Bergamaschi et al. (54). It was noted that these model simulations are for a perpetual year reflecting the climatological wind fields. Model simulations were tuned to a set of previously published stable isotope data which were collected from middle-and high-latitude stations in both hemispheres.
We note that perfect agreement between model and observations is not expected as the Mak et al. (25) samples were taken with nominal conditions of 10 m s −1 wind from the east-south-east, giving an increased chance of sampling air parcels affected by dust, whereas the model data present δ 13 C-CO values predicted for the region regardless of local meteorological conditions. Calculation of Acid Displacement Source Magnitude. Using the reaction rate constant for HCl + OH, k(298 K) = 7.86 × 10 −13 cm 3 s −1 (55) and a relatively high mole fraction of x HCl = 400 ppt HCl (7) and a x OH of 0.14 ppt (double the diurnal average of 0.07 ppt, see SI Appendix, Fig. S6), we calculate a maximum Cl production rate from acid displacement of 3 × 10 4 cm −3 s −1 . This is more than 1 order of magnitude below the Cl production rates by MDSA that were needed to explain the Barbados anomaly. Data, Materials, and Software Availability. Text data have been deposited in University of Copenhagen repository (https://erda.ku.dk/archives/7b38cf06c-c971b21faecca41e8924f6b/published-archive.html) (56).