Dissolved organic matter in the northwestern Pacific marginal seas: Insight into the distribution of its optical properties

Dissolved organic carbon (DOC) and the optical properties of dissolved organic matter (chromophoric- and fluorescent dissolved organic matter; CDOM and FDOM) were measured to determine the distributions and drivers of DOM in the northwestern Pacific marginal seas, including the East China Sea, the Yellow Sea, and the East/Japan Sea, in August 2020. In this study, the concentrations of DOC and CDOM/FDOM in surface water showed good correlations with salinity, indicating a predominant contribution from the Changjiang River. However, significantly high concentrations of DOC and FDOM were also observed in the central Yellow Sea region, which seems to be produced mainly from the continental shelf-water and enriched over the water residence time of the Yellow Sea, while the lowest concentrations of DOC and CDOM/FDOM were found in high-salinity waters near the southern sea of Korea. In addition, the East/Japan Sea showed relatively low DOC concentrations and high FDOM values. To distinguish the physical mixing and biogeochemical processes of DOM, we estimated the water mass fractions using an optimum multi-parameter analysis with hydrological and DOM parameters of the major water masses in this region. Our results showed that five primary external drivers of DOM distribution were intrusions of (1) the Changjiang diluted water and (2) the Kuroshio Current water in the surface water, the mixing of (3) the Yellow Sea bottom cold water and (4) the Kuroshio sub-surface water in the bottom layer, and (5) upwelling of the East/Japan Sea water. We conclude that the distribution of DOM is significantly influenced by the mixing of multiple water masses, and the optical signature of DOM can be an effective proxy for tracing the origins and characteristics of DOM in this region.


Introduction
Dissolved organic matter (DOM), mostly composed of dissolved organic carbon (DOC), is one of the largest reservoirs of reduced carbon on Earth and is the most dynamic organic carbon pool linked with other materials in marine ecosystems (Carlson and Hansell, 2015). The optically sensitive fractions of DOM, which are characterized by light absorption and fluorescence, are termed chromophoric-DOM (CDOM) and fluorescent-DOM (FDOM), respectively (Stedmon and Nelson, 2015). CDOM is the lightabsorbing fraction of DOM, and a fraction of CDOM that fluoresces after absorbing light at a certain wavelength is called FDOM. In general, CDOM accounts for approximately 20-70% of the marine DOC in coastal regions, while it contributes only a minor portion in the open ocean (Coble, 2007). The spectral properties of CDOM and FDOM have been used to decipher the behavior of marine DOM, including its origin, molecular weight, and aromaticity (Coble, 2007;Hansen et al., 2016).
The northwestern Pacific marginal seas, including the East China Sea (ECS), the Yellow Sea (YS), and the East/Japan Sea (EJS), contain one of the largest continental shelves in the world, around the Korean Peninsula, the Japanese archipelago, and the Chinese continent. The productive semi-enclosed shallow water in the northern part is called the YS, while in the southern part, a branch of the Kuroshio Current transports relatively saline, warm, and oligotrophic water into this region (Kim et al., 2004). The EJS resembles the global ocean in terms of its own deep-water circulation (Kim and Kim, 1996). The Changjiang River (also known as the Yangtze River) discharges large amounts of freshwater and terrestrial materials (~80%), contributing significantly to the biogeochemical cycles in this region (Ichikawa and Beardsley, 2002;Wang et al., 2012). Due to the large tidal range, significant contributions of chemical substances were found to originate from organic-rich shelf sediments in the water column Kim et al., 2018). In addition, diffusive fluxes from the submarine groundwater discharge Kim et al., 2018) and atmospheric fluxes (Yu et al., 2018;Ren et al., 2022) to marginal seas play an important role in material budgets in this region. Therefore, the northwestern Pacific marginal seas have complex marine systems. However, a poor understanding of the inflow of multiple external water masses and internal biogeochemical processes still induces less understood behaviors of organic matter in this region.
Recently, the origins and characteristics of DOM have been extensively studied in the northwestern Pacific marginal seas (Ogawa et al., 2003;Bai et al., 2013;Kim et al., 2018;Zhang et al., 2019;Kim et al., 2020;Yamashita et al., 2020;Han et al., 2022). Based on the significant correlation between humic-like FDOM (FDOM H ) and radium isotope and dual carbon isotopic compositions (d 13 C and D 14 C) of shelf-water, this continental margin has been recognized as an important source of both FDOM H and DOC (Kim et al., 2018;Han et al., 2022). These studies further suggested that the shelf-water fluxes of FDOM H and DOC were almost equivalent to 30-40% and comparable to that from the riverine discharge, respectively (Kim et al., 2018;Han et al., 2022). Also, Kim et al. (2020) reported the existence of a 'DOC pulse' in the lower salinity range (24-35) due to the enhanced microbial metabolism in this region. Because of these unique DOM characteristics and behaviors, various approaches to DOM research have been applied in these productive marginal seas.
The purpose of this study was to characterize the distribution of various DOM parameters and to assess the potential of CDOM and FDOM as tracers of water-mass mixing in the northwestern Pacific marginal seas. Here, we hypothesized that DOM properties, including unique absorption and fluorescent properties, are conservative within the source water types (SWTs), which is the intrinsic water mass term of each different property, during the water residence time of the northwestern Pacific marginal seas. Thus, the optical properties of DOM can be used to trace watermass mixing patterns along the entire marginal area to understand the biogeochemical drivers and physical mixing of complex DOM dynamics.

