The air and seawater data in this study were collected during the POLAR-CHANGE cruise, which was conducted from 14 February to 17 March 2023 aboard the B.I.O. Hespérides around the Antarctic Peninsula and in the northern Weddell Sea. Data from this campaign were obtained from two main areas visited on different periods (Fig. 1): (i) near the coast of the Antarctic Peninsula and the Antarctic islands, referred to as coastal, and (ii) farther offshore, referred to as open ocean. The definition of the two areas were arbitrarily set by a 100 km limit from the coast. Two ASIT enclosure experiments were conducted, (A) in coastal waters and (B) in open ocean waters (Fig. 1).

The experimental setup was similar to that described in Sellegri et al. (2023). It consisted of two Air-Sea Interface Tanks (ASITs) installed on deck and connected to a set of air analysis instruments installed in a container close to the tanks (detailed description below) (Fig. 2). The tanks, with a total volume of 1.8 m3, were half-filled with seawater pumped to the ASITs from a 4 m depth (see details below). A continuous flow of clean aerosol-free air at 25 Lmin-1 was introduced into each of the ASITs via a flushing line, resulting in a residence time of 40 min. The outflows from the ASITs (purple and red lines, Fig. 2) were connected to the analytical instruments set via a four-way valve that switched every 20 min, allowing alternate measurements of the air from both ASITs, as well as the incoming air (from the flushing line; green line, Fig. 2), and the ambient air (yellow line, Fig. 2). During the POLAR-CHANGE cruise, one ASIT was covered with a lid transparent to all radiation wavelengths (Control-ASIT), while the other was covered with a lid blocking UV wavelengths below 380 nm i.e. removing UV light exposure inside the ASIT (noUV-ASIT) (Fig. S1 in the Supplement).

The ASIT headspace was continuously flushed with particle-free (filtered) air. In contrast, the lower section containing seawater was operated either in a continuous flow mode or in a static mode without water renewal. The continuous flow mode was used during most of the cruise with seawater being continuously pumped to the ASITs from a 4 m depth and flushed out with an average water residence time of approximately 6 h. This water continuous flow line supply is located 4–5 m below sea level, and consists of an epoxide-free silicone pipe connected to a centrifugal pump (BKMKC-10.11, Tecnium) in which the parts in contact with the fluid are made of polypropylene and glass. The pumping system is designed to be gentle and non-toxic to microorganisms and provide an uncontaminated continuous source of near-surface seawater. Simultaneous measurements of Chl a and DMS concentrations from the underway (UW) pumping system and a dedicated Go-Flo launched overboard, or the uppermost bottle from a conductivity, temperature and depth (CTD) cast have been compared during a previous cruse, and are in very good agreement (Royer et al., 2014). The static flow mode was applied in the two Experiments, A and B, in which the water remained incubated in the ASITs for 2 d.

To characterize the seawater entering the ASITs in continuous flow mode, samples were collected from the ship's underway water outlet on deck every 6–12 h (i.e. 2 to 4 times per day) throughout the entire cruise. This outlet was derived from the same intake water used to fill the ASITs. During experiments A and B, discrete water samples were also collected every 6 h from inside the ASITs (using a tap designed for this purpose), to monitor biological changes in the enclosed natural communities and to identify potential effects from UV light exposure. Cell abundances of microbial groups were estimated by combining flow cytometry and optical and fluorescence microscopy (described below).

The outgoing air line from the ASITs was connected to a set of instruments. An Ozone Analyzer (Thermo Scientific Model 49i) was used to monitor ozone. A Proton Transfer Reaction Quadrupole Mass Spectrometer (PTR-QMS, Ionicon) was used for VOC concentrations including DMS and MeSH. It was operated in standard conditions (Table S1 in the Supplement) and provided a E/N ratio of about 130 Td (with E the electric field strength and N the gas number density, 1 Td = 10⁻¹⁷ Vcm-2). Table S2 gives the m/z monitored during the campaign, with the tentative name of the corresponding compound. As the mass resolution of the PTR-QMS is unit mass, we cannot rule out that interfering compounds were detected at the same m/z. However, the very good correlations found between MeSH concentrations and DMS concentrations, and between sea–air MeSH fluxes and DMS fluxes (Fig. 5), give good confidence to the mass attribution to the target compounds. A Gas Calibration Unit (Ionicon) was used for daily background instrument determination (catalyzer heated at 350 °C) and calibration. Three calibrations with an NPL (National Physics Laboratory, Teddington, UK) standard (including DMS) were made on-board and showed good stability of the instrument. Average instrument sensitivity for the different compounds is presented in Table S2. The Minimum Limit of Detection (MLD) was determined from the standard deviation of blank measurements, following the method described in Rocco et al. (2025b).

