Comparison of Ocean Surface Wind Stress Computed with Different Parameterization Functions of the Drag Coefficient

Depending on the choice of reference wind speed, the quantitative and qualitative properties of the drag coefficient may vary. On the ocean surface, surface waves are the physical roughness at the air-sea interface, and they play an important role in controlling the air-sea exchange processes. The degree of dynamic influence of surface waves scales with wavelength. Drag coefficient computed with the reference wind speed at an elevation proportional to the wavelength (for example, U _λ/2) is fundamentally different from the drag coefficient computed with the wind speed at fixed 10 m elevation (U ₁₀). A comparison has been carried out to quantify the difference in wind stress computation using several different parameterization functions of the drag coefficient. The result indicates that the wind stress computed from U ₁₀ input using a drag coefficient referenced to U _λ/2 is more accurate than that computed with drag coefficient functions referenced to U ₁₀.     

Eddy Diffusion Coefficients of Suspended Particulate Matter: Effects of Wind Energy Transfer and Stratification

Gross sedimentation rates (GSR) have been measured using sediment traps placed at nine different levels above the bed (0·3, 0·5, 0·8, 1·0, 2·0, 4·0, 6·0, 8·0 and 10·0 m). The sediment traps were deployed for 1·25 years and recovered 28 times during the study period. Low average GSR values of 5·5 g m^-2 day^-1 were obtained at 10·0 m, and high average GSR values of 114·8 g m^-2 day^-1 were obtained at 0·3 m. An expression for the eddy diffusion coefficient of suspended particulate matter (K_s), based on the measured GSR is given. The expression has been used for modelling of K_s at the different trap levels above the bed. High values (≈42 cm² s^-1) of K_s were obtained at the upper traps, whereas low values (≈2 cm² s^-1) were obtained near the bed. Comparison between level of turbulent energy in terms of shear stress at the boundaries of the water column, i.e. from the wind and the bed flow, showed that wind energy exceeded that of the bed flow by a factor 16. At 5·0 m K_s was positively correlated (r=0·66) to the eddy diffusion coefficient of momentum (K_m) derived from the wind energy transfer to the water, giving an average β of 0·5 for K_s =βK_m. The density difference between surface and bottom waters has been designated a parameter of stratification, and is discussed in relation to variations of K_s and K_m .     

Short-term variability of fCO2 in seawater and air–sea CO2 fluxes in a coastal upwelling system (Ría de Vigo, NW Spain)

Coastal upwelling systems are regions with highly variable physical processes and very high rates of primary production and very little is known about the effect of these factors on the short-term variations of CO₂ fugacity in seawater (fCO₂^w). This paper presents the effect of short-term variability (<1 week) of upwelling–downwelling events on CO₂ fugacity in seawater (fCO₂^w), oxygen, temperature and salinity fields in the Ría de Vigo (a coastal upwelling ecosystem). The magnitude of fCO₂^w values is physically and biologically modulated and ranges from 285 μatm in July to 615 μatm in October. There is a sharp gradient in fCO₂^w between the inner and the outer zone of the Ría during almost all the sampling dates, with a landward increase in fCO₂^w.CO₂ fluxes calculated from local wind speed and air–sea fCO₂ differences indicate that the inner zone is a sink for atmospheric CO₂ in December only (−0.30 mmol m⁻² day⁻¹). The middle zone absorbs CO₂ in December and July (−0.05 and −0.27 mmol·m⁻² day⁻¹, respectively). The oceanic zone only emits CO₂ in October (0.36 mmol·m⁻² day⁻¹) and absorbs at the highest rate in December (−1.53 mmol·m⁻² day⁻¹).     

One-dimensional modelling of convective CO2 exchange in the Tropical Atlantic