Sampling and hydrographic properties
Sampling was conducted during the summer from August 11 to 21, 2020 aboard the R/V Onnuri of the Korea Institute of Ocean Science and Technology (KIOST). Seawater samples were collected by using the 10 L Niskin bottles at 21 sampling stations in the northwestern Pacific marginal seas ( Figure 1). The hydrographic properties, including temperature and salinity were measured using Map of the sampling stations, bottom topography, and simple schematic patterns of surface currents in the northwestern Pacific marginal seas. The bold color line in yellow indicates the transect used in this study. The arrows indicate the current systems in operation during the sampling periods. a CTD device (SBE911-plus, Sea-Bird Electronics Inc., WA, USA) attached to a rosette sampler.

Analysis of DOC
Seawater samples for DOC analysis were vacuum filtered through pre-combusted GF/F glass microfiber filters (500°C for 5 h; pore size: 0.7 mm; Whatman Inc., Maidstone, UK). The filtrate was acidified (pH~2) using 6 M hydrochloric acid (Sigma-Aldrich, MO, USA) to prevent microbial degradation and stored in precombusted (500°C for 5 h) glass ampoules (Fisher Scientific, NH, USA) at room temperature in a dark cabinet. All samples were measured within two months. The DOC concentration was measured via a high-temperature catalytic oxidation method using a TOC analyzer (TOC-L, Shimadzu, Japan). Prior to the sample analysis, the system blank was reduced until the signal from the DOC-free distilled water showed consistent values with the limit of detection (< 2 mM for DOC) to achieve high accuracy. In addition, the accuracy of the DOC concentration was verified with every sample batch run using deep-sea reference material (DSR: 41-44 mM) obtained from the University of Miami (Hansell, 2005). The precision of these DSR measurements was in good agreement with the consensus values within 2%.