Sea–air fluxes in the ASITs were calculated based on the difference in concentration between the flushing air and the ASIT outflow, using the following equation: 1FASIT-x=Q×CASIT-CflushA where FASIT-x (molm-2s-1) is the net flux of the compound from the water to the air, Q (m3s-1) is the air volume flow rate in the ASITs, CASIT (molm-3) is the concentration of the compound in the air leaving the ASIT, Cflush (molm-3) is the concentration in the flush air entering the ASIT. A (m2) is the water surface area in the ASIT. Using comparisons between measured fluxes and their prediction from empirical relationships using measured seawater dissolved DMS concentrations, SST and windspeed, Rocco et al. (2025b) calculated that the ASIT system generated fluxes equivalent to open ocean 2 ms-1 wind speed fluxes.

Particle concentrations were measured using several Condensation Particle Counters (CPCs), including two from TSI (Models 3010 and 3772) and a water CPC (MAGIC 200/210, Aerosol Devices). A Particle Size Magnifier (PSM, Airmodus A11) coupled to a CPC was also used in stepping mode. The particle counters provided concentrations with four cutoff diameters: 1 and 3 nm using the PSM, 5 nm using the MAGIC CPC, and 10 nm using the TSI CPCs. Loss corrections in the tubing were calculated according to the equation from Willeke and Baron (1993). The correction for the concentration of 1–3 nm particles (N1–3) was made by considering an average diameter of 2 nm, and for 1–10 nm particles (N1–10), by taking an average diameter of 3 nm. According to our analyses of the particle size distribution, this correction results in a slight overestimation of losses, which are generally smaller than estimated here. The formation rate of new particles was calculated with the following equation: 2J=dCASITdt-Q×CflushV where J (particlesm-3s-1) is the rate of formation of new particles, CASIT (particlesm-3) is the concentration of particles in the air exiting the ASIT and Cflush (particlesm-3) is the concentration of particles in the flushair entering the ASIT. Yet, the contribution of the second term, which corresponds to the contribution of particles entering the ASITs, is minimal due to the use of a particles filter on the flushing line. In the calculations of formation rates in the following sections, the concentrations taken will be those of particles greater than 1 nm (N1) corrected for diffusion losses in the lines. t is time in seconds, Q (m3s-1) is the inlet flow into the ASITs in and V (m3) is the volume of the headspace in the ASITs.

The flushing air line for the ASITs was sourced from the bow of the vessel to avoid contamination from engine and kitchen exhausts as the ship moved forward. Data from the air lines were filtered for two hours after contamination events, which were identified based on elevated levels of ozone, benzene, and xylene, as well as changes in particle size distribution (N10 > 0.5 × N1) and variability in particle concentration during a 20 min measurement period.

Seawater concentrations of DMS in the ASIT (DMSw) are calculated from headspace concentrations by using an average departure to the thermodynamical equilibrium, and therefore are more uncertain than direct measurements, as shown by offsets of the calculated values against the measured ones in the underway (Fig. S7). This departure to the thermodynamical equilibrium is based on the excellent correlation between the theoretical seawater concentrations calculated from the equilibrium with the headspace concentrations in the run-through mode of the control ASIT, and the measured concentrations in the underway seawater by Wohl et al. (2026). The slope of the measured versus calculated seawater concentrations was found to be 2.6 for DMS and 7.2 for MeSH. The slope found in the present study for DMS is very similar to the one used in Rocco et al. (2025b) that used the same experimental procedure (MeSH_w concentration was not measured in this previous study). Therefore, DMSw was calculated using: 3DMSw=DMSair×fsal×Hs(T)cc×2.6 where fsal is a dimensionless factor accounting for the effect of salinity on solubility (fsal = 0.848 for DMS, determined experimentally by Dacey et al., 1984) and 2.6 is the deviation from Henry's law equilibrium. Hs(T)cc is the dimensionless Henry's law solubility expressed as concentration in water over concentration in air and is a function of temperature (T) (Sander, 2023).

