High-resolution seafloor mapping using the DSL-120 sonar system: Quantitative assessment of sidescan and phase-bathymetry data from the Lucky Strike segment of the Mid-Atlantic Ridge

High-resolution, side-looking sonar data collected near the seafloor (100 m altitude) provide important structural and topographic information for defining the geological history and current tectonic framework of seafloor terrains. DSL-120 kHz sonar data collected in the rift valley of the Lucky Strike segment of the Mid-Atlantic Ridge near 37° N provide the ability to quantitatively assess the effective resolution limits of both the sidescan imagery and the computed phase-bathymetry of this sonar system. While the theoretical, vertical and horizontal pixel resolutions of the DSL-120 system are <1 m, statistical analysis of DSL-120 sonar data collected from the Lucky Strike segment indicates that the effective spatial resolution of features is 1–2 m for sidescan imagery and 4 m for phase-bathymetry in the seafloor terrain of the Mid-Atlantic Ridge rift valley. Comparison of multibeam bathymetry data collected at the sea-surface with deep-tow DSL-120 bathymetry indicates that depth differences are on the order of the resolution of the multibeam system (10–30 m). Much of this residual can be accounted for by navigational mismatches and the higher resolving ability of the DSL-120 data, which has a bathymetric footprint on the seafloor that is 20 times smaller than that of hull-mounted multibeam at these seafloor depths (2000 m). Comparison of DSL-120 bathymetry with itself on crossing lines indicates that residual depth values are ±20 m, with much of that variation being accounted for by navigational errors. A DSL-120 survey conducted in 1998 on the Juan de Fuca Ridge with better navigation and less complex seafloor terrain had residual depth values half those of the Lucky Strike survey. The quality of the bathymetry data varies as a function of position within the swath, with poorer data directly beneath the tow vehicle and also towards the swath edges.Variations in sidescan amplitude observed across the rift valley and on Lucky Strike Seamount correlate well with changes in seafloor roughness caused by transitions from sedimented seafloor to bare rock outcrops. Distinct changes in sonar backscatter amplitude were also observed between areas covered with hydrothermal pavement that grade into lava flows and the collapsed surface of the lava lake in the summit depression of Lucky Strike Seamount. Small features on the seafloor, including volcanic constructional features (e.g., small cones, haystacks, fissures and collapse features) and hydrothermal vent chimneys or mounds taller than 2 m and greater than 9 m² in surface area, can easily be resolved and mapped using this system. These features at Lucky Strike have been confirmed visually using the submersible Alvin, the remotely operated vehicle Jason, and the towed optical/acoustic mapping system Argo II.     

Radiation balance and heat budget at the ocean/atmosphere interface in the western North Pacific

The downward short- and long-wave radiation fluxes at the sea surface (S, L) were measured aboard the R/VHakuho Maru, University of Tokyo, for the period of 117 days on six cruises from 1981 to 1985 in the western North Pacific near Japan. The upward fluxes of short- and long-wave radiation (S, L) were calculated by Payne's (1972) table and the Stefan-Boltzmann's law, respectively. The sensible and laten heat fluxes (Q _h ,Q _e ) were also estimated from an aerodynamic bulk method.From April to August, the daily mean value ofS varied with the amplitude of 100200 Wm^–2. The value ofS was estimated approximately 6% ofS in all seasons. The difference betweenL andL was so small that the net radiation flux (Q _n ) was dominated byS. In addition, the net heat flux at the sea surface was also dominated byS due to small values ofQ _h andQ _e , and then the ocean was warmed at the rate of 111 Wm^–2 in April and 63 Wm^–2 in August in the Oyashio Area, and 132 Wm^–2 in May and 164 Wm^–2 in June in the Kuroshio Area, respectively.From September to March, a remarkable negative correlation between the day to day variation ofS and that ofL was observed except when an intense cold air outbreak occurred. It was found that the correlation was caused by the cloud climatological feature of the western North Pacific in this period.S was not a dominant factor in the net heat flux. The value ofQ _h +Q _e in the Kuroshio Area ranged from 260 Wm^–2 to 630 Wm^–2, much larger thanQ _n which ranged from –8 Wm^–2 to 92 Wm^–2 in the leg mean values (each leg period was about 10 days). Then the ocean was cooled at the rate of –160–620 Wm^–2 during this period. The net heat flux in the Kuroshio Area averaged over five legs from late November to February was –473 Wm^–2. This value is 50100% larger than the climatological values reported so far.The temporal and spatial variability of radiation fluxes and heat fluxes during each leg was also discussed.     

