Study of the sedimentation and self-weight consolidation behavior of seafloor sediments using a radioisotope densitometer

ABSTRACT

In this study, settling tests were conducted to investigate the sedimentation and self-weight consolidation behavior of seafloor sediments from Isahaya Bay, Ariake Sea, Japan. During the tests, the density variations with depth and time were measured by a gamma-ray transmission radioisotope densitometer. The test results show that the settling process of the seafloor sediments can be classified into the flocculation stage, settling stage, and consolidation stage. The settling rate of the seafloor sediments in the settling stage is dependent on the temperature and initial water content, while the settling rate in the consolidation stage is independent of the temperature and initial water content. The density profile changes from a constant density profile to a linear density profile when the sedimentation process transitions to the self-weight consolidation process. The relations between the void ratio (e) and effective vertical stress (p’) at very low pressures can be calculated from the measured density values, and this can be used for the analysis of the self-weight consolidation of seafloor sediments. For the seafloor sediments tested in this study, the undrained shear strength (s_u) values are almost the same when the density values are less than 1.14?g/cm³, and the s_u values increase linearly with an increase in density when the density values are in the range of 1.14–1.2?g/cm^3.     

Identifying late Quaternary coastal deposits in Kyonggi Bay, Korea, by their geotechnical properties

Based on cone penetration tests with pore pressure measurements (CPTUs) and standard penetration tests (SPTs), the geotechnical properties of five lithostratigraphic units were determined during the construction of Incheon international airport on reclaimed macrotidal flats in Kyonggi Bay, Korea. Two late Pleistocene non-marine units (unit V and unit IV) display largest N values (cf. number of blows required to achieve a standard penetration), reflecting coarse-grained and overconsolidated sediments. Tidal channel and tidal flat facies (unit IIIb) consist of unweathered late Pleistocene tidal sand and mud. The tidal channel facies is characterized by upward-decreasing cone resistance (q _t) and sleeve friction (f _s) with negative pore pressures (u _bt), reflecting a fining-upward textural trend. The tidal flat facies, by contrast, is represented by uniformly low q _t and f _s values with high friction ratios (FRs), suggesting homogeneous muddy deposits. Two overconsolidated units, a weathered late Pleistocene tidal mud (unit IIIa) and an early Holocene organic-rich non-marine mud (unit II), are characterized by high q _t, f _s, FRs and N values, unit IIIa being much more consolidated than unit II. Holocene tidal sands and muds (unit I) show the smallest q _t and f _s values with positive u _bt. These are slightly more consolidated than the tidal flat facies of unit IIIb. Two unconformable boundaries (a sequence boundary and a transgressive surface) have also been identified on some CPTU and SPT profiles. The boundaries are indicated by gradual but sharp increases in q _t, f _s and N values with an abrupt drop of u _bt, which indicates the contact between two units showing contrasting rigidity. The regional pattern produced by the unconformable boundaries indicates the presence of late Pleistocene valleys which pass through the middle of study area. The location of the valleys seems to be controlled by the antecedent basement morphology.     

In situ observations of wave pumping of sediments in the Yellow River Delta with a newly developed benthic chamber

A specially designed benthic chamber for the field observation of sediment resuspension that is caused by the wave-induced oscillatory seepage effect (i.e., the wave pumping of sediments) is newly developed. Observational results from the first sea trial prove that the geometry design and skillful instrumentation of the chamber well realize the goal of monitoring the wave pumping of sediments (WPS) continuously. Based on this field dataset, the quantitative contribution of the WPS to the total sediment resuspension is estimated to be 20–60% merely under the continuous action of normal waves (H_s ≤ 1.5 m) in the subaqueous Yellow River Delta (YRD). Such a large contribution invalidates a commonly held opinion that sediments are purely eroded from the seabed surface by the horizontal “shearing effect” from the wave orbital or current velocities. In fact, a considerable amount of sediments could originate from the shallow subsurface of seabed driven by the vertical “pumping effect” of the wave-generated seepage flows during wavy periods. According to the new findings, an improved conceptual model for the resuspension mechanisms of silty sediments under various hydrodynamics is proposed for the first time.     