Diurnal changes in seawater temperature affect the amount of air–sea gas exchange taking place through changes in solubility and buoyancy-driven nocturnal convection, which enhances the gas transfer velocity. We use a combination of in situ and satellite derived radiometric measurements and a modified version of the General Ocean Turbulence Model (GOTM), which includes the National Oceanic and Atmospheric Administration Coupled-Ocean Atmospheric Response Experiment (NOAA-COARE) air–sea gas transfer parameterization, to investigate heat and carbon dioxide exchange over the diurnal cycle in the Tropical Atlantic. A new term based on a water-side convective velocity scale (w_*w) is included, to improve parameterization of convectively driven gas transfer. Meteorological data from the PIRATA mooring located at 10°S10°W in the Tropical Atlantic are used, in conjunction with cloud cover estimates from Meteosat-7, to calculate fluxes of longwave, latent and sensible heat along with a heat budget and temperature profiles during February 2002. Twin model experiments, representing idealistic and realistic conditions, reveal that over daily time scales the additional contribution to gas exchange from convective overturning is important. Increases in transfer velocity of up to 20% are observed during times of strong insolation and low wind speeds (<6 m s⁻¹); the greatest enhancement from w_*w to the CO₂ flux occurs when diurnal warming is large. Hence, air–sea fluxes of CO₂ calculated using simple parameterizations underestimate the contribution from convective processes. The results support the need for parameterizations of gas transfer that are based on more than wind speed alone and include information about the heat budget.     

The distribution and population dynamics of O-group plaice (Pleuronectes platessa L.) on nursery grounds in the Firth of Forth

Midsummer (1 August) population estimates of about 2 million O-group plaice (Pleuronectes platessa L.) were derived for sandy bays around the Firth of Forth in 1979–1980. This is an order of magnitude less than similar estimates made for the Clyde Sea Area in 1973–1974. Autumn population estimates of 0·4–1·0 million fish were comparable to estimates by the Ministry of Agriculture, Fisheries and Food for the area between the Scottish border and Flamborough Head (2·3 million for 1970 and 1973) which represented 4·8% (1973) to 5·3% (1970) of the total number of O-group fish on the English east coast.Largo Bay was the most important nursery area holding 25% of the total population. It is particularly well situated to receive newly metamorphosed plaice carried in water currents along the north side of the Forth from the spawning ground off Fife Ness. Plaice in the Forth are mainly distributed on fine to medium sandy beaches (186–480 μm), the mean number per haul in midsummer (D) being correlated with the median diameter (m.d. in μm) of the low water sediments by the equation: D=−45·7666+0·2327 m.d. (n=11,r=0·68,P<0·02 but>0·01).The shallow inshore water in sandy bays in the outer Firth was well mixed and more marine than estuarine (27·7–35·0‰). The correlation coefficient between fish density and water temperature was low, while that with salinity (S‰) was: D=6·1618+0·2238S (n=23,r=0·62,P<0·005).Regression analysis demonstrated that the relationship between the instantaneous mortality rate (Z) and the initial population density (D_p) was: Z×100=0·7480+0·0546d_p (n=12,r=0·87,P<0·001).The mean mortality rate for the O-group plaice in the Forth nursery areas was 53% month⁻¹.     

海洋白浪寿命的定义及测量结果

通过对国内外白浪研究和应用的分析,首次提出了有效白浪寿命的定义,给出了计算白浪寿命的公式及测量方法和结果,并报告了以此方法在渤海实测的结果,得到了白浪寿命ＴＬ与海面风速Ｕ１０的关系为ＴＬ＝０．２６Ｕ１０以及白浪寿命概率分布近于瑞利分布等。     

A note on a critical wind speed for air–sea boundary processes

Wind and wind-generated waves were measured in a wind-wave tank. A clear transition was found in the relation between the wind speed U ₁₀ and the wind friction velocity u _* near u _* = 0.2 m/s, where U ₁₀ is the wind speed at 10 m height extrapolated from the measured wind profile in a logarithmic layer, and u _* = 0.2 m/s corresponds roughly to U ₁₀ = 8 m/s in the present measurement. Quite a similar transition was found in the relation between the spectral density of high frequency wind waves and u _*. These results suggest the existence of the critical wind speed for air–sea boundary processes, which was proposed by Munk (J Marine Res 6:203–218, 1947) more than half a century ago. His original idea of the critical wind speed was based on the discontinuities in such phenomena as white caps, wind stress, and evaporation, which commonly appear at a wind speed near 7 m/s. On the basis of the results of our present study and those of earlier studies, we discuss the phenomena which are relevant to the critical wind speed for the air–sea boundary processes. The conclusion is that the critical wind speed exists and it is attributed to the start of wave breaking rather than the Kelvin–Helmholtz instability, but the air–sea boundary processes are not discontinuous at a particular wind speed; because of the stochastic nature of breaking waves, the changes occur over a range of wind speeds. Detailed discussions are presented on the dynamical processes associated with the critical wind speed such as wind-induced change of sea surface roughness and high frequency wave spectrum. Future studies are required, however, to clarify the dynamical processes quantitatively. In particular, there is a need to further examine the gradual change of breaking patterns of wind waves with the increase of wind speed, and the associated change of the structure of the wind over wind waves, such as separation of the airflow at the crest of wind waves, the turbulent stress, and wave-induced stress. Studies on the dynamical structure of the high frequency wave spectrum are also needed.     