Chlorophyll a. Between 200 and 750 mL of seawater were filtered through 25 mm Whatman GF/F glass fibre filters. All filters were stored at -20 °C until analyses on board the R/V Hesperides. Chlorophyll a concentrations were estimated fluorometrically after extraction in 90 % acetone at 4 °C for 24 h (Yentsch and Menzel, 1963). Readings were conducted on a Turner 10AU fluorimeter calibrated with pure chlorophyll extract from spinach (Sigma C5357) using a Beckton-Dickinson spectrophotometer.

Statistical tests were performed using Python (version 3.11) with the SciPy library (version 1.16.2; Virtanen et al., 2020) at a significance level of α = 0.05. For paired measurements, normality of paired differences was first assessed with the Shapiro–Wilk test; when normality was not rejected, the paired t-test was applied, otherwise the non-parametric Wilcoxon signed-rank test was used, with effect size quantified by the rank-biserial correlation coefficient r. For independent, non-concurrent measurements, normality was assessed for each group separately, the Welch's t-test was applied when normality held for both groups and otherwise the Mann–Whitney U test was used, with effect size similarly quantified by the rank-biserial correlation coefficient r. When multiple comparisons were performed on the same variable, p-values were adjusted using the Bonferroni correction to control the family-wise error rate. For the correlation matrices, p-values were instead adjusted using the Benjamini-Hochberg false discovery rate procedure (Benjamini and Hochberg, 1995).

Out of the 22 compounds monitored by the PTR-QMS, only two, DMS and MeSH, were found to be enriched in the ASITs headspace with respect to the entering flushing air. No enrichment of benzene, toluene, ethylbenzene and xylene (BTEX), monoterpenes, isoprene, or iodomethane were detected in the ASITs headspace, as the same molar fractions were observed at the inlet (i.e., in the flushing line) and the outlet (typically 0.13, 0.057, 0.05, and 0.1 ppb, respectively), differing from previous observations using this set-up with different seawaters (Rocco et al., 2021, 2025a). The most plausible explanation for the lack of emissions would be a subsaturation of these compounds in the water phase. For example, isoprene mean ambient atmospheric concentrations were around 50 ppt (Williams, 2026) and mean oceanic concentrations (UW) were 0.07 nM according to the measurement made during the campaign that are presented in Crabeck et al. (2026). They also present concentrations reaching 1.5 nM in top sea-ice, this could therefore be one of the sources of isoprene in the atmosphere that was identified as missing in the study by Ferracci et al. (2024). However, we cannot rule out the possibility of losses or heterogeneous reactivity in the ASIT's headspace that might mask emissions of BTEX or isoprene, but that would mean that emissions of these compounds are particularly low given the measurements taken during a previous campaign using the same set-up (Rocco et al., 2021, 2025b).

Ozone mole fractions in the ambient air varied between 5 and over 20 ppb (Fig. 3) throughout the cruise with an average around 15 ppb. In contrast, ozone concentrations in the flushing line and the ASIT outlets were significantly lower, averaging around 7 ppb, primarily due to losses along the sampling line from the bow of the vessel (see Materials and methods). Notably, ozone levels in the ASITs were consistently lower than those in the flushing air, indicating further chemical reactions within the ASITs headspace or aqueous phase, that were estimated to contribute equally to ozone losses in a previous study (Salignat et al., 2026).