波流作用下海底边界层沉积物再悬浮与影响因素研究

海底沉积物再悬浮及其分布取决于海洋水动力、沉积物类型与床面形态之间复杂的相互作用,准确地理解和确定沉积物再悬浮过程对于沉积物输运的研究具有重要的意义。本文在祥云湾海洋牧场典型海域开展现场原位观测,获取研究区波浪、海流及悬浮沉积物浓度数据;分析了波、流作用下海底边界层悬浮沉积物垂向分布特征,并探究了海洋水动力和床面形态对悬浮沉积物垂向分布的影响。结果表明,研究区波流之间的相互作用不显著,沉积物再悬浮受控于风暴浪作用,风暴浪作用下底床切应力可以达到沉积物临界切应力的10～15倍,沉积物的再悬浮滞后于风暴浪作用2～3 h。在波浪荷载微小的情况下,悬浮沉积物垂向分布呈现"I"型,波浪荷载下,悬浮沉积物垂向分布呈现幂指函数分布,表现为"L"型;床面形态随波、流作用而演化,影响沉积物的再悬浮过程,u_(?w)/u_(?c)=1.00可作为波浪和海流起主导控制作用的床面形态的判别依据,纯波浪荷载作用下的u_(?w)/u_(?c)显著高于波浪主控作用下,但二者之间的界线随着波浪荷载的增加而升高。     

Does suspended matter drained from the Isahaya freshwater reservoir cause organic enrichment in the northern Ariake Bay?

The inner part of Isahaya Bay was converted to a freshwater reservoir following the closure of the land claim dike in 1997. Turbid water drains into Isahaya Bay when water levels increase. We investigated whether particulate organic matter (POM) from the reservoir in Isahaya Bay has caused bottom organic enrichment in the northern part of Ariake Bay. Using potential end-members from before to after the rainy seasons, during which a frequent discharge from the reservoir was expected, stable isotope analyses were performed on sediments collected from Isahaya Bay and northern Ariake Bay. Each end-member was isotopically differentiated by δ¹³C and δ¹⁵N (riverine POM: ?28.5 to ?27.2‰ and 3.3–4.6‰; reservoir POM: ?25.7 to ?25.3‰ and 7.4–8.4‰; marine POM: ?21.8 to ?19.7‰ and 6.7–7.6‰; microphytobenthos estimated from consumers: ?16.1 to ?15.9‰ and 5.2–6.1‰, respectively). Sediment isotopic signatures fell within the mixing space defined by the signatures of the end-members. Marine POM contributed greatly to bottom sediments in both seasons in Isahaya Bay and Ariake Bay, ranging from ca. 60–70 and 40–60%, respectively. Reservoir POM contributed around 10% to bottom sediments. This percentage slightly increased in the sediment of Isahaya Bay after the rainy season, but decreased in the sediment of Ariake Bay. Thus, most of the POM discharged from the reservoir would not reach the northern part of Ariake Bay and would not be a major contributor to organic enrichment. This study is the first to quantitatively describe the contribution of drained reservoir POM outside Isahaya Bay.     

Variability of the partial pressure of CO2 on a daily-to-seasonal time scale in a shallow coastal system affected by intensive aquaculture activities (Bay of Cadiz, SW Iberian Peninsula)

