A Semi-analytical Cdom Absorption Algorithm for Arctic Waters Printer-friendly Version Interactive Discussion Estimating Absorption Coefficients of Colored Dissolved Organic Matter (cdom) Using a Semi-analytical Algorithm for Southern Beaufort Sea (canadian Arctic) Waters: Application to Deriving Co

Biogeosciences Discussions This discussion paper is/has been under review for the journal Biogeosciences (BG). Please refer to the corresponding final paper in BG if available. Abstract A series of papers have suggested that freshwater discharge, including a large amount of dissolved organic matter (DOM), has increased since the middle of the 20th century. In this study, a semi-analytical algorithm for estimating light absorption coefficients of the colored fraction of DOM (CDOM) was developed for Southern Beaufort Sea wa-5 ters using remote sensing reflectance at six wavelengths in the visible spectral domain corresponding to MODIS ocean color sensor. This algorithm allows to separate colored detrital matter (CDM) into CDOM and non-algal particles (NAP) by determining NAP absorption using an empirical relationship between NAP absorption and particle backscattering coefficients. Evaluation using independent datasets, that were not 10 used for developing the algorithm, showed that CDOM absorption can be estimated accurately to within an uncertainty of 35 % and 50 % for oceanic and turbid waters, respectively. In situ measurements showed that dissolved organic carbon (DOC) concentrations were tightly correlated with CDOM absorption (r 2 = 0.97). By combining the CDOM absorption algorithm together with the DOC versus CDOM relationship, it 15 is now possible to estimate DOC concentrations in the near-surface layer of the Southern Beaufort Sea using satellite ocean color data. DOC concentrations in the surface waters were estimated using MODIS ocean color data, and the estimates showed reasonable values compared to in situ measurements. We propose a routine and near real-time method for deriving DOC concentrations from space, which may open the 20 way to an estimate of DOC budgets for Arctic coastal waters.

Estimating absorption coefficients of colored dissolved organic matter (CDOM) using a semi-analytical algorithm for Southern Beaufort Sea (Canadian Arctic) waters: application to deriving concentrations of dissolved organic carbon from space A. Matsuoka 1,2 , S. B. Hooker 3 , A. Bricaud1 , B. Gentili 1 , and M. Babin2

Introduction
The colored fraction of dissolved organic matter (CDOM) plays various roles in physical and biogeochemical processes (e.g.Carder et al., 1989;Miller and Moran, 1997;Nelson et al., 1998;Miller et al., 2002;Matsuoka et al., 2012).The contribution of CDOM light absorption to the total non-water absorption for Arctic waters is significantly higher Introduction

Conclusions References
Tables Figures

Back Close
Full than at lower latitudes (B élanger et al., 2006;Matsuoka et al., 2007Matsuoka et al., , 2009Matsuoka et al., , 2011)), playing a role as a heat absorption mechanism and thereby contributing to sea ice reduction (Hill, 2008;Matsuoka et al., 2011).Recent findings suggest that on land, global warming likely induces thawing of the permafrost containing a significant amount of dissolved organic carbon (DOC) (Freeman et al., 2001;Camill, 2005;Frey and Smith, 2005).Because the freshwater discharge into Arctic rivers has increased since the middle of the 20th century (Peterson et al., 2002;McClelland et al., 2006), a large amount of CDOM, whose absorption is often highly correlated with DOC concentrations at the mouth of a river (e.g.Mannino et al., 2008;Fichot and Benner, 2011;Matsuoka et al., 2012), is expected to be delivered into the Arctic Ocean by the river discharge relative to its autotrophic sources.During ice-free seasons, this material covers the surface ocean and significantly reduces penetration of ultraviolet irradiation that influences marine organisms such as phytoplankton and bacteria (Judd et al., 2007;Leu et al., 2007;Bonilla et al., 2009;Hessen et al., 2012).In the context of the expansion of hole in ozone layer at high northern latitudes (e.g.Wirth and Renger, 1996;Andersen and Knudsen, 2002), variability in CDOM absorption can impact the physiological status or acclimation/adaptation of marine organisms.Quantitative estimation and monitoring of CDOM in Arctic coastal environments are therefore urgently needed to better understand the modification in biogeochemical processes resulting from ongoing global warming.
Because of the similarity in absorption spectra between CDOM and non-algal particles (NAP), the two are often combined together as a single geophysical value, CDM (= CDOM + NAP) absorption (e.g.Garver and Siegel, 1997;Lee et al., 2002;Maritorena et al., 2002).Although basin-scale observations of CDM absorption using ocean color have provided insights into its linkages with phytoplankton biomass or the influence of photobleaching (Siegel et al., 2002;Bricaud et al., 2012), two important issues remain unsolved.First, these studies did not take into account variability in CDM absorption in coastal areas where a large amount of suspended sediments and dissolved organic matter is delivered by rivers.Because of the optical complexity of turbid waters, retrieved inherent optical properties (IOPs) tend to be more variable and thus uncertain Introduction

