Compositional changes in natural gas bubble plumes: observations from the Coal Oil Point marine hydrocarbon seep field

Detailed measurements of bubble composition, dissolved gas concentrations, and plume dynamics were conducted during a 9-month period at a very intense, shallow (22-m water depth) marine hydrocarbon seep in the Santa Barbara Channel, California. Methane, carbon dioxide, and heavier hydrocarbons were lost from rising seep bubbles, while nitrogen and oxygen were gained. Within the rising seawater bubble plume, dissolved methane concentrations were more than 4 orders of magnitude greater than atmospheric equilibrium concentrations. Strong upwelling flows were observed and bubble-rise times were ~40 s, demonstrating the rapid exchange of gases within the bubble plume.     

Time series video analysis of bubble release processes at natural hydrocarbon seeps in the Northern Gulf of Mexico

This research quantifies the rate and volume of oil and gas released from two natural seep sites in the Gulf of Mexico: lease blocks GC600 (1200 m depth) and MC118 (850 m depth). Our objectives were to determine variability in release rates and bubble size at five individual vents and to investigate the effects of tidal fluctuations on bubble release. Observations with autonomous video cameras captured the formation of individual bubbles as they were released through partially exposed deposits of gas hydrate. Image processing techniques determined bubble type (oily, gaseous, and mixed: oily and gaseous), size distribution, release rate, and temporal variations (observation intervals ranged from 3 h to 26 d). A semi-automatic bubble counting algorithm was developed to analyze bubble count and release rates from video data. This method is suitable for discrete vents with small bubble streams commonly seen at seeps and is adaptable to multiple in situ set-ups. Two vents at GC600 (Birthday Candles 1 and Birthday Candles 2) were analyzed. They released oily bubbles with an average diameter of 5.0 mm at a rate of 4.7 bubbles s⁻¹, and 1.3 bubbles s⁻¹, respectively. Approximately 1 km away, within the GC600 seep site, two more vents (Mega Plume 1 and Mega Plume 2) were analyzed. These vents released a mixture of oily and gaseous bubbles with an average diameter of 3.9 mm at a rate of 49 bubbles s⁻¹, and 81 bubbles s⁻¹, respectively. The fifth vent at MC118 (Rudyville) released gaseous bubbles with an average diameter of 3.0 mm at a rate of 127 bubbles s⁻¹. Pressure records at Mega Plume and Rudyville showed a diurnal tidal cycle (24.5 h). Rudyville was the only vent that demonstrated any positive correlation (ρ = 0.60) to the 24.5 h diurnal tidal cycle. However, these observations were not conclusive regarding tidal effects on bubble release.     

基于单波束测深资料的海底冷泉动力学特征反演方法研究

为了准确理解冷泉水体中甲烷气体的分布规律,综合运用单波束测深数据和冷泉水体流场数据,建立了冷泉气泡上升、溶解速率的定量反演方法,利用现场海试资料对反演方法进行了验证,并对海试区域的冷泉气体溢出、溶解通量以及冷泉水体的甲烷浓度进行了估算。计算结果表明,考察船航向与冷泉水体流向的差异会对声学探测结果产生影响,当航向与流向的夹角大于±60°时,声学成像中冷泉倾角的误差将超过50%。同时,冷泉气泡上升速率的衰减与溢出口水深显著相关,相关系数可达0.9,并且冷泉上升流对上升速率的提升效果明显。水合物稳定带的分布对冷泉气泡的收缩速率影响显著,稳定带内、外收缩速率的差异可达3～4倍。冷泉气体通量的计算结果表明,调查区域内冷泉的溢出强度整体较大,同时溶解通量与水深之间呈现明显的分段效应。根据冷泉气体溶解通量估算的冷泉水体甲烷浓度与色谱分析结果具有较好的一致性。相关研究有助于实现对冷泉水体中甲烷气体分布的定量评估,为潜在海底冷泉区的圈划和海域天然气水合物的调查提供技术支撑。     

Megaplume bubble process visualization by 3D multibeam sonar mapping