ASIT DMS and MeSH mole fractions and sea–air fluxes exhibited substantial variability over the cruise (Fig. 3). For the Control-ASIT (effect of UV filtering will be discussed Sect. 3.3) (Fig. 3E, red markers), the FASIT-DMS was on average 3.03 pmolm-2s-1 (i.e., 0.26 µmolm-2d-1), peaking at 8.68 pmolm-2s-1 on 10 March, while the average FASIT-MeSH was 0.64 pmolm-2s-1 (i.e., 0.06 µmolm-2d-1), with a maximum of 2.65 pmolm-2s-1 coinciding with that of FASIT-DMS. Note that fluxes of both compounds were always positive, i.e., these species were always degassing from seawater. Corresponding average seawater concentrations (Fig. S2) estimated for DMS (DMSw) was 1.8 nM, reaching a maximum of 4.5 nM, while estimated MeSHw averaged 0.8 nM, with a peak of 3.2 nM in the Control-ASIT (Fig. S2).

The spatial variability of gaseous DMS molar fractions, estimated aqueous concentrations and sea–air fluxes in the Control-ASIT are shown in Fig. 4. As expected, aqueous DMS concentrations showed the same spatial variability as the gaseous DMS concentrations, with maxima measured in the open Southern Ocean north-east of the Antarctic Peninsula (2.5–3 nM) and minima in the Marginal Ice Zone (MIZ) (1 nM) as well as in the Gerlache Strait (< 1 nM) (location indicated in Fig. 1), which was characterized by lower salinity, reflecting the inflow of fresh water from melting ice. This low salinity trend in the Gerlache Strait was confirmed twice during the ship's outward and return passages through this area. Aqueous DMS concentrations along the west coast of the Antarctic Peninsula (2.5–3 nM) were moderate (around 1.5–2 nM). Ambient air DMS molar fractions (Fig. 4A) were also higher in one open-ocean transect north of the Weddell Sea (0.4–1 ppb) than over most of the rest of the cruise track (< 0.3 ppb). Ambient air concentrations are a function of local fluxes from the ocean, but are also contributed by emissions upwind. Thus, matching observed atmospheric concentrations and emission fluxes is not straightforward. In addition to seawater concentrations, sea–air fluxes in the ASITs are, by definition, a function of incoming ambient air DMS concentrations; therefore, they showed a pattern which is different from that of estimated aqueous DMS concentrations. Fluxes can be low despite high seawater concentrations when incoming gaseous concentrations in the flushing air are high, as occurred over a portion of the transect north of the Weddell Sea. Conversely, fluxes were particularly high in coastal waters to the west of the peninsula, where DMS concentrations in the ambient air were particularly low. Air concentrations, estimated aqueous concentrations and sea–air fluxes of MeSH in the ASIT showed patterns parallel to those of DMS. The ratio of MeSH flux to the sum of sulfur compound fluxes varied from 0.1 to 0.2 for most of the data (Figs. 4F and 6A). Causes of this variability will be further addressed in the following sections.

A strong correlation was observed between sea–air MeSH and DMS fluxes in both ASITs, with R2 values of 0.77 for the Control-ASIT and 0.66 for the noUV-ASIT, across all coastal and open ocean measurement periods (Fig. 5A and B). The median value of the FASIT-MeSH / (FASIT-DMS+FASIT-MeSH) ratio obtained for all data combined was 0.17 (first quartile Q1 = 0.13; third quartile Q3 = 0.22) in the Control-ASIT (Fig. 6A). There do not appear to be major differences between the enclosed experiment periods and the periods with continuous water flow into the tanks, except that the ratio was less variable during the experiments, reflecting the more stable water composition in the ASITs. However, the average FASIT-MeSH/(FASIT-DMS+FASIT-MeSH) in Experiment B (with open ocean waters) was significantly higher than in Experiment A (with coastal water types) in Control-ASIT (Mann–Whitney U = 192, nexpA = 29, nexpB = 33, p = 5.46 × 10⁻⁵, r = 0.599) and in noUV-ASIT (Mann–Whitney U = 47, nexpA = 39, nexpB = 32, p = 2.70 × 10⁻¹¹, r = 0.925).