Conclusions References
Tables Figures

Back Close
Full  IOCCG, 2006;Kostadinov et al., 2007).In contrast to open oceanic waters in which CDOM absorption and DOC concentration are generally not correlated (Siegel et al., 2002;Nelson and Siegel, 2002), the tight correlation between these two variables for coastal waters can be used for deriving DOC concentrations from space when CDOM absorption is estimated using ocean color data.
The second issue related to the use of CDM absorption is that it is difficult to understand the different dynamics of CDOM and NAP separately.For example, particles sink rather quickly, at a vertical speed ranging from a few meters to several hundred meters per day (Ficher and Karakas, 2009), while CDOM can be transported by an ocean current over long distances (> 1000 km) in a few months (Matsuoka et al., 2011(Matsuoka et al., , 2012)).Consequently, partitioning CDM into CDOM and NAP using a semi-analytical algorithm is required, but it is difficult due to their similar spectral shapes (IOCCG, 2006).While recent studies developed local and empirical algorithms for estimating CDOM absorption (D'Sa, 2008;Mannino et al., 2008), the products derived using those algorithms can be significantly affected by spatial and seasonal variations in the empirical relationships.
The objective of this study is to develop a semi-analytical algorithm for estimating CDOM absorption for Arctic waters.The application for deriving DOC concentrations using ocean color data is proposed at the end of this paper.

Datasets for developing the CDOM absorption algorithm
The IOP datasets used for developing our semi-analytical CDOM absorption algorithm are documented in Matsuoka et al. (2007Matsuoka et al. ( , 2011)).Briefly, sampling was conducted in the Western Arctic Ocean in the area from approximately latitudes 65 • to 77 • N and longitudes 125 • to 180 • W (Fig. 1).Data were collected during three cruises from spring to Introduction

Conclusions References
Tables Figures

Back Close
Full  ).Phytoplankton pigments retained on the filters were extracted using methanol (Kishino et al., 1985), absorption coefficients of NAP (a NAP (λ), m −1 ) were measured on the bleached filter, and phytoplankton absorption coefficients (a ϕ (λ), m −1 ) were obtained by subtracting a NAP (λ) from a p (λ).For the absorption measurements of CDOM, water samples were filtered using 0.22 µm Millipore membranes immediately after water sampling and then poured into a 0.1 m quartz cell.Absorption coefficients of CDOM were determined from 280 to 700 nm using a MPS2400 spectrophotometer (Shimazu corp.).All absorption measurements were made following the NASA ocean color protocols (Mitchell et al., 2003).
In this study, chl a concentrations were determined fluorometrically with 90 % acetone (Holm-Hansen et al., 1965) for SBI spr and SBI sum and with dimethylformamide (DMF) (Holm-Hansen et al., 1965;Suzuki and Ishimaru, 1990) for MR aut using a 10-AU field fluorometer (Turner Designs).Suzuki and Ishimaru (1990) showed that chl a concentration obtained using DMF is not significantly different from that determined using 90 % acetone.Introduction

Conclusions References
Tables Figures

Back Close
Full