MultiBeam echosounder data were collected during a surface-ship survey of the 22/4b well site in the North Sea in September 2011 using a Teledyne-Reson 7125. Modern multibeam echosounders are instrumental in providing detection and accurate localization of weak to strong bubble plumes. Two survey profiles effectively insonified the bubble plumes rising from the main crater at the well site, providing snapshot data of bubble plume processes. Additionally, three profiles insonified bubble plumes rising from, in, and to the south of a secondary crater, 1.2 km southeast of the main crater. Data processing included a simple algorithm that muted mislocated echoes from incomplete sidelobe suppression. The data processing produced a Cartesian volume of echo intensity from the water column and seabed.Plume geometry was analyzed to investigate a number of important large-scale plume processes, including plume bubble detrainment due to currents and stratification, downwelling flows, sea surface interaction, plume heterogeneity, and other fluid transport processes. The data showed strong upwelling flows, with bubble vertical motions generally much faster than currents. One important finding was that megaplumes create intrusions above the general thermocline, in part because their extensive upwelling flow lifts the thermocline locally. As a result, the intrusion layer deposits dissolved gases in the upper wave-mixed layer of the water column where it is not isolated from the atmosphere, unlike dissolved gases in the lower water column.The analysis shows that high fidelity multibeam echosounder data can provide a wealth of remote sensing information on bubble plume characteristics and processes, with important applications, including blowout monitoring and response, better understanding of megaplumes such as used in lake destratification, and improved characterization of natural seep emission processes.     

冷泉活动区气泡羽状流数值模拟的进一步研究

以南海气泡羽状流赋存特征为参考,本文继续研究冷泉活动区气泡羽状流的地震响应。为使羽状流建模方案更合理,根据水中气泡在上升过程时其半径的变化情况,修改原建模方案,重新建立了羽状流水体模型。通过正演模拟得到了羽状流炮集地震记录,单炮记录上显示出明显的散射波场特征和模型的周期特点。通过叠前深度偏移处理炮集地震记录,得到边界收敛效果较好的成像剖面,且成像精度较高。以上研究又一次奠定了羽状流地震响应进一步研究的基础,也逐步探索出更适合羽状流地震资料的处理方法,为识别羽状流提供了理论指导。     

Compositional variability and air-sea flux of ethane and propane in the plume of a large, marine seep field near Coal Oil Point, CA

Large quantities of methane (C1), ethane (C2), and propane (C3) emanate from shallow marine seeps near Coal Oil Point (COP), California. Concentrations of these gases were analyzed in the surface water down-current of the seep field over a 15-month period. The variable proportions of C1, C2, and C3 analyzed in gas bubbles emitted from 16 distinct seeps in the COP field encompass much of the variability found in the surface waters down-current. However, waters with disproportionate levels of C1 suggest the presence of additional C1 sources. Based on three spatial surveys, covering areas up to 280 km², C2 and C3 air-sea fluxes were estimated to be in the order of 3.7 and 1.4 μmol day^?1 m^?2, respectively. Only 0.6% of C2 and 0.5% of C3 in the dissolved plume originating from the COP seep field are transferred to the atmosphere in the study area, with the fate of the remainder uncertain.     

Forward modeling and inversion of the relation model between the gas content of plume and its seismic attribute

The methane bubble plume attracts interest because it offers direct evidence of seafloor gas leakage and plays an indirect role in the exploration and identification of natural gas hydrate. In this study, based on established plume models and their migration sections, three amplitude-class attributes were extracted from three formations for the migration sections of five plumes, and the correlation between the gas content and seismic attribute was obtained. As the gas content increases, the amplitude attribute correspondingly increases, and the linear correlation is relatively good. Moreover, correlation coefficients between gas content and amplitude attributes are close to 1.0. By using linear fitting, the relation model between the gas content of the plume and the seismic attribute was obtained. The relation model was subsequently used to invert the gas content from a real data-bearing plume. Comparison of the gas content section of the plume with the attribute section and real seismic section reveals common distribution characteristics, namely, the color of the section in the lower right corner is dark. If the amplitude value is large in the seismic section of the real plume, the amplitude attribute value is also large in the corresponding attribute section, and the inverted value of the gas content is also large (because gas content and amplitude are linearly correlated), which indicates that the plume bubbles of the section in the lower right corner is intensively distributed. Finally, the obtained gas content section of the plume can reflect the distribution of the plume bubble content more simply and intuitively, from which the distribution law of seafloor bubbles can be deduced, and this lays a foundation for the further estimation of the gas content of the plume and hydrate reserves.     

Importance of seep primary production to Lophelia pertusa and associated fauna in the Gulf of Mexico