Merged GLORIA sidescan and Hydrosweep pseudo-sidescan: Processing and creation of digital mosaics

We have replaced the usual band of poor-quality data in the near-nadir region of our GLORIA long-range sidescan-sonar imagery with a shaded-relief image constructed from swath bathymetry data (collected simultaneously with GLORIA) which completely cover the nadir area. We have developed a technique to enhance these pseudo-sidescan images in order to mimic the neighbouring GLORIA backscatter intensities. As a result, the enhanced images greatly facilitate the geologic interpretation of the adjacent GLORIA data, and geologic features evident in the GLORIA data may be correlated with greater confidence across track. Features interpreted from the pseudo-sidescan may be extrapolated from the near-nadir region out into the GLORIA range where they may nt have been recognized otherwise, and therefore the pseudo-sidescan can be used to ground-truth GLORIA interpretations. Creation of digital sidescan mosaics utilized an approach not previously used for GLORIA data. Pixels were correctly placed in cartographic space and the time required to complete a final mosaic was significantly reduced. Computer software for digital mapping and mosaic creation is incorporated into the newly-developed Woods Hole Image Processing System (WHIPS) which can process both low- and high-frequency sidescan, and can interchange data with the Mini Image Processing System (MIPS) most commonly used for GLORIA processing. These techniques are tested by creating digital mosaics of merged GLORIA sidescan and Hydrosweep pseudo-sidescan data from the vicinity of the Juan Fernandez microplate along the East Pacific Rise (EPR).     

Morphology of largeHalimeda bioherms,eastern Java sea (Indonesia): A side-scan sonar study

The most extensive and thickestHalimeda bioherms reported from modern seas occur along both open and protected margins of Kalukalukuang Bank, eastern Java Sea (Indonesia). Seismic profiles suggest these features average 20 to 30 m thick (maximum 50 m) and developed over large areas by coalescence of individual mounds. Bioherm morphologies range from numerous small mounds ( 10–20 m diameter) through hay-stack features ( 100 m diameter) to broad swells. Some modern bank bioherms display surface forms that appear related to onbank flow from Makassar Strait. Upwelling of cold, nutritive water is likely responsible for unusualHalimeda productivity and large bioherm development.     

The evolution of isolated vortices interacted with steep slope

Numerical solutions are examined for isolated, intense vortices as influenced by western bounding bottom topography through the use of a rigid-lid, two-layer primitive -plane numerical model. Systematic studies are made of the sense of rotation (cyclonic/anticyclonic), the consequence of varying the gradient of bottom slope, and the different vertical shear in a two layer ocean. In the basin with a bottom slope, the nearly barotropic anticyclonic vortex forms a modon-like vortex for S with fixedRo ₂<O(1) (where is the ratio between the variation of the Coriolis parameter across the eddy to the Coriolis parameter in the center, S the topographic effect and,Ro ₂ the Rossby number in the lower layer) and its generation is due to a compound effect of the planetary beta, topographic beta, avvection, and mirror image. The formation of the modon-like vortex and the propagation of the original vortex onto the bottom slope depends on the strength of slope gradient and the baroclinicity of the vortex. The nearly barotropic anticyclonic vortex evolves into the stronger upper ocean one with increasing S: the gradient of the bottom slope becomes steeper. Then the original vortex lives longer because the barotropic component of the energy is converted to the baroclinic one and it moves toward southeast in forming a modon-like vortex in the lower layer. The evolution of a vortex in the model results are compared to observational results of a Kuroshio warm core ring (KWCR) obtained from hydrographic data (June, 1985) and from NOAA satellite infrared images (April, 1985 to July, 1985). It is shown that a KWCR (June, 1985) is influenced by the western continental slope/shelf of the East Japan.     

Morphology of a ‘superfast’ Mid-Ocean Ridge crest and flanks: The East Pacific Rise, 7°–9° S

Detailed bathymetric data from a Hydrosweep multibeam sonar survey of a 250 km-long portion of the superfast-spreading southern East Pacific Rise crest and flanks show that the along-axis variation in morphology and axial depth differs significantly from that observed at the fast-spreading northern East Pacific Rise. While the deep mantle upwelling pattern is similar under the northern and southern East Pacific Rise, our observations require that the connectivity of the shallow, subcrestal plumbing system be more efficient beneath the super-fast spreading southern East Pacific Rise than beneath the slower spreading northern East Pacific Rise.     

