Production of fluorescent dissolved organic matter in arctic ocean sediments

ABSTRACT Little is known about the production of fluorescent dissolved organic matter (FDOM) in the anoxic oceanic sediments. In this study, sediment pore waters were sampled from four

different sites in the Chukchi-East Siberian Seas area to examine the bulk dissolved organic carbon (DOC) and their optical properties. The production of FDOM, coupled with the increase of

nutrients, was observed above the sulfate-methane-transition-zone (SMTZ). The presence of FDOM was concurrent with sulfate reduction and increased alkalinity (R2 > 0.96, _p_ < 0.0001),

suggesting a link to organic matter degradation. This inference was supported by the positive correlation (R2 > 0.95, _p_ < 0.0001) between the net production of FDOM and the modeled

degradation rates of particulate organic carbon sulfate reduction. The production of FDOM was more pronounced in a shallow shelf site S1 with a total net production ranging from 17.9 to 62.3

RU for different FDOM components above the SMTZ depth of _ca._ 4.1 mbsf, which presumably underwent more accumulation of particulate organic matter than the other three deeper sites. The

sediments were generally found to be the sources of CDOM and FDOM to the overlying water column, unearthing a channel of generally bio-refractory and pre-aged DOM to the oceans. SIMILAR

CONTENT BEING VIEWED BY OTHERS DOM DEGRADATION BY LIGHT AND MICROBES ALONG THE YUKON RIVER-COASTAL OCEAN CONTINUUM Article Open access 13 May 2021 SUBMARINE MUD VOLCANOES AS A SOURCE OF

CHROMOPHORIC DISSOLVED ORGANIC MATTER TO THE DEEP WATERS OF THE GULF OF CÁDIZ Article Open access 05 February 2021 FACTORS CONTROLLING THE DISTRIBUTIONS OF DISSOLVED ORGANIC MATTER IN THE

EAST CHINA SEA DURING SUMMER Article Open access 16 July 2020 INTRODUCTION The Arctic Ocean is intimately associated with global carbon and hydrological cycling as well as climate change. A

large stock of organic carbon (1100–1500 Pg C), including permafrost, is contained in its drainage basin area, which is sensitive to climate change1. The Arctic Ocean also influences

formation of North Atlantic Deep Water which drives Atlantic thermohaline circulation. The fluvial discharge into the Arctic Ocean is reported to be 25–36 Tg C yr−1 for dissolved organic

carbon (DOC) and around 12 Tg C yr−1 for particulate organic carbon (POC)1,2,3. It was estimated that roughly over 80–90% of the organic carbon, carried by fluvial runoff, is buried in

shallow coastal sediments4. Despite the large fraction of the organic carbon contained in the global oceanic sediments (_ca._ 700 Pg C)5, much is unknown about the dynamics of sediment pore

water dissolved organic matter (PW-DOM). It was previously reported that sediment oxygen consumption rates were much faster in coastal sediments ( > 200 mM O2 m−2d−1) than in deep sea

sediments (_ca._ 0.02 mM O2 m−2d−1)6, implying a higher biological respiration occurring in coastal sediments. The huge difference of more than four orders of magnitude of oxygen consumption

rates may reflect the huge differences of organic matter supply. With the prediction of ice-free summer in the Arctic Ocean over the next few decades7, it is urgent to understand the

dynamics of sediment DOM in this area as increasing organic carbon loading via fluvial runoff can be expected with global warming underway. Whether the organic carbon burial rates will go up

or down, however, is unknown at present. Although the origins of sediment organic matter are diverse including riverine runoff, coastal erosion, sea ice, aeolian input, and autochthonous

primary productivity in the Arctic8, a substantial fraction is known to be of terrigenous origin9,10,11,12. While solar irradiation is usually the major sink for terrigenous chromophoric DOM

(CDOM), autochthonous organic matter produced by surface primary productivity is generally considered susceptible to microbial utilization13,14,15. Meanwhile, the majority (_ca._ 80%) of

the organic matter sinking to seafloor is channeled into benthic respiration16. Once settling into sediments, the remainder organic matter undergoes a complex combination of physical (e.g.,

adsorption), chemical (e.g., redox reactions), and biological (e.g., aerobic and anaerobic respiration) diagenetic processes. The ultimate preservation rate is estimated to be <0.5% of

the gross production of primary productivity17,18,19. Although there have been a number of DOM studies on Arctic Ocean seawaters and arctic rivers12,20,21,22,23,24,25,26, little attention

has been paid to Arctic sediment. A recent report on PW-DOM from surficial sediments in the Eurasian Basin has shown the relationships between the molecular signature and sea-ice coverage,

which illustrated a higher abundance of peptide and aliphatic molecular formula in productive ice margins as compared to the molecular features (e.g., more aromatic) found in multiyear

ice-covered stations27. Other related studies primarily used biomarkers or stable isotopes simply to distinguish between terrigenous and marine sources and found a significant fraction of

terrigenous organic matter in the Arctic Ocean9,10,11,28,29,30. Studies of Arctic PW-DOM are essential in exploring the DOM characteristics and dynamics in the sediments as well as potential

interactions with the overlying water column. Optically active DOM is widely used to trace the dynamics of DOM in aquatic environments. These DOM pools (i.e., CDOM and FDOM) from

terrestrial-derived humic substances in surface waters are generally observed to be bio-refractory while some fresh biomass-derived fractions could be readily utilized13,14. The production

of FDOM in oxic or suboxic ocean interior was previously reported, which was explained by biological oxidation of sinking particulate organic matter (POM) based on a positive relationship

between FDOM and AOU (apparent oxygen utilization)31. These previous findings drive the investigation of CDOM and FDOM in oceanic sediment pore waters with a hope that it could provide

valuable information regarding organic carbon storage and the potential benthic efflux into the ocean water column. Investigation of this potential production and efflux of FDOM in oceanic

sediments can help to reconcile the reactivity of dissolved organic matter and its age in the deep oceans. Excitation-emission matrix coupled with parallel factor analysis (EEM-PARAFAC) is

an ideal method for FDOM studies because it can provide quantitative information on decomposed fluorescent components from three-dimensional fluorescence spectra of DOM32. The optical method

has proven successfully utilized in examining PW-DOM in both oceanic and freshwater ecosystems14,33,34,35,36. The objectives of this study were (1) to characterize the Arctic sediment