CDOM absorption and DOC determination
CDOM absorption datasets for evaluating our CDOM algorithm are documented in Matsuoka et al. (2012).Briefly, measurements were performed in the Southern Beaufort Sea during the France-Canada-USA joint Arctic campaign, MALINA, in the area from approximately latitudes 69 • to 72 • N and longitudes 125 • to 145 • W (Fig. 1).Sampling was conducted from 30 July to 26 August 2009 aboard the Icebreaker CCGS Amundsen, and 37 stations were visited.A barge and/or zodiac were also deployed (32 stations out of 37) during the same time periods of the CTD deployment to obtain surface water samples (see Figure A1 of Matsuoka et al., 2012).The absorption of CDOM in the sample filtrate was measured over the spectral range 200-735 nm using an UltraPath liquid waveguide system (World Precision Instruments, Inc.).A 2 m optical pathlength was used for almost all waters, except for coastal waters in the Mackenzie River mouth where a 0.1 m optical pathlength was used.The detailed methodology is described in Matsuoka et al. (2012).Dissolved organic carbon (DOC) was measured using the high-temperature combustion method and a TOC-V analyzer (Shimadzu) (Benner and Strom, 1993).Water samples using Whatman GF/F filters were stored frozen until analysis in laboratory.In total, 256 samples were collected for DOC analyses at 29 stations (DOC measurements were made by R. Benner).

Calculation of remote sensing reflectance
A detailed methodology for obtaining in-water upwelled radiance (L u ) and downward irradiance (E d ) is documented in Hooker et al. (2012).Briefly, a Compact-Optical profiling System (C-OPS) (Morrow et al., 2010) was deployed at 36 stations during the

Conclusions References
Tables Figures
Subsurface L u at null depth (i.e.L u (0 − , λ)) were obtained from the slope and intercept given by the least-squares linear regression of the log-transformed light parameter versus z.The principal data product used here is the remote sensing reflectance, R rs (λ) = 0.54L u (0 − , λ)/E s (λ), where λ indicates wavelength.In this study, R rs (λ) at six wavelengths (i.e, 412,443, 490, 532, 555, and 670 nm corresponding approximately to MODIS bands) were used.

Semi-analytical algorithms for deriving CDOM absorption and DOC concentration
In this study, we developed a semi-analytical algorithm for estimating CDOM absorption at 443 nm for Arctic waters.The overview of this algorithm is schematically shown in Fig. 2. Basically, we used a mathematical iteration method for minimizing the difference between measured R rs (λ) and R rs (λ) calculated using absorption and backscattering coefficients, as in the Garver-Siegel Maritorena (GSM) model (Garver and Siegel, 1997;Maritorena et al., 2002), but with 5 modifications as follows: (1) chl a-specific absorption coefficients of phytoplankton for Arctic waters (Matsuoka et al., 2011) were here ("the correction factor is generally not necessary" according to S. Maritorena, personal communication, 2011).Finally, DOC concentrations were estimated using an empirical relationship between DOC concentrations and CDOM absorption obtained in the Southern Beaufort Sea, i.e.DOC (µM) = 55 + 357 × a CDOM (443) (Matsuoka et al., 2012).

Ocean color data
Three-weeks composite images of R rs (λ) at 412, 443, 488, 531, 551, and 670 nm in August 2009 were generated using Aqua/MODIS ocean color data.After L1A data were downloaded from the ocean color website http://oceancolor.gsfc.nasa.gov/cgi/browse.pl?sen=am), geometric correction was made using the SeaWiFS Data Analysis System (SEADAS) software.Because the use of MODIS near-infrared bands (748 and 869 nm) for atmospheric correction (i.e.NASA standard atmospheric correction: Gordon and Wang, 1994;Gordon, 1997) often yielded pixels having negative values especially in R rs (412) in coastal areas of the Southern Beaufort Sea (which were flagged), the new method that can be applied to turbid ocean waters, proposed by Wang and Shi (2007), was used to obtain R rs (λ) images in this study.This algorithm basically operates in the same way as that of Gordon and Wang. (1994) but the two NIR bands (i.e.748 and 869 nm) are replaced by two SWIR bands (i.e.1240 and 2130 nm) for atmospheric correction.