Global distribution and atmospheric transport of chlorinated hydrocarbons: HCH (BHC) isomers and DDT compounds in the Western Pacific,Eastern Indian and Antarctic Oceans

Concentrations of chlorinated hydrocarbons such as HCH isomers and DDT compounds were determined in air and surface water samples taken from the Western Pacific, Eastern Indian and Antarctic Oceans. The most interesting finding was their presence in measurable concentrations in the Antarctic Ocean. Chlorinated hydrocarbon pesticides are widely distributed in the open ocean environment over both the northern and southern hemispheres, and some characteristic distribution patterns of pesticide species in different oceanic regions were observed both in air and water samples. HCH residues found in the northern hemisphere were much higher in concentration than those in the southern hemisphere. On the other hand, higher concentrations of DDT residues were found in the tropical regions, but their levels were not so different between both the northern and southern hemispheres. HCH isomers found in the northern hemisphere had the following order of concentrations:-HCH> HCH>-HCH, while in the southern hemisphere-HCH was apparently dominant. DDT compound compositions were rather uniform in all the oceans surveyed, and more than 50% wasp,p-DDT. These facts can be explained by the world wide situation regarding pesticide use and the physicochemical properties of the pesticides such as their vapor pressures and water solubilities. In addition, the meridional circulation of the atmosphere, particularly the mass flows of the Hadley and Ferrel cells in the troposphere, also contributes to the atmospheric transport and global distribution of these pesticides.     

Dissolved Folin Phenol Active Substances (Tannin and Lignin) in the Seawater along the West Coast of India

Information on the distribution of dissolved Folin phenol active substances (FPAS) such as tannin and lignin in the seawater along the west coast of India is provided. Notable amounts of FPAS (surface concentrations: 80 g/l to 147 g/l and bottom concentrations: 80 g/l to 116 g/l) were detected in the seawater along the coast. The distribution pattern brings about a general depth-wise decrease. A seaward decrease was observed in the southern stations whereas reverse was the case in northern stations. A significant negative correlation was observed between FPAS concentration and dissolved oxygen in sub-surface samples. The appreciable amounts of FPAS detected in the coastal waters indicate the presence of organic matter principally originating from terrestrial (upland and coastal marsh) ecosystems in the marine environment. In this context, they may be used as tracers to determine the fate of coastalborn dissolved organic matter in the ocean and to determine directly the relationship between allochthonous and autochthonous organic matter.     

The mechanism of the development of wind-wave spectra

The mechanism of the development of wind-waves will be proposed on the basis of the observed wave spectra in the wind tunnels and at Lake Biwa (Imasato, 1976). It consists of two aspects: One is that the air flow over the wind-waves transfers momentum concentratively to the steepest component waves and the other is that the upper limit of the growth of a wave spectral density is given by the ultimate value in the slope spectral density. The first aspect means that the wave field has the momentum transfer filter on receiving the momentum from the air flow. Wind-waves in the stage of sea-waves receive the necessary amount of momentum by the form drag,e.g. according to the Miles' (1960) inviscid mechanism, through a very narrow frequency region around a dominant spectral peak. On the other hand, wind-waves in the stage of initial-wavelets receive it according to the Miles' (1962a) viscous model through a fairly broad frequency region around the peak. The upper limit ofS _max developing according to viscous mechanism is given byS _max =6.40×10^–4 k _max ^–2cm²s andS _max =2.03C(f _max )^–2cm²s(S _max is the power density of the wave spectral peak with the frequencyf _max ,k _max is the wave number corresponding to the frequencyf _max andC is the phase velocity).From the second aspect, the upper limit of the growth of wave spectral density is given by 33.3f ^–4cm²s in the frequency region of late stage of sea-waves. Therefore, the spectral peak, which has the largest value in the slope spectral density in the component waves of the wave spectrum, rises high over the line 4.15f ^–5cm²s. The energy is transported from the spectral peak to the high frequency part and to the forward face of a wave spectrum by nonlinear wave-wave interaction. This nonlinearity is confirmed by the bispectra calculated from the observed wind-wave data. In the stage of sea-waves, nonlinear rearrangement of the wave energy comes from a narrow momentum transfer filter, and, in the stage of initial-wavelets, it comes mainly from small corrugations and small steepness of the wave field.     

Turbulent structure in water under laboratory wind waves