To investigate the importance of seep primary production to the nutrition of Lophelia pertusa and associated communities and examine local trophic interactions, we analyzed stable carbon, nitrogen, and sulfur compositions in seven quantitative L. pertusa community collections. A significant seep signature was only detected in one of the 35 species tested (Provanna sculpta, a common seep gastropod) despite the presence of seep fauna at the three sample sites. A potential predator of L. pertusa was identified (Coralliophila sp.), and a variety of other trophic interactions among the fauna occupying the coral framework were suggested by the data, including the galatheid crab Munidopsis sp. 2 feeding upon hydroids and the polychaete Eunice sp. feeding upon the sabellid polychaete Euratella sp. Stable carbon abundances were also determined for different sections of L. pertusa skeleton representing different stages in the growth and life of the aggregation. There was no temporal trend detected in the skeleton isotope values, suggesting that L. pertusa settles in these areas only after seepage has largely subsided. Isotope values of individual taxa that were collected from both L. pertusa and vestimentiferan habitats showed decreasing reliance upon seep primary production with average age of the vestimentiferan aggregation, and finally, no seep signature was detected in the coral collections. Together our data suggest that it is the presence of authigenic carbonate substrata, a product of past seep microbial activity, as well as hydrodynamic processes that drive L. pertusa occurrence at seep sites in the Gulf of Mexico, not nutritional dependence upon primary production by seep microbes.     

Geologic control of natural marine hydrocarbon seep emissions, Coal Oil Point seep field, California

High-resolution sonar surveys, and a detailed subsurface model constructed from 3D seismic and well data allowed investigation of the relationship between the subsurface geology and gas-phase (methane) seepage for the Coal Oil Point (COP) seep field, one of the world’s largest and best-studied marine oil and gas seep fields, located over a producing hydrocarbon reservoir near Santa Barbara, California. In general, the relationship between terrestrial gas seepage, migration pathways, and hydrocarbon reservoirs has been difficult to assess, in part because the detection and mapping of gas seepage is problematic. For marine seepage, sonar surveys are an effective tool for mapping seep gas bubbles, and thus spatial distributions. Seepage in the COP seep field occurs in an east–west-trending zone about 3–4 km offshore, and in another zone about 1–2 km from shore. The farthest offshore seeps are mostly located near the crest of a major fold, and also along the trend of major faults. Significantly, because faults observed to cut the fold do not account for all the observed seepage, seepage must occur through fracture and joint systems that are difficult to detect, including intersecting faults and fault damage zones. Inshore seeps are concentrated within the hanging wall of a major reverse fault. The subsurface model lacks the resolution to identify specific structural sources in that area. Although to first order the spatial distribution of seeps generally is related to the major structures, other factors must also control their distribution. The region is known to be critically stressed, which would enhance hydraulic conductivity of favorably oriented faults, joints, and bedding planes. We propose that this process explains much of the remaining spatial distribution.     

Variability of gas composition and flux intensity in natural marine hydrocarbon seeps

The relationship between surface bubble composition and gas flux to the atmosphere was examined at five large seeps from the Coal Oil Point seep field (Santa Barbara Channel, CA, USA). The field research was conducted using a flux buoy designed to simultaneously measure the surface bubbling gas flux and the buoy’s position with differential GPS, and to collect gas samples. Results show that the flux from the five seeps surveyed a total of 11 times ranged from 800–5,500 m³ day^?1. The spatial distribution of flux from the five seeps was well described by two lognormal distributions fitted to two flux ranges. The seafloor and sea surface composition of bubbles differed, with the seafloor bubbles containing significantly more CO₂ (3–25%) and less air (N₂ and O₂). At the sea surface, the mole fraction of N₂ correlated directly with O₂ (R ² = 0.95) and inversely with CH₄ (R ² = 0.97); the CO₂ content was reduced to the detection limit (<0.1%). These data demonstrate that the bubble composition is modified by gas exchange during ascent: dissolved air enters, and CO₂ and hydrocarbon gases leave the bubbles. The mean surface composition at the five seeps varied with water depth and gas flux, with more CH₄ and higher CH₄/N₂ ratios found in shallower seeps with higher flux. It is suggested that the CH₄/N₂ ratio is a good proxy for total or integrated gas loss from the rising bubbles, although additional study is needed before this ratio can be used quantitatively.     

Void fraction measurements and scale effects in breaking waves in freshwater and seawater