Analysis of CDOM and FDOM
Samples for CDOM and FDOM analysis were filtered with the DOC samples and stored in a dark amber vial (Wheaton, NJ, USA) at 4°C. All samples were measured within two weeks. Absorbance spectra and fluorescence excitation-emission matrices (EEMs) were measured by using a spectrofluorometer (Aqualog, HORIBA Jobin Yvon, MA, USA). The CDOM absorbance was blank-corrected, and then a baseline correction was applied at 600 nm, assuming that there is negligible CDOM absorption at that certain wavelength. The CDOM absorbance at 350 nm was converted into the Napierian absorption coefficient (a 350 , unit: m -1 ). The spectral slope coefficients between 275 and 295 nm (S 275-295 , unit: nm -1 ) and between 350 and 400 nm (S 350-400 , unit: nm -1 ) were calculated as the linear regression slope of the log-transformed absorption data of the samples (Helms et al., 2008). The spectral slope ratio (S R ) was determined as the ratio of the slope of the shorter wavelength (S 275-295 ) relative to that of the longer wavelength (S 350-400 ) (Helms et al., 2008). The specific UV absorbance (SUVA 350 , unit: L mg C -1 m -1 ) was derived from the UV absorbance at 350 nm and normalized to the DOC concentration (Helms et al., 2008).
The fluorescence analysis was performed using scanning emission wavelengths of 211.2-617.4 nm in 3.2-nm increments and excitation wavelengths of 240-600 nm in 5-nm increments with an integration time of 3 s. All EEM processing and parallel factor analysis (PARAFAC) modeling were carried out using the MATLAB software R2021a (The MathWorks, MA, USA) with the drEEM toolbox (Murphy et al., 2013). The inner filter effect was corrected using the UV absorbance values of each sample. The nonnegativity constraint was applied to all three modes. The fluorescence intensities of the samples were normalized daily to the Raman peak area of ultrapure Milli-Q water at an excitation wavelength of 350 nm, and these values were stated as Raman Units (R.U.) (Lawaetz and Stedmon, 2009).

Optimum multi-parameter analysis
The optimum multi-parameter (OMP) analysis (Tomczak and Large, 1989) was conducted to determine the fractional contributions of different water masses present in this region. The OMP analysis assumes that the observed water masses can be represented as a mixture of source water types (SWTs). Then, the relative contribution of an SWT to the observed water masses can be determined from the following linear system of equations: where G is the SWT matrix (Table 1), x is the unknown solution vector (water mass fraction of the SWT), d is the vector of the observational data, and R is the residual vector. The system is solved by a least square method to minimize the residuals: where W is the weighting matrix, and W T and R T are the transposes of the matrices W and R, respectively.

Source water types
The salinity and temperature ranged from 26.77 to 34.57 (Avg. 32.34 ± 1.68) and from 4.41 to 29.46°C (Avg. 19.41 ± 6.17°C), respectively, in the entire study region (Figures 2, 3A). The temperature and salinity were found to be strongly stratified in the water column during the sampling campaign ( Figure 3A). In this region, the water masses include the Changjiang diluted water (CDW), the Kuroshio Current water (KCW), the Kuroshio subsurface water (KSSW), the Yellow Sea bottom cold water (YSBCW), and the East/Japan Sea water (EJSW) as previously documented (Hur et al., 1999;Ichikawa and Beardsley, 2002;Chang and Isobe, 2003;Chen, 2009;Wang et al., 2014;Lee et al., 2016;Lee et al., 2020) ( Figure 2). The SWTs, which are the specific properties characterizing the different water masses, were defined based on the temperature versus salinity (T-S) distributions and the literature values (Hur et al., 1999;Ichikawa and Beardsley, 2002;Chang and Isobe, 2003;Chen, 2009;Lee et al., 2014) (Figure 2). The SWT definitions and their specific properties considered in this study are summarized in Table 1.
The CDW is characterized by low salinity (< 31) with high temperatures (> 27°C) in the surface layer (Chang and Isobe, 2003;Lee et al., 2014) (Table 1 and Figure 2). The CDW transports lowsalinity water from the river mouth to the EJS through the surface layer during the summer in this region (Chang and Isobe, 2003). The KCW is a western Pacific boundary current, one of the major ocean currents characterized by high water temperature (19-26°C) and high salinity (34.5-34.7) (Ichikawa and Beardsley, 2002) ( Figures 3A, B). Below the surface layer of the KCW, the KSSW displays little variation in its temperature (slightly lower than KCW; 19-23°C) with slightly higher salinity (34.7-34.9) than the KCW (Zhou et al., 2018) The YS shows relatively large variations in salinity (30.5-34.7) and temperature (8-27°C) with varying contributions of other water masses (the YSBCW, the KCW, and the CDW). In this region, the YSBCW is differentiated from the YS water (YSW) by salinity between 32.5 and 33.5 and temperature near 10°C at depths greater than 50 m (Chen, 2009;Wang et al., 2014;Lee et al., 2016) (Table 1 and Figure 2). Since the water residence time of the YSW is about five years, the biogeochemical characteristics of the YS are significantly altered during its residence time in this region (Nozaki et al., 1991;Hur et al., 1999;Kim et al., 2005;Chen, 2009) (Figure 2). However, the YSW was not considered individually in the water mass mixing model since the hydrographic characteristics of the YS water were determined by mixing between the YSBCW and the CDW.
The EJSW is characterized by a low water temperature (< 10°C) and relatively high salinity (> 34.40), which flows across the Korean strait into the EJS (Kim et al., 2004) (Table 1 and Figure 2). The EJS is a semi-enclosed marginal sea with a maximum depth of over 3500 m. In the EJS, the Changjiang River is known as the major source of DOM in the surface water .