PW-DOM with different depths (up to 10.5 meters below seafloor (mbsf)) in the Chukchi-East Siberia Sea area along a shelf-slope-basin gradient, (2) to infer the potential origins, diagenetic

transformations, and the paleoenvironments associated with the DOM variations, and (3) to explore the potential production of FDOM in the sediments and benthic flux of DOM in the

water-sediment interface. The study area is of interest for this topic in that (1) there is a seasonal sea ice cover in the region so it is very sensitive to climate change, (2) it is

currently affected by nutrient-rich Pacific water inflow through the Bering Strait so it is easily affected by high primary productivity, which is known to be an important source of organic

matter in this region, and (3) the marginal seas are characterized by high sedimentation rates. Sediment cores roughly reflect the geological time spanning up to mid- to late Quaternary

deposit based on the age-depth models at sites in close proximity37,38. The major notable paleoenvironmental events occurring during the geological period include late Quaternary

glacial/interglacial cycles36 and the reopening of the Bering Strait during the Early Holocene at _ca._ 10.5 ka39,40. The opening of the Bering Strait is expected to deliver more nutrients

and thus promote the primary productivity in this region. RESULTS DOWNCORE PROFILES OF WATER CHEMISTRY, DOC, CDOM, AND FDOM The cores from sites S1, S2, S3, and S4 are 6.84, 3.54, 4.92, and

10.50 mbsf, respectively (Fig. 1). The high resolution sub-bottom profiles (SBP) of the sites are shown in the supporting information (Fig. S1). As seen from the SBP, there is no seismic

vent at the study sites. The salinity averaged _ca._ 33.7 to 36.8 psu among the four sites (Tables S1 and S6). A sharp increase of alkalinity in contrast to a sharp decrease of sulfate was

observed at shelf site S1. Alkalinity at site S1 increased from seafloor at 2.6 mM to 4.3 mbsf at 43.3 mM, and it remained a relatively constant below the depth in contrast with the sulfate

concentration decreasing from the seafloor at 26.0 mM to about 4.1 mbsf at 0.0 mM (Fig. S2). The alkalinity is generated and sulfate is consumed via particulate organic matter sulfate

reduction (POCSR: 2CH2O + SO42− → H2S + 2HCO3−) and anaerobic oxidation of methane (AOM: CH4 + SO42− → HS− + HCO3− + H2O) reactions in the sediments. The rapid changes of both alkalinity and

sulfate signify the fast POCSR and AOM reactions which might be due to large amount of decomposable organic matter inputs to the shallow shelf site S1. The sulfate downcore profiles

revealed that the sulfate-methane-transition-zone (SMTZ) is about 4.1 mbsf at site S1, while the other sites had not reached SMTZ yet up to the studied depths (Table S1). The bulk DOC data

for pore waters at sites S1, S2, S3, and S4 exhibited the average concentrations of 6.6 ± 2.4, 4.4 ± 1.9, 2.5 ± 0.8, and 8.4 ± 3.6 mM, respectively. These DOC values were 1–2 orders of

magnitude higher than those reported for the surficial sediments (0–10 cm) from the Amundsen and Nansen Basins located in the Eurasian Basin (0.04–0.35 mM)27. For this study, the DOC values

at the similar surficial sediments with the depths of 0–10 cm were even on the same order of magnitude as those of the deeper cores (Fig. 2). Thus, a much higher DOC appears to be stored in

this study area relative to the Amundsen and Nansen Basins probably due to the high _in situ_ primary productivity (up to 400 g C m−2 yr−1)41 and/or the high input of terrestrial organic

matter. The DOC depth profiles were highly variable among the sites and the core depth, implying the interplay of complex factors (e.g., sea-ice coverage, seasonality of riverine runoff,

ice-rafted debris sediments, nutrients availability, switch of ocean currents, reopening of the Bering Strait) during the late Quaternary glacial/interglacial cycles. The bulk DOC may be

more affected by the above mentioned factors because it encompasses a broader DOM pool than CDOM and FDOM, and thus its reactivity range should also be broader. In contrast, the CDOM,

expressed as absorption coefficient _a_(254) and _a_(350), displayed more obvious trends, especially for site S1, with a general increase followed by a gradual decrease with the core depth.

Meanwhile, an apparent inflection point was observed at the shelf site S1, which coincides with the SMTZ depth (_ca._ 4.1 mbsf) inferred from the sulfate profile. In this study, four

different fluorescent components (C1–C4) were identified from the EEM-PARAFAC modeling (Fig. 3 and Fig. S3). From the spectral features, C1 (Ex/Em maxima: 265/422 nm) was assigned as a

traditionally termed terrestrial humic-like component with Ex maximum in an ultraviolet range. C2 (Ex/Em maxima: 318/410 nm), C3 (Ex/Em maxima: 275(370)/452 nm), and C4 (Ex/Em maxima:

280/308 nm) were referred to marine/microbial humic-like, terrestrial humic-like, and protein-like components, respectively13,14,32. Similar to CDOM, the fluorescent components at site S1

showed an increasing trend with depth down to an obvious turning point around the SMTZ depth (_ca._ 4.1 mbsf) followed by a gradual decrease to the deeper sediments (Fig. 3). VARIATIONS OF

DOM PARAMETERS ALONG THE SHELF-SLOPE-BASIN GRADIENT The average bulk DOC values were 6.6 ± 2.4 mM at the Chukchi Shelf site S1, which may be explained that the shallower shelf area usually

traps more terrestrial organic matter, resulting in a higher organic carbon burial4,8. At the Chukchi Basin site S4, the average DOC values were 8.4 ± 3.6 mM. Since site S4 is a basin site,

the main DOM sources could be associated with _in situ_ marine phytoplankton productivity. This explanation is supported by the low aromaticity proxy (i.e., _a_(254)*) at the site. The

possible causes for the high phytoplankton productivity in this area will be discussed in the following section. The East Siberia Continental Shelf site S3 demonstrated a relatively lower

DOC level (2.5 ± 0.8 mM, _p_ < 0.01) among the four studied sites, implying that other factors such as autochthonous production, sea ice cover, and ocean currents rather than riverine

runoff could be dominantly operative for the carbon preservation. It is notable that even the labile phytoplankton-derived DOM could be preserved in sediments in the case of a fast