This paper follows from the work of Blenkinsopp and Chaplin (2007) and describes detailed measurements of the time-varying distribution of void fractions generated by breaking waves in freshwater, artificial seawater and natural seawater under laboratory conditions, along with flow visualisation of the entrainment process. The measurements were made with highly sensitive optical fibre phase detection probes and the results demonstrate that although an additional population of fine (d < 0.3 mm) bubbles existed in the seawater cases, the total volume and distribution of entrained air, and the spatial and temporal evolution of the bubble plumes were very similar in all three water types. The influence of water type may be relatively insignificant, but a numerical bubble tracking model shows that the effect of scale is an important consideration when modelling the post-entrainment evolution of breaker-entrained bubble plumes. Consequently the results suggest that while the use of freshwater in laboratory models of oceanic processes can be considered valid in most situations, the effect of scale may impact interpretation of the results.     

A microbiological and biogeochemical investigation of the cold seep tubeworm Escarpia southwardae (Annelida: Siboglinidae): Symbiosis and trace element composition of the tube

Tubeworms within the annelid family Siboglinidae rely on sulfur-oxidizing autotrophic bacterial symbionts for their nutrition, and are among the dominant metazoans occurring at deep-sea hydrocarbon seeps. Contrary to their relatives from hydrothermal vents, sulfide uptake for symbionts occurs within the anoxic subsurface sediment, in the posterior ‘root’ region of the animal. This study reports on an integrated microbiological and geochemical investigation of the cold seep tubeworm Escarpia southwardae collected at the Regab pockmark (Gulf of Guinea). Our aim was to further constrain the links between the animal and its symbiotic bacteria, and their environment. We show that E. southwardae harbors abundant sulfur-oxidizing bacterial symbionts in its trophosome. Symbionts are able to fix inorganic carbon using the Calvin-Benson cycle, as reported in most other Siboglinidae, but can also use the reverse Tricarboxilic Acid Cycle. Surprisingly, the observed bacteria appear to be more closely related to symbionts of Escarpia and Lamellibrachia species from very distant sites located in the Gulf of Mexico and eastern Pacific, than to symbionts of a siboglinid occurring at a nearby methane seep site, only a few hundred km away from Regab. Then, by combining scanning electron microscopy and trace element (Mn, Fe, Sr, Zr) analyses of E. southwardae tube, we also show that two distinct oxidation fronts occur along the tube. The first one, near the posterior end of the tube, corresponds to the interface between oxic bottom waters and the underlying anoxic sediment. In contrast, the second redox front is located in the most anterior part of the tube, and could result from active oxygen uptake by the plume of the tubeworm. We speculate that intense oxygen consumption in this region could create favorable conditions for sulfate reduction by specialized bacteria associated with the plume, possibly leading to an additional source of dissolved sulfide that would further enhance the productivity of bacterial symbionts.     

南海北部陆坡气泡羽状流的发现:多波束水体数据

利用2016年在南海西北部陆坡琼东南海域采集的多波束水体数据,发现了海底气体渗漏至海水中形成的羽状流。在多波束数据上, 羽状流成火焰状,直径大约为30~50 m,从1 380多米的海底延伸至大约650 m的深度,高度超过700 m。在经过羽状流的浅剖剖面上,存在显示浅层气存在的声学空白区域,并识别出断裂和裂隙区域,但在水体中并无明显的异常。这可能是由于浅剖数据的分辨率不够未能捕捉到水体异常,或者气体渗漏具有间歇性。该海域存在明显的似海底反射显示,气体渗漏可能与水合物系统之间存在复杂的相互作用。由于缺乏经过羽状流的多道地震数据,难以对羽状流的形成机制进行进一步的推测。南海北部陆坡羽状流的发现对于理解被动大陆边缘的甲烷渗漏机制、水合物的形成与分解具有重要的意义。     

Improved detection and mapping of deepwater hydrocarbon seeps: optimizing multibeam echosounder seafloor backscatter acquisition and processing techniques

Marine seep hunting surveys are a current focus of hydrocarbon exploration surveys due to recent advances in offshore geophysical surveying, geochemical sampling, and analytical technologies. Hydrocarbon seeps are ephemeral, small, discrete, and therefore difficult to sample on the deep seafloor. Multibeam echosounders are an efficient seafloor exploration tool to remotely locate and map seep features. Geophysical signatures from hydrocarbon seeps are acoustically-evident in bathymetric, seafloor backscatter, midwater backscatter datasets. Interpretation of these signatures in backscatter datasets is a fundamental component of commercial seep hunting campaigns. Degradation of backscatter datasets resulting from environmental, geometric, and system noise can interfere with the detection and delineation of seeps. We present a relative backscatter intensity normalization method and an oversampling acquisition technique that can improve the geological resolvability of hydrocarbon seeps. We use Green Canyon (GC) Block 600 in the Northern Gulf of Mexico as a seep calibration site for a Kongsberg EM302 30 kHz MBES prior to the start of the Gigante seep hunting program to analyze these techniques. At GC600, we evaluate the results of a backscatter intensity normalization, assess the effectiveness of 2X seafloor coverage in resolving seep-related features in backscatter data, and determine the off-nadir detection limits of bubble plumes using the EM302. Incorporating these techniques into seep hunting surveys can improve the detectability and sampling of seafloor seeps.     