The present study describes the temporal variability of the water fCO₂ as well as the different driving forces controlling this variability, on time scales from daily to seasonal, in the Rio San Pedro, a tidal creek located in a salt marsh area in the Bay of Cadiz (SW Iberian Peninsula). This shallow tidal creek system is affected by effluents of organic matter and nutrients from the surrounding marine fish farms. Continuous pCO₂, salinity and temperature were recorded for four periods of approximately one month, between February and September in 2004.Major processes controlling the CO₂ variability are related to three different time scales. Daily variations in fCO₂ are controlled by tidal advection and mixing of the water from within the creek and the seawater that enters from the Bay of Cadiz. Significant cyclical variations of the fCO₂ have been observed with the maximum values occurring at low tide. On a fortnightly time scale, the amplitude of the daily variability of fCO₂ is modulated by the variations in the residence time of the water within the creek, which are related to the spring–neap tide sequence.On a third time scale, high seasonal variability is observed for the temperature, salinity and fCO₂. Maximum and minimum values for fCO₂ were 380 µatm and 3760 µatm for February and July respectively. Data suggest that seasonal variability is related to the seasonal variability in discharges from the fish farm and to the increase with temperature of organic matter respiratory processes in the tidal creek. The fCO₂ values observed are in the same range as several highly polluted European estuaries or waters surrounding mangrove forests. From the air–water CO₂ flux computed, it can be concluded that the Rio San Pedro acts as a source of CO₂ to the atmosphere throughout the year, with the summer accounting for the higher average monthly flux.     

Measurements of sediment transport on the sea bed southwest of England

Sand transport rates were measured via bedload traps that were inserted into the sediment surface on the continental shelf southwest of England. Analysis of the trapped mobile sediment shows this to be finer and better sorted than nearby Shipek grab samples. The transport rates are combined with simultaneous, near-bed flow measurements to assess a Bagnold type predictive equation,j = k ₁(u ² ₁₀₀ -u ² 100cr)u ₁₀₀. The measured transport rates vary between 0.41 × 10^–3 and 1.67 × 10^–3 gm/cm s. The data yield a mean calibration coefficient (k ₁) of 0.4 × 10^–6 which is slightly lower than values computed from flume and shallow water data.     

基于Chirp数据和Biot-Stoll模型反演南海北部陆坡海底表层沉积物物理性质

浅地层剖面是基于声学信号(频率在几百至几千赫兹)在沉积物中的传播得到可反映沉积地层结构的数据,海底反射系数与沉积物物理性质密切相关。Biot-Stoll声波传播理论模型可以预测海底沉积物的物理性质,构建反射系数等声学参数与物理参数之间的关系,但在不同的海域采用不同的参数所获得的效果不同。为此,本文基于南海北部陆坡海底表层沉积物的实测物理参数,利用Biot-Stoll模型建立研究区海底反射系数和沉积物物理性质之间的关系,结果表明模型计算值与样品实测值吻合度总体较好,偏差在0.1%～4.9%之间,并建立了频率3.5 kHz时海底反射系数与沉积物孔隙度、密度、平均粒径之间的关系方程,且方程拟合度较高,可决系数R2均大于0.99。在对典型Chirp剖面数据计算其海底反射系数的基础上,反演了海底表层沉积物的孔隙度、密度、颗粒平均粒径等物理性质,其中反演孔隙度、密度、平均粒径与实测孔隙度、密度、平均粒径相对误差均小于5%,结果与实测值基本相符,表明该反演方法在南海北部陆坡区的应用是可行的。     

Technique and Results of Measurements of Turbulent Helicity in a Stratified Surface Layer