Variability in R rs (λ)
Because we used R rs (λ) spectra to estimate CDOM absorption using our semianalytical algorithm (Sect.3.2), those spectra are presented first.There were a variety of R rs (λ) spectral shapes observed in this study from coastal to open ocean waters (Fig. 3a): R rs (λ) spectra in coastal waters had a peak at 555 nm or longer wavelength,

BGD Introduction Conclusions References
Tables Figures

Back Close
Full while those in offshore waters showed a peak at wavelengths ranging from 412 to 490 nm.These spectral shapes look similar to those observed everywhere at lower latitudes (see Werdell and Bailey, 2005 and references therein), but the IOPs for Arctic waters are significantly different from those for lower latitudes (Matsuoka et al., 2007(Matsuoka et al., , 2011)).A parameterization of IOPs specific to Arctic waters is, therefore, required to adequately derive geophysical values from R rs (λ) spectra (Matsuoka et al., 2011).In this study, Arctic-specific parameters are used for estimating CDOM absorption (Sect.3.2).
To generalize R rs (λ) spectral shapes, k-means cluster analysis was performed.The best number of classes (CN) was determined based on a minimum value of the Calinski-Harabasz index, which was found to be equal to 4 (Fig. 3b and Table 1).This index is one of the best methods to objectively obtain an optimal number of classes and may be applied to a variety of datasets (Calinski and Harabasz, 1974).Turbidity of waters tended to increase from cluster 1 to 4, and those clusters were roughly corresponding to those observed in the Chukchi Sea, except for cluster 4 (Matsuoka et al., 2007).Cluster 4 was observed especially in the Mackenzie river mouth.Similar clusters were also identified when a CN of 2 was set (grey solid and dashed lines for oceanic and turbid waters, respectively in Fig. 3b).It is noted that two clusters (oceanic and turbid) are enough for estimating CDOM absorption.Waters can then be classified automatically and systematically before our CDOM absorption algorithm is applied (Sect.4).

Evaluating our algorithm for the retrieval of CDOM absorption
In situ a CDOM (443) values for Southern Beaufort Sea waters varied widely from offshore to the coast (0.018 < a CDOM (443) < 1.08 m −1 ; see Fig. 4 of Matsuoka et al., 2012).In our observations, a CDOM (443) on average accounted for 72 ± 15 % of total non-water absorption, which is consistent with previous findings in the same study area (B élanger et al., 2006;Matsuoka et al., 2007Matsuoka et al., , 2009)).

BGD Introduction Conclusions References
Tables Figures

Back Close
Full The a CDOM (443) values estimated using our CDOM absorption algorithm (a mod CDOM (443), m −1 ) were evaluated against in situ measurements (Fig. 4).In this study, a sensitivity analysis of the algorithm was performed by varying η and S CDM values (i.e. the step 2 in the Fig. 2; see also Table 2), which are often set as constants (Lee et al., 2002;Maritorena et al., 2002); these can actually be parameterized as functions of R rs (see Appendix).Since CDOM absorption values were significantly different between CTD/Niskin and Barge samples (see Fig. A1 of Matsuoka et al., 2012), the performance of the algorithm was evaluated separately for each category (i.e.CTD/Niskin or Barg; CTD/Niskin and Barge correspond to oceanic and coastal waters, respectively).
For oceanic samples, the best performance was obtained when η and S CDM were set to be 1.0 and 0.0185, respectively (Fig. 4a).Excluding only one outlier (dotted circle), all data points fell within the ±35 % errors.The root mean squared error (RMSE) was low and the corresponding coefficient of determination was high (RMSE = 0.06 and r 2 = 0.98).It is important to note that this evaluation was conducted using independent datasets (i.e.MALINA datasets) that were not used for developing the algorithm (i.e. SBI spr, SBI sum, and MR aut datasets were used for this development).While parameterization of IOPs was made for typical oceanic waters (not directly influenced by river discharge) of the Arctic Ocean (Matsuoka et al., 2011), the evaluation results suggest this parameterization can also be applied for Arctic coastal waters directly influenced by the discharge of the Mackenzie River.This idea is partly supported by the fact that phytoplankton absorption properties, which were used for tuning the algorithm (see step 1 in Fig. 2), are similar between oceanic and coastal waters for the same season in the Arctic Ocean (Matsuoka et al., 2009).
A slope slightly lower than 1 was obtained using our algorithm (slope = 0.88; see Table 2 and Fig. 4a).This bias should be examined in a future study by collecting more samples over a wide range of CDOM content.Nonetheless, this semi-analytical algorithm provides estimates of CDOM absorption within a ±35 % accuracy for oceanic waters.Introduction

