The dissolved organic matter in surface waters is composed of fractions of different molecular weight and polarity, characteristics that determine their capacity for complexing different types of pollutants and their environmental impact. In this study, the dissolved organic matter in the surface water of the Bio-Bio River (Central Region of Chile) was characterized chemically and spectroscopically after fractionating by molecular weight and polarity. The technique of direct sample analysis-time of flight-mass spectrometry (DSA-TOF-MS) was used to obtain more information on the composition of dissolved organic matter. It is concluded that dissolved organic matter found in the water of the river from the site of minor human impact (Rucalhue) has a predominantly natural origin, with a high content of aromatic carbon, in contrast to dissolved organic matter found in the waters of the sites that have higher human impact (Laja and Concepción), characterized by a greater molecular size and higher organic carbon content. These results are consistent with those obtained from DSA-TOF-MS, where higher correlation was observed between the mass spectrum of the standard commercial humic acid and dissolved organic matter found in the sectors of Laja and Concepción, unlike the spectrum mass of lignin which is more like dissolved organic matter found in the sector Rucalhue.

1. Introduction

Organic matter is a complex and varied mixture of aromatic and aliphatic hydrocarbons with functional groups such as amides, carboxyl, hydroxyl, and ketones and others in smaller amounts [13]. This organic matter may come from natural sources like the soil or from anthropic sources like organic amendments (compost or biosolids) or spillage from sanitary landfills. Two kinds of organic matter can be identified, usually more abundant, which under normal soil conditions are insoluble and another pool of dissolved organic matter (DOM). The latter can migrate to water from different kinds of sources, either natural or anthropic, such as that added to the ecosystem by spillage from sanitary landfills or through the addition of compost or biosolids to the soil (agricultural land, forests, and pastures) [4].

The substances that compose the DOM have the ability to retain and/or mobilize organic and inorganic pollutants in the environment [1, 5, 6], and it is the kind and nature of the existing functional groups present, like carboxyl (-COOH), phenol (benzene-OH), alcohol (-CH2-OH), and methoxyl (-OCH3), that can influence its interaction ability with the different pollutants [7, 8].

The variety of molecules that constitute the DOM and their low concentration in water systems make their characterization complex [4, 9]. Therefore, to facilitate their study it is often necessary to isolate and concentrate them by methods such as reverse osmosis, evaporation under reduced pressure, or freeze-drying [2, 4, 10] and then fractionate them according to their polarity using fractionating resins like Amberlite XAD or by molecular size through techniques like ultrafiltration, ultracentrifugation, and high performance size exclusion chromatography [6, 11, 12]. When analyzed by H1-NMR, FTIR, fluorescence, UV-Vis, and GC-MS techniques, these fractions provide information related to the components of the DOM, and they are used as its fingerprints [6, 13, 14], characteristics that allow predicting its behavior in the environment and its effect on a variety of pollutants [15].

In recent years research work related with DOM has increased considerably, because of its effect of retention or transport of organic and inorganic pollutants in different natural substrates [5, 16, 17].

In Chile there are several studies aimed at the characterization of basins of great significance to agriculture, such as the Traiguén River [18], the Cachapoal and Coya Rivers in the VI Region [19], and the Laja River in the VIII Region [20]. Such studies have been conducted with the objective of determining the levels of various pollutants in order to assess the possible effects on the soil and the products generated by this activity [18, 21]. However, there is little information about the DOM in the rivers of Chile, even though their characterization is of great importance because it can modify the fate of the pollutants present in the water.

The Bio-Bio basin is part of the VIII Region and it is one of the country’s largest basins (24,264 km2) and water flows, covering part of the Provinces of Malleco and Cautín in the IX Region. At the national level the basin represents an important economic development center, linked to the forestry sector, livestock sector (located mainly in the provinces of Ñuble and Bio-Bio), the industrial sector (represented basically by the metallurgical and chemical industries, oil refineries, the textile industry, and the cellulose industries, among others), and the hydroelectric sector, which is the main source of electric power at the national level. All this confirms the importance of evaluating the presence of DOM and its physicochemical characteristics, as well as the effect that it can have on the pollutants present in the Bio-Bio River.