In August 2014, measurements of the turbulent velocity rotor, turbulent temperature gradient, turbulent helicity, and turbulent potential vortex were performed at the Obukhov Institute of Atmospheric Physics testing ground in Tsimlyansk under different stratification conditions. The measurements were carried out using the technique first used in the Tsimlyansk expedition in 2012 [1]. The measuring facility consisted of four three-component acoustic Gill Windmaster anemometers–thermometers placed at the vertices of a rectangular tetrahedron with a base scale of 0.7 m (in contrast to the experiment in 2012, when the base scale was 5 m). The measuring facility was placed on top of a mast with an adjustable height of 3.5, 5, 13.5, and 25 m and was equipped with a rotator. The temperature profile in the 10–600 m layer was continuously recorded by the Kadygrov microwave profiler [2]. The series of density of instantaneous helicity He = u_i'ω'_i = u₁'ω₁' + u₂'ω'₂ + u₃'ω'₃ and average values of the total and its summands were calculated for 12 daytime and 10 daytime 2-hour intervals. The helicity value averaged over 12 day realizations is about 0.2 m/s², and the average cosine is close to 0.08 ± 0.03. At night, the helicity is estimated as 0.07 ± 0.03 m/s², and the cosine is close to 0.025 ± 0.03. For the abovementioned 12 daytime and 10 daytime 2-hour intervals, the covariance and correlation matrices of temperature components, velocity rotor, velocity, and temperature gradient are calculated. The off-diagonal terms of the covariance matrix exceed by absolute values the diagonal terms several times. Similar characteristics of a potential vortex were estimated in the incompressibility approximation. The systematic error due to spatial averaging of the measured quantities is discussed.     

Evidence of microphytobenthic roles on coastal shallow water of the Seto Inland Sea,Japan

This study describes the temporal variation of microphytobenthic biomass and its controlling factors, as well as the impact of microphytobenthic activities on coastal shallow sediment in the eastern Seto Inland Sea, Japan. The chlorophyll a (Chl a), phaeopigments and sedimentary biophilic element (C, N, P and Si) contents in surface sediments, as well as nutrient concentrations at the sediment-water interface (overlying water and pore water) were observed monthly during November 2003 to May 2005 at one site in Shido Bay (water depth ca. 7 m) and at one site in Harima-Nada (35 m). No correlation was observed between the sedimentary biophilic elements and other parameters. The maximum chlorophyll a content of 48.2 μg g^–1 was found in surface sediments under the photon flux reaching the seafloor of 537 μmol photon m^–2 s^–1 during the winter period when water transparency was the highest at station S (Shido Bay). Our results suggest that higher chlorophyll a content in surface sediment was due to the fresh microphytobenthic biomass (mainly benthic diatom). We also found a significant negative correlation between Chl a and inorganic nutrients in pore water during the low temperature period, especially silicic acid. This result suggests that the silicic acid was assimilated largely during the increase of microphytobenthic biomass in surface sediment. This study suggests that the microphytobenthic roles may have a great effect on nutrient budgets during the large supply of irradiance (winter periods) for their photosynthetic growth in shallow coastal waters.     

Meiobenthos under mariculture conditions of the brown seaweed <Emphasis Type="Italic">Saccharina japonica</Emphasis> in Rifovaya Bay,Peter the Great Gulf,Sea of Japan

The qualitative and quantitative compositions of free-living marine nematodes were studied in Rifovaya Bay (Peter the Great Gulf, Sea of Japan). It was found that the density distribution of nematode populations in bottom sediments of Rifovaya Bay is nonuniform. In total, 72 nematode species were found, including Oncholaimium paraolium, Viscosia epapilosa, and Monoposthia latiannulata on all types of substrates. The dominant trophic group comprised “predators” (2B) and “scrapers” (2A). It was shown that the species composition of nematodes in Rifovaya Bay is very similar to the species composition of nematodes in the other areas of Peter the Great Gulf.     

Horizontal and vertical distributions of <Superscript>137</Superscript>Cs in seabed sediments around the river mouth near Fukushima Daiichi Nuclear Power Plant