sedimentation (typical in inner shelf), in which the fresh organic materials could escape the respiration and preserved in the sediments4,11. For the CDOM, shelf site S1 showed much higher

absorption coefficients (_p_ < 0.0001) of _a_(254) (211 ± 78 m−1) and _a_(350) (53 ± 18 m−1), as compared to the other sites that displayed an order of magnitude lower values of _a_(254)

(12–23 m−1) and _a_(350) (3–9 m−1). These results are interesting in that this site is relatively remote from a large river discharge. Similar to the CDOM patterns, the absolute abundances

of the fluorescent components showed the highest level at site S1 (_p_ < 0.0001). CORRELATIONS AMONG DOM VARIABLES AND WATER CHEMISTRY PARAMETERS Water chemistry could affect the optical

properties of DOM42,43. To assess the influence, Pearson’s correlations were performed among the DOM variables and measured water chemistry parameters (Table S3). No significant correlations

were found between DOC and CDOM or FDOM, suggesting either the uncoupled production mechanisms of DOC vs. optically active DOM or different reactivity between bulk DOC and CDOM or FDOM.

Instead, strong positive correlations (R2 > 0.95, _p_ < 0.0001) were observed for either CDOM or FDOM with nutrients (i.e., NH4+ and PO43−, Figs S4–S5), suggesting the existence of the

potential coupling processes to regulate CDOM, FDOM, and nutrient cycling in the studied sediments. The trends were in line with the positive correlations of FDOM with alkalinity, and the

negative correlations with sulfate (R2 > 0.96, _p_ < 0.0001, n = 39, Fig. 4a,b, Table S3). DISCUSSION The origins of PW-DOM in the study area could be diverse from riverine runoff,

coastal erosion, aeolian inputs, and ice-rafted debris to autochthonous production. The riverine inputs have a strong seasonality although they occur mainly during spring freshet in modern

conditions8. The relative strength of each source likely depends on specific depositional times and locations, and subsequent diagenetic processes. The high variations of DOC observed in

this area may be attributed to the Quaternary glacial/interglacial periods corresponding to the core depth37. Moreover, the Chukchi Borderland and the East-Siberian continental margin were

glaciated by more than 1000 m-thick ice during several Quaternary glaciations prior to the Last Glacial Maximum (Marine Isotope Stage 2)44,45. Grounded ice can significantly contribute to

reworking of older terrestrial and marine sediments. The sea ice coverage regimes are completely different during the glacial versus the interglacial periods, affecting the entire hydrology

and primary productivity. The optical indices of humification index (HIX), biological index (BIX), fluorescence index (FI), specific absorption coefficient _a_(254)*, and slope ratio (SR)

were utilized here to trace the origins and characteristics of the optically active DOM. Detailed definitions and descriptions of these proxies are found elsewhere46,47,48,49,50. The FI

(>1.7) and HIX (<6) suggested a primarily microbial/marine-derived CDOM and FDOM origins with the exception of site S2 (Fig. 2 and Fig. S6). The bulk elemental results of Corg/Norg

were reported to be generally <10 in the Chukchi Sea region38, consistent with the optical data of this study. The Northwind Basin site S2 seems to reflect the optical signals of

potential terrigenous DOM inputs of fulvic acids with relatively low FI (<1.5). This site is close enough to the Canadian Basin to be affected by the discharge of Mackenzie River through

the Beaufort Gyre. Episodic ice-rafted debris inputs as well as reworking by grounded ice on the Northwind Ridge and Chukchi Plateau during glacial times could also be the explanations for

the exceptional results45. However, aromaticity index of _a_(254)* was relatively low (0.6 ± 0.5 L(mgC-m)−1) and the protein-like C4 (could be a mixture of proteinaceous materials and

phenolic moieties)51 was high (2.1 ± 1.9 RU), signifying its mainly marine/microbial-derived PW-DOM source. The highest average DOC values observed at site S4 can be ascribed to the lateral

transport of both nutrients and DOM from the Chukchi Shelf to the Chukchi Basin52, which may stimulate the primary productivity followed by the increase of freshly produced biomass. This

explanation is well supported by our results of the extremely low aromaticity index _a_(254)* (0.1 ± 0.1 L(mgC-m)−1) and the high protein-like C4% (72 ± 9%) at site S4. Although the

potential linkage of the protein-like component with phenolic compounds was previously reported51, the possibility can be excluded here because of the extremely low aromaticity index

_a_(254)*. Isotopic and biomarker analyses indicated a significant fraction of terrigenous organic matter in the Arctic Ocean9,10,28. However, our results suggest that spatial variations

should be taken into account as an important factor. For the areas far away from direct large riverine influences, the DOM sources could be primarily microbial/marine-derived. A study using

polyunsaturated fatty acids and algal sterols has reported relatively unaltered marine-derived organic matter in the Chukchi shelf and slope area11. An organic carbon to nitrogen ratio

(i.e., Corg/Norg) of <10 at nearby sites also supports the dominant source of algal-derived organic matter in the Chukchi Shelf 29,53. Within the oceanic sediment column, sulfates from

overlying seawater or from re-oxidized sulfides operate as important electron acceptors above the SMTZ of sediments (typically about 1 to 100 mbsf), and they are involved in the two chemical

reactions of POCSR and AOM. The POCSR potentially produces DOM. POCSR was reported to dominate organic matter degradation above SMTZ in sediments without seismic vents or faults54. The

study sites here did not display seismic vents in the sub-bottom profiler (Fig. S1). As such, it is reasonable to assume POCSR as the dominant pathway of organic matter degradation at the

sediments above SMTZ for this study. This is the first report on the production of FDOM in the Arctic sediments. The net increases (=the value at the indicated depth minus the value at the

bottom water) of all the fluorescent components were found at all the studied sites, which are summarized in Table 1. The net increases of four EEM-PARAFAC components (C1 to C4) of 62.3,

25.7, 17.9, and 24.1 RU, respectively, were observed for the depth up to the SMTZ at site S1 (_ca._ 4.1 mbsf). For the studied depths at the sites S2, S3, and S4 (3.5 to 10.5 mbsf), the net

increases ranged from 0.2 to 4.5 RU for all EEM-PARAFAC components. The FDOM increases with depth should not be due to the diffusive loss to the water column giving the huge differences of