FDOM components
In this study, four FDOM components were determined using the PARAFAC model ( Figure S1). To characterize the properties and the origins of each FDOM component, the fluorescence spectra were compared with the OpenFluor database (Murphy et al., 2014). The C1, C2, C3, and C4 components yielded statistical matches with FDOM components from a total of 21, 9, 5, and 6 previous studies, respectively (Table S1). The C1, C3, and C4 components presented similar fluorescent spectra at relatively longer wavelength ranges to the FDOM H component. The maximum peak location of the C1 component (Max ex/em = 320/398.4 nm) was similar to that of the marine humic-like component, identified as the M peak (Coble, 2007). The C2 component (Max ex/em = 290/323.3 nm) was a relatively sharp spectrum located at a shorter wavelength, described as a protein-like (tryptophan-like) component associated with autochthonous organic substances (Coble, 2007). The C3 component (Max ex/em = 305(425)/517.4 nm) seemed to be a mixture of two humic-like components, peaks A and C, which are classified as having terrestrial origins (Coble, 2007). However, both components have also been reported in deep-ocean and porewater samples (Burdige et al., 2004;Kim and Kim, 2016). The C4 component (Max ex/em = 365/444.4 nm) was also observed at relatively longer wavelengths. Based on the results from the OpenFluor database, the C4 component was characterized as the   (Chang and Isobe, 2003;Lee et al., 2014). KCW, Kuroshio Current water (Zhang et al., 2007;Zhou et al., 2018). KSSW, Kuroshio sub-surface water (Zhang et al., 2007;Zhou et al., 2018). YSBCW, Yellow Sea bottom cold water (Wang et al., 2003;Li et al., 2015). EJSW, East/Japan Sea water (Kim et al., 2004). YSW, Yellow Sea water (Hur et al., 1999;Chen, 2009). terrestrial source and the best predictor of lignin phenol concentration by Walker et al. (2013).