Investigations including a bathymetric survey, sonic prospecting, and vibrocoring were performed to understand the horizontal and vertical distribution of ¹³⁷Cs in seabed sediments in shallow seas with depths less than 30 m near the Fukushima Daiichi Nuclear Power Plant. Especially, features of ¹³⁷Cs distributions in deeper sections of the seabed sediments were studied to evaluate the vertical heterogeneity of ¹³⁷Cs distribution in the seabed sediments in shallow seas. The distribution area of the seabed sediments was less than half of the investigation area, and the locations of the seabed sediments were divided into flat and terrace-like seafloors based on their topographical features. The thicknesses of the seabed sediment layers were mostly <2 m. The ¹³⁷Cs inventories in the seabed sediments varied from 13 ± 1 to 3,510 ± 26 kBq m^?2, and continuous distributions of ¹³⁷Cs at depths greater than 81 cm were observed. The ¹³⁷Cs distributions were not uniform; however, the ¹³⁷Cs inventories tended to be larger near the base of the steeper ascending slopes than in the terrace-like seafloors themselves. Based on the relationship between the ¹³⁷Cs inventories and mean shear stress, features of the seafloor topography were inferred to be significant control factors governing the horizontal and vertical distribution of ¹³⁷Cs in the seabed sediments. Rapid changes and multiple peaks in the vertical profile of the ¹³⁷Cs distributions suggest that they are related to pulse input caused by heavy-rain events. Change in the ¹³⁷Cs inventories with depth in this study are larger than those reported in previous studies, indicating earlier results of ¹³⁷Cs inventories per unit in seabed sediments in shallow seas, especially near the river mouth, which drains a radiologically highly-contaminated basin, were underestimated.     

Tidal boundary layer measurements in the presence of waves

Measurements of tidal current and wave velocity made at 0.69 and 1.85 m above a rough seafloor exhibit large current gradients (boundary layer) in the water column. The logarithmic boundary layer flow model was fitted to the measurements, and thus roughness (z₀) and friction velocity (u^*) parameters were derived. The roughness parameter values were generally consistent with the observed upstream physical roughness. The values of both parameters for conditions in the rough turbulence flow regime are generally larger (much larger for ebb) than earlier published values for similar measurements of currents in the absence of significant waves but are comparable to values from recent measurements of currents in the presence of storm waves. The high parameter values here appear to relate more to the magnitude of the current and to the upstream physical bottom roughness than to the magnitude of the seastate. Large boundary layers in the flow at the seabed have a profound effect on the design of offshore structures such as offshore pipelines.     

Relations between Sea Surface Temperature and Air-Sea Heat Flux at Periods from 1 Day to 1 Year Observed at Ocean Buoy Stations around Japan

Relations between sea-surface temperature (T _s) and heat flux at the sea surface (F) have been investigated using data from ocean observation buoys located off Shikoku in the Sea of Japan and in the East China Sea. Wavelet transformation decomposed F and T _s to wavelet coefficients (WLC) in the period-time domain. Assuming one-dimensional heat transfer by eddy diffusion in the upper ocean, the phase difference (δθ) defined as the difference between the phase of the temporal change rate of T _s, and the phase of F ranges statistically from 0 to +π/4 when F changes T _s, and is around −π/2 when heat convergence in the sea (Av) forces T _s. The δθ values are distributed from 0 to +π/4 at one-day and one-year periods at all buoys. WLC amplitude (WLCA) of F at periods from 16 to 32 day periods, which may be caused by the atmospheric ridge-trough systems, maintains energy longer than WLCA at periods from 2 to 16 days, which may be caused by monsoonal surges. At periods from 2 to 64 days, δθ values distribute from 0 to +π/4 or around −π/2 at each event, reflecting the surroundings of each ocean, i.e., Kuroshio recirculation in the off-Shikoku area, water-temperature front in the Sea of Japan, and water exchange in the continental shelf edge in the East China Sea. We demonstrate that the wavelet analysis can characterize the correspondence between irregular signals of F and T _s in various time scales and locations. This revised version was published online in August 2006 with corrections to the Cover Date.     

The international at-sea intercomparison of fCO2 systems during the R/V Meteor Cruise 36/1 in the North Atlantic Ocean