depth profiles between the shelf and basin sites. The quantity of the FDOM production was two orders of magnitude higher at site S1 than at the other sites. The shallower shelf site has a

much shallower water depth of 100 m in comparison with the other slope and basin sites with much deeper water depths (S2: 2077 m; S3: 715 m; S4: 2240 m, Table S1). Furthermore, the C1 FDOM

component at site S1 showed the highest production (62.3 RU) among all the fluorescent components and all the sites. This component has been reported to be a photo-refractory and/or

photoproduct in a coastal wetland ecosystem14,55,56, and it was also found in a high abundance in canal water draining agricultural areas with “old” historic peat-derived DOM57. It is

notable that the above-mentioned environments could be similar to the case of the “old” permafrost-derived organic matter immobilized from the surrounding continents via the coastal erosion

and fluvial runoff to the Arctic sediments. Although DOM optical signatures were primarily marine/microbial-derived in this study area as discussed above, the biogenic particles may carry

terrigenous organic matter, which could be transformed to DOM during diagenesis through biological oxidation and re-working. In this study, the production of FDOM was accompanied by the

increases of nutrients of ammonium and phosphate (i.e., NH4+ and PO43−), concurrent with increased alkalinity and reduced sulfate (SO42−, R2 > 0.96, _p_ < 0.0001, n = 39, Fig. 4a–c,

Figs S4–S5), suggesting the presence of a biological pathway for the production via anaerobic respiration rather than physical dissolution or fragmentation. If FDOM was produced simply by

fragmentation or physical dissolution, such strong correlations with the consumption of sulfate, releases of nutrients, and increase of alkalinity would not be observed. Benefiting from the

advantages of EEM-PARAFAC to decompose fluorescence EEMs into the four individual components, we found that the most dominant FDOM component produced was the terrestrial humic-like C1, known

as photoproducts and/or photo-refractory substances, which is potentially associated with the permafrost-derived organic matter. Decreasing HIX with core depth, concomitant with the

increases of BIX and SR values, indicate increased biological activities and a lower molecular weight of CDOM with core depth (Fig. S6), supporting our inference on the origins of FDOM and

the subsequent diagenesis in the sediments. Moreover, positive correlations between EEM-PARAFAC components and ammonium concentration and between net FDOM production and the modeled POCSR

rates (see supporting information for the details of POCSR rates modeling) were observed (R2 > 0.95, _p_ < 0.0001, n = 4, Fig. 4d), corroborating the FDOM production via anaerobic

degradation of particulate organic matter. Although n is only four here, the correlation is strong as seen from the high R2 (>0.95) as well as significant test value of _p_ < 0.0001

for all four EEM-PARAFAC components C1–C4. The correlation remains strong (R2 > 0.87, _p_ < 0.001, n = 3) for most of the components (C1 and C3) at the relatively clustered sites S2,

S3, and S4. Our observation is also in line with a previous report on the assimilation of particulate organic carbon by archaea above SMTZ in sediments58. As the rates of POCSR were

simulated with two extreme scenarios of input parameters (i.e., the highest and lowest estimations, see supporting information), we expect that the potential uncertainties involved in POCSR

rates modeling would be marginal. The apparent production of FDOM above the SMTZ in the Arctic sediments was observed, which was associated with the organic matter degradation in the

sediments as evidenced by the positive correlation with nutrients as well as the POCSR rates, and by the inverse correlation with sulfates. The subsequent flux estimation based on the Fick’s

first law of diffusion (Diffusion flux J = sediment porosity at the water-sediment interface × diffusion coefficient × concentration gradient, i.e., J = ɸoDs(əC/əz)o) showed that most of

the sites served as the sources of CDOM and FDOM (primarily for shelf site S1) to the overlying water column (see Supporting Information, Table S5). The continental shelf receiving a large

terrestrial organic matter loading seems to be a hotspot of FDOM production in the sediments and subsequent efflux to the water column. In contrast, the slope and basin sites S2, S3, and S4

illustrated much lower FDOM production and therefore lower FDOM benthic flux. Some exceptions did occur with potential influx of the protein-like component, presumably due to autumn

phytoplankton bloom in the water column. Our results provide further insight into the refractory nature of deep oceanic DOM, partially explaining the discrepancy between the average age of

4,000–6,000 yrs of oceanic DOM and the oceanic mixing time of about 1,000 yrs, in that humic-like FDOM are generally bio-refractory in nature13,14. The continental shelf appears to be an

important area for the benthic efflux of FDOM. Since the Arctic Ocean has the widest continental margin, occupying about one third of the global ocean shelf areas, our major findings

including the production of FDOM, and the benthic fluxes of FDOM, are worthy of being treated as important factors in constructing conceptual models of global carbon cycling. In this study,

EEM-PARAFAC revealed the potential co-existence of the efflux and the influx of different DOM moieties in Arctic sediments. In addition, there was high variability of DOM quantity and

quality with core depth and among different sites. METHODS SITE DESCRIPTION Sampling sites are located in the shallow Chukchi Shelf (site ARA06C-JPC-1a), Northwind Basin (site ARA06C-JPC-2),

East Siberia Continental Slope (site ARA06C-JPC-3), and Chukchi Basin (site ARA06C-JPC-4) of the western Arctic Ocean (Fig. 1). For simplicity, the four sampling sites are denoted here as

the sites S1, S2, S3, and S4, respectively. The locations and the site description are summarized in Table S1. The sub-bottom profiler (SBP) data (see supporting information and Fig. S1)

showed that the sediments are Holocene marine transgressive sediments (S1), hemi-pelagic sediments (S3), and mostly fine-grained homogeneous pelagic sediment (S2 and S4). The modern Chukchi

Sea is well known as seasonally ice-covered and having an expansive shallow continental shelf and extremely high primary productivity41 driven by nutrient-enriched Bering Shelf Anadyr Water

entraining through the Bering Strait. The Northwind Basin is separated from the Canadian Basin by the Northwind Ridge, in which Pacific Summer Water (PSW) was reported to cause the

disproportionally large recent ice retreat in this area59. The East Siberia Sea is supposed to receive a larger riverine discharge as compared to the Chukchi Sea, due to two large Siberian

rivers (i.e., Kolyma and Indigirka). The Chukchi Sea circulation is dominated by the Bering Sea Water (BSW, >80% of inflow, about 0.7 Sv) and the Alaskan Coastal Water (ACW)60. The BSW is