Distributions and origins of DOM
In this study, the concentrations of DOC and FDOM showed good negative correlations with salinity, indicating the predominant contribution of the Changjiang River (Figures 2, 4). However, there was an enhanced concentration of DOC found above the endmember (the Changjiang River source) mixing line, which is marked as the CDW source in Figure 4A. This additional input of DOC has been attributed to the microbial activity in the water column near the estuarine mixing zone . Also, almost similar or even higher concentrations of DOC than those of the CDW were found in the YSW ( Figure 4A). In this region, this additional DOC source was suggested to be produced from the shelf-water based on the dual carbon isotope analysis (marine d 13 C signature and young 14 C age) (Kim et al., 2018;Han et al., 2022). These enhanced concentrations and distributions of DOC were found to be highly dependent on the age of the shelf-water, which is also consistent with those of enhanced FDOM H concentrations in the shelf region (Kim et al., 2018;Han et al., 2022). Although photochemical degradation is one of the major removal mechanisms of FDOM in the surface layer of aquatic environments (Coble, 2007), the C3 (r 2 = 0.49; p<0.05) and C4 (r 2 = 0.40; p<0.05) components showed good correlations with A B D C  salinity in this study ( Figures 4B, C). Such a negative correlation may suggest that the water mixing process is more distinct than the other processes, perhaps owing to a lower degree of photochemical degradation of DOM than the water mass mixing during the longrange transport in this region (Kim et al., 2018;Kim et al., 2020). In the YSW, both DOC and FDOM concentrations showed similar mixing patterns, and the water masses were mainly consisted of the CDW and the YSBCW (Figure 2). The KCW was differentiated by the lowest concentrations of DOC and FDOM in the high-salinity waters (Figures 2, 4). However, the EJSW showed low concentrations of DOC with high concentrations of FDOM in high-salinity waters compared with that of the KCW (Figure 2). This high concentration of FDOM could be caused by microbial activity in the ocean interior (> 200 m) and transported to the surface layer by own circulation system and upwelling .

Evaluation of DOM parameters for OMP analysis
The multiple water masses can be clearly differentiated based on three end-member models using the classical T-S method: the CDW, the YSW, and the KCW. However, it is difficult to estimate the distributions and proportions of water masses where various water masses (more than three) are present, with varying contributions. Thus, in this study, the OMP analysis was applied together with various DOM parameters (DOC, C3, C4, and SUVA 350 ) to determine the relative mixing fractions of the different water masses and to determine the driving force of the distributions of DOM parameters in this shelf region.
There could be several limitations with applying DOM parameters to the OMP model in the shallow marginal seas. In general, the OMP method is not applicable in such productive marginal seas, where the oxygen and nutrients are not conservative (Chen, 1988;Zhou et al., 2018). However, these methods using radium isotope and dissolved inorganic iodine have been successfully applied as conservative tracers to determine the distribution of different water masses in the ECS shelf region in previous studies (Zhang et al., 2007;Zhou et al., 2018).
In addition, some chemical parameters such as oxygen, nutrients, and isotope values may not be useful for water source discrimination due to their large involvement in the air-sea exchange and biological activities (Chen, 1988;Zhang et al., 2007;Amaral et al., 2020). In particular, Amaral et al. (2020) suggested that the distribution of DOC was not dependent on the water mass mixing due to its multiple sources and complicated biogeochemical processes. However, DOM parameters used in this study, including DOC (r 2 = 0.55; p<0.05), C3 (r 2 = 0.49; p<0.05), C4 (r 2 = 0.40; p<0.05), and SUVA 350 (r 2 = 0.19; p<0.05) (Figure 4), showed no remarkable internal changes during the water mass mixing in this region. This indicates that the water mass mixing is a more dominant process than the other biogeochemical processes and/or external source inputs in this region.
Both CDOM and FDOM are known as reliable proxies for DOM, and the correlations between CDOM/FDOM and salinity showed consistent behavior with those of DOC. Although DOM components (DOC, FDOM, and CDOM) showed good correlations with salinity, each component had a different slope in this study, indicating the different end-member values and mixing processes of each water mass in this region (Figure 4). This shows the individual characteristics of SWTs and validates whether these DOM parameters can be used as conservative tracers for each water mass in this region. However, although FDOM is a powerful tracer for determining the DOM origins in the ocean, different FDOM components can be selectively degraded during their transport and/or mixing (Romera-Castillo et al., 2011). In this study, the C1 component was classified as the marine humic-like M peak, which could be preferentially degraded by marine bacteria, and the terrestrial humic-like peaks seem more resistant to bioavailable degradation (Romera-Castillo et al., 2011). Thus, the C3 and C4 components (terrestrial humic-like peaks) that showed consistent behavior with salinity were applied to the OMP analysis ( Figures 4B, C). The optical signature of DOM, including FDOM, a 350 , and SUVA 350 , showed distinct characteristics in each water mass because different external origins contributed to the unique DOM behaviors in the water masses in this region ( Figure 3D-J). For example, SUVA 350 is known to be highly associated with the molecular weight and aromaticity of DOM (Helms et al., 2008). In this study, SUVA 350 showed no remarkable internal variations during the water mass mixing and the concentrations of DOC and a 350 values also showed a positive correlation (r 2 = 0.40; p<0.05) among different water masses ( Figures 4D, S2), which may derive such a conservative mixing pattern from the source to the receptor region, which is similar to other DOM parameters used in this study. Thus, SUVA 350 may allow us to identify different aromaticbiased carbon contents over the overall organic carbon content in the water masses (Figures 3, 4D). These characteristics may also allow for the use of DOM parameters as conservative tracers in this region.