The structure of the turbulent boundary layer underneath laboratory wind waves was studied by using a combination of a high-sensitivity thermometer array with a two-component sonic flowmeter. The temperature fluctuations are used to detect movements of water parcels, with temperature as a passive quantity. The turbulence energy was dominant in the frequency range (0.01 0.1 Hz), which was much smaller than the wind-wave frequency (2 5 Hz), and in which the turbulence was anisotropic. There was a frequency range (0.2 2 Hz for velocity, 0.2 5 Hz for temperature fluctuation) where the turbulence was isotropic and had a –5/3 slope in the energy spectrum. These points are the same as those in previous works. However, by analyses of the time series by using a variable-interval time-averaging technique (VITA), it has been found that conspicuous events in this main turbulence energy band are the downward bursting from the vicinity of the water surface. Thus the structure of the water layer underneath the wind waves has characters which are similar to the familiar turbulent boundary layer over a rough solid wall, as already conceived. It has been found that, at the same time, the turbulence energy can be related to quantities of the wind waves (the root mean squared water level fluctuation and the wave peak frequency), for different wind and wave conditions. That is, the turbulence underneath the wind waves develops under a close coupling with the wind waves.     

Some characteristics of the development process of the wind-wave spectrum

The development process of wind-waves of which spectral peak distributes from 0.6 cps to 9.3 cps will be discussed on the basis of the wind tunnel experiments and of the field observations performed at Lake Biwa. The characteristics of power and slope spectra are here presented. The development process of these wind-waves is characterized by three stages;i.e. initial-wavelets, transition stage and sea-waves. In the wind tunnel experiments, the transition from the stage of the initial-wavelets to the transition stage occurs when the wave spectral peak arrives at the line 6.40×10^–4 k ^–2cm²·sec (wherek is wave number) or when the slope spectral density at the frequencyf _max becomes larger than 6.40×10^–4 sec. In the stage of sea-waves, the component wave of a wave-spectral peak is steepest in the component waves. And the wave spectral peak develops along the line 1.02×10² f ^–6 cm²·sec (wheref is the frequency corresponding to the wave numberk) untill it reaches the line 33.3f ^–4cm²·sec, and thereafter develops along the latter line, which indicates the constant density of slope spectrum. It is suggested that the nonlinearity of wind-waves must become stronger as wind-waves develop. The effective momentum flux _ws from the air flow to wind-waves in this stage is evaluated to be about 49% of the total stress ₀.     

Seamount gravity anomaly modelling with variably thick sediment cover

Inversion modelling of marine gravity anomalies to derive predicted seafloor topography has provided significant advance in delineating deep-ocean bathymetry where the seafloor both conforms to the half-space cooling model of seafloor spreading, and largely sediment-free. Similar modelling for elevated ridges and seamounts, that are formed by processes other than seafloor spreading and/or have proximal sediment sources (e.g., continental margins and volcanic arcs), have significantly higher errors when validated against modern shipborne echo-sounding data. A three-dimensional, five-layer gravity model is emulated for the cases of both synthetic and real seamounts, with varying degrees of sediment burial, to establish the sensitivity of variable sediment cover as a source of error. A simple `Gaussian' seamount with base radius of 30 km, 2000 m of relief, has a maximum 140–160 mGal anomaly, that decreases to 50 mGal with the addition of 1 km of sediment cover with simple `flood' geometry. Complete burial, with a typical sediment density of 2300 kg m^–3, results in a 120 mGal difference from a sediment-free seamount model. Increasing sediment density results in an exponential decay of the seamount anomaly. More complex synthetic geometries of varying basement relief and sediment thickness show that the anomaly amplitude remains significant, especially where the latter is >700–800 m thick. For the real case, seamounts of the Three Kings Ridge (northern New Zealand) imaged with seismic reflection data, with varying degrees of sediment cover of up to 1 km, when modelled both with and with-out the inclusion of a sediment layer, typically have rms differences of 30 mGal between observed and modelled gravity anomalies. Significantly, the rms errors are reduced by 50% with the inclusion of a sediment layer that corresponds to a reduction of predicted seafloor topography rms errors of 192–684 m to 78–360 m.     

Intrusive Stratification of the Gulf-Stream Frontal Zone According to the Results of Investigations Performed during Cruise 43 of the R/V Akademik Vernadskii

We analyze the data of investigation of the intrusive structure of the Gulf-Stream frontal zone obtained in making frequent drift sections with the help of an MGI-8102 probing complex, study the regularities of variation of temperature, salinity, and density along separate intrusions, and present a series of results connected with the specific features of initiation and development of intrusions and the types of exchange processes determining their transformations. It is shown that the T-diagrams of all intrusions are well separated into segments with different slopes. Moreover, by comparing the slopes and locations of these segments with each other and with the T-diagrams of the Gulf Stream and slope waters, we can fairly reliably attribute the corresponding segments of intrusions to one of the following four types: initial dynamic folds of the frontal zone, layers of domestic water separating the intrusive segments of foreign water from each other, intrusive segments characterized by the penetration of ambient water, and segments not intrusive initially but getting the required slope as a result of interaction with upper and lower intrusive segments. For segments of the last two types, it is possible to specify the predominant type of exchange.     

Variation of water system distribution in the Sanriku Coastal Area

We analyzed hydrographic data obtained monthly by the Iwate Prefectural Fisheries Laboratory during 1977–1981. Our attention was focussed on the classification of waters based on the Temperature-Salinity (T-S) scatter diagram analyses, and on the extraction of the mean state of their distribution for each season. The water system was defined as all waters occupying some specific region on the T-S plane, and the mean state for each season was expressed by the distribution of the mode water systems (i.e., the most commonly occuring water system).Because of the coexistence of the Tsugaru Warm Current, the Oyashio and the Kuroshio, and their large seasonal variability, the sea waters had a wide variety of T-S combinations. The detailed T-S scatter diagram analyses allowed us to classify the waters into six water systems. A time-longitude diagram of water systems and the distribution of mode water systems are presented. The latter showed clear seasonal variations of the Tsugaru Warm Current and the Oyashio.Though the water system analysis was successfully able to present the distribution of the various waters in the Sanriku Coastal Area and to clearly extract the mean state of sea conditions, it was found that the classification of water systems on the T-S plane was not complete in all cases for this area.     

Spatial Patterns in Nutrient and In Vivo Fluorescence Distributions in the Marginal Ice Zone and the Seasonally Open Oceanic Zone in the Indian Sector of the Antarctic Ocean, in Austral Summer

Surface temperature, salinity, concentrations of silicate (Si) and nitrate + nitrite (N), and in vivo fluorescence (Fluor) were investigated in the marginal ice zone (MIZ) and the seasonally open oceanic zone (SOOZ) (32–40°E, 64–69°S) from February 23 to 28 1992. In the MIZ the mean Si and N were 67.8 ± 2.2 M and 32.5 ± 1.7 M, respectively. There was a trend that low N values coincided with high Fluor values. Observation conducted at one point (64°S, 38°E) revealed a diel variation pattern in Fluor. Applying this pattern of deviation from noon value, all Fluor data were normalized to value at local noon. In the MIZ a significant negative correlation was observed between the normalized Fluor and N but not Si. On the other hand, Si decreased continuously from south to north in the SOOZ and was negatively correlated with the normalized Fluor. Difference in Si concentration was about 30 M between the sea around 64°S and the MIZ, while the difference in N concentration was estimated as less than 10 M. If diatoms take up silicate and nitrogen at an approximate ratio of 1:1, additional nitrogenous nutrients other than nitrate and nitrite (e.g. ammonia, urea etc.) would be required. In this case, an f-ratio of lower than 33% is obtained. It is suggested that in the MIZ abundance of phytoplankton community dominated by non-diatom increases utilizing nitrate while in the SOOZ abundance of phytoplankton community dominated by diatoms increases consuming Si and regenerated nitrogen.     

Chlorinated hydrocarbons in the North Pacific and Indian Oceans

PCBs, DDT compounds and HCH isomers were detected in the air and surface waters of the North Pacific and Indian Oceans, including the Bering Sea, East China Sea, South China Sea, Bay of Bengal and the Arabian Sea. The general concentrations of each chlorinated hydrocarbon were as follows: water PCBs 0.1 to 1.0, DDT 0.01 to 1.0, HCH 1.0 to 10 ngl ^–1; air DDT 0.01 to 1.0, HCH 0.1 to 10 ng m^–3. PCB concentrations in surface waters were slightly lower than those of the North Atlantic and North Sea previously reported, while DDT concentrations in the air and water were higher. Remarkably high concentrations of DDT and HCH were found in the air off the western coast of India. Also in the Pacific site off Central America, a fairly high concentration of DDT was observed in an air sample. These data suggest that large amounts of DDT and HCH are being used in the tropical zone, especially in southern Asia. Furthermore, high concentrations were observed both in the air and water of the Northwest Pacific between 30°N and 40°N latitude. There is a possibility that both pesticides are not only still being used in lower latitude countries but also in the mid-latitude ones of the Asian continent excluding Japan. In addition to this atmospheric circulation may also contribute to the concentration of these pesticides in the mid-latitudinal zone.     