Anomalously low alkenone temperatures caused by lateral particle and sediment transport in the Malvinas Current region, western Argentine Basin

We analysed the alkenone unsaturation ratio (U^K′₃₇) in 87 surface sediment samples from the western South Atlantic (5°N–50°S) in order to evaluate its applicability as a paleotemperature tool for this part of the ocean. The measured U^K′₃₇ ratios were converted into temperature using the global core-top calibration of Müller et al. (1998) and compared with annual mean atlas sea-surface temperatures (SSTs) of overlying surface waters. The results reveal a close correspondence (<1.5°C) between atlas and alkenone temperatures for the Western Tropical Atlantic and the Brazil Current region north of 32°S, but deviating low alkenone temperatures by −2° to −6°C are found in the regions of the Brazil–Malvinas Confluence (35–39°S) and the Malvinas Current (41–48°S). From the oceanographic evidence these low U^K′₃₇ values cannot be explained by preferential alkenone production below the mixed layer or during the cold season. Higher nutrient availability and algal growth rates are also unlikely causes. Instead, our results imply that lateral displacement of suspended particles and sediments, caused by strong surface and bottom currents, benthic storms, and downslope processes is responsible for the deviating U^K′₃₇ temperatures. In this way, particles and sediments carrying a cold water U^K′₃₇ signal of coastal or southern origin are transported northward and offshore into areas with warmer surface waters. In the northern Argentine Basin the depth between displaced and unaffected sediments appears to coincide with the boundary between the northward flowing Lower Circumpolar Deep Water (LCDW) and the southward flowing North Atlantic Deep Water (NADW) at about 4000 m.     

Effects of alkaline aluminate waste dumping on seawater chemistry

The alkaline aluminate waste, of which 1000–2000 tonnes are dumped a few times a year off the Belgian coast in the Southern Bight of the North Sea, contains 5·4% NaCl, 1·8% dissolved Al and 7·4% NaOH, in addition to traces of heavy metals and some aniline- and phenol-derivatives. The pH rises locally to 8-5 and the total Al-concentration reaches 120niglitre⁻¹(corresponding to an initial waste dilution factor of only 150) in the 10-m wide track just beyond the discharging barge, but these decay quickly to pH 8·1 and 1 mg litre^-1 in the 30-m wide track, 500m behind the barge. The relation between the waste concentration and seawater pH was studied. The white precipitate that forms immediately in the sea was identified as Mg₆---Al₂---CO₃---(OH)₁₆---4H₂0 (hydrotalcite-manasseite like). No trace of it was found in the local sediments.     

基于一维综合孔径微波辐射计的海面温度反演研究

Due to the low spatial resolution of sea surface temperature(T_S) retrieval by real aperture microwave radiometers,in this study, an iterative retrieval method that minimizes the differences between brightness temperature(T_B)measured and modeled was used to retrieve sea surface temperature with a one-dimensional synthetic aperture microwave radiometer, temporarily named 1 D-SAMR. Regarding the configuration of the radiometer, an angular resolution of 0.43° was reached by theoretical calculation. Experiments on sea surface temperature retrieval were carried out with ideal parameters; the results show that the main factors affecting the retrieval accuracy of sea surface temperature are the accuracy of radiometer calibration and the precision of auxiliary geophysical parameters. In the case of no auxiliary parameter errors, the greatest error in retrieved sea surface temperature is obtained at low T_S scene(i.e., 0.710 6 K for the incidence angle of 35° under the radiometer calibration accuracy of0.5 K). While errors on auxiliary parameters are assumed to follow a Gaussian distribution, the greatest error on retrieved sea surface temperature was 1.330 5 K at an incidence angle of 65° in poorly known sea surface wind speed(W)(the error on W of 1.0 m/s) over high W scene, for the radiometer calibration accuracy of 0.5 K.     

Carbon dioxide flux techniques performed during GasEx-98