Distributions of SWTs
Based on the OMP calculation, we determined the mixing proportions of five different water masses in this region: the CDW, the YSBCW, the KCW, the KSSW, and the EJSW ( Figure 5). The shelf-water, which flows in the YS and the ECS regions, was mainly consisted of the CDW (10-90%), the YSBCW (40-70%), and the KCW (10-40%), which is also consistent with previous studies (Chen, 2009;Zhou et al., 2018) (Figure 5). As expected, the contribution of the YSBCW was up to 70% of the total water in the YS region ( Figure 5). The contribution of the CDW was greater in the surface layer near the discharge region and gradually decreased with increasing distance from the river mouth (Figure 4). The contributions of the KCW (~80%) and the KSSW (~90%) were greater near the southern sea of Korea where the Kuroshio Current mainly passes through ( Figure 5). The contribution of the EJSW was significant only at the stations near the EJS area ( Figure 5).
Since the water masses flow from the innermost region of the YS to the EJS passing the ECS and the southern sea of Korea following the pathway of the water current, the water mass distributions with the mixing proportions found from the OMP analysis clearly reflected the pattern of the water current in this region ( Figure 5). In this study, the YSBCW and the KCW showed the most abundant distributions, with the average contributions of 61 ± 18% and 54 ± 30%, respectively. The contribution of the CDW was estimated to be approximately 41 ± 32% in the surface water of the shelf region. These results were also consistent with the previous studies in this region (Zhang et al., 2007;Lee et al., 2014;Zhou et al., 2018). Lee et al. (2014) reported that the contributions of the YSW, the KCW, and the CDW were approximately 30-50 ± 4%, 65 ± 15%, and 30 ± 7%, respectively. On the other hand, Zhang et al. (2007) summarized that the water mass contributions were about 14% from the KSW, 43% from the KSSW, 22% from the CDW, and 22% from the Taiwan Strait Warm Water (TSWW) in this region. Zhou et al. (2018) also found that the surface water in the ECS is mainly mixed by the CDW, the KCW, and TSWW, and the contribution of each water mass was approximately 30-90%, 50%, and 70%, respectively. Some discrepancies or large ranges could be found due to the spatial and interannual variabilities among different studies.
In addition, the relationships between the value calculated by the OMP analysis versus the observed values of DOC (r 2 = 0.63; p < 0.01), FDOM H (r 2 = 0.52; p < 0.01), and SUVA 350 (r 2 = 0.54; p < 0.01) showed good correlations. Although the coefficients of determination were not significantly high (close to 1), this also suggests a possibility and more evidence for using DOM parameters with this model in such productive marginal seas. In a previous study, the results of the OMP analysis were inconclusive in ocean systems with multiple sources and complicated internal biogeochemical cycles (Amaral et al., 2020). Although our study region has a complex mixing trend between various water masses, the optical properties of DOM in the water masses were distinct and influenced by the conservative mixing in the entire northwestern Pacific marginal seas. Also, CDOM and FDOM, which is the sufficiently easy and simple methods of estimating DOM, could be more applicable and efficient in such productive marginal seas. A large deviation or scatter found from the linear regression may be associated with the additional inputs (atmospheric depositions and diffusive groundwater fluxes) and outputs or internal cycles (biogeochemical processes, biodegradation, and photodegradation) during water mixing.