DMS and MeSH molar fractions both have a diel variation, stronger for MeSH than for DMS (Fig. 7). The diel variation of the DMS mole fraction appear clearly when the fractions are normalized to the daily maxima molar fraction (Fig. S5A) with a decrease between 12:00 and 19:00 LT and an increase of the molar fractions in the hours following the decline of radiation level, with significantly stronger diel contrast in the Control-ASIT (paired t-test, n = 193, t = -4.95, p = 1.6 × 10⁻⁶), most notably around midnight (00:00 LT: mean difference = 0.13 ± 0.12, r = 0.94, p = 0.016). The diel variation of MeSH molar fraction was even clearer with high fraction at night and decrease through the day (Figs. 7B and S5B). This pattern was marked in the noUV-ASIT both during day and night in molar fraction but is clearer for Control-ASIT when the fractions are normalized with significant differences (Wilcoxon signed-rank test, n = 193, W = 6118, p = 6.7 × 10⁻⁵, r = 0.35), particularly around 14:00 and 16:00 LT (r = 0.944 and 1.00, p < 0.05) (Fig. S5B).

The concentrations of 1–3 nm (N1--3), 1–10 nm (N1–10), > 1 nm (N1) and > 10 nm (N10) particles are presented in Fig. 8. The formation rates of 1 nm particles (J1) are obtained using the method presented in Sect. 2. In the Control-ASIT, very few and weak new particle formation events were observed (Fig. 8) with N1 particle concentrations reaching up to 38 particlescm-3, mainly composed of 1 to 3 nm particles, and formation rates peaked at 0.02 particlescm-3s-1. Few of these particles reached higher size than 10 nm which is consistent with the short residence time (40 min) in the ASITs. Most of the particles in the Control-ASIT were formed with coastal waters and new particles formation occurs only in the middle of the day, around maximum radiations periods (Fig. 9). The effect of filtering UV light on new particle formation will be addressed in the next section.

DMS and MeSH positive sea–air fluxes indicate emission from seawater linked to marine biological processes (Fig. 3). The equivalent wind speed inside the ASIT headspace was estimated as 2–3 ms-1 (Rocco et al., 2025b), i.e., significantly lower than typical conditions at the ocean surface. Consequently, a comparison to ambient fluxes in the literature is not useful but a comparison to similar experiments performed in the Southwestern Pacific Ocean during the Sea2Cloud voyage can be done (Sellegri et al., 2023). In Sea2Cloud, the average DMS fluxes in the Control-ASIT were 23.2, 8.2, and 2.9 pmolm-2s-1 (i.e. 2.0, 0.7 and 0.3 µmolm-2d-1) for frontal, subantarctic, and subtropical waters, respectively. The fluxes observed during POLAR-CHANGE (3.03 pmolm-2s-1, peaking at 8.68 pmolm-2s-1) were similar to Sea2Cloud fluxes in subantarctic and subtropical waters. Consequently, calculated aqueous DMS concentrations (0.3–4 and 1–4.3 nM for Control and noUV ASITs, respectively) were similar to those measured during the Sea2Cloud in subantarctic seawaters (3–5 nM; Rocco et al., 2025b).

During Experiments A and B, the estimated lower concentrations of DMSw in the Control-ASIT (Fig. S2) are likely due to an aqueous-phase reactivity (loss) of DMS, which is accumulating over several days in static flow mode (Salignat et al., 2026) as opposed to the continuous flow mode. Indeed, the gas-phase photochemical loss of DMS was shown to be negligible within the 40 min residence time of our experimental set-up (Rocco et al., 2025b; Salignat et al., 2026). Salignat et al. (2026) also report that the diel profile observed for DMS in the ASIT incubated seawater during Sea2Cloud, with maxima during the night and minima during the day was explained by aqueous-phase ozonolysis when a constant (day and night) production was prescribed. Earlier, Brimblecombe and Shooter (1986) demonstrated the photooxidation of DMS in seawater under UV irradiation. Various studies have reported that photochemical degradation can account for 7 % to 40 % of dissolved DMS loss in the ocean (Galí et al., 2023), depending on regional conditions, and have linked these losses to oxidants produced by chromophoric dissolved organic matter (CDOM) and NO3- photolysis when NO3- concentrations are higher than 1 µM, as is the case in Antarctic waters (Bouillon and Miller, 2005). The hypothesis of aqueous-phase oxidation also provides a plausible explanation for the persistence of emission differences at night, when UV radiation is absent. Although the average residence time of the water in the ASITs was approximately six hours, the residence time in a continuously stirred reactor is a distribution of residence times with some parcels being expelled from the reactor before the six hours and some parcels having much longer residence time. This means that daytime chemical transformations in the water produce compounds that persist during nocturnal periods. Additionally, ozone-mediated oxidation of DMS in seawater, as discussed by Hoffmann et al. (2016) and Salignat et al. (2026), could contribute to the observed differences between the ASITs. Since ozone concentrations were generally higher in the Control-ASIT for most of the campaign, and the ozone differential between ASITs showed a weak anti-correlation with the differences in DMS and MeSH fluxes, this mechanism may play a role. Although aqueous oxidation of MeSH could also be enhanced under these conditions, this process remains poorly characterized in the literature. A different DMS production induced by UV stress such as DMSP production inhibition could also explain part of the differences between the ASITs.

