{'fraction_non_alphanumeric': 0.03399114244597896, 'fraction_numerical': 0.031023750292834752, 'mean_word_length': 4.654576420866713, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A global end-member approach to derive aCDOM(440) from near-surface optical measurements

This study establishes an optical inversion scheme for deriving the absorption coefficient of colored (or chromophoric, depending on the literature) dissolved organic material (CDOM) at the 440 nm wavelength, which can be applied to global water masses with near-equal efficacy. The approach uses a ratio of diffuse attenuation coefficient spectral end-members, i.e., a shortand long-wavelength pair. The global perspective is established by sampling “extremely” clear water plus a generalized extent in turbidity and optical properties that each span 3 decades of dynamic range. A unique data set was collected in oceanic, coastal, and inland waters (as shallow as 0.6 m) from the North Pacific Ocean, the Arctic Ocean, Hawaii, Japan, Puerto Rico, and the western coast of the United States. The data were partitioned using subjective categorizations to define a validation quality subset of conservative water masses (i.e., the inflow and outflow of properties constrain the range in the gradient of a constituent) plus 15 subcategories of more complex water masses that were not necessarily evolving conservatively. The dependence on optical complexity was confirmed with an objective methodology based on a cluster analysis technique. The latter defined five distinct classes with validation quality data present in all classes, but which also decreased in percent composition as a function of increasing class number and optical complexity. Four algorithms based on different validation quality end-members were validated with accuracies of 1.2 %–6.2 %, wherein pairs with the largest spectral span were most accurate. Although algorithm accuracy decreased with the inclusion of more subcategories containing nonconservative water masses, changes to the algorithm fit were small when a preponderance of subcategories were included. The high accuracy for all end-member algorithms was the result of data acquisition and data processing improvements, e.g., increased vertical sampling resolution to less than 1 mm (with pressure transducer precision of 0.03–0.08 mm) and a boundary constraint to mitigate wave-focusing effects, respectively. An independent evaluation with a historical database confirmed the consistency of the algorithmic approach and its application to quality assurance, e.g., to flag data outside expected ranges, identify suspect spectra, and objectively determine the in-water extrapolation interval by converging agreement for all applicable end-member algorithms. The legacy data exhibit degraded performance (as 44 % uncertainty) due to a lack of high-quality near-surface observations, especially for clear waters wherein wave-focusing effects are problematic. The novel optical approach allows the in situ estimation of an inwater constituent in keeping with the accuracy obtained in the laboratory.


Introduction
The colored (or chromophoric, depending on the literature) dissolved organic matter (CDOM) spectral absorption 1 coefficient, a CDOM (λ), where λ is wavelength, is widely used to investigate terrestrial and oceanic biogeochemical 2 processes, as summarized in the review by Nelson and Siegel (2013). The selection of a CDOM (440) as the parameter 3 of interest is a consequence of the relationships between CDOM and the solar illumination of aquatic ecosystems, 4 as follows: a) CDOM protects microorganisms from harmful ultraviolet (UV) radiation, albeit while reducing pho-5 tosynthetically available radiation (Nelson and Siegel 2013); b) CDOM affects the heat content of a water mass, 6 e.g., causing stratification for brown lakes (Houser 2006); c) CDOM supplies inorganic nutrients, i.e., ammonium, 7 (Bushaw et al. 1996) and can be a source of labile organic substances (Mopper et al. 1991) through photochemical 8 degradation and mineralization processes; and d) CDOM is a potentially useful parameter to distinguish and trace 9 water masses in the coastal zone and open ocean (Nelson et al. 2007 andTanaka et al. 2016). 10 This study evaluates whether a proposed algorithm (Hooker et al. 2013) for deriving a CDOM (440) from a ratio 11 of diffuse attenuation coefficient spectral end members, K d (λ 1 )/K d (λ 2 ), with the shortest wavelength denoted λ 1 12 and the longest denoted λ 2 , can be applied to global water masses with equal efficacy. Typically, a CDOM (440) is 13 determined in the laboratory using an optical instrument and a water sample obtained in situ. In this study, the 14 water sample is collected in temporal close proximity to in-water optical sampling used to derive K d (λ) from vertical 15 profiles starting sufficiently close to the water surface to accurately derive the widest range of K d wavelengths. 16 The Hooker et al. (2013) a CDOM (440) algorithm is based on a straightforward principal, as follows: if a water 17 mass is studied optically in a homogeneous near-surface interval of the water column, optical data products can be 18 derived for all wavelengths and the most sensitive parts are the spectral end members. The end members exhibit 19 the greatest range in values as a function of the absorption and scattering processes responsible for the attenuation 20 of light and can be inverted to derive typical constituents as a function of changes in attenuation properties. 21 Although an ability to derive an in-water constituent from optical measurements provides a follow-on opportunity 22 for airborne or satellite synoptic surveys, this is not a principal objective of this study. The reason for de-emphasizing 23 remote sensing is the remote sensing instruments typically available do not provide the spectral range used herein, 24 although a legacy pair of wavelengths are considered below. In addition, the principal parameter used here is 25 not measured directly by a remote sensor. Although K d (λ) can be derived from remote sensing data for part of 26 the needed spectrum (Cao et al. 2014), the inversion is incomplete and the introduced inaccuracies compromise 27 a principal objective, which is to determine a CDOM (λ) in situ with an accuracy commensurate with laboratory 28 analyses. Finally, despite the success with high-spatial-resolution remote sensing platforms for studying coastal and 1 inland waters (Palmer et al. 2015 andMouw et al. 2015), this study includes water bodies too small to be studied by 2 such platforms, and the reduced spectral resolution or range, coupled with the methods for characterizing a CDOM (λ), 3 have proved inadequate even within relatively large, lacustrine water bodies (Kutser et al. 2015). 4 The global perspective refers to a generalized concept of sampling a multitude of geographical areas and watersheds 5 wherein three broad categories are sampled: open ocean, coastal zone (e.g., shelf waters, bays, estuaries, lagoons, 6 etc.), and inland water bodies (e.g., rivers, lakes, reservoirs, wetlands and marshes, etc.). The near-surface viewpoint 7 is not driven exclusively by the desire to produce data products at all wavelengths. The other reasons for sampling 8 and deriving data products close to the water surface are as follows: a) establish a technique that can ultimately 9 support remote sensing objectives as the technologies advance, wherein the spaceborne and airborne approaches 10 obtain data products directly from the sea surface signal; b) use the same protocols for sampling and deriving data 11 products for all water masses, so the widest dynamic range in water properties can be considered (the shallowest 12 water depth sampled was 0.6 m); and c) improve the use of the global solar irradiance observations (obtained with a 13 separate solar reference) in setting a constraint for the fitting procedures used to derive the in-water data products 14 (Antoine et al. 2013), the effectiveness of which is related to how close the extrapolation interval is to the surface. 15 Using a homogeneous near-surface interval to derive all data products ensures the spectral interrelationships 16 coincide with the same water used to determine the in-water constituents by laboratory analysis. The perspectives of 17 natural changes and typical properties are also important, because some water bodies are not automatically assumed 18 to have typical water properties. For example, endorheic lakes are enclosed, so ground seepage and evaporation are 19 the principal outflow mechanisms with evaporation continuously concentrating constituents (Yapiyev et al. 2017). 20 Over time, a narrowly defined ecosystem evolves to withstand the increasingly extreme conditions, and in some cases, 21 higher-order life ceases to exist, e.g., the Dead Sea (Bodaker et al. 2010). 22 Endorheic lakes are an end point in the expression of water masses, because the range in the temporal gradient 23 of a constituent, e.g., salt, is somewhat unbounded and the water body does not evolve conservatively due to the 24 significant outflow versus inflow imbalance. For purposes exclusive to this discussion, a conservative water body is 25 defined to have an inflow and outflow of properties that constrain the range in the gradient of a constituent. This 26 natural range is usually established by seasonal factors, although unnatural or atypical stressors can add optical 27 complexities, which may or may not be seasonal. Examples of the latter are anthropogenic sources (e.g., pollution or 28 agricultural water diversion) and severe weather (e.g., typhoon-induced bottom resuspension in coastal ecosystems). 29 Consequently, other water bodies subjected to an unexpected stressor that allows an unbounded gradient in 1 a constituent, e.g, long-term drought, are anticipated to not evolve conservatively and the constituents might be 2 expressed as extreme values as a function of time. Once the stressor is no longer applied, the water mass evolves 3 semi-conservatively, wherein the atypical properties are diluted or removed, and at some point in time the water 4 body reverts to a conservative evolution, i.e., the gradient in the constituent is within an expected or natural range. 5 A global perspective is constructed with overlapping ranges in the natural gradient of the constituent and the 6 optical inversion parameters within conservative water masses (Lee and Hu 2006). If the assembled dynamic range 7 extends across a sufficiently dense sampling of clear to turbid water masses, an explicit sampling of every possible 8 global water mass is not deemed necessary. The turbid-water endpoint of the dynamic range is somewhat undefined, 9 because of present limitations in obtaining in-water optical measurements in extremely shallow or turbid waters (a 10 case in point, White Lake, is presented below), but the clear-water endpoint is defined by the pure-water limit, i.e., 11 no constituents (Morel 1974). Consequently, the dynamic range herein can only be extended in one direction and any 12 turbid additions involve necessarily small volumes, so the global perspective is at most only marginally incomplete. 13 Based on the degree of complexity for a water mass not evolving conservatively, it is anticipated that such water 14 masses are unsuitable for validating a global algorithm established to invert the optical properties of conservative 15 waters. Whether a validation site is inappropriate is a function of the severity of the stressor creating the non-16 conservative evolution. For example, a short-term water diversion from a lake is expected to create a short-term 17 complexity, whereas a long-term drought likely creates a time series of increasingly extreme values (Vazquez et al. 2011 18 andGuarch-Ribot andButturini 2016). Consequently, the sampling objective used here was to obtain measurements 19 in conservative water masses plus water bodies subjected to one or more stressors. To assess performance, both 20 subjective and objective classifications of water mass complexity are included in the algorithm evaluation process. 21