Because of the complexity of the DOM, in a first stage it is necessary to extract it and then fractionate it either according to molecular weight by means of ultrafiltration or according to polarity by the method described by Leenheer [22]. Then the structural characterization of each fraction is made by techniques like elemental analysis, FTIR and UV-Vis analysis, and H1-NMR spectrometry.

Currently there is a lack of a standard analytical methodology for the structural characterization of the DOM, and this can be described only partially from various analyzes. In recent decades, new techniques have been used for the characterization of DOM, such as Fourier Transform Ion Cyclotron Resonance Mass Spectrometry [3], Electrospray Ionization Fourier Transform Ion Cyclotron Resonance Mass Spectrometry [23], and Matrix Assisted Laser Desorption Ionization Time-of-Flight [24, 25].

Direct sample analysis-time of flight-mass spectrometry (DSA-TOF-MS) has proven to be a useful technique in the analysis of natural DOM, with advantages such as speed of analysis, fast screening of samples in seconds, without sample preparation [26]. Direct sample analysis (DSA) is an open atmosphere ambient ionization technique which utilizes power of field-free atmospheric pressure chemical ionization (APCI) source for rapid and direct analysis of samples (solid, liquid, or gas). In relation to DSA ionization mechanism, the heated nitrogen gas volatilizes the analytes present in the sample and hydronium ions donate protons to the analyte forming the ion [M+H]+. The use of TOF-MS also allows greater resolution mass spectrometer with greater accuracy of the measured mass.

In soils and sediments we can identify organic matter of different sizes; within this, carbonaceous geosorbents (CG) can be identified, corresponding to organic material in the 63–500 µm diameter range, with high hydrophobicity; this property allows CG to intervene in the processes of sorption of different organic compounds such as terpenes [27] or contaminants such as polycyclic aromatics hydrocarbons (PAHs) [28]. The origin of the CG is varied; it may be of an anthropic origin, such as those by-products of the combustion of hydrocarbons, or natural, such as pollen or fungi spores; it is possible that part of this material can be transferred to rivers or lakes but due to their size, it becomes part of the sediment of these ecosystems [29]. The DOM, operationally defined as less than 0.45 μm, resembles CG in relation to interaction with organic pollutants. It has been reported that CG, due to the long hydrological residence, can become transformed by natural geochemical processes, contributing to the DOM [30].

The purpose of this paper was to isolate and fractionate the DOM and then compare the structural and chemical characteristics of the isolated fractions of DOM from the same source, but from different sectors of the basin that differ in their anthropic activity. The study involves the use of DSA-TOF-MS technique, which has not been described in the characterization of DOM, hoping to provide new information on the characteristics of DOM from different sources. Such knowledge can assist in our understanding of the composition and chemical properties of DOM samples in aqueous environments.

2. Materials and Methods

2.1. Studied Area

The Bio-Bio River receives industrial and household discharges, and its course was followed from its source to its mouth in the Pacific Ocean, considering three sampling zones that differ in their associated anthropic impact. The first sampling point was the Rucalhue sector, an area with a low water volume characterized by low anthropic impact. The second sampling point was in Laja, between Nacimiento and San Rosendo, with a river width of approximately 2 km, and the samples were taken from a boat; this sector receives wastes from the nearby cellulose plants. The third sampling point was in Concepción, a zone with a large water volume and a high anthropic impact (Figure 1).

2.2. Water Sample Collection and Characterization

The samples were collected every 10 min in high-density polyethylene bottles that had been decontaminated with 10% HNO3, completing a total of ten samples with a volume of approximately 2.5 L each. The samples were transported to the laboratory and were stored at 4°C for their immediate analysis, with a part of them being frozen for later use.

The pH, conductivity, organic carbon content, and suspended solids of the samples were determined. The method of Walkley-Black modified for water was used to determine organic carbon in the water samples, prior correlation of the data with the quantification of total organic carbon (TOC) determined with a Vario TOC Select analyzer (Elementar). Reproducible TOC values with an accuracy of 1% were always found by injecting 500 µL aliquots in the analyzer. Suspended solids were determined by the ASTM method [31].

2.3. Characterization of the DOM

An initial volume of 2 L of sample was used, and it was concentrated in a rotavapor down to a volume of 50 mL, followed by drying in a vacuum oven (WTB Winder), getting a representative solid that was studied by FTIR, H1-NMR, and elemental analysis.

2.4. Collection and Fractionation of DOM