A strong positive correlation was found in DMS and MeSH fluxes (Fig. 5), as also previously reported earlier by Novak et al. (2022) over non-enclosed coastal waters (R2 = 0.65) and for molar ratio in previous ASIT measurements by Rocco et al. (2025b) (R2 = 0.81 in the Control-ASIT). The literature, mainly recent, reports ratios of MeSH to total volatile methylated sulfur pools and fluxes close to those reported here, generally between 0.10 and 0.40. Gros et al. (2023) reported a ratio of MeSH/(DMS+MeSH) concentrations in water averaging 0.20 in the North Atlantic but varying greatly with latitude and regularly reaching 0.40 above 70° N. Leck and Rodhe (1991) estimated the average ratio of MeSH/ (sulfur species) in the Baltic and North Seas at around 0.07, based on concentrations in water and including CS2 and dimethyl disulfide (DMDS) among the sulfur species. Lawson et al. (2020) calculated flux ratios of 0.14 to 0.24 in the Southwest Pacific, neglecting oxidative losses at night. A ratio value of 0.16 in dissolved concentration and 0.17 in fluxes was found by Kettle et al. (2001) in the South and North Atlantic. Finally, during the Sea2Cloud campaign, the FASIT-MeSH/(FASIT-DMS+FASIT-MeSH) ratio ranged between 0.11 and 0.18 in the Control-ASIT (Rocco et al., 2025b). Thus, the values obtained in the present study are consistent with previous measurements, with values mainly between 0.10 and 0.20, but reaching 0.30 to 0.45 mostly in open ocean waters. Air concentration ratios within the Control-ASIT headspace (median of 0.11; Q1 = 0.8; Q3 = 0.14) were in the range of those reported by Mynard et al. (2025), between their Antarctic Ice-Edge (median of 0.07) and the subtropical front region (0.11–0.15).

The diel pattern of MeSH concentrations showed a decrease in the afternoon, in both ASITs, that may mainly be resulting from the higher gas-phase OH-reactivity (or photochemical loss) of MeSH compared to DMS, within the 40 min residence time in the ASIT's headspace. Greater chemical loss for MeSH than for DMS would then be behind the observed diel variation of the air MeSH/(MeSH+DMS) concentration ratio, as also suggested by Mynard et al. (2025), and could explain the significant differences in MeSH concentration between the two ASITs and a lower MeSH/(MeSH+DMS) ratio during the day compared to night, even outside the experimental periods. These observations closely resemble the diel pattern reported by Novak et al. (2022). Lawson et al. (2020) reported broadly similar diel profiles for DMS and MeSH, with a production minimum around 4 p.m. (local time) for both compounds, a trend more evident for MeSH in our experiments. However, a closer look at our measurements reveals that the MeSH concentration differs from one ASIT to the other most at night, which points to the oxidation of MeSH in the Control-ASIT being from another oxidant than OH (i.e., a nighttime aqueous-phase oxidant at larger concentrations in the presence of daytime UV), or to an impact of UV on the biology processes leading to a change in the production of MeSH in the seawater at night.