Contributions of SWTs
To investigate the DOM drivers and their relationship to water mass mixing in this region, we compared the water mass fractions (f) of the SWTs together with the DOM parameters using the Pearson correlation coefficient value (Table 2). Here, the positive regression coefficient (r) values between the f value and DOM parameters indicate the production and/or input of DOM by the water mass f, while the negative values would indicate the dilution effect and/or degradation process by the additional water mass inputs ( Table 2).
The fraction of the CDW (f CDW ) was significantly correlated with the concentrations of DOC (r=0.88), C1 (r=0.78), C3 (r=0.81), C4 (r=0.74), and a 350 (r=0.89) ( Table 2). This indicates that the Changjiang River primarily contributes to the DOM budget, which mostly contains humic organic matter in this region. This is supported by the fact that the Changjiang River is one of the largest rivers in the world and discharges approximately 80% of the freshwater and terrestrial materials in the northwestern Pacific marginal seas (Wang et al., 2012). The river-driven DOC budget is believed to be produced by direct discharge from the land and indirect biogeochemical processes in the estuary mixing zone in this region . The negative r value of S R (r= −0.48) with high positive r value of SUVA 350 (r= 0.66) values may represent more aromatic and less irradiated DOM originating from the CDW. This was attributed to the fact that river water contains a significant amount of aromatic terrestrial organic matter from the soils and sediments (Hopkinson et al., 1998;Miller et al., 2009).
The fraction of the YSBCW (f YSBCW ) also showed relatively high coefficient values for DOC (r=0.49), C3 (r=0.43), and a 350 (r=0.47) ( Table 2). As mentioned above, higher concentrations of DOC and FDOM were found in the YS region (high salinity water) than those in the CDW (low salinity water). It seems once the river-driven DOC is enriched in the estuarine mixing zone , transported to the shelf region, and continuously accumulated in the shelf-water over the residence time as previously documented (Han et al., 2022) (Figure 4A). Here, the source of DOC in the YS shelf-water was found to be from the fresh marine DOM sources based on the dual carbon isotopic compositions (marine d 13 C signature and young 14 C age) (Han et al., 2022). Similarly, Kim et al. (2018) suggested that FDOM H (both M peak and C peak) was mainly produced in the bottom sediment in the central YS region and enriched over the residence time (~5 years) in the shelf-water, as previously explained. In this study, the C1 component slightly increased in the central YSW ( Figure 3F) and correlated with DOC concentrations (r 2 = 0.24; p<0.05) ( Figure S3), which can be also evidenced by the enhanced marine origin DOM in the shelf waters. Also, lower SUVA 350 values in the YSBCW than that in the CDW were obtained from the lower a 350 values caused by the aged water mass (~5 years) found in the deep water of the YS region. Thus, the fresh DOM source or DOM production seems to be mainly contributed by the CDW and the YSBCW in the YS region ( Figure 5).
In the case of the KCW (f KCW ) and the KSSW (f KSSW ) fractions, the DOM parameters showed negative coefficient values, indicating a dilution effect of those DOM parameters by these water masses. However, the f KCW and the f KSSW showed negative correlation values with absorbance measurements including, a 350 (a proxy of the abundance of CDOM) (r=−0.39 and r=−0.86) and S R (r=0.38 and r=0.60) ( Table 2). This indicates that these water masses were mainly dependent on the input of the highly irradiated DOM (high S R value) with a low absorptivity and less aromatic content derived from the old Pacific Ocean water. The distinct characteristics of the CDOM parameters (a 350 , S R , and SUVA 350 ) showed that these parameters were mainly derived from the river water source, whose coefficient values were positive ( Figures 3D, E, J). Thus, the Kuroshio-derived water masses (the KCW and the KSSW) were characterized by low biological productivity and a relatively old Pacific Ocean water that intrudes into the continental shelf area (Chen et al., 2021) (Figure 3). The fraction of the EJSW (f EJSW ) showed a significant relationship with S (r=−0.43) and SUVA 350 (r=0.36). Although the f EJSW showed weak coefficient values compared with other water masses, the f EJSW also showed significantly different characteristics from the EJSW itself due to its own deep-water circulation system (Kim and Kim, 1996) (Table 2). The concentration of DOC showed relatively low values with a negative r value (r=−0.11), while the FDOM components showed positive r values (Table 2). In previous studies, FDOM H in the deep layer was found to be produced mainly by the degradation of sinking organic particles (Nelson and Siegel, 2013;Kim and Kim, 2016).
Thus, it seems that the distributions of DOM production and degradation showed a similar pattern to the water mass mixing pattern in this region ( Figure 5). Our results suggest that the distributions of DOM parameters are significantly influenced by the mixing of multiple water masses in the northwestern Pacific marginal seas (Figure 6). In Figure 5, the area where the contribution of each SWT is over 50% was considered to be fully composed of that particular SWT, while the area where the contribution of each SWT is less than 50% was considered to be a mixed area with various water masses. In the YS region, the water masses were mostly composed of the YSBCW (Figure 6). In the ECS area, the surface water (< 25 m) was mainly contributed by the CDW, whereas the bottom water was mostly contributed by the KCW (Figure 6). In the southern sea of Korea, the water masses were mainly mixed by the KCW and the KSSW ( Figure 6). The EJSW enters to the southern sea area as upwelling water from the EJS coastal region ( Figure 6). Because the spatial distribution of dominant SWTs was appropriate, the optical properties of DOM can be used as a powerful proxy for the origins and characteristics of DOM and water-mass distributions in this region.