Conclusions References
Tables Figures

Back Close
Full For coastal samples, on the other hand, the CDOM algorithm with the fixed η and S CDM values did not perform well especially for high CDOM absorption values (RMSE = 0.18 and r 2 = 0.65; see Fig. 4b and Table 2).The slope of a mod CDOM (443) versus a obs CDOM (443) was high (i.e.1.51).The best performance was obtained for coastal samples when using fixed S CDM and varing η as a function of R rs (λ) (Fig. 4d; see also appendix A1).Using this algorithm, we were able to estimate CDOM absorption for turbid waters within ±50 % accuracy except for a few data points.This suggests that unlike S CDM , the choice of the η value impacts the performance of the CDOM algorithm for turbid waters (Fig. 4c).This result is consistent with the fact that optical properties in a thin layer at the top of the surface waters were characterized by high backscattering of various types of particles especially in the Mackenzie river mouth (Hooker et al., 2012;Doxaran et al., 2012), altering the η value.Due to scattering of particles as well as high CDOM absorption, the light penetration depth was very shallow in the Mackenzie River mouth (1.38 ± 1.82 m), but was deep in offshore areas (16.65 ± 9.84 m).Thus, η and S CDM parameters need to be set separately for oceanic and turbid waters.S CDM varies widely in global ocean basins (Bricaud et al., 2012), while this value tends to be constant in coastal areas (Babin et al., 2003;Magnuson et al., 2004).Our result may imply that waters in the Southern Beaufort Sea are strongly influenced by one given end-member, the Mackenzie river discharge, making S CDM stable (Matsuoka et al., 2012).The variability in this parameter needs to be examined further in Arctic waters.

Deriving DOC concentration using ocean color
In this section, we applied our CDOM absorption algorithm for Aqua/MODIS ocean color images.The three-week composite images in August 2009 were generated in order to correspond to the period of the MALINA cruise.
Overall, a sat CDOM (443) values were high in the western channel compared to those in the eastern channel of the Mackenzie delta.This result is consistent with our field observation (Fig. 5a; see Fig. 6 in Matsuoka et al., 2012).Geometric mean and geometric

Conclusions References
Tables Figures

Back Close
Full , respectively) were of the same order of magnitude compared to our observations for the MALINA cruise (0.055 and 2.265 m −1 , respectively) (Fig. 5b; Matsuoka et al., 2012).
In the Mackenzie river mouth, DOC concentrations (ranging 55 µM < DOC < 500 µM) in the surface layer were tightly correlated with a CDOM (443) (ranging 0.018 m −1 < a CDOM (443) < 1.08 m −1 ) (r 2 = 0.97; p < 0.0001; Matsuoka et al., 2012).Combining the CDOM absorption algorithm developed in this study with the DOC versus CDOM relationship published by Matsuoka et al. (2012), it is now possible to derive DOC concentrations in the surface layer of the Southern Beaufort Sea using satellite ocean color data.We applied this method to a CDOM (443) images for deriving surface DOC concentrations.Note that the DOC versus a CDOM (443) relationship we use is valid for waters from the Mackenzie river mouth to offshore waters near the ice edge area.
Consistently with a sat CDOM (443) values, high DOC concentrations (> 100 µM) were widely distributed in the Mackenzie shelf area (Fig. 6a).The geometric mean of DOC sat was 76 µM, which was similar to that observed for the MALINA cruise (Matsuoka et al., 2012).This result suggests that our estimation of DOC concentrations performs reasonably well, highlighting the potential utility of the algorithm for deriving DOC concentrations in the surface layer and further estimating the DOC budget for Arctic coastal waters.