Methods
The Hooker et al. (2013) study excluded lacustrine water bodies; the largest inland water masses were tidal 22 estuaries. Consequently, the new validation data set includes a large variety of lakes and reservoirs, wherein some 23 were selected precisely because compliance with the original (Hooker et al. 2013) algorithm was not anticipated. 24 These nonconservative water bodies provide an important test of the algorithmic approach, because if they do not 25 appear as outliers with respect to the original algorithm, the principles behind the algorithm are challenged.

26
To improve the quality of optical measurements obtained in near-surface waters, which is essential for studying 27 shallow ecosystems, methodological advancements were included for this study. Consequently, the methods described 28 herein are distinguished with respect to the original Hooker et al. (2013) research as follows: a) a significantly enlarged 1 study area with new water body types (e.g., lakes and reservoirs, more numerous rivers, the marginal ice zone, etc.); 2 and b) the use of more advanced optical technologies to improve sampling efficiency and data quality. 3

Optical Instrumentation
The optical instrument suite deployed for this study is a handheld, free-falling Compact-Optical Profiling System 4 (C -OPS), as first described by Morrow et al. (2010), that measures the downward irradiance and upwelled radiance, 5 E d (z, λ) and L u (z, λ), respectively, where z is depth. An above-water reference, sited to avoid shadows and reflections, 6 simultaneously measures the global solar irradiance, E d (0 + , λ), where 0 + indicates above the water surface. This 7 configuration was deployed by Hooker et al. (2013), except the study documented herein used advanced radiometers 8 with three gain stages rather than two. The majority of profiles were obtained with the Compact-Propulsion Option 9 for Profiling Systems (C -PrOPS) plus a conductivity sensor for improved water mass characterization. The former 10 uses two small digital thrusters to maneuver the backplane (Hooker et al. 2018a) beyond the influence of platform 11 perturbations, which does not remove the self-shading effect, thus, necessitating an L u (λ) correction (Gordon and 12 Ding 1992).  provides the negative consequences of common positioning alternatives. 13 A transparent drawing of the next-generation C -OPS with C -PrOPS is presented in Fig. 1. When the weak thrust 14 holding the (slightly) negatively buoyant profiler at the surface is removed, one or more bladders in the hydrobaric 15 buoyancy chamber slowly compress and increase the near-surface loitering of the profiler, which results in a vertical 16 sampling resolution (VSR) to within 1 cm or less. The VSR is defined as the vertical extent of the extrapolation 17 interval used to derive the data products, e.g., K d (λ), divided by the number of retained data points in the interval, 18 wherein retention requires planar orientation to within 5 • of vertical. For coastal and inland waters, the average 19 VSR was 6.0 mm, but for very shallow or turbid waters, the average VSR was 0.9 mm. In comparison, the Hooker et 20 al. (2013)  open-ocean profiles were in a more turbulent wave field, so the profiler was ballasted to sink faster and descend 22 deeper. The open ocean deep mixed layers mean a slightly coarser vertical resolution is not a limitation. 23 The two digital thrusters have the same cant angle with respect to the vertical, which directs the weak turbulence 24 from the thrusters downward and below the irradiance instrument, thereby ensuring both light apertures are observing 25 undisturbed water; the opposite occurs if thrust is reversed. To steer the backplane like a remotely operated 26 vehicle, differential thrust is applied to the two thrusters (Hooker et al. 2018a) and allows for real-time positioning 27 adjustments, which is a significant advantage in shallow waters, e.g., away from a shoreline or within a wetland. 28 Once the profiler reaches the desired position for obtaining a vertical profile of measurements, or cast, it is kept 1 in position by maintaining weak forward thrust. While at the surface, the pressure transducer measures atmospheric 2 pressure right before a profile commences, which allows a pressure tare for every cast and improves the accuracy 3 of depth measurements . The pressure transducers used herein have a depth resolution, in terms of 4 precision, of 0.03-0.08 mm in all water masses. When thrust is removed, the thruster-induced bias in the roll axis 5 relaxes (the pitch angle is negligible, because prior thrust aligns the backplane with almost no pitch angle), and 6 the profiler descends with stable tilts (Hooker et al. 2018a). Unlike rocket-shaped profilers (Hooker et al. 2001), the 7 profiler has no significant righting moment and the planar orientation of the radiometers is maintained from the start 8 of data acquisition, which significantly improves the VSR. In deep waters, the thrusters are used at the bottom of 9 the cast to reduce the time between casts; in very shallow waters, the profiler was hauled closer to the surface before 10 the thrusters were used to prevent resuspension of bottom material. All profiles were obtained in waters wherein the 11 10% light level was above the bottom depth to ensure all data products were uncontaminated by bottom reflections. 12 Seven different instrument suites including a next-generation hybridspectral profiler with fixed-wavelength and 13 hyperspectral detector components plus a C -PrOPS (Hooker et al. 2018b) were used for this study. The fixed-14 wavelengths of the radiometers had similar configurations such that all measured the same nine spectral end members 15 from 320-412 nm and 670-780 nm, plus six common wavelengths (Table 1). All optical instruments were calibrated 16 at the same manufacturer facility with traceability to the National Institute of Standards and Technology (NIST) 17 as described by Hooker et al. (2018b). NIST traceability is a requirement of the NASA Ocean Optics Protocols 18 (hereafter, the Protocols). The Protocols set the standards for calibration and validation activities (Mueller and 19 Austin 1992), which were revised (Mueller and Austin 1995) and updated over time (Mueller 2000(Mueller , 2002(Mueller , and 2003. 20 A comparison of C -OPS acquisition with and without thrusters (Hooker et al. 2018a) verified the former improved 21 efficiency in all waters by a factor of two or more, plus either no or minor adjustments to the extrapolation interval used 22 to derive data products for replicate casts were needed (thrusters minimize the negative influences of heterogeneity 23 across all wavelengths). The improved efficiency yields a closer temporal matchup between the collection of optical 24 profiles and the water sample. Approximately 98% of water samples were obtained from the surface using a bucket. 25 For some inland waters, the profiler was launched from a shoreline or dock e.g., when the boat ramp was out of service 26 (due to drought, flooding, invasive species regulations, etc.). If a water sample could not be otherwise retrieved from 27 the profiling location, the Compact-Profiler Underway Measurement Pumping System (C -PUMPS) was used (Hooker 28 et al. 2018a). C -PUMPS provides a 20 ml s −1 flow rate from the profiler and fills a 1 l container in less than 1 min.

Field Sampling
The Hooker et al. (2013) study area was the Beaufort Sea in proximity to the Mackenzie River outflow, the 1 Gulf of Maine and vicinity, including major portions of its inland watershed plus minor watershed drainage from 2 smaller rivers and a saltwater marsh. A validation data set from observations made in U.S. coastal waters within the 3 southern Mid-Atlantic Bight were also used. Neither of the two data sets included typical lacustrine water masses. 4 The new sampling area for the study herein included the western U.S. (i.e., California,Oregon,Washington,Nevada,5 Utah, and Idaho), Hawaii, Puerto Rico, Japan, the western North Pacific Ocean (e.g., the Kuroshio and Oyashio 6 Currents), the central North Pacific Ocean, the Bering Sea, the Chukchi Sea, and the Beaufort Sea (Fig. 2). The 7 latter is the only region that slightly overlaps the Hooker et al. (2013) study. The new data set includes sampling in 8 a wide diversity of inland rivers, lakes, and reservoirs, including hypersaline and alkaline lakes. 9 The new field data are divided into the aforementioned three primary categories according to whether or not the  waters, a single seawater sample was usually collected. A selected volume of each water sample was filtered through 22 a 0.22 µm filter under a gentle vacuum and collected in appropriate (e.g., pre-combusted) clean glass vials or bottles. 23 Typically, CDOM absorption was measured within a few hours after sampling. In some cases when this was not 24 possible, samples were stored at −20 • C or less until subsequent laboratory analysis at a shore facility (Sect. 2.4). 25 While Fellman et al. (2008) reported that freezing changed the chemical composition of DOC for freshwater (stream) 26 samples, Hancke et al. (2014) found no such effect for Arctic marine samples, nor was any systematic bias observed 27 between the subset of frozen samples compared to the full suite of samples used for this analysis. 28 The surface water sample was obtained as quickly as possible after three optical casts were performed. In 1 some cases when the heterogeneity or turbulence of the water mass was considered to be excessive, an additional 2 three optical casts were executed immediately after the water sample was collected. The determination of excessive 3 conditions was based on the stability achieved in optical variables (e.g., the average vertical tilt in the upper 2 m of 4 the water column, changes in the 10% light level depth, etc.) during data acquisition for the first three casts. 5

Optical Data Processing
All optical data products discussed herein, e.g., K d (λ), were estimated in a near-surface interval of the water 6 column with homogeneous properties as confirmed with physical data and analysis of the linearity of extrapolations 7 provided by the processing scheme. The processor used here is based on a well-established methodology (Smith and 8 Baker 1984) that Hooker et al. (2001) showed is capable of agreement at the 1% level within an international round 9 robin, when the processing options are as similar as possible and both data acquisition and processing strictly adhere 10 to the Protocols. Summary details of the data acquisition and processing capabilities are provided in Antoine et al. 11 (2013), , and Hooker et al. (2018aand Hooker et al. ( , 2018band Hooker et al. ( , and 2018c, so only brief overviews are presented herein. In-water radiometric parameters in physical units are normalized with respect to separate, but simultaneous, 13 E d (0 + , λ, t) measurements, with t expressing time dependence. After solar normalization and ±5 • tilt filtering, a 14 near-surface homogeneous portion of E d (z, λ) centered at depth z 0 and extending from z 1 = z 0 +∆z and z 2 = z 0 −∆z 15 is established separately for the blue-green and red wavelengths; the UV and near infrared (NIR) wavelengths are 16 included in the blue-green and red intervals, respectively. Both intervals begin at the same shallowest depth, but the 17 blue-green interval is allowed to extend deeper if the extrapolation linearity, as determined statistically, is thereby 18 improved (this only occurs in oligotrophic, optically simple, waters with deep mixed layers). The negative value of 19 the regression slope yields K d (λ), which is used to extrapolate the fitted portion of the E d profile to z = 0 -.

20
A principal benefit of profile data with a high VSR is that the aliasing caused by wave-focusing effects (Zaneveld 21 et al. 2001) can be significantly reduced during data processing. The separately obtained above-and in-water E d 22 values at z = 0 + and z = 0 -, respectively, can be compared using The appropriateness of the E d extrapolation interval, initially established by z 1 and z 2 , is evaluated by determining 2 if (1) is satisfied to within the 2.3-2.7% uncertainty (k = 2 coverage factor) of the optical calibrations (Hooker et al. 3 2018b); if not, z 1 and z 2 are redetermined-while keeping the selected depths within the shallowest homogeneous layer 4 possible-until the disagreement is minimized (usually to within 5% to include some inevitable variance from natural 5 processes to the calibration uncertainty). In this procedure, selection of the near-surface extrapolation interval uses 6 a boundary condition or constraint (Antoine et al. 2013), wherein the central tendency of the distribution of data 7 within the extrapolation interval, which are typically subjected to wave-focusing effects, satisfies (1). 8 The linear decay of all light parameters in the selected near-surface layer are then evaluated, and if linearity is 9 acceptable, the entire process is repeated on a cast-by-cast basis. Subsurface quantities at null depth are obtained 10 from the slope and intercept given by the least-squares linear regression versus z within the extrapolation interval 11 specified by z 1 and z 2 . A secondary benefit of profile data with a high VSR is that the extrapolation interval can have 12 a restricted vertical extent, but still have sufficient data to satisfy (1) and produce data products at all wavelengths. 13 This is an important advantage in optically complex water masses, which are usually turbid and shallow. 14

Water Sample Processing
For the Pacific Ocean samples and approximately half of the Arctic samples, the absorption spectrum of CDOM 15 was determined using a spectrophotometer (Shimadzu UV-1800) according to Yamashita et al. (2013). Briefly, after 16 the water sample was thawed and reached room temperature, the spectral absorbance was measured from 200-800 nm 17 at 0.5 nm intervals with a 10 cm quartz-windowed cell. Absorbance spectra of a blank (Milli-Q water) and samples 18 were obtained against air, and a blank spectrum was subtracted from each sample spectrum. The blank-corrected 19 absorbance spectrum was baseline-corrected by subtracting average values ranging from 590-600 nm (Yamashita 20 and Tanoue 2009), and then converted to the absorption coefficient (Green and Blough 1994). A single absorbance 21 analysis was generally carried out for the open ocean samples, with an average accuracy from replicates of 2.5%.

22
For the other Arctic samples, CDOM absorption coefficients were determined using an UltraPath liquid waveguide 23 (Matsuoka et al. 2012 and2017). Briefly, CDOM absorbance for a filtrate (less than 0.2 µm) was measured relative to a 24 reference water within a few hours after sampling. The reference water was prepared in advance using pre-combusted 25 pure salt (450 • C for 4 h) with Milli-Q water to adjust salinity within ±2 between a sample and a reference to correct 26 for the refractive index effect. A 2 m optical path was used for all waters except some coastal sites wherein a 0.1 m 27 path was used (Matsuoka et al. 2012 and. After baseline correction, absorbance was converted into absorption 28 coefficients by including the optical path length. The detection limit was within 0.001 m −1 (Matsuoka et al. 2017).

1
For the western U.S. coastal and inland waters, water samples were passed through a 0.2 µm syringe filter 2 (Whatman GD/X) and absorbance of CDOM was measured on either a Cary Varian 50 spectrophotometer using a 3 10 cm quartz cell or an UltraPath liquid waveguide spectrometer with 2 m path length. The syringe filter was rinsed 4 with sample prior to collection, with the sample stored in an amber, acid-washed and combusted (450 • C for 4 h) glass 5 vial with Teflon septa, and kept in the dark at 4 • C until analysis. Absorption spectra of the filtered samples were 6 measured using ultrapure water from a Millipore Milli-Q A10 pure water system with UV to reduce total organic 7 carbon to less than 10 ppb. The absorption coefficient (calculated as absorbance divided by path length, multiplied 8 by 2.303 to convert to natural log units) at 440 nm to represent CDOM abundance was estimated using the Single 9 Exponential Model (SEM) for absorption from 300-700 nm as described by Twardowski et al. (2004). The three 10 different methods used in this study for determining CDOM absorption do not influence the results (see a sensitivity 11 analysis in Sect. 4). For all the CDOM data-regardless of the collection and processing method-the spectral slope 12 for the wavelength range 350-500 nm was calculated by fitting a nonlinear least-squares model to a CDOM (λ). 13

Data Subcategories
Algorithm validation requires an assessment prior to data collection whether or not the sampling violates the 14 assumptions made to create the algorithm. For this study, water masses wherein the evolution of a CDOM (440) or 15 K d (λ) was considered conservative, i.e., within the likely range in the gradient of such properties because no stressors 16 to challenge that perspective were evident, are considered validation quality with primary categories of open ocean, 17 coastal zone, or inland waters. Those water masses thought to contradict the algorithmic approach, because one or 18 more stressors challenging the conservative evolution perspective were evident prior to sampling, are subcategorized 19 to exclusively assess algorithm performance as more complex water masses are included in validation, as follows: 20 1. Waters closer to an ice field contain meltwater properties, which freshen the neighboring water body and can 21 result in additional particles or compounds not usually found in the parent water mass. 22 2. Waters farther from an ice field, but within proximity, containing lesser amounts of meltwater properties. 23 3. Resuspension occurs naturally when a sufficient flow (e.g., an ebb or flood tide) or turbulent wave field 24 (e.g., created by sufficiently strong winds) interacts with shallow bottom sediment to create concentrations 25 of constituents that would otherwise not be present; it occurs unnaturally when a boat propeller (or other 26 mechanical device) churns up shallow bottom sediment (e.g., in a harbor, marina, or navigation channel).

4.
A refilled lake experiences a rapid inflow of alluvium (e.g., gravel, sand, silt, and clay) from riverbeds and 1 eroding banks, plus floating and partially submerged debris, that can also resuspend bottom sediment. If the 2 refilled lake is a controlled reservoir and exceeds the normally maintained fill level, new lake bottom is added, 3 which can be a source for additional, perhaps atypical, water constituents in terms of type or concentration. 4

5.
A drought-stricken lake has a longer residence time (the amount of time for the time-elapsed outflow to equal 5 the lake volume) than normal, because once the water level remains below the overflow elevation, evaporation 6 and ground seepage are the primary outflows. Increased residence time can concentrate constituents, plus 7 dried and exposed bottom material can be resuspended into the shrinking lake volume due to wind or rain. can be a source of pollutants and bottom resuspension, and typically include structures (breakwaters, jetties, 10 piers, etc.) for shelter from severe weather, which can alter residence times by restricting water exchange. 11

A harmful algal bloom (HAB) is a toxic or noxious algae in a concentration producing a deleterious effect on 12
humans or the environment. HABs are usually influenced by chemical, physical, and biological factors. 13 8. A wetland (plus marsh or mangrove) filters dissolved and suspended water constituents (e.g., from tidal cycles, 14 weather events, etc.) through settling and plant consumption, but might not completely remove them. 15 9. A polluted water mass is contaminated from an anthropogenic source that alters the natural water properties. 16 10. An alkaline (or soda) lake has limited biodiversity due to an elevated pH of 9-12 with high carbonate and 17 complex salt concentrations affecting the solubility and toxicity of chemicals and heavy metals (Grant 2006). 18 11. A hypersaline lake contains high concentrations of sodium chloride (or other salts) surpassing seawater, which 19 limits biodiversity to organisms tolerating high saline levels (e.g., Mono Lake had a salinity of about 50).

12.
A river mouth is where significant amounts of alluvium are deposited into a larger water body (e.g., a delta). 21 13. An atypical bloom is a generic case of high biomass based on local reports evaluated with respect to typical 22 temporal and spatial conditions, which may involve weather effects (e.g., wind) concentrating algae advectively. 23 14. An invasive species is an introduced plant, fungus, or animal that is not native to a water body and is 24 anticipated to alter the heretofore established properties and perhaps with a significantly negative outcome 25 (e.g., damage to the environment, economy, or health of organisms, including humans).

26
15. A parent water mass modifier is a localized alteration of water properties, e.g., a creek inflow into a lake, and 27 demonstrates the sensitivity of the methods used herein to distinguish small changes. 28 The above 15 subcategories plus the original validation quality category results in 16 categorizations, and as the 1 former more complex waters are incrementally added to the latter, quantitative cause-and-effect validation scenarios 2 are created (Sect. 4). If a sample was applicable to more than one subcategory, e.g., a wetland can experience 3 resuspension from tidal currents, a dominant subcategory was selected based on observations prior to sampling. 4 Categorization ambiguities are not worrisome, because complex water masses are a small fraction of global ecosystems. 5 The proximity to ice (farther and closer) subcategories are based on the relative position of safely operating a 6 small vessel in and around an ice field. Sampling was usually as close to the ice as possible, and then as far from the 7 ice within line of sight of the larger ship the small boat was launched from. In comparison, categorizing a refilled lake 8 is a straightforward comparison of the water level datum available from local authorities with respect to the outflow 9 elevation and historical norms. All refilled lakes were at 100% capacity or more, e.g., Washoe Lake was overfilled.

10
The categorization of bottom resuspension is primarily based on visual evidence, wherein resuspended particles 11 are visible and produce a significant change in water color (e.g., Akkeshi Bay after the passage of typhoon Vongfong). 12 Although a subset of sampling obtained in harbors could be classified as resuspension stations, a harbor (or marina) 13 is identified based on local identification of such facilities. Similarly, wetlands (plus marshes and mangroves) are 14 identified based on navigation charts, maps, and local descriptions. Hypersaline (endorheic) lakes are similarly 15 categorized by state and local authorities (e.g., Mono Lake, Great Salt Lake, and Salton Sea), as are alkaline lakes 16 (e.g., Mono Lake, with dual classification, plus Borax Lake and Soda Lake).

17
Categorizing drought-stricken lakes relies on local authorities reporting lake elevations and inflow water volumes 18 with respect to historical norms. Examples from 2015 are as follows: a) Shasta Lake water storage was 56% below 19 normal storage; b) Lake Almanor was 118.5 ft (36.1 m) below normal elevation; c) the Truckee River flow into Pyramid 20 Lake (Nevada) was near historical lows and was dry for three days prior to sampling; and d) Eagle Lake had a water 21 level of 5,091.5 ft (1,551.9 m), which was within 0.5 ft (0.2 m) of the lowest level recorded in 1935.

22
The categorization of a river mouth is a combination of the geographical location (e.g., the Columbia River) and 23 evidence of the presence of the water mass the river flows into (e.g., salt water intrusion from the near-shore bay). 24 The inflow of smaller rivers, streams, and creeks into a larger water body (e.g., Ward Creek flowing into Lake Tahoe) 25 are not classified as river mouths, but rather as parent water mass modification from a creek inflow. The creek 26 designation is to ensure the understanding that the inflow volume is small, but the expectation is changes in water 27 properties are nonetheless discernible, because of the enhanced sensitivity (e.g., VSR) of the methods used. 28 The categorization of a water mass subjected to pollution, an atypical bloom, HAB, or an invasive species relied 29 principally on eutrophic chlorophyll concentrations plus historical reporting or on-site local representatives. The 1 latter are frequently present at boat ramps to oversee measures to mitigate health concerns or prevent the spread of 2 invasive species, which is presently a significant and escalating problem throughout the western U.S. 3

NOMAD Archive
The NASA bio-Optical Marine Algorithm Dataset (NOMAD) v2.a (Werdell and Bailey 2005) is a small, qual-4 ity controlled subset of a larger data repository established early in the Sea-viewing Wide Field-of-view Sensor 5 (SeaWiFS) satellite mission (Hooker and Esaias 1993) called the SeaWiFS Bio-optical Archive and Storage System 6 (SeaBASS) and is described by Hooker et al. (1994). The NOMAD database does not include applicable a CDOM (440) 7 measurements with contemporaneous UV and NIR spectral end members as used in this study. Consequently, the 8 Hooker et al. (2013) algorithmic approach, which is based on UV and NIR end members, cannot be evaluated. 9 The NOMAD database, however, does include K d (λ) with legacy visible (VIS) wavelengths plus matching dis-10 solved (Gelbstoff) spectral absorption coefficient at 443 nm, a g (443), which is functionally equivalent to a CDOM (443) 11 following Röttgers and Doerffer (2007). The consequences of the 3 nm shift in a g (443) with respect to a CDOM (440) 12 are considered negligible for a generalized inquiry involving legacy optical data, because the fixed wavelengths in-13 volved have 10 nm bandwidths and there are multiple sources of uncertainties in the derived data products of equal 14 or greater importance (Hooker et al. 2013), e.g., pressure tares, aperture offsets, dark currents, wave focusing, etc.

Alternative Classification Scheme
The subcategory scheme (Sect. 2.5) included factual knowledge combined with careful inspection regarding one or 16 more significant constituents or stressors influencing water mass complexity. Previous studies applied fuzzy c-means 17 (FCM) classification to ocean color algorithm development (Moore et al. 2001) or uncertainty estimation (Moore 18 et al. 2015), thereby demonstrating the usefulness of the FCM over a crisp or hard (e.g., k-means) classification. 19 Application of the FCM classification is to log-transformed K d (λ) data based on the Calinski and Harabasz (1974) 20 index (Matsuoka et al. 2013

Results
All data not categorized as one of the 15 subcategories prior to data collection (Sect. distinctions: a) the addition of lacustrine water bodies (including Lake Shikotsu in Hokkaido, Japan), which almost 5 span the entire three decades of dynamic range (e.g., Crater Lake to Pinto Lake); and b) the expansion of the 6 dynamic range to include clearer and more turbid water masses, which also means deeper and shallower waters. 7 Given the diversity of sampling in this study, it is unlikely that the more than three decades in a CDOM (440)  S presented in this study is comparable to similar recent global analyses spanning the majority of marine waters. 16 The new validation quality data significantly adhere to the original algorithm, as evidenced by how the red, 17 green, and blue circles in Fig. 3 are well contained within the approximately ±15% gray boundaries that denote the 18 dispersion in the original algorithm. The category that spans the largest percentage of the dynamic range for both 19 axes is the inland waters data, although the coastal zone is somewhat similar because of the clear Hawaiian coastal 20 waters that were sampled. The open ocean category has the smallest dynamic range, but this does not diminish the 21 importance of this category because it represents the greatest surface area and volume of water on the planet.

22
The new data values, V n , are compared to the corresponding algorithm value, with the latter being the reference 23 values, V r , in the comparison calculation. The relative percent difference (RPD) between the new data and the 24 algorithm is computed as RPD = 100 (V n − V r )/V r, and is expressed as a percent. The average RPD for all the 25 new data is 0.02%, i.e., the new data show a negligible bias with respect to the original algorithm. The absolute 26 percent difference (APD), which provides an estimate of the dispersion of differences between the new data and the 27 algorithm, is the absolute value of the RPD. The average APD value for all the new data is 3.9%, i.e., the new 28 validation quality data are usually to within 5% of the original algorithm (as visually confirmed by Fig. 3).

Drought-Stricken, Alkaline, Hypersaline, and Refilled Lakes
The new lacustrine data are presented in Fig. 4. Data from the hypersaline and alkaline (endorheic) lakes do not 1 conform with the algorithm. Drought-stricken lakes exhibit a wider range of departure, with the most significant 2 occurring for the most depleted water bodies, e.g., Lake Almanor and Shasta Lake. Endorheic drought-stricken lakes, 3 e.g., Eagle Lake and Pyramid Lake, are the most extreme. Refilled lakes also do not conform with the algorithm, 4 and refilled drought-stricken lakes exhibit an increase in CDOM and turbidity, e.g., Shasta Lake and Pyramid Lake. 5 The refilled lakes in Fig. 4 are frequently more different with respect to the algorithm than hypersaline or alkaline 6 lakes, especially in terms of turbidity as determined by the K d ratio. This is because some of the refilled lakes are 7 overfilled, wherein the shore of the lake extends beyond the normal acreage of the lake (e.g., Washoe Lake and Little 8 Washoe Lake). In overfilled lakes, land that is not normally flooded is added as new lake bottom, and the new 9 acreage is a source of atypical constituents, either in composition or concentration (e.g., atypical constituents from 10 land-use activities can be added when a lake overfills, because these activities are not possible in the water mass). 11 The refilling of a normally dry endorheic basin, e.g., White Lake, wherein the flood waters and the reclaimed lake 12 bottom provide the maximum areal and volumetric source of dissolved and suspended constituents results in some 13 of the most extreme results, both in terms of turbidity and with respect to the algorithm.

14
The discharge from overfilled reservoirs also has significant deviations with respect to the algorithm. Thermalito 15 Afterbay receives the discharge from Lake Oroville and was sampled after the overfilling of the parent water mass 16 during the drought-breaking California 2016-2017 winter. Data were obtained in two locations, with higher CDOM 17 data obtained in a shallow marsh. The ability to distinguish small localized differences establishes the sensitivity of 18 the methods used herein. For example, the three refilled Shasta Lake samplings in Fig. 4 were conducted in different 19 locations subjected to the inflow of a creek, as well as a large floating and partially submerged debris field.

River Mouth, Resuspension, and Ice Edge Proximity
The inflow of dissolved and suspended constituents to a parent water mass is explored further by considering a 21 variety of sources that can add to water mass complexity. The new data are shown in Fig. 5 and were obtained 22 in river mouths, water bodies with known suspension or visible resuspension, plus samples obtained closest to or 23 farthest from the ice edge within an oceanic ice field. Water bodies with known suspension or visible resuspension 24 are primarily from tidal and riverine flows, which are shown in Fig. 5 as triangles. Almost all of the resuspension 25 data were obtained at peak tidal flow to ensure safe navigation in the necessarily shallow waters. The Akkeshi Bay 26 data were obtained the day after the passage of typhoon Vongfong, wherein the shallow bay waters were a distinctly 1 different color than normal. The Sacramento River data were obtained after heavy rains, wherein the boat ramp to 2 be used was closed due to flood waters. The difference between the flooded Sacramento River with respect to the 3 inland riverine data in Fig. 3 not in flood conditions (i.e., as conservative water masses), shows the classification of 4 the Sacramento River is appropriate and the subjective classification approach has merit.

5
The resuspension data in Fig. 5 also include Bear Lake (which straddles the Utah-Idaho border) plus the effects 6 of a large ship docking in the shallow RWC Channel with the aid of a tugboat. The latter involved the churning up 7 of bottom material that significantly changed the color of the water. The resuspension sampling occurred shortly 8 after the ship was docked. The Bear Lake scattering anomaly is created primarily through ground water seepage, 9 which is rich in calcium carbonate particles (Davis and Milligan 2011). The ground water is nutrient poor and the 10 small amount of riverine input to the lake is through a swamp and wetlands, wherein plants consume the nutrients 11 and sediments settle out. Consequently, the Bear Lake data represent a significant clear-water scattering anomaly.

12
The other clear-water data in Fig. 5 were obtained principally in proximity to Arctic oceanic ice fields, and are 13 distinguished as being closer to, or farther from, the ice. These data are displaced above or below the algorithm, 14 respectively, even in the turbid waters of Kotzebue Sound. The majority of these data were obtained using a small 15 boat launched from a larger ice breaker, so the data obtained closer to the ice are as safely close to the ice as possible 16 while beyond the shading of the water mass by the ice field. The classification of closer to and farther from is 17 qualitative and in complicated ice fields misclassification is possible. The Fig. 5 data show only two farther from 18 points that are likely not classified correctly, and this category is the most vulnerable to a qualitative error. 19 In regards to the resuspension data, which all cluster below the algorithm in Fig. 5, river mouth data are the 20 opposite-the data cluster above, but the number of such observations is much smaller. The reduced number is 21 due to the difficulty of operating a trailered small boat in a shallow river, and then safely navigating the vessel out 22 into the river mouth through a frequently narrow channel, wherein the higher sea state of the coastal ocean can 23 be significantly amplified and boat traffic can make station work hazardous. Within the plume of a river mouth, 24 two usually rather different water masses meet and mix over short time scales. Under those conditions, short-term 25 deviations, with respect to the algorithm, can emerge and that is what is shown in the Fig. 5 data.

Atypical Blooms, Invasive Species, and Harbors
The presence of an atypical bloom, particularly a HAB, was anticipated to create additional optical complexity, 27 because one or more significant stressors are frequently involved, e.g., an overabundance of nutrients, which can be 28 anthropogenic in origin (Heisler et al. 2008). In a generic context, an atypical bloom includes the concentration of 1 biomass to artificial levels (Kudela et al. 2015), perhaps due to local weather, (e.g., advective processes from winds 2 and waves). Invasive species and harbors were also expected to increase optical complexity. 3 The new data obtained in harbors and water bodies experiencing an invasive species, or atypical bloom, including 4 a HAB, are shown in Fig. 6. Some of these data could have had two classifications. For example, the Tahoe  few extreme values, which is expected because harbors exchange water with the parent water mass. San Leandro 12 and El Granada have the largest expression, but San Leandro is in a heavily urbanized area immediately south of 13 Oakland International Airport in the San Francisco Bay area, so significant anthropogenic sources are anticipated.
14 Like some coastal harbors, El Granada vessels are moored in an inner shallow harbor protected by an outer 15 deeper area, both with perimeter breakwaters and narrow channels. The two harbor areas cannot exchange water 16 completely, i.e., a portion of the water volume is trapped during each tidal cycle, and are more turbid than the parent 17 water mass, Half Moon Bay. The inner harbor is a likely and persistent anthropogenic source with a longer residence 18 time, so it is anticipated to have an a CDOM (440) value exceeding the neighboring bay. The increased residence time 19 and reduced exchange rates through the narrow channels are a possible mechanism to increase a CDOM (440). Other 20 harbors wherein a protected moorage has elevated a CDOM (440) include Las Vegas (Lake Mead) and Crescent City.

21
The HAB data in Fig. 6 were frequently obtained opportunistically and, thus, were not necessarily from the peak 22 of the phenomenon. Also, a bloom is heterogeneous and navigation within the bloom is mostly based on visual 23 observations, so the relationship with respect to the algorithm is not always extreme. The Monterey Bay HAB data 24 are the most extensive, because there was the opportunity for scheduling some of the data collection during a time 25 period when a HAB was likely to occur. In all cases, a HAB observation has a larger K d ratio than the algorithm 26 predicts, and this is principally caused by an increase in the K d (320) value, i.e., increased attenuation in the UV, 27 which might indicate the presence of mycosporine-like amino acids (Jessup et al. 2009 andKwon et al. 2018). 28 An atypical bloom is primarily a combination of local reporting, and a heterogeneous eutrophic water mass, i.e., 29 chlorophyll concentration exceeds 1 mg m −3 , with some water bodies having concentrations greater than 10 mg m −3 . 1 Consequently, the lack of sophistication and specificity related to explaining this subcategory does not exclude a 2 simpler explanation. For example, local wind conditions could elevate the values associated with a typical bloom 3 into atypical concentrations. This phenomenon was observed in more than one lake, e.g., Pyramid Lake and Upper 4 Klamath Lake. The majority of the atypical blooms are in rather close agreement with the algorithm. 5

Wetlands, Pollution, and Water Mass Modifiers
The new data obtained in wetlands or polluted waters are presented in Fig. 7. The former are almost all marsh 6 grass except two, which are labeled as to their types. The two unlabeled at the top of the plot are from Cutoff 7 Slough in California and are marsh grass. All wetlands exhibit the same relationship, that is, they are all displaced 8 above the algorithm, although four are in rather close agreement with the algorithm. The polluted water masses are 9 associated with agricultural (Upper Klamath Lake and Upper Elkhorn Slough) or mining (Clear Lake) runoff, with 10 the latter being the most severe. For both Upper Klamath Lake and Clear Lake, blue-green algae were plainly visible 11 with extreme maximum chlorophyll concentrations of 1.117 g m −3 and 1.420 g m −3 . The chlorophyll concentrations 12 in Upper Elkhorn Slough were less, but are still extreme with a maximum value exceeding 100 mg m −3 .
13 Figure 7 also includes examples of a small inflow from a creek or another source modifying the neighboring 14 parent water mass. These data provide a measure of the sensitivity of the data acquisition, processing, and analysis 15 techniques used herein. Although other sensitivity examples are documented above, e.g., the distinction between 16 sampling closer to, or farther from, the ice edge (Fig. 5), the Fig. 7 examples span diverse spatial scales, e.g., creek 17 inflows, a fish kill in the Salton Sea, and a large floating and partially submerged debris field in Shasta Lake. In all 18 cases, the anticipated algorithmic relationships appear different than the parent water mass. The water properties 19 of the creek inflow are not known, because access to the source from a small boat was problematic.

20
The generalized properties of the inflowing creek waters, determined visually, are as follows: a) the Lake Tahoe 21 inflow was turbid, milky meltwater from snow and ice melting on shorelands; b) the Shasta Lake inflow was from 22 rocky, tree-covered terrain and was significantly clearer than the lake water (the water pooled into a small pond before 23 flowing into the lake and was easily observed); c) the Donner Lake inflow was from a rocky, tree-lined canyon; d) the 24 Mono Lake inflow was across a mostly barren, rock-strewn shore with loose soil and was notably brown compared 25 to the green lake; and e) the Pinto Lake inflow was from a densely vegetated buffer zone adjacent to farmland. The 26 displacement of the modified waters with respect to the parent water mass are in keeping with these observations, 27 i.e., the waters subjected to turbid or clear inflows had larger or smaller K d (320)/K d (780) ratios, respectively.