Conclusions
The different distributions of water masses and DOM components (DOC, FDOM, and CDOM) were determined based on field observations and OMP analysis in the northwestern Pacific marginal seas in August 2020. This study aims to solve the complex mixing processes and detailed proportions of various water masses in productive marginal seas. Thus, a new approach was applied to determine the driving force of the distributions of various DOM parameters by using an OMP analysis with DOM parameters in these productive marginal seas and to validate whether these DOM parameters could be an effective proxy for the water mass tracers. In this region, five different water masses were found, including the CDW, the YSBCW, the KCW, the KSSW, and the EJSW, with varying contributions and mixing proportions. In the surface layer, the external drivers of DOM distribution were observed to be intrusions of (1) the CDW and (2) the KCW, while in the bottom layer, the seawater was mainly contributed by the mixing of (3) the YSBCW and (4) the KSSW together with (5) an upwelling of the EJSW. Our results showed that fresh DOM was mainly derived from the CDW and the YSBCW, while the old DOM was mainly derived from the KCW and the KSSW in this region. In addition, the distributions of the optical properties of DOM were significantly influenced by the mixing of multiple water masses. Our results suggest that the optical signature of DOM can be a good proxy for tracing the origins, characteristics, and distribution of DOM in the northwestern Pacific marginal seas. Further studies are required to investigate the seasonal and interannual variability of the DOM parameters and water masses mixings in this region.

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions
JK and HH conceived and designed the study. HH collected the samples and performed the sample analyses. HK performed the OMP and water mass analyses. JK and HH conducted the statistical data analyses. HH and JK wrote the manuscript. GK, JH, and SN reviewed and commented on the manuscript. All authors contributed to the article and approved the submitted version.

Funding
This research was supported by Korea Institute of Marine Science and Technology Promotion (KIMST) funded by the Ministry of Oceans and Fisheries (20220533). This work was also supported by the National Research Foundation of Korea (NRF) funded by the Korean government (NRF-2021R1C1C1004733). HH was partially supported by Basic Science Research Program through the NRF funded by the Ministry of Education (2021R1A6A3A01086864).