Conclusions
In this study, we developed a semi-analytical algorithm for estimating CDOM absorption for Arctic waters.Comparison with in situ measurements showed that a CDOM (443) can be estimated accurately to within an uncertainty of 35 % and 50 % for oceanic and turbid waters of the Arctic Ocean, respectively.While we acknowledge that the estimation of a CDOM (443) in highly turbid waters is still a challenging task, to our knowledge, this is the first semi-analytical algorithm for estimating CDOM absorption, especially for Introduction

Conclusions References
Tables Figures

Back Close
Full Arctic waters.This approach allows a partitioning of CDM (i.e.CDOM plus NAP) into CDOM and NAP absorption contributions, thereby promoting an understanding of the dynamics for each component.As a result of combining the CDOM absorption algorithm together with the DOC versus CDOM relationship, we propose a routine and near real-time method for deriving DOC concentrations for Arctic coastal waters using satellite ocean color data (Fig. 7).First, satellite R rs (λ) spectra are categorized into two classes according to their spectral shapes.This can be made using supervised classification referring to clusters 1 and 2 obtained in this study (Fig. 3b).Second, parameterization of η and S CDM is set to be either constant or as a function of R rs (λ) depending on a R rs (λ) spectral class.Third, our CDOM absorption algorithm is applied.Finally, DOC concentrations are estimated using a relationship between DOC and a CDOM (443) derived from in situ measurements.Although this kind of method is potentially useful for estimating and monitoring change in DOC budgets for Arctic coastal waters as result of ongoing global warming, the DOC versus a CDOM (443) relationship was obtained empirically during summer in the Southern Beaufort Sea and might vary seasonally; a peak DOC flux occurs during ice melt in the early spring, significantly contributing to the DOC budget (Holmes et al., 2008).Further field observations are thus required to address this issue.where r rs (λ) represents below-surface R rs (λ) as calculated by R rs (λ)/0.5238(Gordon et al., 1988).Consequently, we decided to apply these equations to our CDOM algorithm for turbid waters.Detailed performance of the algorithm using those is shown in Fig. 4 and Table 2.  Full Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | in coastal areas compared to those in open ocean areas (e.g.
Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

(
used; (2) spectral slopes of CDM absorption, S CDM (m −1 ) and of b bp (λ), η (m −1 ) were set either to be constants or as a function of R rs (λ) (those parameters were determined based on a sensitivity analysis conducted in this study; see Sect.3.2 for details; (3) b bp (555) was calculated using b bp (443) and η; (4) a NAP (443) was calculated using the empirical relationship between a NAP (443) and b bp (555) obtained for Arctic waters (i.e. a NAP (443) = b bp (555)/0.2393;Matsuoka et al., 2007); and (5) a CDOM (443) was obtained by subtracting a NAP (443) from a CDM (443).Note that the correction factor for the a CDM (443) estimate of 0.754188 applied by Maritorena et al. (2002) was not applied Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Werdell, P. J. and Bailey, S. W.: An improved in-situ bio-optical data set for ocean color algorithm development and satellite data product validation, Remote Sens. Environ., 98, 122-140, 2005.Wirth, M. and Renger, W.: Evidence of large scale ozone depletion within the arctic polar vortex 94/95 based on airborne LIDAR measurements, Geophys.Res.Lett., 23, 813-816, 1996Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 1 .
Fig. 1.Locations of sampling stations for the SBI spr (red open circles), SBI sum (blue diamonds), MR aut (green crosses), and MALINA (red closed circles) in the Western Arctic Ocean.

Fig. 5 .Fig. 7 .
Fig. 5. (a) A composite image of CDOM absorption, a sat CDOM (443) for August 2009 using three weeks of MODIS R rs (λ) data, with sea ice concentration shown in gray.(b) Histogram of a sat CDOM (443) values.

Table 2 .
Comparison of a CDOM (443) estimates with measured a CDOM (443) for either oceanic or coastal water samples.