a mixture of Bering Shelf Anadyr Water (BSAW) enriched with nutrients (NO3− ≥ 20 μM) and Bering Shelf Water with a low level of nutrients (NO3− ≤ 1 μM)60. A stable stratification is

generally maintained in Arctic seawater due to a large riverine freshwater discharge and sea ice melting61. Water masses usually exist in several layers, including the low-salinity surface

mixed layer, a complex halocline intermediate water (50–200 m, Pacific-derived on the Bering Strait side while Atlantic-derived on the Fram Strait side), a warm and saline Atlantic Deep

Water, and a layer of bottom water (>900 m)8,62. The water residence times (WRTs) for seawaters are variable depending on the layers. In general, relatively shorter (1 to10 yrs) WRTs were

found for the surface and the halocline layers, while longer WRTs (25 to 30 yrs), for the Atlantic deep layer, and up to about 300 yrs for the bottom water8. Two major current systems

(i.e., the Beaufort Gyre and the Transpolar Drift) dominate surface water circulation in the Arctic Ocean. The Beaufort Gyre occupies most of the Amerasian Basin, in which our study sites

are located, while the Transpolar Drift is dominant in the Eurasian Basin. A switching between cyclonic and anticyclonic Beaufort Gyre was also observed63. SAMPLING Sediment cores were

collected at the sites using both Jumbo Piston Corer (JPC) and Box Corer from August to September 2015, during the RV ARAON Arctic Expedition (ARA06C). The surficial sediment porosity (Table

S2) was calculated from Wet Bulk Density (WBD) determined by whole-core logging64 using the undisturbed box-cores (maximum length of 50 cm) to avoid the potential alteration of the surface

sediments (usually centimeters) sampled by JPC. The overlying bottom water was carefully collected from undisturbed surface sediment cores. The pore waters in the sediment cores were slowly

extracted by Rhizon samplers. The bottom and pore waters were both collected in acid-prewashed syringes and filtered through an in-line 0.20 μm disposable polytetrafluoroethylene filter. The

Rhizon sampler was reported to extract relatively more humified FDOM, and lower DOC, CDOM, and protein-like FDOM as compared to the extraction based on centrifugation35. The aliquots were

transferred into acid-prewashed plastic bottles (Nalgene® high density polyethylene). The pore and bottom water samples for DOM analyses were both immediately stored in a freezer. The

samples for the onboard analyses and ion measurements were stored in a refrigerator (about 4 °C). ANALYTICAL MEASUREMENTS Salinity was determined by a Reflectometer, and alkalinity was

analyzed by titration using 0.1 N HCl onboard. The reproducibility of alkalinity was <2% based on the repeated analyses of the IAPSO (i.e., International Association for the Physical

Sciences of the Oceans) standard seawater. The PO43− and NH4+ were measured using colorimetric methods using a spectrophotometer (Shimadzu UV-2450) at 885 nm and 640 nm, respectively65.

Details regarding the measurements of other ions are described elsewhere66. The sediment porosity was measured from box-cores. DOC concentrations of pore waters were measured using a total

organic carbon analyzer (Shimazu TOC-VCPH) as non-purgeable organic carbon (NPOC) with an analytical reproducibility of <2%67. DOC of bottom water was not available due to limited sample

volume. Absorption spectra were scanned from 240 to 800 nm on a Shimadzu 1800 ultraviolet-visible (UV-Vis) spectrophotometer (Shimadzu Inc., Japan). Fluorescence EEMs were obtained using a

Hitachi F-7000 luminescence spectrometer (Hitachi Inc., Japan) at the excitation/emission (Ex/Em) wavelengths of 250–500/280–550 nm. The excitation and the emission scans were set at 5 nm

and 1 nm steps, respectively. All samples were transferred into a special capped cuvette in a nitrogen-filled glove box for the UV-Vis and EEMs measurements. Further details for the EEM

measurements and Raman Unit (RU) normalization are described elsewhere68. The procedures of post-acquisition corrections are available in a previous report56. PARAFAC modeling was performed

using the MATLAB7.0.4 with the DOMFluor toolbox69. All the corrected pore water (n = 39) and bottom water (n = 4) EEMs were used for modeling. The number of fluorescent components was

determined based on the split-half validation and an example of measured, modeled, and residual EEMs are shown in supporting information (Fig. S3). ADDITIONAL INFORMATION HOW TO CITE THIS

ARTICLE: Chen, M. _et al_. Production of fluorescent dissolved organic matter in Arctic Ocean sediments. _Sci. Rep._ 6, 39213; doi: 10.1038/srep39213 (2016). PUBLISHER'S NOTE: Springer

Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. REFERENCES * Rachold, V. et al. In The Organic Carbon Cycle in the Arctic Ocean

(eds Ruediger Stein & Robie W. MacDonald ) 33–55 (Springer Berlin Heidelberg, 2004). * Holmes, R. M. et al. Seasonal and Annual Fluxes of Nutrients and Organic Matter from Large Rivers

to the Arctic Ocean and Surrounding Seas. Estuaries and Coasts 35, 369–382 (2012). CAS  Google Scholar  * Raymond, P. A. et al. Flux and age of dissolved organic carbon exported to the

Arctic Ocean: A carbon isotopic study of the five largest arctic rivers. Global Biogeochemical Cycles 21, GB4011, doi: 10.1029/2007gb002934 (2007). Article  ADS  Google Scholar  * Hedges, J.

I. & Keil, R. G. Sedimentary organic matter preservation: an assessment and speculative synthesis. Marine Chemistry 49, 81–115 (1995). CAS  Google Scholar  * Houghton, R. A. Balancing

the Global Carbon Budget. Annual Review of Earth and Planetary Sciences 35, 313–347, doi: 10.1146/annurev.earth.35.031306.140057 (2007). Article  ADS  CAS  Google Scholar  * Andersson, J. H.

et al. Respiration patterns in the deep ocean. Geophysical Research Letters 31, L03304, doi: 10.1029/2003gl018756 (2004). Article  ADS  Google Scholar  * Wang, M. & Overland, J. E. A

sea ice free summer Arctic within 30 years? Geophysical Research Letters 36, L07502, doi: 10.1029/2009gl037820 (2009). Article  ADS  Google Scholar  * Stein, R. In Developments in marine