Alternative Spectral End Members
The end-member wavelengths used in alternative K d (λ 1 )/K d (λ 2 ) ratios, hereafter Λ λ1 λ2 , follow the combinations 1 first used by Hooker et al. (2013), i.e., the UV-NIR Λ 340 710 pair, as well as the VIS Λ 412 670 pair. Shortly after the start 2 of this study, C -OPS system 021 was upgraded (Table 1), so the Λ 313 875 pair is also available and provides the widest 3 spectral span (562 nm) between end members. A plot of the end-member combinations is presented in Fig. 8, which 4 also includes the linear fits and the root mean square error (RMSE) of the data with respect to the fits. The data 5 in Fig. 8 are only those observations provided in Fig. 3 The fits in Fig. 8 show the end-member pair with the best accuracy is Λ 320 780 , although the Λ 313 875 and Λ 340 710 fits are 8 to within the calibration uncertainty of the radiometers plus inevitable environmental variance, i.e., to within 5%. 9 The slope of the Λ 320 780 fit is to within 1.1% of the original Hooker et al. (2013)  data creates gaps in the data distribution, which partially explains why these data do not yield the lowest RMSE.
14 The Λ 313 875 Fig. 8 data show the variance also increases after the transition from more turbid to clear waters, 15 i.e., a CDOM (440) = 0.02 m −1 , and continues to increase with increasing water clarity. The larger variance as a 16 function of water clarity is caused by the increasing importance of wave-focusing effects coupled with increasing 17 NIR attenuation. Both problems are tractable for Λ 320 780 , but contribute to the difficulty of deriving data products 18 and ultimately producing a stable Λ 313 875 . The increased Λ 340 710 and Λ 412 670 variances are not restricted to the problems 19 described for Λ 313 875 . As end members are brought spectrally closer together, the range of expression available to 20 distinguish two similar but optically different water masses decreases. Consequently, choosing the extrapolation 21 interval is more sensitive to small changes in the parameters that ultimately determine the fit for the extrapolation 22 interval. For legacy end members, clear waters have a lesser range of expression and turbid waters have the greatest, 23 so this problem decreases as turbidity increases, which is seen in the Λ 340 710 and Λ 412 670 Fig. 8 data. Mid-Atlantic Bight and parts of Massachusetts Bay were part of the data used by Hooker et al. (2013), so data from 6 these areas are anticipated to be compliant with the end-member algorithm. The average depth of Chesapeake Bay 7 is relatively shallow (6.4 m) with a signifiant portion (over 24%) less than 2 m deep. Given the extensive contribution 8 of rivers, tributaries, and tides to bay dynamics, resuspension of material is anticipated to be a source of bias in 9 optical properties with respect to end-member algorithms for some bay stations (as shown in Fig. 5).