Geophysical and geochemical observations on actively seeping hydrocarbon gases on the south-eastern Yellow Sea continental shelf

In the southeastern Yellow Sea, active seepage of hydrocarbon gases has been observed by high-resolution (3.5 kHz) seismic profiling both in 1987 and 2001, occurring through a large number of plumes from the topmost pre-Holocene sedimentary layer. It is strong enough to compensate for current speed, extending vertically up to the sea surface. The gas seepage often appears to be explosive to form craters and diapirs, although pockmarks are rare due to the redistribution of mobile palimpsest sands. In core-top seawater and sediments, the gases are characterized by high amounts of C₂, homogenous ¹³C₁ values and a large difference (19.7 on average) between ¹³C₁ (–55.2 to –53.6 PDB) and ¹³C₂ (–36.8 to –32.5 PDB) values. The gases are considered to be generated with a smaller amount of C₁ at the early thermal cracking stage of labile source materials, after which the C₂ gas is enriched in ¹³C by diffusion or biological alternation at the generation or accumulation site. The homogenous ¹³C₁ values may be one of the geochemical characteristics of gases acquired at depth which are less altered in the case of rapid diffusive gas migration to the seafloor.     

A stratigraphic and tectonofacies framework of the “calcari aLucina” in the Apennine Chain,Italy

The Apennine Chain provides the first example of stratigraphic (time) and synsedimentary tectonic (space) distribution of the calcari aLucina Miocene equivalents of modern cold-vent carbonates. Chemosynthetic faunal assemblages and related carbonate deposits are found at different stratigraphic levels, with peaks during Langhian-Serravallian and late Tortonian-early Messinian times. A general increase in frequency and volume occurs with time. A genetic link between venting and the Messinian Evaporite event is difficult to demonstrate. However,Lucina limestones are limited to preevaporitic times, and their maximum abundance is reached just before the onset of the Messinian Evaporite accumulation.Lucina limestones occur in almost all tectofacies of the orogen, from backland to foreland.     

Petrography and stable isotope aspects of cold-vent activity imprinted on Miocene-age “calcari aLucina” from Tuscan and Romagna Apennines,Italy

Over 20 occurrences of discontinuous limestone blocks, locally called calcari aLucina, were mapped in the Tuscan—Romagna region of the northern Italian Apennines. The limestones, consisting of a variable mixture of authigenic carbonates (calcite, dolomite, and aragonite), sulfides (primarily pyrite), and allogenic silicates, occur in association with turbidite and hemipelagite units that were deposited in foredeep basins during early to late Miocene times. The limestone blocks are interpreted to represent relicts of carbonate buildups formed around methane-rich fluid vents on the basis of their (1) striking petrographic similarities to carbonates from cold vents in the modern oceans; (2) unique chemosynthetic-like fauna, and (3) anomalously negative ¹³C values ( ¹³C = – 16 to – 58 PDB). The contemporaneous tectonism of the Apennine orogeny is likely to be the primary cause for the expulsion of the methane-rich fluids to the seabed in a manner analogous to the fluid-flow processes occurring at modern accretionary prisms.     

Volume transport of the western boundary current penetrating into a marginal sea

The amount of penetration of a western boundary current into a marginal sea which is connected to an open ocean by two narrow straits is estimated from a linear, steady and barotropic theoretical model. In this model the western boundary current in the open ocean is driven by a wind stress imposed at the sea surface. The inflow of the water of the open ocean into the marginal sea is caused by the pressure difference between two straits produced by the wind-driven circulation in the open ocean.Main external parameters are combined into two non-dimensional parameters; and (the ratio of the depth of the marginal sea to that of the open ocean), whereb is the distance between north and south boundaries of the ocean,D ₀ is the depth of the open ocean, is the latitudinal variation of the Coriolis parameter andR is the coefficient of friction. The friction is assumed to be proportional to the flow velocity.In the limit of infinite the volume transport into the marginal sea is not affected by the width of two straits and . It is mainly controlled by the wind stress and the positions of two straits. For finite values of , however, the volume transport depends considerably on and the width of the straits.Guided by both this model and physical considerations, we obtained a relation between the volume transport into the marginal sea and the external parameters. This relation predicts that about 2 % of the volume transport of the Kuroshio penetrates into the Japan Sea.