geology (ed. Chamley, H. ) (Elsevier, 2008). * Boucsein, B. & Stein, R. Particulate organic matter in surface sediments of the Laptev Sea (Arctic Ocean): application of maceral analysis

as organic-carbon-source indicator. Marine Geology 162, 573–586 (2000). ADS  CAS  Google Scholar  * Schubert, C. J. & Calvert, S. E. Nitrogen and carbon isotopic composition of marine

and terrestrial organic matter in Arctic Ocean sediments: implications for nutrient utilization and organic matter composition. Deep-Sea Research Part I-Oceanographic Research Papers 48,

789–810 (2001). ADS  CAS  Google Scholar  * Belicka, L. L., Macdonald, R. W., Yunker, M. B. & Harvey, H. R. The role of depositional regime on carbon transport and preservation in Arctic

Ocean sediments. Marine Chemistry 86, 65–88 (2004). CAS  Google Scholar  * Benner, R., Louchouarn, P. & Amon, R. M. W. Terrigenous dissolved organic matter in the Arctic Ocean and its

transport to surface and deep waters of the North Atlantic. Global Biogeochemical Cycles 19, GB2025, doi: 10.1029/2004gb002398 (2005). Article  ADS  Google Scholar  * Chen, M. & Jaffé,

R. Photo- and bio-reactivity patterns of dissolved organic matter from biomass and soil leachates and surface waters in a subtropical wetland. Water Research 61, 181–190 (2014). CAS  PubMed

  Google Scholar  * Chen, M. & Jaffé, R. Quantitative assessment of photo- and bio-reactivity of chromophoric and fluorescent dissolved organic matter from biomass and soil leachates and

from surface waters in a subtropical wetland. Biogeochemistry, 129, 273–289, doi: 10.1007/s10533-016-0231-7 (2016). Article  CAS  Google Scholar  * Vähätalo, A. V. & Wetzel, R. G.

Photochemical and microbial decomposition of chromophoric dissolved organic matter during long (months–years) exposures. Marine Chemistry 89, 313–326 (2004). Google Scholar  * Deming, J. W.

& Yager, P. L. In Deep-Sea Food Chains and the Global Carbon Cycle eds Gilbert T. Rowe & Vita Pariente ) 11–27 (Springer Netherlands, 1992). * Burdige, D. J. Preservation of Organic

Matter in Marine Sediments: Controls, Mechanisms, and an Imbalance in Sediment Organic Carbon Budgets? Chemical Reviews 107, 467–485 (2007). CAS  PubMed  Google Scholar  * Gudasz, C. et al.

Temperature-controlled organic carbon mineralization in lake sediments. Nature 466, 478–481 (2010). ADS  CAS  PubMed  Google Scholar  * Shimotori, K., Watanabe, K. & Hama, T.

Fluorescence characteristics of humic-like fluorescent dissolved organic matter produced by various taxa of marine bacteria. Aquatic Microbial Ecology 65, 249–260 (2012). Google Scholar  *

Opsahl, S., Benner, R. & Amon, R. M. W. Major flux of terrigenous dissolved organic matter through the Arctic Ocean. Limnology and Oceanography 44, 2017–2023 (1999). ADS  CAS  Google

Scholar  * Amon, R. M. W. & Benner, R. Combined neutral sugars as indicators of the diagenetic state of dissolved organic matter in the Arctic Ocean. Deep-Sea Research Part

I-Oceanographic Research Papers 50, 151–169 (2003). ADS  CAS  Google Scholar  * Stedmon, C. A., Amon, R. M. W., Rinehart, A. J. & Walker, S. A. The supply and characteristics of colored

dissolved organic matter (CDOM) in the Arctic Ocean: Pan Arctic trends and differences. Marine Chemistry 124, 108–118 (2011). CAS  Google Scholar  * Gueguen, C. et al. The nature of colored

dissolved organic matter in the southern Canada Basin and East Siberian Sea. Deep-Sea Research Part II-Topical Studies in Oceanography 81–84, 102–113 (2012). ADS  Google Scholar  * Walker,

S. A., Amon, R. M. W. & Stedmon, C. A. Variations in high-latitude riverine fluorescent dissolved organic matter: A comparison of large Arctic rivers. Journal of Geophysical

Research-Biogeosciences 118, 1689–1702 (2013). ADS  CAS  Google Scholar  * Cory, R. M., Ward, C. P., Crump, B. C. & Kling, G. W. Sunlight controls water column processing of carbon in

arctic fresh waters. Science 345, 925–928 (2014). ADS  CAS  PubMed  Google Scholar  * Mann, P. J. et al. Pan-arctic trends in terrestrial dissolved organic matter from optical measurements.

Frontiers in Earth Science 4, Article 25, doi: 10.3389/feart.2016. 00025 (2016). * Rossel, P. E., Bienhold, C., Boetius, A. & Dittmar, T. Dissolved organic matter in pore water of Arctic

Ocean sediments: Environmental influence on molecular composition. Organic Geochemistry 97, 41–52 (2016). CAS  Google Scholar  * Yunker, M. B., Macdonald, R. W., Cretney, W. J., Fowler, B.

R. & McLaughlin, F. A. Alkane, terpene and polycyclic aromatic hydrocarbon geochemistry of the Mackenzie River and Mackenzie shelf: Riverine contributions to Beaufort Sea coastal

sediment. Geochimica et Cosmochimica Acta 57, 3041–3061 (1993). ADS  CAS  Google Scholar  * Naidu, S. A. et al. Organic carbon isotope ratios (δ13 C) of Arctic Amerasian Continental shelf

sediments. International Journal of Earth Sciences 89, 522–532 (2000). ADS  CAS  Google Scholar  * Park, Y.-H. et al. Distribution, source and transportation of glycerol dialkyl glycerol

tetraethers in surface sediments from the western Arctic Ocean and the northern Bering Sea. Marine Chemistry 165, 10–24 (2014). ADS  CAS  Google Scholar  * Yamashita, Y. & Tanoue, E.

Production of bio-refractory fluorescent dissolved organic matter in the ocean interior. Nature Geoscience 1, 579–582 (2008). ADS  CAS  Google Scholar  * Stedmon, C. A., Markager, S. &