The difference in ozone concentration between ASITs may be attributed to several factors, including chemical reactions with other gaseous species, interfacial reactions during the air residence time (∼ 40 min), or aqueous-phase reactions occurring over the 6 h water residence time. These ozone losses, which occasionally reached 100 % between the flushing air and ASIT headspace, were variable over time, suggesting that chemical reactivity, rather than physical deposition alone, played a significant role. If physical losses to surfaces were dominant, they would likely remain relatively constant or vary in relation to temperature. The noUV-ASIT lid blocked nearly all radiation below 375 nm (Fig. S1), preventing most photochemical processes. In the presence of NO, the photolysis of NO2 and NO3 (which requires wavelengths < 410 nm; Bower et al., 1993; Andersen et al., 2022) would be inhibited, potentially leading to a reduction in ozone production while ozone consumption continues unabated. One other potential explanation for the increased ozone losses in the noUV-ASIT is a higher ozone consumption, that could be due to enhanced water-to-air net fluxes of reactive compounds, particularly DMS and MeSH, whose fluxes were indeed higher in the noUV-ASIT during both experiments (see Fig. 3). It could also be the ozonolysis of organic matter, accumulated over time especially during experiments, potentially leading to more VOC emissions (Kilgour et al., 2024). The increased difference between ASITs in Experiments A and B suggests that UV removal had a significant impact on the aqueous-phase chemistry (as discussed earlier) that in turn impacts ozone losses to the aqueous phase during experiments. The temporal pattern of these changes also provides insight into the relevant timescales: the differences between the two ASITs emerged within the first 12 h of the experiments, reaching 37 % higher ozone in the Control-ASIT compared to the noUV-ASIT in Experiment A and 38 % in Experiment B at that time point. This indicates that reactions in the aqueous phase and associated emission processes occur over timescales longer than six hours, but less than 12 h in this experimental setup.

We explored the potential links among the different biological elements and chemical compounds measured in the ASIT experiments by considering the temporal patterns and the correlations between the different parameters (Figs. 12 and S8). Note that significant correlations should not always be interpreted as a direct cause-effect relationship, nor a short-term response occurring in the ASITs. Furthermore, the correlations should be taken individually, as indicators of the diverse, dynamic and complex biological processes undergoing in the ASIT experiments, somewhat apparently contradictory in some cases. In Experiment A, which involved incubating coastal water in the ASITs for 2 d, DMS and MeSH seawater concentrations slightly decreased over time, particularly in the Control-ASIT, where MeSH concentrations became almost zero, before rising up again (Fig. 3). In Experiment B, which was similar to Experiment A but with a different phytoplankton composition from open ocean waters, concentrations were higher, and the difference between the two ASITs was greater than in Experiment A, with MeSH concentrations initially 50 % higher in the noUV-ASIT than in the Control-ASIT. As Experiment B progressed, the differences in DMS and MeSH concentrations between the two ASITs diminished, until a point where concentrations were even higher in the Control-ASIT.

The diel profile of Fv/Fm, typically observed in the surface ocean, where diurnal exposure to high light levels (including UV) leads to some degree of stress on the photosynthetic machinery, showed differences between the Control-ASIT and the noUV-ASIT and therefore could partially explain differences in production of MeSH during the experiments. Notably, during enclosed-water experiments, a diurnal increase in HNA bacterial abundance coincided with a severe decrease in phytoplankton's photosynthetic efficiency (Fv/Fm), and both were more acute in the presence of UV. It can be speculated that solar radiation stress in phytoplankton cells results in the release of organic compounds, including DMSP, that would have fostered bacterial growth as reflected in the increase of (active) HNA bacteria. In turns, this could have led to higher gross MeSH production by bacterial DMSP catabolism. The fact that air MeSH concentrations in the ASITs were higher during the night and higher in the noUV-ASIT points to aqueous, solar radiation dependent losses overcoming potential increases in gross MeSH production.