The water samples were filtered through paper (Advantec, 70 mm, number 2) and then the DOM was separated by filtration through a 0.45 μm cellulose nitrate membrane (Sartorius Stedim Biotech), and this fraction was called DOM.

For the fractionation by molecular weight 4.0 L of DOM solution and a continuous ultrafiltration cell (Amicon, Millipore) were used, using an ultrafiltration membrane with 1 kDa cutoff. The fractions greater than 1 kDa (F1) and lower than 1 kDa (F2) were concentrated to a volume of 200 mL in a rotavapor, and this volume was then divided into two fractions: one 150 mL fraction was dried in a vacuum drying oven (WTB Winder) and the remaining 50 mL was treated with sodium azide at a concentration of 100 mgL−1 (100 μL of sodium azide in 50 mL of solution) to avoid the degradation of the DOM by microorganisms, and it was stored at 4°C for its later analysis. The dry residue obtained, representative of the DOM fractions, was studied by FTIR and by the organic carbon analysis described in Section 2.2 (Figure 2).

2.5. Fractionation of the DOM by Polarity

For the fractionation of the DOM according to polarity, use was made of the analytical procedure described by Leenheer [22], which consists of a separation of organic solutes in three fractions operationally defined as hydrophobic base (F3), hydrophobic acid (F4), and hydrophilic fraction (F5). This separation was done according to their adsorption by the ion exchange resin XAD-8 (Supelite, Supelco) (Figure 2).

The procedure was applied in triplicate to the F1 and F2 DOM fractions from each of the sampling points (Rucalhue, Concepción, and Laja), generating 18 fractions in total. In each of these fractions, the organic carbon content was determined in duplicate as described in Section 2.2.

2.6. Analysis of the DOM and Its Fractions by UV-Vis Spectroscopy

A solution with a carbon concentration of 10 mgL−1 was prepared at pH 6.5 (10−3 molL−1 phosphate buffer) with each of the DOM fractions obtained, and their absorbance was then determined by UV-VIS spectroscopy (UV/Vis spectrometer UV2, Unicam) at wavelengths of 204, 254, 250, 365, and 436 nm. The specific UV-absorbance (SUVA) for each sample was determined as the quotient between the absorbance at 254 nm and the DOC concentration in the sample [32].

2.7. Analysis of the DOM and Its Fractions by DSA-TOF-MS

DOM representative of each sampling point and its fractions was used for this experiment. A solution of 5 mgmL−1 of DOM in HPLC grade water containing 0.1 molL−1 NaOH was prepared. Furthermore, standard solutions of commercial humic acid (HA) and alkali lignin (LG) at the same concentration were considered, both standards supplied by Sigma Aldrich. The liquid sample (10 μL) was placed in a mesh holder for analysis.

The AxION direct sample analysis-time of flight-mass spectrometry (AxION-DSA-TOF-MS, Perkin Elmer) conditions were as follows: corona current of 30 μA, heater temperature of 300°C, auxiliary gas (N2) flow of 4 Lmin−1, nebulizer gas (N2) pressure of 552 kPa, drying gas (N2) flow of 3 Lmin−1, and drying gas (N2) temperature of 25°C. The AxION-DSA-TOF-MS was run in positive ionization mode with a flight tube voltage of −10000 V. The capillary exit voltage was set to 100 V for normal MS analysis and 155 V for collision induced dissociation analysis. Mass spectra were acquired with a mass range of 100–3000 m/z and an acquisition rate of one spectrum per second. All samples were analyzed for 10 sec. To maintain mass accuracy, five lock mass ions were used (m/z 121.0509, m/z 622.0299, m/z 922.0119, m/z 1521.9771, and m/z 2121.9405).

3. Results and Discussion

3.1. Characterization of the Water

The amount and composition of DOM vary spatially and in time, not only for the surrounding environment, but also by environmental processes to which it is subjected. The main phenomena affecting the DOM are biodegradation and photodegradation, which at the end of the process would convert inorganic species, for example, CO2 [33]; this implies that one of the determinant factors for studying the DOM is the time of analysis, which must be controlled to minimize variability in the amount and composition of the DOM; for this reason we opted to collect 10 samples per sampling point. This allowed the analysis of samples in parallel, not avoiding degradation of the DOM, but if obtaining DOM under similar conditions for the three sampling points. While the information published recommends applying chemicals such as biocides or work at low temperatures to avoid degradation of the DOM, we check in the laboratory that these methods were not effective.