Bro, R. Tracing dissolved organic matter in aquatic environments using a new approach to fluorescence spectroscopy. Marine Chemistry 82, 239–254 (2003). CAS  Google Scholar  * Burdige, D. J.

& Komada, T. Sediment pore waters. 2nd edn, 536–569 (Academic Press, 2014). * Hur, J., Lee, B. M. & Shin, K. H. Spectroscopic characterization of dissolved organic matter isolates

from sediments and the association with phenanthrene binding affinity. Chemosphere 111, 450–457 (2014). ADS  CAS  PubMed  Google Scholar  * Chen, M., Lee, J.-H. & Hur, J. Effects of

sampling methods on the quantity and quality of dissolved organic matter in sediment pore waters as revealed by absorption and fluorescence spectroscopy. Environmental Science and Pollution

Research 19, 14841–14851 (2015). Google Scholar  * Chen, M. & Hur, J. Pre-treatments, characteristics, and biogeochemical dynamics of dissolved organic matter in sediments: A review.

Water Research 79, 10–25 (2015). CAS  PubMed  Google Scholar  * Stein, R. et al. Towards a better (litho-) stratigraphy and reconstruction of Quaternary paleoenvironment in the Amerasian

Basin (Arctic Ocean). Polarforschung 79, 97–121 (2010). Google Scholar  * Stein, R., Fahl, K., Schade, I., Manerung, A., Wassmuth, S., Niessen, F. & Nam, S.-I. Holocene variability in

sea ice cover, primary production, and Pacific-Water inflow and climate change in the Chukchi and East Siberian Seas (Arctic Ocean). _Journal of Quaternary Science_ (accepted manuscript in

press). * Elias, S. A., Short, S. K., Nelson, C. H. & Birks, H. H. Life and times of the Bering land bridge. Nature 382, 60–63 (1996). ADS  CAS  Google Scholar  * Elias, S. A., Short, S.

K. & Birks, H. H. Late Wisconsin environments of the Bering Land Bridge. Palaeogeography, Palaeoclimatology, Palaeoecology 136, 293–308 (1997). ADS  Google Scholar  * Sakshaug, E. In

The organic carbon cycle in the Arctic Ocean (eds Stein, R. & Macdonald, R. W. ) 57–82 (Springer-Verlag, 2004). * Pace, M. et al. pH change induces shifts in the size and light

absorption of dissolved organic matter. Biogeochemistry 108, 109–118 (2012). CAS  Google Scholar  * Romera-Castillo, C., Chen, M. L., Yamashita, Y. & Jaffe, R. Fluorescence

characteristics of size-fractionated dissolved organic matter: Implications for a molecular assembly based structure? Water Research 55, 40–51 (2014). CAS  PubMed  Google Scholar  * Niessen,

F. et al. Repeated Pleistocene glaciation of the East Siberian continental margin. Nature Geoscience 6, 842–846 (2013). ADS  CAS  Google Scholar  * Dove, D., Polyak, L. & Coakley, B.

Widespread, multi-source glacial erosion on the Chukchi margin, Arctic Ocean. Quaternary Science Reviews 92, 112–122 (2014). Google Scholar  * Zsolnay, A., Baigar, E., Jimenez, M., Steinweg,

B. & Saccomandi, F. Differentiating with fluorescence spectroscopy the sources of dissolved organic matter in soils subjected to drying. Chemosphere 38, 45–50 (1999). ADS  CAS  PubMed 

Google Scholar  * Huguet, A. et al. Properties of fluorescent dissolved organic matter in the Gironde Estuary. Organic Geochemistry 40, 706–719 (2009). CAS  Google Scholar  * McKnight, D. M.

et al. Spectrofluorometric characterization of dissolved organic matter for indication of precursor organic material and aromaticity. Limnology and Oceanography 46, 38–48 (2001). ADS  CAS 

Google Scholar  * Weishaar, J. L. et al. Evaluation of Specific Ultraviolet Absorbance as an Indicator of the Chemical Composition and Reactivity of Dissolved Organic Carbon. Environmental

Science & Technology 37, 4702–4708 (2003). ADS  CAS  Google Scholar  * Helms, J. R. et al. Absorption spectral slopes and slope ratios as indicators of molecular weight, source, and

photobleaching of chromophoric dissolved organic matter. Limnology and Oceanography 53, 955–969 (2008). ADS  Google Scholar  * Maie, N., Scully, N. M., Pisani, O. & Jaffé, R. Composition

of a protein-like fluorophore of dissolved organic matter in coastal wetland and estuarine ecosystems. Water Research 41, 563–570 (2007). CAS  PubMed  Google Scholar  * Hioki, N. et al.

Laterally spreading iron, humic-like dissolved organic matter and nutrients in cold, dense subsurface water of the Arctic Ocean. Scientific Reports 4, Article number 6775, doi:

677510.1038/srep06775 (2014). * Frigstad, H., Andersen, T., Bellerby, R. G. J., Silyakova, A. & Hessen, D. O. Variation in the seston C:N ratio of the Arctic Ocean and pan-Arctic

shelves. Journal of Marine Systems 129, 214–223 (2014). ADS  Google Scholar  * Hong, W.-L., Torres, M. E., Kim, J.-H., Choi, J. & Bahk, J.-J. Towards quantifying the reaction network

around the sulfate–methane-transition-zone in the Ulleung Basin, East Sea, with a kinetic modeling approach. Geochimica et Cosmochimica Acta 140, 127–141 (2014). ADS  CAS  Google Scholar  *

Stedmon, C. A. et al. Photochemical production of ammonium and transformation of dissolved organic matter in the Baltic Sea. Marine Chemistry 104, 227–240 (2007). CAS  Google Scholar  *

Chen, M. L., Price, R. M., Yamashita, Y. & Jaffe, R. Comparative study of dissolved organic matter from groundwater and surface water in the Florida coastal Everglades using

multi-dimensional spectrofluorometry combined with multivariate statistics. Applied Geochemistry 25, 872–880 (2010). Google Scholar  * Yamashita, Y., Scinto, L. J., Maie, N. & Jaffe, R.

Dissolved Organic Matter Characteristics Across a Subtropical Wetland’s Landscape: Application of Optical Properties in the Assessment of Environmental Dynamics. Ecosystems 13, 1006–1019

(2010). CAS  Google Scholar  * Biddle, J. F. et al. Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru. Proceedings of the National Academy of Sciences of the United