Furthermore, a significant positive relationship between estimated aqueous concentrations of MeSH and DMS and HNA bacterial abundances was also observed in both ASITs (Spearman ρ = 0.43 and 0.49 in Control-ASIT; ρ = 0.57 and 0.63 for noUV-ASIT, respectively for MeSH and DMS). This relationship was particularly noticeable in open ocean waters (Fig. 12). A relationship of sulfur compounds with V1 and V3 virus abundance was also observed in both ASITs during experiments but not with V2 virus (Fig. S8). As explained above, dynamic and diverse processes occurred in the ASITs experiments which differently affected the particular biogeochemical elements at play, in this case targeting the virus groups.

The positive relationship between bacteria and sulfur compounds had been previously observed during the Sea2Cloud campaign (Rocco et al., 2025b). This relationship has been documented in several studies and is explained by the ability of bacteria to cleave and demethylate-demethiolate DMSP, giving rise to DMS and MeSH, respectively (Kiene, 1996; Kilgour et al., 2022). The proportion of each compound in the production, as well as in the concentrations, depends on several factors, including (1) DMSP release by exudation or mortality of the algal producers, (2) whether or not the main algal DMSP producers harbor DMSP lyases and are therefore capable of releasing DMS as well, (3) the relative proportion of bacterial genes (and the corresponding enzymes) for cleavage and demethylation (which mostly depends on the taxonomic composition of the bacterial assemblage), and (4) the relative activity of these enzymes (which may depend on bacterial sulfur and carbon demands relative to DMSP availability) (Simó 2001; Sun et al., 2016; Hopkins et al., 2023).

None of the biological variables available showed a strong correlation with aqueous DMS and MeSH concentrations throughout the entire campaign when considering the continuous flow mode of our ASITs set-up, except with Chl a and bacterial abundances during the open ocean period. The non-systematic measurement of diatoms and dinoflagellates, known emitters of DMSP (Wittek et al., 2020; Sheehan and Petrou, 2020; Shaw et al., 2022), could explain this lack of correlations. During the experiments, Chl a was inversely correlated with both DMSw (Spearman ρ = -0.81) and MeSHw (ρ = -0.67), as were 1–2 and 15–20 µm SED nanophytoplankton size classes (Fig. S8A). This contrasts with the positive associations between Chl a and DMSP, the main DMS precursor, reported in Antarctic coastal and glacially influenced waters (Li et al., 2024; Petrou et al., 2026), and points to a decoupling between phytoplankton biomass and dissolved volatile sulfur compounds, whose concentrations depend on the balance between DMSP release, cleavage and subsequent removal as discussed earlier. Given the absence of DMSP measurements in the ASITs, this interpretation remains hypothetical but it must be noted that significant positive correlation was observed between DMSP and nanophytoplankton (particularly cryptophytes) in sea surface samples of the same cruise (Wohl et al., 2026).

The maximum particle concentrations observed during this campaign were slightly higher than those observed during the previous campaign using the same ASIT setup (Chamba et al., 2023), where N1 particles reached concentrations of 180 particlescm-3, and nucleation rates were up to 0.08 particlescm-3s-1. The present nucleation rates in the noUV-ASIT, reaching 0.15 particlescm-3s-1, are slightly above those reported for Antarctic areas (0.05 particlescm-3s-1) and below those reported for Arctic areas in the ambient open atmosphere (0.51 particlescm-3s-1) (Kerminen et al., 2018).

Air MeSH concentrations showed the strongest positive correlation to N1 and N10 particle concentrations in the noUV-ASIT in continuous flow mode with coastal waters (R2 = 0.56 for N10; yellow in Fig. 13). This correlation degraded in the case of DMS (R2 = 0.31). No significant relationship of MeSH and DMS to particle concentrations was observed in the noUV-ASIT with open ocean waters or during Experiments A and B, nor in the Control-ASIT in any observational period when particle concentrations remained low.

It is important to note that, like during the Sea2Cloud campaign (Chamba et al., 2023), no biological variable or gaseous compound showed differences between the ASITs as pronounced as the concentrations of particles, which were up to a factor of 10 higher in the noUV-ASIT. As concluded from the Sea2Cloud study, the presence of DMS and MeSH alone did not explain the formation of new particles in the enclosures. However, a significant positive correlation was found between N1 and the 7–15 µm nanophytoplankton cell abundance during experiments (Fig. S8). This calls for further studies identifying the new particle precursors related to nanophytoplankton.