The ‘International Intercomparison Exercise of fCO₂ Systems’ was carried out in 1996 during the R/V Meteor Cruise 36/1 from Bermuda/UK to Gran Canaria/Spain. Nine groups from six countries (Australia, Denmark, France, Germany, Japan, USA) participated in this exercise, bringing together 15 participants with seven underway fugacity of carbon dioxide (fCO₂) systems, one discrete fCO₂ system, and two underway pH systems, as well as systems for discrete measurement of total alkalinity and total dissolved inorganic carbon. Here, we compare surface seawater fCO₂ measured synchronously by all participating instruments. A common infrastructure (seawater and calibration gas supply), different quality checks (performance of calibration procedures for CO₂, temperature measurements) and a common procedure for calculation of final fCO₂ were provided to reduce the largest possible amount of controllable sources of error. The results show that under such conditions underway measurements of the fCO₂ in surface seawater and overlying air can be made to a high degree of agreement (±1 μatm) with a variety of possible equilibrator and system designs. Also, discrete fCO₂ measurements can be made in good agreement (±3 μatm) with underway fCO₂ data sets. However, even well-designed systems, which are operated without any obvious sign of malfunction, can show significant differences of the order of 10 μatm. Based on our results, no “best choice” for the type of the equilibrator nor specifics on its dimensions and flow rates of seawater and air can be made in regard to the achievable accuracy of the fCO₂ system. Measurements of equilibrator temperature do not seem to be made with the required accuracy resulting in significant errors in fCO₂ results. Calculation of fCO₂ from high-quality total dissolved inorganic carbon (C_T) and total alkalinity (A_T) measurements does not yield results comparable in accuracy and precision to fCO₂ measurements.     

Quantitative distribution of the starfish in some bays of the far east

The quantitative distribution of the starfish was analyzed in Vostok Bay, Nakhodka Bay, and Kievka Bay (the Sea of Japan). Nine species were found. Three species of Asteroidea dominated (Asterina pectinifera, Asterias amurensis, and Leptasterias sp.). The population density of A. pectinifera and A. amurensis was higher in the shallow waters compared to the sublittoral zone. The frequency of the occurrence, the biomass, and the population density of Leptasterias sp. were the highest for the shallow waters and decreased accordingly when the water temperature decline on the transect from Vostok Bay to Kievka Bay. In Kievka Bay, which was characterized by the lowest water temperatures, Leptasterias sp. dominated in the littoral and sublittoral zones followed by A. pectinifera and A. amurensis.     

Sediment distribution and transport across the continental shelf and slope under idealized wind forcing

Resuspension, transport, and deposition of sediments over the continental shelf and slope are complex processes and there is still a need to understand the underlying spatial and temporal dynamical scales. As a step towards this goal, a two-dimensional slice model (zero gradients in the alongshore direction) based on the primitive flow equations and a range of sediment classes has been developed. The circulation is forced from rest by upwelling or downwelling winds, which are spatially uniform. Results are presented for a range of wind speeds and sediment settling speeds. Upwelling flows carry fine sediments (low settling speeds) far offshore within the surface Ekman layer, and significant deposition eventually occurs beyond the shelf break. However, coarser sediments quickly settle out of the deeper onshore component of the circulation, which can lead to accumulation of bottom sediments within the coastal zone. Downwelling flows are more effective at transporting coarse sediments off the shelf. However, strong vertical mixing at the shelf break ensures that some material is also carried into the surface Ekman layer and returned onshore. The concentrations and settling fluxes of coarse sediments decrease offshore and increase with depth under both upwelling and downwelling conditions, consistent with trends observed in sediment trap data. However, finer sediments decrease with depth (upwelling) or reach a maximum around the depth of the shelf break (downwelling). It is shown that under uniform wind conditions, suspended sediment concentrations and settling fluxes decay offshore over a length scale of order τ_s/ρf|w_s|, where τ_s is the wind stress, ρ the water density, f the Coriolis parameter, and w_s is the sediment settling velocity. This scaling applies to both upwelling and downwelling conditions, provided offshore transport is dominated by wind-driven advection, rather than horizontal diffusion.     