Molecular signals for anaerobic methane oxidation in Black Sea seep carbonates and a microbial mat

Linked to gas seeps on the Ukrainian shelf (northwestern Black Sea), massive authigenic carbonates form as a result of anaerobic methane oxidation. Lipid distributions in these ‘cold seep’ carbonates and an associated microbial mat were investigated for process markers reflecting the presence and metabolic activity of distinctive methane-related biota. The samples contain free, irregular isoprenoid hydrocarbons, namely the tail-to-tail linked acyclic C₂₀-isoprenoid 2,6,11,15-tetramethylhexadecane (crocetane), its C₂₅-homologue 2,6,10,15,19-pentamethylicosane (PMI), and several unsaturated derivatives thereof. Furthermore, specific acyclic and cyclic C₄₀-isoprenoids were released upon ether cleavage of the polar fraction from the carbonate. The abundance of these compounds indicates a pronounced role of particular Archaea in the biogeochemical cycling of carbon at methane seeps. Stable carbon isotopic analyses of these lipids reveal extraordinary depletions in ¹³C corresponding to δ-values in the range of −100±30‰ PDB, whereas other compounds show isotopic compositions normally observed for marine lipids (around −30‰ PDB). The isotope data imply that the biosynthesis of the archaeal isoprenoids occurred in situ and involved the utilization of isotopically depleted, i.e. methane-derived, carbon. Apart from archaeal markers, the carbonate and the mat contain authigenic, framboidal pyrite and isotopically depleted fatty acids, namely iso-, and anteiso-branched compounds most likely derived from sulphate-reducing bacteria (SRB). The indications for a tight association of these normally competitive organisms support a model invoking a syntrophic relationship of SRB with Archaea responsible for the anaerobic oxidation of methane. The biomarker patterns obtained from the Black Sea samples were further compared to those from a Oligocene seep carbonate (Lincoln Creek Formation, WA, USA) in order to evaluate their biomarker potential for ancient settings. The prominent occurrence of isotopically light crocetane (−112‰) and PMI (−120‰) meets the findings for the contemporary materials. Thus, isotopically depleted isoprenoids provide diagenetically stable fingerprints for the reconstruction of carbon cycling in both, modern and ancient methane seep systems.     

Dissolved and particulate organic carbon in hydrothermal plumes from the East Pacific Rise, 9°50′N

Chemoautotrophic production in seafloor hydrothermal systems has the potential to provide an important source of organic carbon that is exported to the surrounding deep-ocean. While hydrothermal plumes may export carbon, entrained from chimney walls and biologically rich diffuse flow areas, away from sites of venting they also have the potential to provide an environment for in-situ carbon fixation. In this study, we have followed the fate of dissolved and particulate organic carbon (DOC and POC) as it is dispersed through and settles beneath a hydrothermal plume system at 9°50′N on the East Pacific Rise. Concentrations of both DOC and POC are elevated in buoyant plume samples that were collected directly above sites of active venting using both DSV Alvin and a CTD-rosette. Similar levels of POC enrichment are also observed in the dispersing non-buoyant plume, ∼500 m downstream from the vent-site. Further, sediment-trap samples collected beneath the same dispersing plume system, show evidence for a close coupling between organic carbon and Fe oxyhydroxide fluxes. We propose, therefore, a process that concentrates POC into hydrothermal plumes as they disperse through the deep-ocean. This is most probably the result of some combination of preferential adsorption of organic carbon onto Fe-oxyhydroxides and/or microbial activity that preferentially concentrates organic carbon in association with Fe-oxyhydroxides (e.g. through the microbial oxidation of Fe(II) and Fe sulfides). This potential for biological production and consumption within hydrothermal plumes highlights the importance of a multidisciplinary approach to understanding the role of the carbon cycle in deep-sea hydrothermal systems as well as the role that hydrothermal systems may play in regulating global deep-ocean carbon budgets.     