States of America 103, 3846–3851 (2006). ADS  CAS  PubMed  PubMed Central  Google Scholar  * Sumata, H. & Shimada, K. Northward transport of pacific summer water along the Northwind

Ridge in the Western Arctic Ocean. Journal of Oceanography 63, 363–378 (2007). Google Scholar  * Grebmeier, J. M., McRoy, C. P. & Feder, H. M. Pelagic-benthic coupling on the shelf of

the northern Bering and Chukchi Seas. I. Food supply source and benthic biomass. MARINE ECOLOGY - PROGRESS SERIES 48, 57–67 (1988). ADS  Google Scholar  * Cauwet, G. & Sidorov, I. The

biogeochemistry of Lena River: organic carbon and nutrients distribution. Marine Chemistry 53, 211–227 (1996). CAS  Google Scholar  * Anderson, L. G. & Amon, R. M. W. In Biogeochemistry

of marine dissolved organic matter (eds Hansell, D. A. & Carlson, C. ) 609–630 (Elsevier, 2014). * Morison, J. et al. Changing Arctic Ocean freshwater pathways. Nature 481, 66–70 (2012).

ADS  CAS  PubMed  Google Scholar  * Best, A. I. & Gunn, D. E. Calibration of marine sediment core loggers for quantitative acoustic impedance studies. Marine Geology 160, 137–146

(1999). ADS  Google Scholar  * Gieskes, J. M., Gamo, T. & Brumsack, H. Chemical Methods for Interstitial Water Analysis Aboard JOIDES Resolution. Ocean Drilling Program Technical Note

15, 1–60, http://www-odp.tamu.edu/publications/tnotes/tn15/technote15.pdf (Date of access: 19/10/2016) (1991). * Kim, J.-H. et al. Pore fluid chemistry from the Second Gas Hydrate Drilling

Expedition in the Ulleung Basin (UBGH2): Source, mechanisms and consequences of fluid freshening in the central part of the Ulleung Basin, East Sea. Marine and Petroleum Geology 47, 99–112

(2013). Google Scholar  * Hur, J. & Cho, J. Prediction of BOD, COD, and Total Nitrogen Concentrations in a Typical Urban River Using a Fluorescence Excitation-Emission Matrix with

PARAFAC and UV Absorption Indices. Sensors 12, 972–986 (2012). CAS  PubMed  Google Scholar  * Chen, M., Kim, S., Parka, J.-E., Kim, H. & Hur, J. Effects of dissolved organic matter (DOM)

sources and nature of solid extraction sorbent on recoverable DOM composition: Implication into potential lability of different compound groups. Anal. Bioanal. Chem. 408, 4809–4819 (2016).

CAS  PubMed  Google Scholar  * Stedmon, C. A. & Bro, R. Characterizing dissolved organic matter fluorescence with parallel factor analysis: a tutorial. Limnology and Oceanography-Methods

6, 572–579 (2008). CAS  Google Scholar  * Spreen, G., Kaleschke, L. & Heygster, G. Sea ice remote sensing using AMSR-E 89-GHz channels. Journal of Geophysical Research: Oceans 113,

C02S03, doi: 10.1029/2005jc003384 (2008). Article  Google Scholar  * Schlitzer, R. Ocean Data View, http://odv.awi.de (2016). Download references ACKNOWLEDGEMENTS This work was supported by

a National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (No. 2014R1A2A2A09049496), 2015 Arctic Expedition with RV Araon conducted through the Basic Research

Program (PE16062) of the Korea Polar Research Institute, and research program PACES-II (T3, WP3.2) of the Alfred Wegener Institute, and it was based on a preliminary survey for marine

energy resources in the Arctic region (GP2016-041). The authors thank the captain and crews of RV Araon for help with coring processing during the Arctic expedition in 2015 and Ms. Yun-Kyung

Lee for help with analytical measurements. W.L. Hong would like to acknowledge the support from the Research Council of Norway through its Centres of Excellence funding scheme, project

number 223259. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Department of Environment & Energy, Sejong University, 05006, Seoul, South Korea Meilian Chen & Jin Hur * Petroleum &

Marine Division, Korea Institute of Geoscience and Mineral Resources, 124 Gwahang-no, Yuseong-gu, 34132, Daejeon, South Korea Ji-Hoon Kim & Moo-Hee Kang * Division of Polar

Paleoenvironment Korea Polar Research Institute, Incheon, 21990, South Korea Seung-Il Nam * Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI), Am Alten Hafen 26,

Bremerhaven, 27568, Germany Frank Niessen * Department of Geology, CAGE - Centre for Arctic Gas Hydrate, Environment and Climate, UiT The Arctic University of Norway, Tromsø, Norway Wei-Li

Hong Authors * Meilian Chen View author publications You can also search for this author inPubMed Google Scholar * Ji-Hoon Kim View author publications You can also search for this author

inPubMed Google Scholar * Seung-Il Nam View author publications You can also search for this author inPubMed Google Scholar * Frank Niessen View author publications You can also search for

this author inPubMed Google Scholar * Wei-Li Hong View author publications You can also search for this author inPubMed Google Scholar * Moo-Hee Kang View author publications You can also

search for this author inPubMed Google Scholar * Jin Hur View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS J.H. contributed to the project

planning. J.H.K. and S.I.N. collected the samples and did the onboard analysis. M.C. analyzed the data and wrote the draft of the manuscript. F.N. carried out core logging and

sediment-porosity calculation. W.L.H. did the POCSR rates modeling. M.H.K. performed the sub-bottom profiler. All authors contribute to the final version of the manuscript. ETHICS

DECLARATIONS COMPETING INTERESTS The authors declare no competing financial interests. ELECTRONIC SUPPLEMENTARY MATERIAL SUPPLEMENTARY INFORMATION RIGHTS AND PERMISSIONS This work is

licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license,

unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce

the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Chen, M., Kim, JH., Nam, SI. _et

al._ Production of fluorescent dissolved organic matter in Arctic Ocean sediments. _Sci Rep_ 6, 39213 (2016). https://doi.org/10.1038/srep39213 Download citation * Received: 27 July 2016 *

Accepted: 21 November 2016 * Published: 16 December 2016 * DOI: https://doi.org/10.1038/srep39213 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this

content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative