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.     

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.     

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

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

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⁻¹.     

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.     

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.     

Seasonal variations of alkenones and U37 in the Chesapeake Bay water column

Alkenone unsaturation indices (U^K₃₇ and U^K′₃₇) have long been used as proxies for surface water temperature in the open ocean. Recent studies have suggested that in other marine environments, variables other than temperature may affect both the production of alkenones and the values of the indices. Here, we present the results of a reconnaissance field study in which alkenones were extracted from particulate matter filtered from the water column in Chesapeake Bay during 2000 and 2001. A multivariate analysis shows a strong positive correlation between U^K₃₇ (and U^K′₃₇) values and temperature, and a significant negative correlation between U^K₃₇ (and U^K′₃₇) values and nitrate concentrations. However, temperature and nitrate concentrations also co-vary significantly. The temperature vs. U^K₃₇ relationships (U^K₃₇=0.018 (T)−0.162, R²=0.84, U^K′₃₇=0.013 (T)−0.04, R²=0.80) have lower slopes than the open-ocean equations of Prahl et al. [1988. Further evaluation of long-chain alkenones as indicators of paleoceanographic conditions. Geochimica et Cosmochimica Acta 52, 2303–2310] and Müller et al. [1998. Calibration of the alkenone paleotemperature index U^K′₃₇ based on core-tops from the eastern South Atlantic and the global ocean (60°N–60°S). Geochimica et Cosmochimica Acta 62, 1757–1772], but are similar to the relationships found in controlled studies with elevated nutrient levels and higher nitrate:phosphate (N:P) ratios. This implies that high nutrient levels in Chesapeake Bay have either lowered the U^K₃₇ vs. temperature slope, or nutrient levels are the main controller of the U^K₃₇ index. In addition, particularly high abundances (>5% of total C₃₇ alkenones) of the tetra-unsaturated ketone, C_37:4, were found when water temperatures reached 25 °C or higher, thus posing further questions about the controls on alkenone production as well as the biochemical roles of alkenones.     

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.     

Distribution of the CO2 partial pressure along an Atlantic meridional transect

Atmospheric and oceanic pCO₂ were measured continuously along an Atlantic Meridional transect (50°N–50°S) in September–October 1995 and 1996 (U.K. to the Falklands Islands) and in April–May 1996 (Falklands Islands to the UK). The Atlantic ocean was a net sink for atmospheric CO₂ for all 3 transects. The largest sinks were located at high latitudes, in regions of high wind speed, where strong CO₂ undersaturations, associated with high biological activity, were observed. In these regions the partial pressure difference between the ocean and the atmosphere reached −110 μatm. A CO₂ source occurred in the equatorial region between 0° and 10°S, where ΔpCO₂ of up to 40 μatm was found. Another source was in the northern subtropical gyre where its extension varied according to the season. Along the whole transect the October cruises exhibited similar pCO₂ distributions suggesting a dominance of the seasonal variability and small year to year changes.     

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.     

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.     

Elemental mercury concentrations and formation rates in the Scheldt estuary and the North Sea

Concentrations of Hg⁰ in surface waters and atmosphere of the Scheldt estuary and the North Sea are presented and their relationship with biological processes is discussed. Hg⁰ concentrations in the Scheldt estuary range from 0.1 to 0.38 pmol·l⁻¹ in the winter and from 0.24 to 0.65 pmol·l⁻¹ in the summer and show a positive relationship with phytoplankton pigments. In the North Sea Hg⁰ concentrations range from 0.06 to 0.8 pmol·l⁻¹ and are higher in coastal stations. Transfer velocities across the air–sea interface were calculated using a classical shear turbulence model. Volatilization fluxes of Hg⁰ were calculated for the Scheldt estuary and the North Sea. For the Scheldt estuary the fluxes range from 226–284 pmol·m⁻²·d⁻¹ in winter and 500–701 pmol·m⁻²·d⁻¹ in summer and for the North Sea the fluxes range from 59–1110 pmol·m⁻²·d⁻¹ for an average windspeed of 8.1 m·s⁻¹. These fluxes are comparable to the wet and dry depositional fluxes to the North Sea. Hg⁰ formation rates necessary to balance the volatilization fluxes vary from 0.2 to 4% d⁻¹.     

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.     

Thorium isotopes as tracers of particles dynamics and deep water circulation in the Indian sector of the Southern Ocean (ANTARES IV)

Dissolved and particulate samples were collected to study the distribution of thorium isotopes (²³⁴Th, ²³²Th and ²³⁰Th) in the water column of the Indian sector of the Southern Ocean (from 42°S to 47°S and from 60°E to 66°E, north of the Polar Front) during Austral summer 1999. Vertical profiles of excess ²³⁰Th (²³⁰Th_xs) increases linearly with depth in surface water (0–100 m) and a model was applied to estimate a residence time relative to the thorium scavenging (τ_scav). Low τ_scav in the Polar Front Zone (PFZ) are found, compared to those estimated in the Subtropical Front Zone (STZ). Changes in particle composition between the PFZ and STZ could influence the ²³⁰Th_xs scavenging efficiency and explain this difference. An innovative coupling between ²³⁴Th and ²³⁰Th_xs was then used to simultaneously constrain the settling velocities of small (0.6–60 μm) and large (above 60 μm) particles. Although the different hydrological and biogeochemical regimes visited during the ANTARES IV cruise did not explain the spatial variation of sinking velocity estimates, our results indicate that less particles may reach the seafloor north (60 ± 2 m d^− 1, station 8) than south of the Agulhas Return Current (119 ± 23 and 130 ± 5 m d^− 1 at stations 3 and 7, respectively). This information is essential for understanding particle transport and by extension, carbon export. In the deep water column, the ²³⁰Th_xs concentrations did not increase linearly with depth, probably due to lateral transport of North Atlantic Deep Water (NADW) from the Atlantic to the Indian sector, which renews the deep waters and decreases the ²³⁰Th_xs concentrations. A specific ²³⁰Th_xs transport model is applied in the deep water column and allows us to assess a “travel time” of NADW ranging from 2 to 15 years.     

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.