10
The retained NOMAD data are separated into two regimes: north Chesapeake Bay (NCB) and all other water 11 masses, which consist of 106 stations for each. The dividing line for the NCB is the latitude of the Wicomico River 12 in the Maryland Eastern Shore (slightly north of the Potomac River mouth). The separation is arbitrary and is 13 used to compare the 106 NCB observations from NOMAD with 174 C -OPS K d ratios and a CDOM (440) data pairs 14 obtained in the NCB (not shown in Fig. 2), albeit at different times and locations than the NOMAD data. During 15 data collection, the C -OPS sampling was with system 021 (Table 1) and included notations about in situ conditions 16 useful for establishing a resuspension subcategory, but the procedures predated and were not as rigorous as Sect. 2.5. 17 The C -OPS and NOMAD data plotted in Fig. 9 show general agreement (linearity) of the NOMAD data with 18 respect to the algorithm, which independently confirms the Hooker et al. (2013) algorithmic approach (and as 19 evaluated in more detail herein). Within the narrower turbidity range of the NOMAD and C -OPS NCB data 20 without likely resuspension, there is improved agreement. The C -OPS NCB resuspension data appear properly 21 categorized, because of their relationship to the algorithm (Fig. 5). There is evidence the C -OPS data considered 22 free of resuspension effects nonetheless include some resuspension (e.g., some solid circles in Fig. 9  compared to the 6.2% value determined with C -OPS data (Fig. 8). The more extreme NOMAD values suggest a 27 subcategorization methodology that could be applied to archival data would improve agreement with the algorithm 28 (already demonstrated with the removal of 13 observations using the Λ 412 670 algorithm in Fig. 8).
If the NOMAD data are partitioned into turbid and clear subsets, using a CDOM (440) > 0.2 and a CDOM (440) ≤ 0.2 1 as thresholds, respectively, the fit equation for the turbid Λ 412 670 data is y = 0.3437x − 0.2404. The slope of this turbid 2 NOMAD fit is similar to the corresponding end-member fit presented in Fig. 8 for which y = 0.3504x − 0.1033, and 3 agrees to within 1.9%. The fit for the clear NOMAD data, however, is y = 0.0758x + 0.0648, which is significantly 4 different at the 78.4% level. The 13 NOMAD stations that were not retained out of the original 225 NOMAD stations 5 were in clear waters, which is another indicator that the NOMAD data in clearer waters are problematic. 6 With respect to the algorithm, the increased bias, variance, and 13 negative derived values obtained with NOMAD 7 data in clearer waters suggests the legacy data are degraded by sampling artifacts. Example degradations in legacy 8 data include wave-focusing effects (because of the slower sampling rates), coarser VSR (because of faster descent 9 rates), and deeper extrapolation intervals (because of near-surface data loss from large vertical tilts and large aperture 10 depth offsets). Although some legacy data problems are absent from C -OPS data (e.g., because there is no righting 11 moment when C -OPS sampling begins and C -PrOPS stabilizes the planar orientation of all apertures, some aspects 12 of these limitations are present in the Fig. 8 data, but they are not significant, i.e., they result in a small increase in 13 variance, which slightly degrades algorithm performance. Inclusion of the C -OPS NCB data without resuspension to 14 deriving the Λ 412 670 algorithm (Fig. 8) results in rather small changes to the fit coefficients. The slope is to within 4.3% 15 and the intercept is to within 4.6% (both within the net 5% uncertainty for calibration and environmental variance). 16

Objective versus Subjective Classification
The data set established herein has an extensive number of observations directly suitable for validation exercises 17 however, showed archival data provided a significantly less global data set in terms of spectral expanse and dynamic 22 range as used herein. Archival data usually do not include a subcategory parameter for the observations, e.g.,  MAD has no applicable keyword. Although some subcategories could be determined from geolocation, temporal, and 24 survey information (e.g., a harbor, wetland, alkaline lake, etc.), other influences are not usually established without 25 an observer (e.g., atypical bloom or resuspension caused by a vessel). Consequently, a subcategorization based on 26 the optical measurements alone might be advantageous to the validation process, particularly for archival data.

27
The subcategory approach is evaluated using K d (λ) spectra for the aforementioned 16 subcategories of data, 28 which can be described objectively based on spectral shapes and magnitudes. A small number of observations are 1 excluded to ensure consistency in the determination of all K d (λ) values, e.g., the White Lake data had estimated 2 values in the UV domain, Bear Lake is a unique scattering anomaly created by calcium carbonate particles, ship-3 induced resuspension is anthropogenic in origin, etc. With the additional restriction of wavelength commonality 4 spanning 320-780 nm (Table 1), a total of 1,171 spectra are used for the objective classification analysis.

5
Application of the FCM classification to log-transformed K d (λ) data for the 16 subcategorizations successfully 6 classifies all K d (λ) spectra into five classes based on the Calinski and Harabasz (1974) index (Sect. 2.7). Spectral 7 shapes, as well as their magnitudes, uniquely vary between the five classes and span a continuum of water masses 8 from oceanic and lacustrine case-1 to extreme case-2 inland waters. The continuum of water mass composition is 9 summarized by the centered spectra for the five FCM classes shown in Fig. 10, where N i is the class number set 10 by index i. The diversity achieved in sampling lacustrine water masses is revealed in Fig. 10 by the range of K d (λ) 11 spectra for the example drought-stricken and refilled lakes shown, which are compared with Crater Lake and K w . 12 The refilled lakes are shown to emphasize the complex relationships presented in Figs. 4-7 are not detectable by 13 K d (λ) alone, but require an understanding of the applicable end-member ratio and the a CDOM (440) value. 14 The two lakes at the bottom and top of the dynamic range in Fig. 10, are Crater Lake and White Lake, respectively, 15 with the latter having the largest displacement with respect to the algorithm in Figs. 4-7. Using the inverse of K d (λ) 16 as a proxy for the vertical scale that must be properly sampled (i.e., a sufficient number of observations must 17 be obtained within the vertical scale to derive data products), the vertical scale in Fig. 10 ranges from meters 18 (Crater Lake) to millimeters (White Lake), with the latter only being possible with unprecedented VSR. The shorter 19 wavelengths for White Lake are shown as "estimated" because these wavelengths had to be processed with individual 20 extrapolation intervals to provide the K d estimates, which is not in keeping with all the other data products.

21
The average K d (313) value obtained for deep water (586 m) sampling in Crater Lake was 0.072 m −1 and the 22 coefficient of variation from the six casts was 1.1%. According to Morel et al. (2007) Crater Lake is "extremely" clear 23 water, and the CV shows exceptional reproducibility, i.e., 1% radiometry. Note that as shown in Fig. 3, Crater Lake 24 was sampled twice. A station with higher a CDOM (440) values was conducted in shallower water above submerged 25 moss that grows in large, dense mats. The constituent properties of the water were anticipated to be influenced by 26 the moss matts, and the a CDOM (440) values are elevated with respect to the deeper station.

27
The proportionate composition of each FCM class in Fig. 10 as a function of the original subjective subcategories, 28 wherein contributions less than 5% are reported but not considered significant, is presented in Table 2. Using Fig. 29 1 N 1 The smallest contributor to the dynamic ranges, although arguably accounting for most of the pixels in a global 2 CDOM image, with the K d (λ) maximum in the NIR and the minimum in the blue domain (400-490 nm). 3 The spectral shape is consistent with typical case-1 waters (Morel and Maritorena 2001). The proportional 4 makeup is dominated by the validation quality subcategory (Figs. 3 and 8) at 83% with farther from ice data 5 (Fig. 5), which are almost exclusively from slightly modified case-1 waters, contributing an additional 14%. 6 N 2 The K d (λ) maximum in the NIR is similar to case-1 waters, and the minimum is shifted to longer wavelengths 7 (490-565 nm) from case-1 modifications, principally from proximity to ice effects. The validation quality 8 proportion decreases to 66% and is supplanted with proximity to ice subcategories (Fig. 5). (710 nm), and maximum values compared to N 1 -N 4 that peak in the UV. The resuspension subcategory is 20 dominant at 38%, followed by drought-stricken lakes at 14%, and the validation quality subcategory is reduced 21 to 13%. The remaining principal contributors are wetlands, refilled lakes, and polluted water bodies.

22
The decrease in the percent composition of the validation quality data as a function of increasing class number (N 1 -23 N 5 ) is an indicator of the difficulty of validating an algorithm within increasingly complex waters. The recurring 24 contribution of a relatively small number of principal subjective subcategories to the gradient in optical complexity 25 confirms the subcategory approach has merit and reveals the cause-and-effect relationships of the subcategories. Table 2 also indicates the resolution or granularity for the 15 subcategories of more complex water masses was 27 more nuanced than required. For example, alkaline and hypersaline lakes could be one subcategory, as could refilled 28 and drought-stricken lakes. Combining subcategories does not ensure an eventual convergence with the objective 1 FCM classification, because the latter partitions the processes present in the subcategories into varying degrees of 2 contribution for each identified class. This partitioning involves both direct and indirect evidence, which is perhaps 3 best realized with the original granularity of the 15 subcategories as revealed by considering resuspension processes. 4 The direct evidence of resuspension is provided by the resuspension subcategory (created for this phenomenon). 5 Indirect resuspension is present or likely in multiple subcategories, however. For example, melt water releases particles 6 at the ice edge, bottom deposits in harbors are stirred up by boat propellers, wind and rain redeposit exposed bottom 7 material into drought-stricken lakes, refilled lakes contain suspended material from riverine inflow, tidal and wave 8 action suspend material in shallow wetlands, etc. Table 2 data show the percent contribution of direct and indirect 9 suspension increases with increasing class number N 1 (17%) to N 3 (46%) to N 5 (77%), so a significant part of the 10 cascade towards complexity is correlated with water masses not evolving conservatively, e.g., due to resuspension.

11
A small number of principal subcategories (N p ) significantly determine the K d (λ) classification spectra, although 12 the paucity of observations for some subcategories, in part opportunistic or planned depending on the subcategory, is 13 a largely unknown mitigating factor. Nonetheless, if refilled and drought-stricken lakes plus hypersaline and alkaline 14 lakes, are respectively considered one subcategory rather than two, which appears reasonable based on Fig. 4, five 15 or less principal subcategories determine 97%, 92%, 85%, 80%, and 96% of the composition for N 1 -N 5 , respectively. 16 Without the subcategory scheme, the importance of direct or indirect resuspension would not have emerged with 17 the clarity provided in Table 2. The composition of the FCM classes confirm the applicability of subcategories while 18 revealing which ones are important to each class and which ones are marginally or not important. For the latter, the 19 parent water mass modifier is not significant, as expected, because this small subcategory was created to demonstrate 20 the sensitivity of the in situ optical and laboratory methods. Other subcategories might appear insignificant because 21 of the difficulty of obtaining the sample (e.g., river mouth data are underrepresented) or the paucity of opportunities 22 to collect a sample (e.g., local restrictions to prevent the spread of invasive species).

23
Because metadata to subcategorize NOMAD data are unavailable, the objective classification scheme was applied 24 to NOMAD and to determine the number of data in each class, as follows: N 1 6, N 2 13, N 3 135, N 4 49, N 5 0, plus 25 9 were unclassified. All the turbid data are in classes N 3 and N 4 ; the clear data are in classes N 1 and N 2 , plus 9 are 26 unclassified. This means the slope of the aforementioned clear partition (Sect. 3.6) had 19 points that were classified 27 and 9 that were not, and likely accounts for the poor performance with respect to the algorithm. It also suggests 28 the measurements were the issue with the NOMAD data, because some spectra could not be classified.

Discussion
The optical data herein had a near-surface VSR less than 1 mm, which allowed data products spanning 313-875 nm 1 while encompassing a global perspective of water masses as described by approximately three decades of generalized 2 water properties. The validation approach was based on the concept that water masses evolving conservatively (i.e., 3 free from stressors that might alter the natural range in the gradient of a constituent) are suitable for validating the 4 original Hooker et al. (2013) inversion algorithm for deriving a CDOM (440) from K d (λ) spectral end members.

5
The identification of 15 subcategories that were likely not evolving conservatively yielded 609 validation quality 6 data points spanning extremely clear to highly turbid waters sampled within the open ocean, coastal zone, and inland 7 waters. The new data adhered to the original Λ 320 780 algorithm (Fig. 3) with an RPD of 0.02% and an APD of 3.86%. 8 Alternative spectral end members (e.g., Λ 313 875 , Λ 340 710 , and Λ 412 670 ) had increasingly larger RMSE values, but were to 9 within the calibration uncertainty of the radiometers plus inevitable environmental variance (i.e., a net uncertainty 10 to within 5%) except for the narrowest spectral span of legacy end members, which was 6.2% for Λ 412 670 (Fig. 8).

11
Although no data archive exists with the spectral and spatial coverage used herein, NOMAD Λ 412 670 data showed 12 general agreement (linearity) with respect to the original Hooker et al. (2013) algorithm, and independently confirmed 13 the algorithmic approach. There was also general agreement between the NOMAD NCB data and the correspond-14 ing C -OPS NCB data without likely resuspension, wherein the C -OPS NCB resuspension data appeared properly 15 categorized, because of their relationship with the algorithm (Fig. 5).

16
The near-surface VSR is quantitative evidence of the successful mitigation of a variety of sampling difficulties 17 (e.g., large aperture offsets, righting moment instabilities, wave focusing effects, etc.). The high VSR achieved with 18 C -OPS resulted in an increased sensitivity for deriving data products in turbid water masses with vertical scales 19 on the order of millimeters (Fig. 10). This sensitivity allowed small localized changes in a parent water mass to 20 be distinguished (Figs. 5 and 7). The laboratory analyses were similarly sensitive, so the optical determinations of 21 in-water constituents using field measurements were commensurately as effective as laboratory analyses. The ability 22 to distinguish small differences in water properties ensured the discrimination of 15 subcategories, wherein each 23 represented a more complex water mass not evolving conservatively and not automatically used for validation.

24
Plots of the a CDOM (440) and Λ 320 780 relationships for the 15 subcategories of more complex water bodies revealed 25 some data were significantly different with respect to the original algorithm and others that were not . The wide range in complexity is largely the result of the substantial effort that was made to adopt a global 27 perspective and sample the greatest diversity of water bodies possible, some of which were very difficult to access 28 and sometimes required shoreline launches of the optical instrumentation that C -PrOPS made possible, e.g., severely 1 drought-stricken and refilled lakes (high and dry, flooded, or debris-blocked boat ramps) plus hypersaline and soda 2 lakes (Mono Lake, Salton Sea, Borax Lake, and Soda Lake had no boat ramps). 3 The accuracy of the algorithm as a function of including increasing proportions of the 15 subcategories of more 4 complex water masses not necessarily appropriate for validation exercises , because they were likely not 5 evolving conservatively, is explored by expanding the 609 validation quality observations of Λ 320 780 end members (Figs. 6 3 and 8) to include the following subcategories (hereafter referred to as the second data set): inflows to a parent water 7 mass that are not hypersaline or drought-stricken lakes, closer to or farther from ice edge proximity, river mouth, 8 resuspension (but not including Bear Lake and the ship-induced RWC Channel resuspension), atypical blooms, HABs, 9 and wetlands. This second data set has 930 observations, the linear fit is y = 0.2511x − 0.0046, the RMSE is 5.7%, 10 and the new slope is to within 1.7% of the original value presented by Hooker et al. (2013), i.e., y = 0.2556x − 0.0030. 11 The reason the slope of the second data set is not significantly different than the original fit coefficients is the 12 data that were added are situated above and below the distribution of the validation quality data set, as shown in 13 Figs. 4 and 5. The C -OPS NCB data without resuspension cluster on or below the algorithmic relationship (Fig. 9). 14 If these data are added to the second data set used to derive the Λ 320 780 algorithm (the 320 nm and 780 nm wavelengths 15 were always part of system 021 as shown in Table 1), the resulting slope and intercept is y = 0.2561x − 0.0076, which 16 is to within 2.0% of the slope determined for the second data set (y = 0.2511x − 0.0046).

17
If a third data set is created by adding drought-stricken and refilled lakes to the second data set, but not including 18 White Lake (which had some estimated data products), this third data set has 1,044 observations. The linear fit 19 is y = 0.2249x + 0.0044, the RMSE is 6.8%, and the new slope is reduced (as expected, because the added data 20 are all below the algorithmic relationship and primarily turbid (Fig. 4), but still within 10.4% of the original value 21 presented by Hooker et al. (2013). The exclusion of White Lake, as well as hypersaline, alkaline, and polluted water 22 bodies, from the third data set is for a practical reason: they are, or have significant characteristics of, extreme water 23 masses and White Lake rarely exists.

24
There are other lakes presented herein that can be considered extreme, e.g., other endorheic lakes (e.g., Pyramid 25 Lake, Eagle Lake, etc.), plus shallow lakes in high wind areas wherein bottom material is resuspended on a near-26 continuous basis (e.g., Washoe Lake and Little Washoe Lake). If all subcategories, except extreme lacustrine water 27 bodies, are used to create a fourth data set it has 1,086 observations, i.e., almost 90% of the 1,230 maximum and 28 93% of the data used in Table 2 to create the five objective FCM classifications. The linear fit of this fourth more 29 comprehensive data set is y = 0.2379x − 0.0049, the RMSE is 6.2%, and the new slope is not reduced as much as in 1 the third data set (as expected) and is to within 6.9% of the original value presented by Hooker et al. (2013).

2
A fifth data set, using classes N 1 -N 4 spans the entire a CDOM (440) dynamic range (Table 2). This alternative 3 for an arguably global algorithm yields a linear fit of y = 0.2317x − 0.0053, the RMSE is 5.3%, and the slope is to 4 within 9.3% of the original Hooker et al. (2013) value. This result is significantly similar to the subjective results 5 of the fourth data set (the slopes agree to within 2.7%). Consequently, the robustness of the algorithm is directly 6 supported by the combination of subjective and objective classifications, with the latter using FCM. 7 With the exception of the third data set that added primarily turbid data exclusively below the algorithmic 8 relationship (drought-stricken and refilled lakes), all of the results from the expanded data sets are rather indistin-9 guishable from the original Λ 320 780 fit provided by Hooker et al. (2013), wherein y = 0.2556x − 0.0030, or the separate 10 validation quality data set presented in Fig. 8, for which y = 0.2583x − 0.0053. Ignoring the third data set, the 11 x-intercept for the expanded data sets is approximately equal to what can be expected for pure water, i.e., 0.02. 12 The close agreement of the various expanded data sets with adherence to the pure-water limit is another significant 13 confirmation of the algorithmic approach using spectral end members. 14 The robustness is further established by creating a so-called universal algorithm, which is assumed to mean that 15 any water mass wherein an optical profiler can be deployed is expected to be part of the evaluation of the end-16 member approach. In this case, the universal algorithm is constructed from all the data from all subcategories, and 17 is distinguished from the global perspective in that the universal dataset includes a far greater proportion of complex 18 water types than exist globally. The linear fit of this universal data set is y = 0.2206x + 0.0088, the RMSE is 7.5%, 19 and the new slope is to within 13.7% of the original value presented by Hooker et al. (2013). The removal of data that 20 are, or have significant characteristics of, extreme water masses improves performance. For example, if hypersaline, 21 alkaline, and polluted lakes are removed, the linear fit of this sixth comprehensive data set is y = 0.2250x + 0.0024, 22 the RMSE is 6.8%, and the new slope is to within 12.0% of the original value presented by Hooker et al. (2013). 23 This study used three different laboratory methods to determine a CDOM (440) from water samples, and seven 24 different optical instrument suites to determine K d (λ). Despite the agreement between the Λ 320 780 fits and their x-25 intercepts, the possibility the results are the result of an unidentified stochastic process has not been addressed. The 26 latter is not likely for the optical data, because the radiometers were calibrated at one facility and deployed with 27 strict adherence to the Protocols using the same acquisition software. Furthermore, data products were derived using 28 the same processing software with one operator, and the variance in optical data products is shown in Figs. 3-8. 29 A Global End-Member Approach to Derive a CDOM (440) from Near-Surface Optical Measurements The capabilities of the above-and in-water radiometers for C -OPS systems 021 (with C -PrOPS) and 039 (without 1 C -PrOPS) were intercompared to next-generation hybridspectral instruments (Hooker et al. 2018c) containing a 2 hyperspectral detector system plus 18 fixed wavelengths (system 038 in Table 1). The comparisons showed an 3 agreement of 4.2-4.8%, which is to within the calibration uncertainty (2.3-2.7%) plus natural variability (i.e., a net 4 uncertainty to within 5%), as long as the stability threshold for a backplane without thrusters or the noise limit of 5 the hyperspectral sensor was not exceeded (Hooker et al. 2018b and2018c, respectively). 6 Although field data demonstrate small changes in parent water mass modifications can be discriminated (Figs. 7 5 and 7), the laboratory methods and instruments were not systematically intercalibrated to establish an overall 8 uncertainty. In addition, the laboratory methods included different temperature controls and storage procedures for 9 the water samples, as well as application of null-point corrections in different spectral ranges. A lack of systematic 10 intercalibration is not a significant detraction, because it is a common difficulty when combining observations from 11 databases (e.g., NOMAD data), and methods exist to nonetheless determine the efficacy of the combined data. 8. Briefly, for each optical and water sampling, normality of distribution for a CDOM (440) was created using the 15 measured value as the mean (µ) and 7% of the measured value as the standard deviation (σ). Similarly, normality 16 of distribution for Λ 320 780 was created using the measured value as µ and 5% as σ. A lower σ percentage was used for 17 the latter, because the original data set retained each optical cast whereas only one water sample was obtained, so 18 the variability of the optical data was already significantly represented . 19 For each original data pair, 10 5 variations were prepared for both Λ 320 780 and a CDOM (440). Of these data, 10 3 were 20 randomly selected and the mean computed for each Λ 320 780 and a CDOM (440) pair. This exercise was repeated 10 3 times 21 (bootstrap) and an overall µ, and σ were obtained. For the different combinations of applying µ ± σ to bound the 22 dispersion of the original data, the resulting linear fits showed the new algorithm slopes changed by 0.3-1.1% with 23 respect to the Λ 320 780 validation quality data set in Fig. 8. For all combinations of µ±σ, the maximum RMSE was 1.1%, 24 which is similar to the 1.2% value for Λ 320 780 in Fig. 8. Consequently, the use of three different laboratory methods to 25 determine a CDOM (440), which were not intercalibrated, does not significantly influence the results presented here.

26
The state-of-the-art (approximately 1%) accuracy achieved with the Λ 320 780 end members (Fig. 8) is, therefore, due 27 to the strict adherence to sampling and processing protocols coupled with an unprecedented VSR for the optical data. 28 This combination also resulted in all end-member pairs-including the legacy (VIS) Λ 412 670 pair-having a superior 29 accuracy compared to many common global inversion algorithms. For example, a CDOM (λ) algorithms based on the 1 water-leaving radiance, L W (λ), and its normalized forms have an RMSE exceeding 10% or more (Mannino et al. 2 2008 and 2014), and the ocean color (OC) chlorophyll a variants that use band ratios (O'Reilly et al. 1998 and2000), 3 typically exceed 20% or more while generally excluding complex and inland waters as used herein. 4 The robustness of the end-member approach, is further confirmed by how the expanded data sets (up to six), 5 which had respectively increased amounts of water masses that were not evolving conservatively, nonetheless yielded 6 fits with RMSE values outperforming the aforementioned L W (λ) algorithms. In regards to which of the variants best 7 represents a global perspective, the Λ λ1 λ2 pairs in Fig. 8 are considered appropriate. Screening of newly collected or 8 archival data, e.g., NOMAD data (Sect. 3.7), with respect to a selected algorithm can be accomplished by initially 9 flagging data points more than 12% from the expected relationship, and then more carefully examining those points 10 using both objective and subjective criteria (based on available metadata) to determine whether the results are 11 expected, or are more likely to indicate a problem with data collection procedures.
If an algorithm is to be applied to a water mass that is not evolving conservatively, an individualized relationship 13 between a CDOM (λ) and Λ λ1 λ2 should be established, especially for an extreme water mass, e.g., drought-stricken. This 14 maintains the accuracy of the global relationship while providing a mechanism for improving the study of noncon-15 servative waters. Identification of water masses not evolving conservatively can likely be determined using K d (λ) 16 values determined from in-water optical (e.g, C -OPS) data without a need for laboratory analyses (Fig. 9). With 17 additional research to produce K d (λ) for all wavelengths, e.g., expanding upon Cao et al. (2014), the identification 18 of water masses evolving conservatively or not can be made from above-water observations. This determination 19 provides a sensitive indicator of water masses subjected to stresses influencing water quality, e.g., drought. The 20 onset of next-generation satellites, e.g., the Japanese Second-generation Global Imager (SGLI) mission (Honda et al. 21 2012), offers unique opportunities for such monitoring because of the expanded spatial and spectral domain, with the 22 latter allowing improved Λ λ1 λ2 algorithm accuracy with respect to legacy missions, e.g., Moderate Resolution Imaging 23 Spectroradiometer (MODIS). While planned high spectral resolution sensors, such as the Plankton, Aerosol, Cloud, 24 ocean Ecosystem (PACE) and Surface Biology and Geology (SBG) missions, may support more sophisticated re-25 trievals of parameters like CDOM, the simplified approach provided by end-member analysis can be used with both 26 legacy and next-generation sensors, thereby providing continuity in space and time, as well as a capability to generate 27 high quality in-water data with a simplified measurement approach (assuming strict adherence to the Protocols).  quality (e.g., White Lake, Bear Lake, ship-induced resuspension, etc.). The five classes N i , where i is the class 2 number index, are shown with the number of K d (λ) spectra, N s , within each class in parentheses, as well as the 3 percent composition of the original 16 subjective subcategories equalling or exceeding a 1% contribution threshold 4 for each class. The principal subcategories in each class, i.e., the most numerous (approximately 5% contribution 5 or more), are shown in bold typeface. After the subcategories, the number of principal subcategories (N p ) and the 6 number of all subcategories (N a ) with a 1% composition or more are summarized in slanted typeface followed by the 7 percent composition from all the principal subcategories. The extent of the dynamic range in percent for the optical 8 and biogeochemical data are shown in the last two lines as a function of applying successive class numbers.  and support electronics; e) rotating V-block for pitch adjustment; f) two-point harness attachment; g) hydrobaric 3 buoyancy chamber, which accommodates up to 3 compressible bladders; h) slotted flotation and i) bronze weights 4 for buoyancy and roll adjustment; j) water temperature probe and k) pressure transducer port on the radiance end 5 cap; l) conductivity sensor; m) electronics module; n) digital thruster (one of two, on each side); and o) thruster 6 guard. The aggregator and support electronics control the 19 microradiometers as a single device. The side bumpers 7 and thruster guards protect the radiometers and digital thrusters from unanticipated side impacts, respectively.  15.0% dispersion. The headwaters of the San Francisco Bay Redwood Creek (RWC) Channel, which is surrounded 16 by wetlands, is the most turbulent coastal water mass and has the highest a CDOM (440) value.
17 Fig. 4. The new data from lacustrine water bodies that were drought-stricken (magenta diamonds), refilled after 18 drought (blue triangles), alkaline (red squares), or hypersaline (green circles) in relation to the original algorithm 19 (gray circles). Each data point represents a single water sample with multiple optical casts for each, which results in 20 a series of results (typically 3-6) along the x-axis. Columbia and Umqua River data that are upstream of the river mouth. diamonds, respectively). The NOMAD NCB data are compared to C -OPS NCB data obtained at different times 9 and locations with the latter separated into two categories, wherein one is likely subjected to bottom resuspension 10 (light blue open circles) and the other is not (dark blue solid circles).
11 Fig. 10. The centered K d (λ) spectra of the five classes (N i ) determined from an objective FCM classification of the 12 data presented in Figs. 3-8 with a few omissions for data consistency (e.g., White Lake, Bear Lake, ship-induced 13 resuspension, etc.) and shown with respect to K w . Example K d spectra from drought-stricken and refilled lakes plus 14 Crater Lake, obtained by averaging the results from multiple optical casts, are also shown to demonstrate the more 15 than three decades of dynamic range in turbidity that was sampled. The shorter wavelengths for White Lake (open 16 dark red squares) required individual wavelength processing to provide the estimated K d values, whereas all other 17 data were obtained with a single processing.   Figure 10