Mixed layer depth front and subduction of low potential vorticity water in an idealized ocean GCM

It is known that there is a front-like structure at the mixed layer depth (MLD) distribution in the subtropical gyre, which is called the MLD front, and is associated with the formation region of mode water. In the present article, the generation mechanism of the MLD front is studied using an idealized ocean general circulation model with no seasonal forcing. First, it is shown that the MLD front occurs along a curve where u _g ·∇T _s = 0 is satisfied (u _g is the upper ocean geostrophic velocity vector, T _s is the sea surface temperature and ∇ is the horizontal gradient operator). In other words, the front is the boundary between the subduction region (u _g ·∇T _s > 0) and the region where subduction does not occur (u _g ·∇T _s < 0). Second, we have investigated subduction of low potential vorticity water at the MLD front, which has been pointed out by past studies. Since u _g ·∇T _s = 0 at the MLD front, the water particles do not cross the outcrop at the MLD front. The water that is subducted at the MLD front has come from the deep mixed layer region where the sea surface temperature is higher than that at the MLD front. The temperature of the water in the deep mixed layer region decreases as it is advected eastward, attains its minimum at the MLD front where u _g ·∇T _s = 0, and then subducts under the warmer surface layer. Since the deep mixed layer water subducts beneath a thin stratified surface layer, maintaining its thickness, the mixed layer depth changes abruptly at the subduction location.     

Factors influencing the magnitude of phytoplankton primary production in a high-energy surf zone

The main factors influencing phytoplankton primary production in the surf zone of the Sundays River Beach, Algoa Bay have been characterized. These factors include cell concentration, chlorophyll concentration, irradiance, temperature and salinity. Good relationships have been obtained between cell concentration, chlorophyll concentration and primary production. The P-I curves showed dependence on temperature with a linear regression between temperature and I_k values. Light saturation was shown to occur between 300 and 510 μmol m⁻² s⁻¹ at normal field temperatures. T_max and T_min were found to be 34°C and 0°C, respectively; P_max was 25°C. Salinity had a marked effect on primary production with S_max occurring at 60 ppt and an extrapolated S_min at 0 ppt. P_max was found to occur at 30 ppt.     

The spatial distribution of surface <Emphasis Type="Italic">f</Emphasis>CO<Subscript>2</Subscript> in the Southwestern East Sea/Japan sea during summer 2005

Fugacity of CO₂ (fCO₂), temperature, salinity, nutrients, and chlorophyll-a were measured in the surface waters of southwestern East Sea/Japan Sea in July 2005. Surface waters were divided into three waters based on hydrographic characteristics: the water with moderate sea surface temperature (SST) and high sea surface salinity (SSS) located east of the front (East water); the water with high SST and moderate SSS located west of the front (West water); and the water with low SST and SSS located in the middle part of the study area (Middle water). High fCO₂ larger than 420 μatm were found in the West water. In the Middle water, CO₂ was undersaturated with respect to the atmosphere, with values between 246 and 380 μatm. Moderate fCO₂ values ranging from 370 to 420 μatm were observed in the East water. For the East and West waters, estimates of temperature dependency of fCO₂ (12.6 and 15.1 μatm °C⁻¹, respectively) were rather similar to a theoretical value, indicating that SST is likely to be a major factor controlling the surface fCO₂ distribution in these two regions. In the Middle water, however, the estimated temperature dependence was somewhat lower than the theoretical value, and relatively high concentrations of surface chlorophyll-a coincided with the low surface fCO₂, implying that biological uptake may considerably affect the fCO₂ distribution. The net sea-to-air CO₂ flux of the study area was estimated to be 0.30±4.81 mmol m⁻² day⁻¹ in summer, 2005.