A comparison of black smoker hydrothermal plume behavior at Monolith Vent and at Clam Acres Vent Field: Dependence on source configuration

Quantitative visualization of acoustic images is used to compare the properties and behavior of high temperature hydrothermal plumes at two sites with different source configurations, increasing our understanding of how plume behavior reflects source configuration. Acoustic imaging experiments were conducted at the Clam Acres area of the Southwest Vent Field, 21°N East Pacific Rise and at Monolith Vent, North Cleft segment, Juan de Fuca Ridge. At Clam Acres, black smokers discharge from two adjacent chimneys which act as point sources, whereas multiple vents at Monolith Vent define a distributed elliptical source. Both plumes exhibit consistent dilution patterns, reasonable fits to the expected power law increase in centerline dilution with height, and simple bending of plume centerlines in response to ambient currents. Our data suggest that point source vents are associated with ordered plume structure, normal entrainment rates, and initial expansion of isosurfaces while distributed source vents are associated with disorganized plume structure, variable entrainment rates, and initial contraction of isosurfaces.     

A new optical instrument for the study of bubbles at high voidfractions within breaking waves

An optical imaging system (BubbleCam) has been tested for the quantification of bubble distributions at high void fractions formed beneath breaking waves. The instrument consists of a CCD video camera, stroboscopic light source, and optics allowing adjustable magnification, a fixed imaging volume, and the resolution of bubbles 3 pixels in radius and larger (equivalent to a minimum bubble radius of about 200 μm in the test configuration). BubbleCam has been deployed in a shore-based configuration (data and power supplied via shore-connected cables) as well as an autonomous device in the open sea with its own power supply and data storage. The resulting images are processed using a variant of the Hough transform which allows computer-automated counting and measurement of the bubbles within the video frames. In addition, images can be qualitatively examined to provide insights into bubble plume evolution and creation mechanisms     

热液和冷泉10种无脊椎动物基因组大小测定及比较

基因组大小(或称C值)作为生物单倍体细胞中全套染色体的DNA总量,在一定程度上是恒定的,因而C值可以作为生物物种的一个特定参数。深海热液和冷泉为更好地理解C值与不同环境之间的关系提供了一个特征性的模型。本文采用流式细胞术,测定了来自热液和冷泉环境中的10种深海无脊椎动物的C值,其分布范围从0.87 pg到12.28 pg,其中,相比于软体动物和多毛类,甲壳生物基因组大小及变异均较大。对比热液和冷泉两个群落中共有种(深海偏顶蛤Bathymodiolus platifrons、柯氏潜铠虾Shinkaia crosnieri以及长角阿尔文虾Alvinocaris longirostris)的基因组大小,发现C值差异并不显著。同时,综合已有的数据,对深海化能极端环境与其他环境条件下的物种C值进行对比分析,结果显示深海化能极端环境下生物的基因组大小并没有发现明确的变化趋势。     

Satellite-measured temporal and spatial variability of the Tokachi River plume

River plumes have important effects on marine ecosystems. Variation in the extent and dispersal of river plumes is often associated with river discharge, wind characteristics and ocean circulation. The objectives of this study were to identify the Tokachi River plume by satellite, determine its relationship with river discharge and clarify its temporal and spatial dynamics. SeaWiFS multispectral satellite data (normalized water-leaving radiance: nL_w) with 1.1 km spatial resolution were used to determine the spatial and temporal variability of the plume during 1998–2002. Supervised maximum likelihood classification using six channels of nL_w at 412, 443, 490, 510, 555 and 670 nm with each band's spectral signature statistic was used to define classes of surface water and to estimate the plume area. Supervised maximum likelihood classification separated three to four classes of coastal water based on optical characteristics as a result of wind stress events. The satellite-observed plume area was correlated with the amount of river discharge from April to October. The plume distribution patterns were influenced by wind direction and magnitude, the occurrences of a near-shore eddy field and surface currents. Empirical orthogonal function (EOF) was used to express the spatial and temporal variability of the plume using anomalies of nL_w(555) monthly averaged images. The first mode (44% of variance) showed the turbid plume distribution resulting from re-suspension by strong wind mixing along the coast during winter. This mode also showed the plume was distributed along-shelf direction in spring to early autumn. The second mode (17% of variance) showed spring pattern across-shelf direction. EOF analysis also explained the interannual variability of the plume signature, which might have been affected by the flow of the Oyashio Current and the occurrence of a near-shore eddy field.