Nitric oxide (NO) in the Bohai Sea and the Yellow Sea

. Nitric oxide (NO) is a short-lived compound of the marine nitrogen cycle; however, our knowledge about its oceanic distribution and turnover is rudimentary. Here we present the measurements of dissolved NO in the surface and bottom layers at 75 stations in the Bohai Sea (BS) and the Yellow Sea (YS) in June 2011. Moreover, NO photoproduction rates were determined at 27 stations in both seas. The NO concentrations in the surface and bottom layers were highly variable and ranged from below the limit of detection (i.e., 32 pmol L − 1 ) to 616 pmol L − 1 in the surface layer and 482 pmol L − 1 in the bottom layer.


Introduction
Nitric oxide (NO) is a short-lived intermediate of the oceanic nitrogen cycle (Bange, 2008). It received limited attention so far because its determination in seawater is challenging (Zafiriou and McFarland, 1980;Lutterbeck and Bange, 2015;Liu et al., 2017). NO in surface seawater can be produced via the photolysis of nitrite (NO − 2 ) (Zafiriou and McFarland, 1981;Olasehinde et al., 2009Olasehinde et al., , 2010Liu et al., 2017): This reaction may account for 10 % of nitrite loss in surface waters of the central equatorial Pacific (Zafiriou et al., 1980). Zafiriou and McFarland (1981) suggested that photochemically produced NO is a potential source of atmospheric NO during daylight, which could further lead to ozone hole, acid precipitation, and photochemical smog. Apart from the photochemical production, various microbial pathways of NO have been identified, including denitrification, nitrification, and anammox (Schreiber et al., 2014;Martens-Habbena et al., 2015;Caranto and Lancaster, 2017;Kuypers et al., 2018). Additionally, NO is a messenger molecule in marine organisms: phytoplankton not only responds to exogenous NO  but also produces NO dur-ing their growth (Zhang et al., 2006a, b;Kim et al., 2006Kim et al., , 2008. Chen et al. (2015) reported that calmodulin (a messenger protein expressed in eukaryotic cells) of the tropical sea cucumber participates in the production of NO during immune response. Morrall et al. (1998Morrall et al. ( , 2000 characterized the NO synthase activity in the tropical sea anemone Aiptasia pallida, suggesting that NO and NO synthase can act as ecotoxicological biomarkers in the tropical marine environment. Moreover, the characterization of NO synthase gene in the intertidal copepod Tigriopus japonicas has been found, and the intracellular production of NO in shrimp haemocytes has been observed (Xian et al., 2013;Jeong et al., 2016). Thus, NO seems to be widespread, with different functions in marine organisms.
The current understanding of the oceanic NO distribution is mainly limited to the ocean surface (Zafiriou and McFarland, 1981;Olasehinde at al., 2009Olasehinde at al., , 2010Liu et al., 2017) and oxygen minimum zones (OMZs; Ward and Zafiriou, 1988;Lutterbeck et al., 2018). Only recently have the distribution of NO as well as its seasonal variation in the Jiaozhou Bay and adjacent waters been studied (Feng et al., 2011;Xue et al., 2012;Tian et al., 2016).
In this study, we present the first measurements of dissolved NO in the Yellow Sea (YS) and the Bohai Sea (BS). The overarching objective of our study was to decipher the biogeochemical fluxes of NO in the BS and YS. The specific objectives were (i) to determine the spatial variation in dissolved NO concentrations in the water column, (ii) to determine NO photoproduction rates, and (iii) to estimate the sea-to-air gas exchange fluxes of NO.

Study area
The BS and the YS are marginal seas of the western Pacific Ocean. The BS is a shallow and almost enclosed sea with a surface area of 77×10 3 km 2 and a volume of 1.39×10 3 km 3 . Its maximum depth is 83 m, with an average depth of only 18 m. The Yellow River is a major source of freshwater to the BS whereas the YS is the source of salt for the BS through water exchange via the Bohai Strait. The YS has a surface area of 380 × 10 3 km 2 and total volume of 16.7 × 10 3 km 3 . Maximum depth is 140 m with an average depth of 44 m, and like the BS, it is a continental shelf sea. The Yangtze River at the southwest corner of the YS is the major source of freshwater for the southern Yellow Sea (SYS) and the East China Sea. The hydrographic properties of this region are mainly influenced by the Yellow Sea Cold Water Mass (YSCWM) on the shelf (Lü et al., 2010;Li et al., 2016) and the Yellow Sea Coastal Currents on the western side of the basin (Su, 1998;Lee et al., 2002;Zhang et al., 2004). The BS and the YS are regions surrounded by areas of high popula- tion growth and economic development in China and Korea (Zhan et al., 2010;Jiang et al., 2014).

Sampling
Samples were collected from 13 to 28 June 2011 on board the R/V Dong Fang Hong 2 in the BS and the YS. Sampling for NO from both the surface (at 1 m) and the bottom layer (= 1 m above the ocean bottom) were performed at 75 sampling stations (including one 24 h anchor station: B65) shown in Fig. 1 and listed in Table 1 (Tian, 2019). Water samples were collected using 8 L Niskin bottles equipped with silicon O-rings and Teflon-coated springs and mounted on a Sea-Bird CTD (conductivity, temperature, depth) instrument (Sea-Bird Electronics, Inc., USA). A 500 mL Wheaton glass serum bottle was rinsed with in situ seawater three times and was then filled with seawater quickly through a siphon. When the overflowed sample reached the half volume of the bottle, the siphon was withdrawn rapidly, 0.5 mL saturated HgCl 2 (aq) solution was added, and the bottle was sealed quickly. The surface water samples were immediately analyzed after collection, and samples from other depths were temporarily placed in the dark in a water bath with a surface seawater circulation system and analyzed within 1 h (Liu et al., 2017). Samples for the analysis of NH + 4 -N, NO − 2 -N, and NO − 3 -N were immediately filtered through 0.45 µm Millipore membrane (pre-ignited at 450 • C for 6 h). The filtrates were stored in polyethylene bottles (pre-marinated with 1 : 10 HCl for 24 h) at −20 • C. Photoproduction rates were determined in surface water samples from 22 stations in the YS and 4 stations in the BS (shown in Fig. 1), respectively.

Analytical procedures
Analysis for dissolved NO was conducted with the improved method of Liu et al. (2017) by a combined purge-and-trap and fluorometric detection method. The precision of the analytical method was better than ±7 %, and the limit of detection (LOD) was 32 pmol L −1 in the field experiment (Liu et al., 2017).
The photolysis experiments were conducted under natural light conditions on deck. Surface water samples were filtered with 0.45 µm Millipore membrane and transferred into 10 mL cleaned quartz vials wrapped with aluminum foil with no headspace. Then 200 µL NaN 3 solutions were added to remove the microbial influence. After adding 10 µL of 1 × 10 −3 mol L −1 2, 3-diaminonaphthalene (DAN, a NO probe) solution and gentle mixing, the fluorescence of the mixed solution was measured before irradiation (Liu et al., 2017). Then the duplicate sample made without the addition of DAN was capped with a Teflon-lined silicone septum and without aluminum foil; following this, the vial was placed in a shallow circulating seawater bath. After 0.5 h irradiation time, 10 µL of 1 × 10 −3 mol L −1 DAN was added into the duplicate sample without addition of DAN before and the fluorescence of the solution was analyzed. The NO concentrations were measured with the method described in Liu et al. (2017). The NO apparent photolysis production rates were computed based on the time-dependent difference between the NO concentrations before and after irradiation. For dark controls, vials were wrapped in aluminum foil. The quartz vials and syringe used in the experiment were soaked in a 10 % (v/v) HCl bath for 24 h, rinsed with Milli-Q water, and baked at 500 • C for 4 h.
Chlorophyll a (Chl a) was fluorometrically measured by an F-4500 fluorescence spectrophotometer after filtration of 200 mL seawater through a Whatman glass fiber filter and extraction in 90 % acetone according to Strickland and Parsons (1968). The wind speeds were measured at a height of 10 m above the sea surface using a Model 27600-4X shipborne weather instrument (Young, USA). Illuminance was measured by a digital illuminometer (TES-1330A, Shenzhen, China). Dissolved oxygen (DO) was measured in discrete water samples by the Winkler method (Grasshoff et al., 2009). The concentrations of dissolved inorganic nitrogen (nitrate, nitrite, and ammonium) were analyzed using a nutrient automatic analyzer (Auto Analyzer 3, SEAL Analytical, USA) in the laboratory. The limits of the detection were 0.14 µmol L −1 for nitrate, nitrite, and ammonium, with the precision of the method found to be better than 3 % .

Calculation of NO flux and NO saturation
Fluxes of NO across the sea-to-air interface were estimated following the approach of McGillis et al. (2000) for a sparingly soluble gas that is also moderately reactive in the atmo-sphere: where F is the flux density (mass area −1 time −1 ) across the air-sea interface, k sea is the gas transfer velocity (length per time), and c sea is the measured concentration of NO in the surface seawater (mol per volume). The partial pressure of the atmospheric NO (pNO air ) was calculated as follows: where x NO air is the mole fraction of atmosphere NO (dimensionless). We used the value of 2.13 ppb for x NO air , which is the average atmospheric NO mole fraction over the YS (Hu Min, Peking University, personal communication, 2018). p ss is the barometric pressure at sea surface, which was set to 101.325 kPa as the average pressure, and p w is the water vapor pressure at sea surface, which was calculated following Weiss and Price (1980) H cp is the Henry's law constant, which is calculated as follows:  (2015). k sea was calculated as follows: where k w is the water side air-sea gas transfer coefficient for sparingly soluble gases (length per time) calculated according to Wanninkhof (2014), γ a is the fraction of the entire gas concentration gradient across the air side boundary layer, as a fraction of the entire gradient from the bulk water to the bulk air (dimensionless) (McGillis et al., 2000), k a is the air side air-sea gas transfer coefficient (length per time) according to McGillis et al. (2000), H cc is the Henry coefficient (dimensionless) (Sander, 2015), M NO and M H 2 O are the relative molecular mass of NO and H 2 O (dimensionless), and u is the wind speed at 10 m height under neutral boundary conditions (length per time). The Schmidt number (Sc) is the kinematic viscosity of water divided by the molecular diffusion coefficient of the gas in (sea)water (Jähne et al., 1987;Wanninkhof, 2014). Seawater 4490 Y. Tian et al.: Nitric oxide (NO) in the Bohai Sea and the Yellow Sea dynamic viscosity (µ sw ) is a function of temperature (T ) and salinity (S) and was estimated using the following equations (Sharqawy et al., 2010): A = 1.541 + 1.998 × 10 −2 T − 9.52 × 10 −5 T 2 , Seawater density was estimated using Millero's empirical equation (Millero et al., 1976), and NO diffusion coefficient D L (×10 −5 cm 2 s −1 ) in water was calculated according to Wise and Houghton (1968): The saturation factor (α) is defined as α = c sea /(pNO air × H cp ), α>1 represents that NO was supersaturated and that the flux was from sea to air.

NO in the surface and bottom layers
The NO concentrations from the surface and bottom layers of the BS and the YS, as well as the local sampling time, bottom depth (D), temperature (T ), salinity (S), Chl a, wind speed (u), and DO are listed in Table 1. In the study area, temperature varied from 8.8 to 21.7 • C and salinity varied from 29.57 to 32.73 ‰ in the surface water. DO fluctuated from 178 to 271 µmol L −1 . The average concentrations of NH + 4 -N, NO − 2 -N, and NO − 3 -N were 2.11, 0.20, and 2.59 µmol L −1 , respectively. While in the bottom layer, temperature varied from 3.8 to 21.6 • C, salinity varied from 30.00 to 33.71 ‰, and DO varied from 147 to 266 µmol L −1 .
The surface concentrations of NO ranged from below the LOD to 616 pmol L −1 , with an overall average of 186 ± 108 pmol L −1 , and exhibited a considerable spatial variability (Fig. 2a). The mean NO surface concentrations in the BS (203 ± 107 pmol L −1 ) and northern YS (NYS) (212 ± 130 pmol L −1 ) were higher than the mean NO concentration in the SYS (159 ± 84 pmol L −1 ). The NO concentrations in the bottom layer ranged from below the LOD to 482 pmol L −1 , with an overall average of 174 ± 123 pmol L −1 . The mean concentrations of NO in the bottom layers of the BS and NYS were 228 ± 116 and 210 ± 138, respectively, and were higher than the mean (127 ± 98 pmol L −1 ) for the SYS. The maximum NO surface and bottom concentrations were measured at stations B21 and B28 in the NYS, respectively (Fig. 2b), whereas the lowest NO surface and bottom concentrations were measured in the center of the SYS (Table 1). Overall, there were no statistically significant differences (p>0.05) between the mean NO concentrations in the surface and bottom layers.
NO surface concentrations did not show any statistically significant relationship (p>0.05) with sea surface temperature, salinity, DO, Chl a, and illuminance. Trends of NO concentrations with salinity and DIN were only found for the stations affected by the outflow of the Yellow River in the southern BS (Fig. 3), where we found an inverse relationship between salinity and NO surface concentrations. Here high NO concentrations were associated with high DIN concentrations (data not shown; see Liu et al., 2015;Yang et al., 2015) indicating that DIN, especially NO − 2 , was a prerequisite for enhanced NO concentrations.
The water columns of the BS and YS were welloxygenated during our study, and thus no suboxic or anoxic conditions were detected. Therefore, we did not find any enhancement of NO concentrations in the bottom layers with low DO or NO − 2 concentrations as observed in the OMZ of the eastern tropical South Pacific Ocean off Peru (Lutterbeck et al., 2018). However, NO concentrations in the bottom layer showed negative correlations with salinity (p<0.05, r = −0.272, n = 60) and bottom depth (p<0.05, r = −0.298, n = 61) indicating a decrease in NO concentrations from the coast toward offshore waters. It is known that both the BS and the YS are affected by submarine ground water discharge (Kim et al., 2005;Taniguchi et al., 2008;Liu et al., 2017a, b). Moreover, NO has been detected in groundwater and aquifers (Smith et al., 2004;Smith and Yoshinari, 2008). Therefore, we suggest that input of NO from submarine groundwater discharge contributed to the distribution of NO bottom concentrations as well.
An overview of published NO surface concentrations is given Table 3. The average NO surface and bottom concentrations from this study are comparable to the concentrations measured in the Jiaozhou Bay, in the waters off Qingdao, and in the Seto Inland Sea. However, our mean concentrations are considerably higher than the NO concentrations reported from the central equatorial Pacific and the eastern tropical North Pacific Ocean. The maximum concentration reported here is at the lower end of the NO concentrations recently reported from the anoxic OMZ off Peru. Overall, NO surface concentrations seem to be generally higher in coastal waters compared to those found in offshore waters.

Diurnal variability
The diurnal variability of surface NO concentrations, illuminance, DO, and Chl a were investigated at anchor station B65 (Fig. 4). NO concentrations varied from 64 to 424 pmol L −1 , exhibiting a significant diurnal variation with the maximum concentration 8-fold higher than the minimum concentration. The NO concentrations reached the maximum concentration in the early afternoon (about 13:00 LT) and then decreased to the minimum concentration at 22:00 LT. A less pronounced second maximum (209 pmol L −1 ) was reached at 04:00 LT. DO showed a similar diurnal cycle but shifted by 3 h with maxima at 16:00 LT and 07:00 LT. Chl a concentra-  tions peaked at 19:00 LT and 07:00 LT. The illuminance had its maximum at 13:00 LT, coinciding with the NO maximum indicating that the first NO maximum was indeed resulting from a photochemical production during daytime (Zafiriou and McFarland, 1981). However, the second maximum of NO at 04:00 LT, when it was still dark, must have resulted from an alternative chemical and/or biological production.

Photoproduction rates
The results of the NO photoproduction experiments are listed in Table 2. The photoproduction rates of NO in the BS, the NYS, and the SYS varied from not detectable to 5.07 × 10 −11 , 0.09 to 0.69 × 10 −11 , and 0.32 to 1.54 × 10 −11 mol L −1 s −1 , respectively. The average photo- production rate of the whole study area was 1.14 ± 1.37 × 10 −11 mol L −1 s −1 . The photoproduction rates from the BS and YS are in good agreement with the rates reported from the Seto Inland Sea (Olasehinde et al., 2009(Olasehinde et al., , 2010. However, the mean NO photoproduction rates are higher than those from the central equatorial Pacific Ocean (Zafiriou and McFarland, 1981) and lower than those from Kurose River, Japan (Olasehinde et al., 2009) (Table 2), which mirrors the available NO − 2 concentrations that are low in the open ocean and but high in a river.
The average photoproduction rate in the NYS was obviously lower than those in the BS and the SYS, consistent with illuminances of 22 450, 20 433, 27 852 lux and mean surface downward UV radiation flux 34.7, 32.1, 40.6 W m −2 (European Centre for Medium-Range Weather Forecasts (ECMWF), reanalysis data) for the BS, the NYS, and the SYS, respectively. Enhanced photoproduction rates (≥ 2.00 × 10 −11 mol L −1 s −1 ) occurred in the SYS, especially in the central part of the SYS and stations influenced by the Yellow Sea Cold Current (see Fig. 1). This is apparently in contrast to the distribution of NO, which showed the lowest concentrations in the central SYS (see above). The high illumination observed in the SYS (Table 2) not only leads to enhanced NO photoproduction (see Table 2) but also generates reactive oxygen species like O q − 2 , ROO q , and other OH-derived radicals, which in turn would efficiently scavenge NO (Olasehinde et al., 2010). Overall, the NO photoproduction rates showed a positive relationship between illuminance (p<0.01, r = 0.884, n = 26) and the mean surface downward UV radiation flux (p<0.01, r = 0.865, n = 26), indicating that the NO concentrations in the surface layer during daylight were dominated by photochemical production. However, we did not find a significant difference between the mean NO concentrations sampled during the day (179±80 pmol L −1 ) and at night (195±140 pmol L −1 ). This suggests that there was also a non-photochemical NO dark production in the surface layer, such as the nitrification process, or other chemical process, such as the process of ammonium (NH + 4 /NH 3 ) oxidation into NO − 2 and NO − 3 (Joussotdubien and Kadiri, 1970). Caranto and Lancaster (2017) found that NO is an obligate bacterial nitrification intermediate produced by hydroxylamine oxidoreductase; Ward and Zafiriou (1988) also found that NO might be an intermediate for soluble byproducts of nitrification such as hydroxylamine.
Besides, Olasehinde et al. (2010) found that filtered and unfiltered seawater samples collected from the Seto Inland Sea showed no significant difference in NO photoformation rates, which suggested a negligible contribution of NO produced by photobiological processes from particle matter in seawater. However, Liu et al. (2017) reported that the rate difference between filtered (0.45 µm, 1.46×10 −12 mol L −1 s −1 ) and unfiltered (1.52 × 10 −12 mol L −1 s −1 ) seawater samples from coastal waters of the YS indicated that particles in seawater could increase the NO production rate. The difference might be due to the composition of sample, filter membrane, etc. Thus, further research is needed.

Sea-to-air fluxes of NO
In the present study, the NO supersaturation was ubiquitous at all investigated sites. The supersaturation factors varied from 8 to 154, with an average of 47. Together with the fact that NO was supersaturated during daytime and nighttime at the 24 h station B65, which indicates that the BS and YS were a source of NO to the atmosphere. The sea-to-air flux densities ranged from 5.8 × 10 −19 to 3.6 × 10 −15 mol cm −2 s −1 , with an average value of 4.5 × 10 −16 mol cm −2 s −1 (Fig. 5). The comparison of wind speeds and flux densities reveal that the flux densities are mainly driven by the wind speed (Fig. 5) and that they showed significant positive relationship (p<0.01, r = 0.938, n = 20). Anifowose and Sakugawa (2017) also found that wind speed was an important factor governing NO flux at the air-sea interface. Our flux densities from the BS and YS were similar to those computed for the central equatorial Pacific Ocean and the Seto Inland Sea, while they are slightly lower than those computed for of Jiaozhou Bay waters (Table 2). Based on the YS area of 380 × 10 3 km 2 and the BS area of 77 × 10 3 km 2 , the emission of NO to the atmosphere was estimated to be 9.0 × 10 8 g N yr −1 or 7.6 × 10 7 g N month −1 . Ding et al. (2018) reported a satellite-derived average NO x emission estimate for June over the BS and the YS in the period from 2007 to 2016 of about 1.3 × 10 10 g N month −1 . The obvious very large discrepancy between the satellite-derived emission estimate and the one presented here results from the fact that Ding et al.'s (2018) estimate is dominated by the NO x emissions from ships' diesel engines. This indicates that oceanic NO emissions to the atmosphere only account for a negligible fraction (∼ 0.6 %) of the NO x emissions observed over the BS and YS.

NO mixed-layer budget
In order to estimate the contribution of different sources and sinks of NO in the surface layer of the BS and YS, we applied a simple box model. We assume that the surface layer is represented by the mixed layer, with a mean water depth of 15 m (Qiao et al., 2004). At steady state the loss of NO by air-sea exchange (F ase ) must be equal to the sum of the photoproduction rate (F pp ) in the mixed layer, the input from below into the mixed layer by diapycnal diffusion (F dia ) and other production or consumption pathways (F poc ), and advection into or out of the BS and YS mixed layer (F adv ). To this end, F ase is given by The mean F ase was 4.5×10 −16 mol cm −2 s −1 (see above). We assume that the NO photoproduction decreases linearly from the average rate of 1.14 ± 1.37 × 10 −11 mol L −1 s −1 at the surface (see above) to 0 at 15 m. The mean F pp was calculated to be 8.6×10 −12 mol cm −2 s −1 in mixed layer. Since the mean surface NO concentration in the surface layer was statistically not different from the mean NO concentration in the bottom layer (see above), it is reasonable to assume that F dia = 0. Moreover, it seems reasonable to assume that advection of NO into or out of the BS and YS surface layer is zero. Since F ase < F pp , we conclude that F poc should be negative, indicating that the bulk of the produced NO was rapidly consumed in the surface layer before its release to the atmosphere. Please note that we have some indications that there is also a dark production of NO (see above), meaning that the "true" NO surface production might be even higher. Chemical reactions with DO, OH, or ROO q etc. are potential sinks for NO in the surface layer of the BS and YS (Ford et al., 1993;Olasehinde et al., 2010;Carpenter and Nightingale, 2015).

Conclusions
This study reports the distribution and photoproduction rates of dissolved NO measured during a cruise in June 2011 in the Bohai and Yellow seas. The NO concentration distribution in both the surface and bottom layers were highly variable; however, there was no significant difference between the mean NO concentrations in the surface and bottom layers. NO concentrations in the bottom layer showed a significant decrease with salinity, which indicates NO input by submarine groundwater discharge. NO in the surface layer was supersaturated at all stations during day and night. The accumulation of NO during daytime was a result of photoproduction, and the measured NO photoproduction rates were correlated to illuminance. The persistent nighttime NO supersaturation pointed to a non-photochemical (so far unidentified chemical and/or biological) NO dark production. On the basis of a simple box model calculation, we conclude that the bulk of the NO produced in the surface layer was rapidly consumed before its release to the atmosphere. Overall, the BS and YS were a persistent source of NO to the atmosphere at the time of our measurements. However, the oceanic NO emissions were negligible compared to the NO x emissions from ships' diesel engines.