The organic carbon content (Table 1) of sampling point 3 was greater by one unit than that of points 1 and 2, whose values were 2.4 and 2.9 mg·L−1, respectively. The conductivity was 89.4, 114, and 125 μScm−1 for points 1, 2, and 3, respectively. The results of the determination of the pH of the water were similar, with an average of 7.6 for the three sampling points. With respect to suspended solids, sampling points 1 and 2 had similar values of 12.9 gL−1, in contrast with sampling point 3, whose value was almost three units lower, in agreement with the flow volume differences between the sampling points, with the flow at point 3 greater than that at the other points.

DOM concentrations measured as organic carbon are reported as 1–10 mgL−1 in rivers, 1–50 mgL−1 in lakes, 8.6 mgL−1 in wastewater, and 0.5–52 mgL−1 in soil [34]. This organic matter can interact with the different pollutants present in the soil, and in the specific case of heavy metals, the organic matter can favor their stabilization in the soil solution [35], while in the specific case of DOM from residues such as compost, it can mobilize heavy metals through the soil profile [36].

3.2. Characterization of the DOM

One of the main problems of the study of DOM is the variability in the amount and composition of the material obtained. Therefore, analyzes were performed continuously, verifying the content of organic carbon in the subfractions obtained in each stage, to ensure the presence of DOM throughout the fractionation process by either polarity or molecular weight; such procedures were performed in triplicate, with relative standard deviation (RSD) of organic carbon content lower than 15%. In this sense the DSA-TOF-MS technique is a great alternative, because for each fraction of DOM it is possible to obtain a highly reproducible mass profile, in short analysis times; however, as these profiles are highly complex it is impossible to associate a numerical variability with these results.

With respect to the determination of SUVA, the RSD of the determination of organic carbon and absorbance determined in the DOM were lower than 10%; however, the values are not shown since the results of SUVA and relationships absorbance correspond to determinations obtained as quotients between the individual parameters.

The mass of DOM from river water at sampling points 1, 2, and 3 was 70.7, 63.4, and 56.4 mgL−1, respectively. With respect to elemental microanalysis of the solid fraction of the DOM, the percentages for sampling point 1 were N: 0.1, C: 4.9, and H: 1.02; for sampling point 2, N: 0.5, C: 4.6, and H: 1.04; and for sampling point 3, N: 0.7, C: 5.3, and H: 1.00.