A comprehensive study of air–sea interactions focused on improving the quantification of CO₂ fluxes and gas transfer velocities was performed within a large open ocean CO₂ sink region in the North Atlantic. This study, GasEx-98, included shipboard measurements of direct covariance CO₂ fluxes, atmospheric CO₂ profiles, atmospheric DMS profiles, water column mass balances of CO₂, and measurements of deliberate SF₆–³He tracers, along with air–sea momentum, heat, and water vapor fluxes. The large air–sea differences in partial pressure of CO₂ caused by a springtime algal bloom provided high signals for accurate CO₂ flux measurements. Measurements were performed over a wind speed range of 1–16 m s⁻¹ during the three-week process study. This first comparison between the novel air-side and more conventional water column measurements of air–sea gas transfer show a general agreement between independent air–sea gas flux techniques. These new advances in open ocean air–sea gas flux measurements demonstrate the progress in the ability to quantify air–sea CO₂ fluxes on short time scales. This capability will help improve the understanding of processes controlling the air–sea fluxes, which in turn will improve our ability to make regional and global CO₂ flux estimates.     

A global wave parameter database for geophysical applications. Part 1: Wave-current–turbulence interaction parameters for the open ocean based on traditional parameterizations

Ocean surface mixing and drift are influenced by the mixed layer depth, buoyancy fluxes and currents below the mixed layer. Drift and mixing are also functions of the surface Stokes drift U_ss, volume Stokes transport T_S, a wave breaking height scale H_swg, and the flux of energy from waves to ocean turbulence Φ_oc. Here we describe a global database of these parameters, estimated from a well-validated numerical wave model, that uses traditional forms of the wave generation and dissipation parameterizations, and covers the years 2003–2007. Compared to previous studies, the present work has the advantage of being consistent with the known physical processes that regulate the wave field and the air–sea fluxes, and also consistent with a very large number of in situ and satellite observations of wave parameters. Consequently, some of our estimates differ significantly from previous estimates. In particular, we find that the mean global integral of Φ_oc is 68 TW, and the yearly mean value of T_S is typically 10–30% of the Ekman transport, except in well-defined regions where it can reach 60%. We also have refined our previous estimates of U_ss by using a better treatment of the high frequency part of the wave spectrum. In the open ocean, U_ss  0.013U₁₀, where U₁₀ is the wind speed at 10 m height.     

Time-series sediment trap record of alkenones from the western Sea of Okhotsk

C₃₇–C₃₉ alkenones were measured in time-series sediment trap samples collected from August 1998 to June 2000 at two depths in the seasonal sea ice region of the western Sea of Okhotsk, off Sakhalin, in order to investigate alkenone production and water-column processes in the region. Measurable export fluxes of alkenones are ranged from < 0.1 to 5.8 μg/m²/day and clearly showed that the alkenone production was restricted to autumn. In 1998, maximum export flux of alkenones occurred in September when surface water column was well stratified with low nutrients in the surface mixing layer. In the next year, the maximum flux is observed in October. Comparison between alkenone temperature and satellite based sea surface temperature (SST) shows that the estimated alkenone temperatures in August 1998 were found to be  10 °C lower than the temporal satellite SST, suggesting that alkenones are produced in surface to subsurface thermocline layers during the period. Annual mean flux of alkenones is lower in the lower traps than that of the upper traps, suggesting rapid degradation of alkenones in water column, but the U^K₃₇′ value is not significantly altered. This study indicates that U^K₃₇′ values preserved in the surface sediments off Sakhalin reflect the seasonal temperature signal of near surface water, rather than annual mean surface temperature.     

Modelling the monthly sea surface fCO2 fields in the Indian Ocean

In order to construct monthly fields of sea surface fugacity of carbon dioxide (f_CO2) on a large scale in the Indian Ocean, we use a one-dimensional model which takes into account the main physical and biogeochemical processes controlling f_CO2 variations in the ocean. Physical and biogeochemical processes are constrained by the monthly variations of sea surface temperature, salinity, chlorophyll concentration, wind speed and mixed-layer depth. The model is applied to four locations in the Indian Ocean and it well predicts observed temporal variations in f_CO2 at these locations. Regarding to monthly f_CO2 observations, the model also well simulates the f_CO2 distribution and its temporal variations along a track located between 20 ° and 50 °S with a maximal error of + 10 μatm. The model is also used to predict f_CO2 for 2 ° × 2 ° grids over the entire Indian Ocean and simulates seasonal cycles that are consistent with observations. The monthly f_CO2 fields derived from the model are used to estimate a global air-sea CO₂ flux over the Indian Ocean basin. We estimate a net sink of 0.5 Gt/yr C for the Indian Ocean (20 °N-50 °S), with the main sink located between 20 ° and 50 °S.     

Seasonal and interannual variability of the air–sea CO2 flux in the Atlantic sector of the Barents Sea

The seasonal and interannual variability of the air–sea CO₂ flux (F) in the Atlantic sector of the Barents Sea have been investigated. Data for seawater fugacity of CO₂ (fCO₂^sw) acquired during five cruises in the region were used to identify and validate an empirical procedure to compute fCO₂^sw from phosphate (PO₄), seawater temperature (T), and salinity (S). This procedure was then applied to time series data of T, S, and PO₄ collected in the Barents Sea Opening during the period 1990–1999, and the resulting fCO₂^sw estimates were combined with data for the atmospheric mole fraction of CO₂, sea level pressure, and wind speed to evaluate F.The results show that the Atlantic sector of the Barents Sea is an annual sink of atmospheric CO₂. The monthly mean uptake increases nearly monotonically from 0.101 mol C m^− 2 in midwinter to 0.656 mol C m^− 2 in midfall before it gradually decreases to the winter value. Interannual variability in the monthly mean flux was evaluated for the winter, summer, and fall seasons and was found to be ± 0.071 mol C m^− 2 month^− 1. The variability is controlled mainly through combined variation of fCO₂^sw and wind speed. The annual mean uptake of atmospheric CO₂ in the region was estimated to 4.27 ± 0.68 mol C m^− 2.     

Evaluation of the Diurnal Warming of Sea Surface Temperature Using Satellite-Derived Marine Meteorological Data

In order to produce a high-quality sea surface temperature (SST) data set, the daily amplitude of SST (ΔSST) should be accurately known. The purpose of this study was to evaluate the diurnal variation of sea surface temperature in a simple manner. The authors first simulated ΔSST with a one-dimensional numerical model using buoy-observed meteorological data and satellite-derived solar radiation data. When insolation is strong, the model-simulated 1-m-depth ΔSST becomes much smaller than the in situ value as wind speed decreases. By forcibly mixing the sea surface layer, the model ΔSST becomes closer to the in situ value. It can be considered that part of this difference is due to the turbulence induced by the buoy hull. Then, on the assumption that the model results were reliable, the authors derived a regression equation to evaluate ΔSST at the skin and 1-m depth from daily mean wind speed (U) and daily peak solar radiation (PS). ΔSST is approximately proportional to In(U) and (PS)², and the skin ΔSST estimated by the equation is not inconsistent with in situ observation results reported in past studies. The authors prepared maps of PS and U using only satellite data, and demonstrated the ΔSST evaluation over a wide area. The result showed that some wide patchy areas where the skin ΔSST exceeds 3.0 K can appear in the tropics and the mid-latitudes in summer. This revised version was published online in July 2006 with corrections to the Cover Date.     

Nitrogen cycling in deep-sea sediments of the Porcupine Abyssal Plain, NE Atlantic

Rates of transformation, recycling and burial of nitrogen and their temporal and spatial variability were investigated in deep-sea sediments of the Porcupine Abyssal Plain (PAP), NE Atlantic during eight cruises from 1996 to 2000. Benthic fluxes of ammonium (NH₄) and nitrate (NO₃) were measured in situ using a benthic lander. Fluxes of dissolved organic nitrogen (DON) and denitrification rates were calculated from pore water profiles of DON and NO₃, respectively. Burial of nitrogen was calculated from down core profiles of nitrogen in the solid phase together with ¹⁴C-based sediment accumulation rates and dry bulk density. Average NH₄ and NO₃-effluxes were 7.4 ± 19 μmol m⁻² d⁻¹ (n = 7) and 52 ± 30 μmol m⁻² d⁻¹ (n = 14), respectively, during the period 1996–2000. During the same period, the DON-flux was 11 ± 5.6 μmol m⁻² d⁻¹ (n = 5) and the denitrification rate was 5.1 ± 3.0 μmol m⁻² d⁻¹ (n = 22). Temporal and spatial variations were only found in the benthic NO₃ fluxes. The average burial rate was 4.6 ± 0.9 μmol m⁻² d⁻¹. On average over the sampling period, the recycling efficiency of the PON input to the sediment was 94% and the burial efficiency hence 6%. The DON flux constituted 14% of the nitrogen recycled, and it was of similar magnitude as the sum of burial and denitrification. By assuming the PAP is representative of all deep-sea areas, rates of denitrification, burial and DON efflux were extrapolated to the total area of the deep-sea floor (>2000 m) and integrated values of denitrification and burial of 8 ± 5 and 7 ± 1 Tg N year⁻¹, respectively, were obtained. This value of total deep-sea sediment denitrification corresponds to 3–12% of the global ocean benthic denitrification. Burial in deep-sea sediments makes up at least 25% of the global ocean nitrogen burial. The integrated DON flux from the deep-sea floor is comparable in magnitude to a reported global riverine input of DON suggesting that deep-sea sediments constitute an important source of DON to the world ocean.     

Some effects of certain metals on development and mortality within the moult cycle of crangon crangon (L.)

Crangon crangon (L.) (Crustacea: Decapoda) were subjected to three metals (Cd, Cu and Zn) at respective concentrations (0·35, 1·6 and 14·4mg/litre) that would induce 50% mortality in 5–6 days. Significantly different (P < 0·05) rates of mortality occurred in the three groups, although their ET₅₀'s were almost identical. No moultstage-dependent mortality was observed for Cd, but in Cu and Zn the post-moult stages (A and B) were more susceptible than either the intermoult (C_a and C_b) or premoult (D₀–D₄) stages.No significant (P> 0·05) sex- or size-dependent mortality was evident at the tested concentrations of these metals.Feeding and ecdysis were inhibited in all metals after day 3 and development from one moult-stage to the next was greatly reduced in Cd.Possible modes of toxicity have been discussed in terms of the physiological condition inferred from the moult-stage at the time of death.     

Scyphomedusae in surface waters near the Oregon coast, May–August 1981

The scyphomedusae of the surface waters off Oregon and southern Washington were collected with commercial purse seines from May–August 1981. Twelve east-west transects, located from north of the Columbia River to south of Coos Bay were sampled from the 37 m isobath to distances up to 48 km from shore. Chrysaora fuscescens was the dominant species collected in each month. Maximum sampled abundances reached 18001 of medusae per 10⁵ m³. Using an estimated carbon content of 0·280% of wet weight, this medusa density was calculated to contain approximately 50 mg Cm⁻³. Seven of the 263 samples contained so many medusae that they exceeded the capacity of the sampling gear. In all months but May, when medusa densities were relatively low, the density of C. fuscescens was greatest closest to shore and decreased rapidly offshore. Mean umbrella diameter increased from 8·6 cm in May to 18·5 cm in August, while the largest specimens increased from 19 cm in May to 37 cm in August. Aurelia aurtia, Cyanea capillata and Phacellophora camtschatica were also collected, but were much less abundant than C. fuscescens. The relative abundance of C. fuscescens was compared with the maximum abundance of copepods off the Oregon coast, and the hydrographic features influencing medusa distribution patterns are discussed.     

Reconstruction of sea surface dimethylsulfide in the North Pacific during 1970s to 2000s

We proposed an empirical equation of sea surface dimethylsulfide (DMS, nM) using sea surface temperature (SST, K), sea surface nitrate (SSN, μM) and latitude (L, °N) to reconstruct the sea surface flux of DMS over the North Pacific between 25°N and 55°N: ln DMS = 0.06346 · SST − 0.1210 · SSN − 14.11 · cos(L) − 6.278 (R² = 0.63, p < 0.0001). Applying our algorithm to climatological hydrographic data in the North Pacific, we reconstructed the climatological distributions of DMS and its flux between 25 °N and 55 °N. DMS generally increased eastward and northward, and DMS in the northeastern region became to 2–5 times as large as that in the southwestern region. DMS in the later half of the year was 2–4 times as large as that in the first half of the year. Moreover, applying our algorithm to hydrographic time series datasets in the western North Pacific from 1971 to 2000, we found that DMS in the last three decades has shown linear increasing trends of 0.03 ± 0.01 nM year^− 1 in the subpolar region, and 0.01 ± 0.001 nM year^− 1 in the subtropical region, indicating that the annual flux of DMS from sea to air has increased by 1.9–4.8 μmol m^− 2 year^− 1. The linear increase was consistent with the annual rate of increase of 1% of the climatological averaged flux in the western North Pacific in the last three decades.