Because of the complexity of the structure of the DOM, the H1-NMR spectra were considered as a pattern of signals rather than as a means of identifying specific compounds (Supplementary Figure in Supplementary Material available online at http://dx.doi.org/10.1155/2016/1537370). It must be noted that the distribution found resembles the spectra reported by [37], for fulvic acid, and that corresponding to the hydrophilic fraction with a signal at 0–4 ppm for the aliphatic region and 6–8.6 ppm for the aromatic region.

The results of the FTIR analyses of the DOM are presented in Supplementary Figure for the three sampling points, showing agreement between the spectroscopic profiles of the different DOM, indicating that the structures and functional groups are similar. The main signals are seen at 3422.3 cm−1, corresponding to the H-bonded OH stretching of carboxyl, phenol, and alcohol; a lower intensity signal at 1637 cm−1 is corresponding to the C=O stretching of COO, ketonic C=O, and aromatic C=C conjugated with the COO of pyridine to amino acid; a signal at 1425,1 cm−1 is corresponding to -C-H of carbohydrates; and a strong signal at 1096.6 cm−1 is corresponding to the C-O group of polysaccharides; only for the DOM from points 2 and 3 a signal is seen at 1384.6 cm−1 that corresponds to polyhydroxylated phenols. The spectroscopic profiles obtained for the DOM from the three sampling points agree largely with the FTIR spectrum of a humic acid, with the difference that the DOM does not present the signal at 2900 cm−1 corresponding to aliphatic C-H stretching [38].

The fractionating of DOM by molecular weight was done using a 1 kDa cutoff, because it was determined that the amount of organic carbon present in the fractions obtained with molecular weights greater than 30 kDa was ten times smaller than that obtained using the 1 kDa cutoff size. Table 2 shows the results of the determination of the organic carbon content in fractions F1 and F2. Upon carrying out the mass balance it is seen that part of the DOM is lost during the ultrafiltration process, which in general takes several days for its performance (120 continuous hours), so it is possible that DOM degradation processes as well as adsorption of it on the ultrafiltration membrane used may continue taking place [4]. The results showed that the predominant fraction in the DOM for the three sampling points is the one greater than 1 kDa, with a higher organic carbon content.

In relation to the FTIR of the F1 of the DOM (Supplementary Figure ), a stronger signal is seen for sampling point 2 at 3500 cm−1 corresponding to a H-bonded OH stretching of carboxyl, phenol, and alcohol; at 1630 cm−1 corresponding to C=O stretching of COO, ketonic C=O, and aromatic C=C conjugated with COO of pyridine to amino acid; at 1400 cm−1 corresponding to -C-H of carbohydrates; and at 1225 cm−1 corresponding to C-O stretching and O-H bending mainly of carboxyl groups. The signal at 1512 cm−1, corresponding to aromatic C=C, N-H bending, and C=N stretching, is seen only in fraction F1 from sampling point 1.

In relation to the FTIR of F2 (Supplementary Figure ), it is noted that the signal at 2920 cm−1, corresponding to -OH groups, is not seen, and the signals at 1630 and 1400 cm−1 are weaker. For fractions with the same molecular weight, a concordance is observed in the spectroscopic profile of the DOM from the three sampling points.

3.3. DOM Fractionation by Polarity

From the fractionation of the DOM by polarity two fractions were recognized: a hydrophilic fraction composed of simple organic acids, polyhydroxylated phenols, carbohydrates, amino acids, and amino sugars and a hydrophobic fraction that may contain phenols, hydrocarbons, fats, nucleic acids, and quinones.

The organic carbon content in the extracts fractionated by polarity (Figure 3) showed that the predominant fraction is hydrophobic, corresponding specifically to the hydrophobic base (F3), with a proportion of organic carbon between 55 and 84% with respect to the sum of the organic carbon of the three fractions, for each sampling point. The highest organic carbon content was that of F3 of the fraction smaller than 1 kDa from point 1, with a mainly natural origin, followed by F5 and finally F4.

When the organic carbon results of the different polarity fractions and the FTIR spectra of F1 and F2 are compared, it is seen that for F1 from point 2 the signals show a DOM with a highly hydrophobic character, because of the signal corresponding to hydrocarbons at 2930–2860 and 1455 cm−1, in agreement with the higher organic carbon content present in F3 [22].

3.4. Analyzes of DOM and Its Fractions by UV-Vis

The SUVA parameter (Table 2), which allows describing the hydrophobic and hydrophilic character of the DOM obtained from natural water, was determined from the results of UV-Vis spectroscopy and the organic carbon content of each of the fractions. The results obtained for the DOM fraction were 2.19, 3.41, and 6.97 for points 1, 2, and 3, respectively, where a SUVA > 4 indicates mainly hydrophobic and especially aromatic matter, while a SUVA < 3 refers to mainly hydrophilic matter [32].

The results obtained for the extract fractionated by molecular weight showed that F2 has a strongly hydrophobic character for the three sampling points, but F1 does not, where a hydrophobic character is found only for sampling points 1 and 2 (Table 2).

The 254/436 absorbance ratios (Table 2) have been described in the range of 4.37 to 11.34 in the river water, associated with a terrestrial/allocthonous DOM that has a higher aromatic carbon content, associated with the presence of substances such as tannins and/or humic compounds derived from higher plants and organic matter from the soil [39, 40]. The result obtained for the DOM from point 3 is within this range, but not so for points 1 and 2, whose average value was 25, indicating that the DOM is derived from a microbial source [41]. In relation to the extracts fractionated by molecular weight, F2 had higher 254/436 ratios than F1 for the three sampling points.

The 254/204 and 250/365 absorbance ratios have been reported to be useful for the characterization of the DOM, whose 254/204 ratio allows the correlation with the amount of aromatic carbon present in the DOM [42, 43], and the 250/365 ratio is used as an indicator of molecular size. The results show that the 254/204 values (Table 2) were similar for the DOM from the three sampling points, but not so for the different molecular weight fractions, where the ratio is greater in F1 and F2 from sampling point 2.

As to the 250/365 ratio (Table 2), the values obtained for the different molecular weight fractions coincide, pointing to the efficiency of the fractionating procedure with the ultrafiltration cell. The 250/365 results of the DOM indicate that the one from sampling point 1 has a lower molecular weight, followed by the DOM from sampling points 3 and 2.

3.5. Analysis of the DOM and Its Fractions by DSA-TOF-MS

The results of the DOM analysis by DSA-TOF- MS are presented in Figures 4 to 8. The calibration results are presented in Table 3, where the residues obtained were lower than 0.0006. The volume of calibrant solution used in the determinations was 20 μL to increase the reproducibility of the results.

Studies have been conducted to characterize natural organic matter using electrospray ionization mass spectrometry; with this technique only the polar compounds are ionized, while purely aliphatic moiety in humic and fulvic acids is not detected, and Kujawinski et al. [23] have suggested that compounds with high carboxylic acid content may be susceptible to fragmentation during electrospray ionization.

The DSA is a type of ambient ionization source. Ambient mass spectrometry is classified by the ability to sample and ionize analytes directly from surfaces with little to no sample preparation.

The results of the standards of HA and LG, analyzed by DSA-TOF-MS, showed differences in the distribution of signals (Figures 4 and 5), which can be attributed to the higher number, size, and structural complexity of the masses of DOM obtained from HA.

In relation to the mass spectra obtained for the DOM of Bio-Bio River, the fragmentation pattern varies according to the origin of the DOM, where for sampling point 1 (Figure 6), of mainly natural origin, a spectrum similar to the pattern obtained for the masses of LG can be observed unlike the mass spectra of the sample points 2 and 3 (Figures 7 and 8), in which the distribution of the masses is more similar to that observed for the HA. This change in the DOM would reflect the greatest human activity that occurs in these sampling points. The importance of the humic fraction is because it has a high influence on the physicochemical characteristics of natural water systems because it increases the light attenuation and maintains the pH constant by the presence of organic acids that act as buffer and because it has the ability to form complexes with some elements, affecting the redox chemistry of trace metals [44].

A water system is a complex system, where the presence of organic matter can even modify some physicochemical properties such as electrical conductivity or pH, but the content of organic matter does not necessarily correlate with the amount of DOM, since it will depend on the origin of organic matter. There is evidence that the concentration of dissolved organic carbon and pH are closely correlated with the molecular composition of the DOM in different types of aquatic systems. In this regard, Roth et al. [45] indicate that a higher pH of the molecular formulae of DOM had a lower degree of unsaturation and oxygenation, with a lower molecular size and greater abundance of nitrogen containing compounds. This may be because, at a higher pH, a greater biological activity related to the greater abundance of nitrogen containing compounds; however at a lower pH there would be a greater abundance of tannins in the system, which can inhibit the biological breakdown of the DOM.

4. Conclusions

From the results it can be concluded that the DOM present in natural water is composed of macromolecules with a molecular weight range, chemical composition, and a specified spectroscopic profile that are similar for the three sampling points.

From the fractionation by molecular weight it was possible to determine for the three sampling points that the DOM has mainly a molecular weight greater than 1 kDa, with a hydrophobic character, where the molecular size is greater for point 2, a sector with a high anthropic impact, followed by point 3, and finally point 1.

The effect of human activity is reflected in increased molecular size of the DOM, as well as increased organic carbon content and conductivity.

From the spectroscopic results it can be concluded that the DOM of the Rucalhue sector (point 1) has a terrestrial origin, with a high aromatic carbon content, distinguishing it from the DOM of the Laja and Concepción sectors (points 2 and 3), which have a greater anthropic impact.

The results of the DSA-TOF-MS show mass spectra of DOM points 2 and 3 of mainly anthropogenic origin were similar to that obtained for the deviation AH, unlike that DOM of mainly natural origin of point 1 whose mass spectrum is consistent with that obtained for the standard LG. This shows that this technique offers advantages in structural characterization of DOM mainly because of its short time of analysis and ease of sample preparation.

Competing Interests

The authors declare that they have no competing interests.


This study was supported by National Fund for Scientific and Technological Development (Fondecyt), through the Project Fondecyt no. 11110223 and the Project Fondequip no. EQM-130119.

Supplementary Materials

In the section Supplementary Material H1-NMR and FTIR spectra of the DOM obtained from the three sampling points are included.

  1. Supplementary Material