Generalized ocean color inversion model for retrieving marine inherent optical properties

Inherent Optical Properties in Lake Taihu Derived from VIIRS Satellite Observations

Using in situ remote sensing reflectance and inherent optical property (IOP) measurements, a near-infrared (NIR)-based IOP algorithm is developed and tuned for Lake Taihu, in order to derive the particle backscattering coefficient bbp(λ), total absorption coefficient at(λ), dissolved and detrital absorption coefficient adg(λ), and phytoplankton absorption coefficient aph(λ), with satellite observations from the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP). The IOP algorithm for Lake Taihu has a reasonably good accuracy. In fact, the determination coefficients between the retrieved and in situ IOPs are 0.772, 0.638, and 0.487 for at(λ), adg(λ), and aph(λ), respectively. The IOP products in Lake Taihu that have been derived from VIIRS-SNPP observations show significant spatial and temporal variations. Southern Lake Taihu features enhanced bbp(λ) and adg(λ), while northern Lake Taihu shows higher aph(λ). The seasonal and interannual variability of adg(λ) and bbp(λ) in Lake Taihu is quantified and characterized with the highest bbp(λ) and adg(λ) in the winter, and the lowest in the summer. In the winter, bbp(443) and adg(443) can reach over ~1.5 and ~5.0 m−1, respectively, while they are ~0.5–1.0 and ~2.0 m−1 in the summer. This study shows that in Lake Taihu adg(λ) is the most significant IOP, while aph(λ) is the least in terms of the IOP values and contributions to remote sensing reflectance. The highest bbp(λ) and adg(λ) occurred in the winter between 2017–2018, and the lowest bbp(λ) and adg(λ) occurred in the summer of 2014. In comparison, the seasonal and interannual variability of mean aph(λ) for Lake Taihu is less significant, even though enhanced seasonal and interannual variability can be found in some parts of Lake Taihu, such as in the northern Lake Taihu region.

Real Time HABs Mapping Using NASA Glenn Hyperspectral Imager

The hyperspectral imaging system (HSI) developed by the NASA Glenn Research Center was used from 2015-2017 to collect high spatial resolution data over Lake Erie and the Ohio River. Paired with a vicarious correction approach implemented by the Michigan Tech Research Institute, radiance data collected by the HSI system can be converted to high quality reflectance data which can be used to generate near-real time (within 24 hours) products for the monitoring of harmful algal blooms using existing algorithms. The vicarious correction method relies on imaging a spectrally constant target to normalize HSI data for atmospheric and instrument calibration signals. A large asphalt parking lot near the Western Basin of Lake Erie was spectrally characterized and was determined to be a suitable correction target. Due to the HSI deployment aboard an aircraft, it is able to provide unique insights into water quality conditions not offered by space-based solutions. Aircraft can operate under cloud cover and flight paths can be chosen and changed on-demand, allowing for far more flexibility than space-based platforms. The HSI is also able to collect data at a high spatial resolution (^~1 m), allowing for the monitoring of small water bodies, the ability to detect small patches of surface scum, and the capability to monitor the proximity of blooms to targets of interest such as water intakes. With this new rapid turnaround time, airborne data can serve as a complementary monitoring tool to existing satellite platforms, targeting critical areas and responding to bloom events on-demand.

Effects of forward models on the semi-analytical retrieval of inherent optical properties from remote sensing reflectance

Inherent optical properties (IOPs) play a key role in modulating an aquatic light field; they are the core link for remotely sensing water constituents based on ocean color remote sensing. Many semi-analytical algorithms (SAAs) have been developed to obtain IOPs from remote sensing reflectance (R_rs) data; these algorithms require a forward model (FM) to link the IOPs to R_rs. Most currently available SAAs use the FM presented by Gordon et al.[J. Geophys. Res.93, 10909 (1988)JGREA20148-022710.1029/JD093iD09p10909] (G88 hereafter) without knowledge of how other models would impact the retrieval of IOPs from R_rs. This study evaluates the effects of two popular SAAs, namely, the quasi-analytical algorithm (QAA) and the generalized IOP algorithm (GIOP), combined with six different FMs on the retrieval of IOPs from a synthetic data set generated with Hydrolight software. The results indicated that different FMs can have quite different effects on the computed R_rs(λ), and the effects were not uniform across the R_rs spectrum. Of the six FMs tested, G88 and P05 [Appl. Opt.44, 1236 (2005)APOPAI0003-693510.1364/AO.44.001236] produced the best estimates of R_rs(λ) at 350, 440, and 550 nm in both oceanic and coastal sub-datasets; they also were less impacted by changes in the particle phase function. M02 also produced a good estimation of R_rs but only at 440 nm, and L04 performed well only in the oceanic condition. When the two SAAs were combined with the six FMs, in the oceanic condition, QAA and GIOP combined with M02 (QM02 and GM02) provided better quality for the absorption coefficient [a(λ)] at 350, 440, and 550 nm when compared with the SAAs combined with the other models. However, for the retrieval of the particle backscattering coefficient [b_bp(λ)] in the oceanic condition, QAA and GIOP combined with L04 (QL04 and GL04) performed better than the others, and GL04 always provided a better estimation of b_bp(λ) than QL04. In the coastal condition, QAA and GIOP combined with G88 or P05 produced slightly better quality of IOPs compared with the other four FMs. Compared with GIOP in the coastal condition, QAA combined with G88 or P05 always showed better quality of retrieval of a(λ) but weaker quality of retrieval of b_bp(λ).

Evaluation of Weighting Average Functions as a Simplification of the Radiative Transfer Simulation in Vertically Inhomogeneous Eutrophic Waters

Current water color remote sensing algorithms typically do not consider the vertical variations of phytoplankton. Ecolight with a radiative transfer program was used to simulate the underwater light field of vertical inhomogeneous waters based on the optical properties of a eutrophic lake (i.e., Lake Chaohu, China). Results showed that the vertical distribution of chlorophyll-a (Chla(z)) can considerably affect spectrum shape and magnitude of apparent optical properties (AOPs), including subsurface remote sensing reflectance in water (rrs(λ, z)) and the diffuse attenuation coefficient (Kx(λ, z)). The vertical variations of Chla(z) changed the spectrum shapes of rrs(λ, z) at the green and red wavelengths with a maximum value at approximately 590 nm, and changed the Kx(λ, z) from blue to red wavelength range with no obvious spectral variation. The difference between rrs(λ, z) at depth z m and its asymptotic value (Δrrs(λ, z)) could reach to ~78% in highly stratified waters. Diffuse attenuation coefficient of downwelling plane irradiance (Kd(λ, z)) had larger vertical variations, especially near water surface, in highly stratified waters. Three weighting average functions performed well in less stratified waters, and the weighting average function proposed by Zaneveld et al., (2005) performed best in highly stratified waters. The total contribution of the first three layers to rrs(λ, 0−) was approximately 90%, but the contribution of each layer in the water column to rrs(λ, 0−) varied with wavelength, vertical distribution of Chla(z) profiles, concentration of suspended particulate inorganic matter (SPIM), and colored dissolved organic matter (CDOM). A simple stratified remote sensing reflectance model considering the vertical distribution of phytoplankton was built based on the contribution of each layer to rrs(λ, 0−).

Calibration of an airborne oceanographic lidar using ocean backscattering measurements from space

We used satellite measurements of the optical backscattering coefficient to calibrate the signal from an airborne oceanographic lidar. This technique provided the radiometric calibration for the lidar signal and a local estimate of the ratio of the particulate backscattering coefficient, bbp, to the volume scattering function at the scattering angle of 180°, βp(180). Results using an ordinary regression, a reduced major axis regression, and a least squares bisector suggest that either of the latter two provided a better result than an ordinary regression. The statistical errors in the two recommended regressions, the difference in calibrations factors between them, and the difference between these and a laboratory calibration were all less than 5%.

Sensitivity of inherent optical properties from ocean reflectance inversion models to satellite instrument wavelength suites

The Earth science community seeks to develop climate data records (CDRs) from satellite measurements of ocean color, a continuous data record which now exceeds twenty years. Space agencies will launch additional instruments in the coming decade that will continue this data record, including the NASA PACE spectrometer. Inherent optical properties (IOPs) quantitatively describe the absorbing and scattering constituents of seawater and can be estimated from satellite-observed spectroradiometric data using semi-analytical algorithms (SAAs). SAAs exploit the contrasting optical signatures of constituent matter at spectral bands observed by satellite sensors. SAA performance, therefore, depends on the spectral resolution of the satellite spectroradiometer. A CDR spanning SeaWiFS, MODIS, OLCI, and PACE, for example, would include IOPs derived using varied wavelength suites if all available wavelengths were considered. Here, we explored differences in derived IOPs that stem simply from the use of (eight) different wavelength suites of input radiometric measurements. Using synthesized data and SeaWiFS Level-3 mission-long composites, we demonstrated equivalent SAA performance for all wavelength suites, but that IOP retrievals vary by several percent across wavelength suites and as a function of water type. The differences equate to roughly ≤ 6, 12, and 7% for a_dg (443), a_ph (443), and b_bp (443), respectively, for waters with C_a ≤ 1 mg m^-3. These values shrink for sensors with similar wavelength suites (e.g., SeaWiFS, MODIS, and MERIS) and rise to substantially larger values for higher C_a waters. Our results also indicate that including 400 nm (in the case of OLCI) influences the derived IOPs, using longer wavelengths (>600 nm) influences the derived IOPs when there is a red signal, and, including additional spectral information shows potential for improved IOP estimation, but not without revisiting SAA parameterizations and execution. While modest in scope, we believe this study contributes to the knowledge base for CDR development. The implication of ignoring such an analysis as CDRs continue to be developed is a prolonged inability to distinguish between algorithmic and environmental contributions to trends and anomalies in the IOP time-series.

Evidence of Environmental Changes Caused by Chinese Island-Building

This paper quantifies environmental effects of island-building operations in the South China Sea, which result from dredging and can negatively impact marine flora, fauna, and ecosystems. The extent of the damage caused by island-creation is believed to be large, as the South China Sea reefs support the largest concentration of marine biodiversity on Earth. Through use of satellite imagery, we investigate the island-construction on Mischief Reef in the South China Sea, showing backscatter increases of up to 350% in waters surrounding the reef, with plumes of excess sediment exceeding 250 km² at times during island-construction, and the cumulative area impacted by dredging exceeding 1,200 km². Comparison of satellite-derived chlorophyll-a, backscatter, absorption and remote sensing reflectance at 412 nm suggest that dredging activities led to a decrease in biological health of the region resulting from the smothering of natural benthic habitats and reef complexes with sediment. We anticipate this ex post facto quantification of the connectivity between island-construction, large particulate plumes and a decrease in absorption related to marine life in the water column to establish a starting point for further study into ecosystem impact. The potential associations between these damages and a long-term reduction in ocean life and resources could serve inter-governmental bodies with a baseline metric for evaluating the level of damage caused. This may result in both forward-looking deterrent policies that limit island-building as well as backward-looking compensation.

Deriving inherent optical properties from classical water color measurements: Forel-Ule index and Secchi disk depth

Secchi disk depth (ZSD) and Forel-Ule index (FUI) are the two oldest and easiest measurements of water optical properties based on visual determination. With an overarching objective to obtain water inherent optical properties (IOPs) using these historical measurements, this study presents a model for associating remote-sensing reflectance (Rrs) with FUI and ZSD. Based upon this, a scheme (FZ2ab) for converting FUI and ZSD to absorption (a) and backscattering coefficients (bb) is developed and evaluated. For a data set from HydroLight simulations, the difference is <11% between FZ2ab-derived a and known a, and <28% between FZ2ab-derived bb and known bb. Further, for a data set from field measurements, the difference is < 30% between FZ2ab-derived a and measured a. These results indicate that FZ2ab can bridge the gap between historical measurements and the focus of IOP measurements in modern marine optics, and potentially extend our knowledge on the bio-optical properties of global seas to the past century through the historical measurements of FUI and ZSD.

Optical Classification of the Remote Sensing Reflectance and Its Application in Deriving the Specific Phytoplankton Absorption in Optically Complex Lakes

Optical water types (OWTs) were identified from remote sensing reflectance (Rrs(λ)) values in a field-measured dataset of several large lakes in the lower reaches of the Yangtze and Huai River (LYHR) Basin. Four OWTs were determined from normalized remote sensing reflectance spectra (NRrs(λ)) using the k-means clustering approach, and were identified in the Sentinel 3A OLCI (Ocean Land Color Instrument) image data over lakes in the LYHR Basin. The results showed that 1) Each OWT is associated with different bio-optical properties, such as the concentration of chlorophyll-a (Chla), suspended particulate matter (SPM), proportion of suspended particulate inorganic matter (SPIM), and absorption coefficient of each component. One optical water type showed an obvious characteristic with a high contribution of mineral particles, while one type was mostly determined by a high content of phytoplankton. The other types belonged to the optically mixed water types. 2) Class-specific Chla inversion algorithms performed better for all water types, except type 4, compared to the overall dataset. In addition, class-specific inversion algorithms for estimating the Chla-specific absorption coefficient of phytoplankton at 443 nm (a*ph(443)) were developed based on the relationship between a*ph(443) and Chla of each OWT. The spatial variations in the class-specific model-derived a*ph(443) values were illustrated for 2 March 2017, and 24 October 2017. 3) The dominant water type and the Shannon index (H) were used to characterize the optical variability or similarity of the lakes in the LYHR Basin using cloud-free OLCI images in 2017. A high optical variation was located in the western and southern parts of Lake Taihu, the southern part of Lake Hongze, Lake Chaohu, and several small lakes near the Yangtze River, while the northern part of Lake Hongze had a low optical diversity. This work demonstrates the potential and necessity of optical classification in estimating bio-optical parameters using class-specific inversion algorithms and monitoring of the optical variations in optically complex and dynamic lake waters.

The Assessment of Landsat-8 OLI Atmospheric Correction Algorithms for Inland Waters

The OLI (Operational Land Imager) sensor on Landsat-8 has the potential to meet the requirements of remote sensing of water color. However, the optical properties of inland waters are more complex than those of oceanic waters, and inland atmospheric correction presents additional challenges. We examined the performance of atmospheric correction (AC) methods for remote sensing over three highly turbid or hypereutrophic inland waters in China: Lake Hongze, Lake Chaohu, and Lake Taihu. Four water-AC algorithms (SWIR (Short Wave Infrared), EXP (Exponential Extrapolation), DSF (Dark Spectrum Fitting), and MUMM (Management Unit Mathematics Models)) and three land-AC algorithms (FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes), 6SV (a version of Second Simulation of the Satellite Signal in the Solar Spectrum), and QUAC (Quick Atmospheric Correction)) were assessed using Landsat-8 OLI data and concurrent in situ data. The results showed that the EXP (and DSF) together with 6SV algorithms provided the best estimates of the remote sensing reflectance (Rrs) and band ratios in water-AC algorithms and land-AC algorithms, respectively. AC algorithms showed a discriminating accuracy for different water types (turbid waters, in-water algae waters, and floating bloom waters). For turbid waters, EXP gave the best Rrs in visible bands. For the in-water algae and floating bloom waters, however, all water-algorithms failed due to an inappropriate aerosol model and non-zero reflectance at 1609 nm. The results of the study show the improvements that can be achieved considering SWIR bands and using band ratios, and the need for further development of AC algorithms for complex aquatic and atmospheric conditions, typical of inland waters.

Atmospheric Correction of Satellite Ocean-Color Imagery During the PACE Era

The Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission will carry into space the Ocean Color Instrument (OCI), a spectrometer measuring at 5nm spectral resolution in the ultraviolet (UV) to near infrared (NIR) with additional spectral bands in the shortwave infrared (SWIR), and two multi-angle polarimeters that will overlap the OCI spectral range and spatial coverage, i. e., the Spectrometer for Planetary Exploration (SPEXone) and the Hyper-Angular Rainbow Polarimeter (HARP2). These instruments, especially when used in synergy, have great potential for improving estimates of water reflectance in the post Earth Observing System (EOS) era. Extending the top-of-atmosphere (TOA) observations to the UV, where aerosol absorption is effective, adding spectral bands in the SWIR, where even the most turbid waters are black and sensitivity to the aerosol coarse mode is higher than at shorter wavelengths, and measuring in the oxygen A-band to estimate aerosol altitude will enable greater accuracy in atmospheric correction for ocean color science. The multi-angular and polarized measurements, sensitive to aerosol properties (e.g., size distribution, index of refraction), can further help to identify or constrain the aerosol model, or to retrieve directly water reflectance. Algorithms that exploit the new capabilities are presented, and their ability to improve accuracy is discussed. They embrace a modern, adapted heritage two-step algorithm and alternative schemes (deterministic, statistical) that aim at inverting the TOA signal in a single step. These schemes, by the nature of their construction, their robustness, their generalization properties, and their ability to associate uncertainties, are expected to become the new standard in the future. A strategy for atmospheric correction is presented that ensures continuity and consistency with past and present ocean-color missions while enabling full exploitation of the new dimensions and possibilities. Despite the major improvements anticipated with the PACE instruments, gaps/issues remain to be filled/tackled. They include dealing properly with whitecaps, taking into account Earth-curvature effects, correcting for adjacency effects, accounting for the coupling between scattering and absorption, modeling accurately water reflectance, and acquiring a sufficiently representative dataset of water reflectance in the UV to SWIR. Dedicated efforts, experimental and theoretical, are in order to gather the necessary information and rectify inadequacies. Ideas and solutions are put forward to address the unresolved issues. Thanks to its design and characteristics, the PACE mission will mark the beginning of a new era of unprecedented accuracy in ocean-color radiometry from space.

Progress in Forward-Inverse Modeling Based on Radiative Transfer Tools for Coupled Atmosphere-Snow/Ice-Ocean Systems: A Review and Description of the AccuRT Model

A tutorial review is provided of forward and inverse radiative transfer in coupled atmosphere-snow/ice-water systems. The coupled system is assumed to consist of two adjacent horizontal slabs separated by an interface across which the refractive index changes abruptly from its value in air to that in ice/water. A comprehensive review is provided of the inherent optical properties of air and water (including snow and ice). The radiative transfer equation for unpolarized as well as polarized radiation is described and solutions are outlined. Several examples of how to formulate and solve inverse problems encountered in environmental optics involving coupled atmosphere-water systems are discussed in some detail to illustrate how the solutions to the radiative transfer equation can be used as a forward model to solve practical inverse problems.

Ocean Color Analytical Model Explicitly Dependent on the Volume Scattering Function

An analytical radiative transfer (RT) model for remote sensing reflectance that includes the bidirectional reflectance distribution function (BRDF) is described. The model, called ZTT (Zaneveld-Twardowski-Tonizzo), is based on the restatement of the RT equation by Zaneveld (1995) in terms of light field shape factors. Besides remote sensing geometry considerations (solar zenith angle, viewing angle, and relative azimuth), the inputs are Inherent Optical Properties (IOPs) absorption a and backscattering bb coefficients, the shape of the particulate volume scattering function (VSF) in the backward direction, and the particulate backscattering ratio. Model performance (absolute error) is equivalent to full RT simulations for available high quality validation data sets, indicating almost all residual errors are inherent to the data sets themselves, i.e., from the measurements of IOPs and radiometry used as model input and in match up assessments, respectively. Best performance was observed when a constant backward phase function shape based on the findings of Sullivan and Twardowski (2009) was assumed in the model. Critically, using a constant phase function in the backward direction eliminates a key unknown, providing a path toward inversion to solve for a and bb. Performance degraded when using other phase function shapes. With available data sets, the model shows stronger performance than current state-of-the-art look-up table (LUT) based BRDF models used to normalize reflectance data, formulated on simpler first order RT approximations between rrs and bb/a or bb/(a + bb) (Morel et al., 2002; Lee et al., 2011). Stronger performance of ZTT relative to LUT-based models is attributed to using a more representative phase function shape, as well as the additional degrees of freedom achieved with several physically meaningful terms in the model. Since the model is fully described with analytical expressions, errors for terms can be individually assessed, and refinements can be readily made without carrying out the gamut of full RT computations required for LUT-based models. The ZTT model is invertible to solve for a and bb from remote sensing reflectance, and inversion approaches are being pursued in ongoing work. The focus here is with development and testing of the in-water forward model, but current ocean color remote sensing approaches to cope with an air-sea interface and atmospheric effects would appear to be transferable. In summary, this new analytical model shows good potential for future ocean color inversion with low bias, well-constrained uncertainties (including the VSF), and explicit terms that can be readily tuned. Emphasis is put on application to the future NASA Plankton, Aerosol, Cloud, and ocean Ecosystem (PACE) mission.

Concentrations of Multiple Phytoplankton Pigments in the Global Oceans Obtained from Satellite Ocean Color Measurements with MERIS

The remote sensing of chlorophyll a concentration from ocean color satellites has been an essential variable quantifying phytoplankton in the past decades, yet estimation of accessory pigments from ocean color remote sensing data has remained largely elusive. In this study, we validated the concentrations of multiple pigments (Cpigs) retrieved from in situ and MEdium Resolution Imaging Spectrometer (MERIS) measured remote sensing reflectance (Rrs(λ)) in the global oceans. A multi-pigment inversion model (MuPI) was used to semi-analytically retrieve Cpigs from Rrs(λ). With a set of globally optimized parameters, the accuracy of the retrievals obtained with MuPI is quite promising. Compared with High-Performance Liquid Chromatography (HPLC) measurements near Bermuda, the concentrations of chlorophyll a, b, c ([Chl-a], [Chl-b], [Chl-c]), photoprotective carotenoids ([PPC]), and photosynthetic carotenoids ([PSC]) can be retrieved from MERIS data with a mean unbiased absolute percentage difference of 38%, 78%, 65%, 36%, and 47%, respectively. The advantage of the MuPI approach is the simultaneous retrievals of [Chl-a] and the accessory pigments [Chl-b], [Chl-c], [PPC], [PSC] from MERIS Rrs(λ) based on a closure between the input and output Rrs(λ) spectra. These results can greatly expand scientific studies of ocean biology and biogeochemistry of the global oceans that are not possible when the only available information is [Chl-a].

Analytic Network Process-Based Multi-Criteria Decision Approach and Sensitivity Analysis for Temporary Facility Layout Planning in Construction Projects

In construction projects, the planning objectives include the safety and labor productivity of the activities, along with the cost influence affecting the overall performance of construction. Temporary facilities are critical in supporting structures and equipment that impact the direct task executions and resource transportation during various construction stages. The positioning of temporary facilities and the relative spatial relation between multiple facilities are critical factors that affect the success of construction delivery. To ensure the balance among safety, labor productivity and construction cost performance, all influencing factors associated with the construction objectives should be taken into account in temporary facility layout planning. This paper proposes a novel multi-criteria temporary facility layout planning model that integrates Analytic Network Process (ANP) modeling and simulation-based sensitivity evaluation, which effectively transforms the spatial layout planning problem into a mathematical decision problem. An application example is analyzed to demonstrate its capabilities of optimizing temporary facility layout planning. The proposed framework provides the construction managers and layout planners with a useful tool for selecting the optimal alternative of temporary facility layout plan on the construction site.

Approach for identifying optically shallow pixels when processing ocean-color imagery

Ocean-color remote sensing is routinely used to derive marine geophysical parameters from sensor-observed spectral water-leaving radiances. However, in clear geometrically shallow regions, traditional ocean-color algorithms can be confounded by light reflected from the seafloor. Such regions are typically referred to as "optically shallow". When performing spatiotemporal analyses of ocean color datasets, optically shallow features such as submerged sand banks and coral reefs can lead to unexpected regional biases. Most contemporary approaches mask or flag suspected optically shallow pixels based on ancillary bathymetric data. However, the extent to which seafloor reflectance contaminates the water-leaving radiance is dependent on bathymetry, water clarity and seafloor albedo. In this paper, an approach for flagging optically shallow pixels has been developed that considers all three of these variables. In the method, the optical depth of the water column at 547 nm, ζ(547), is predicted from bathymetric data and estimated water-column optical properties. If ζ(547) is less then the pre-defined threshold, a pixel is flagged as potentially optically shallow. Radiative transfer modeling was used to identify a conservative threshold value of ζ(547) = 20 for a bright sand seafloor. In addition, pixels in waters shallower than 5 m are also flagged. We also examined how varying bathymetric datasets may affect the optically shallow flag using MODIS data. It is anticipated that the optically shallow flag will benefit end-users when quality controlling derived ocean color products. Further, the flag may prove useful as a mechanism for switching between optically deep and shallow algorithms during ocean color processing.

Evaluation of Semi-Analytical Algorithms to Retrieve Particulate and Dissolved Absorption Coefficients in Gulf of California Optically Complex Waters

Two semi-analytical algorithms, Generalized Inherent Optical Property (GIOP) and Garver-Siegel-Maritorena (GSM), were evaluated in terms of how well they reproduced the absorption coefficient of phytoplankton (aph(λ)) and dissolved and detrital organic matter (adg(λ)) at three wavelengths (λ of 412, 443, and 488 nm) in a zone with optically complex waters, the Upper Gulf of California (UGC) and the Northern Gulf of California (NGC). In the UGC, detritus determines most of the total light absorption, whereas, in the NGC, chromophoric dissolved organic material (CDOM) and phytoplankton dominate. Upon comparing the results of each model with a database assembled from four cruises done from spring to summer (March through September) between 2011 and 2013, it was found that GIOP is a better estimator for aph(λ) than GSM, independently of the region. However, both algorithms underestimate in situ values in the NGC, whereas they overestimate them in the UGC. Errors are associated with the following: (a) the constant a*ph(λ) value used by GSM and GIOP (0.055 m2 mgChla−1) is higher than the most frequent value observed in this study’s data (0.03 m2 mgChla−1), and (b) satellite-derived chlorophyll a concentration (Chla) is biased high compared with in situ Chla. GIOP gave also better results for the adg(λ) estimation than GSM, especially in the NGC. The spectral slope Sdg was identified as an important parameter for estimating adg(λ), and this study’s results indicated that the use of a fixed input value in models was not adequate. The evaluation confirms the lack of generality of algorithms like GIOP and GSM, whose reflectance model is too simplified to capture expected variability. Finally, a greater monitoring effort is suggested in the study area regarding the collection of in situ reflectance data, which would allow explaining the effects that detritus and CDOM may have on the semi-analytical reflectance inversions, as well as isolating the possible influence of the atmosphere on the satellite-derived water reflectance and Chla.

Enhance field water-color measurements with a Secchi disk and its implication for fusion of active and passive ocean-color remote sensing

Inversion of the total absorption (a) and backscattering coefficients of bulk water through a fusion of remote sensing reflectance (R_rs) and Secchi disk depth (Z_SD) is developed. An application of such a system to a synthesized wide-range dataset shows a reduction of ∼3 folds in the uncertainties of inverted a(λ) (in a range of ∼0.01-6.8  m^-1) from R_rs(λ) for the 350-560 nm range. Such a fusion is further proposed to process concurrent active (ocean LiDAR) and passive (ocean-color) measurements, which can lead to nearly "exact" analytical inversion of an R_rs spectrum. With such a fusion, it is found that the uncertainty in the inverted total a in the 350-560 nm range could be reduced to ∼2% for the synthesized data, which can thus significantly improve the derivation of a coefficients of other varying components. Although the inclusion of Z_SD places an extra constraint in the inversion of R_rs, no apparent improvement over the quasi-analytical algorithm (QAA) was found when the fusion of Z_SD and R_rs was applied to a field dataset, which calls for more accurate determination of the absorption coefficients from water samples.

Retrieval of aerosol properties and water-leaving reflectance from multi-angular polarimetric measurements over coastal waters

Ocean color remote sensing is an important tool to monitor water quality and biogeochemical conditions of ocean. Atmospheric correction, which obtains water-leaving radiance from the total radiance measured by satellite-borne or airborne sensors, remains a challenging task for coastal waters due to the complex optical properties of aerosols and ocean waters. In this paper, we report a research algorithm on aerosol and ocean color retrieval with emphasis on coastal waters, which uses coupled atmosphere and ocean radiative transfer model to fit polarized radiance measurements at multiple viewing angles and multiple wavelengths. Ocean optical properties are characterized by a generalized bio-optical model with direct accounting for the absorption and scattering of phytoplankton, colored dissolved organic matter (CDOM) and non-algal particles (NAP). Our retrieval algorithm can accurately determine the water-leaving radiance and aerosol properties for coastal waters, and may be used to improve the atmospheric correction when apply to a hyperspectral ocean color instrument.

Performance metrics for the assessment of satellite data products: an ocean color case study

Performance assessment of ocean color satellite data has generally relied on statistical metrics chosen for their common usage and the rationale for selecting certain metrics is infrequently explained. Commonly reported statistics based on mean squared errors, such as the coefficient of determination (r²), root mean square error, and regression slopes, are most appropriate for Gaussian distributions without outliers and, therefore, are often not ideal for ocean color algorithm performance assessment, which is often limited by sample availability. In contrast, metrics based on simple deviations, such as bias and mean absolute error, as well as pair-wise comparisons, often provide more robust and straightforward quantities for evaluating ocean color algorithms with non-Gaussian distributions and outliers. This study uses a SeaWiFS chlorophyll-a validation data set to demonstrate a framework for satellite data product assessment and recommends a multi-metric and user-dependent approach that can be applied within science, modeling, and resource management communities.

Vertical variations in optical properties of the waters in the Yellow Sea and Bohai Sea at seasonal scales and their influencing mechanisms

This research used the profile data measured extensively in the Yellow Sea and Bohai Sea (YSBS) to explain the temporal and spatial distribution characteristics of optical properties and systematically analyzed the influencing mechanisms of the seasonal variations of optical properties in the YSBS in conjunction with synchronously measured hydrological and biogeochemical data in vertical profiles. The main conclusions obtained are as follows: the vertical variations in the optical properties in the YSBS are mainly influenced by the stratification effect, vertical mixing function, and phytoplankton growth process; and the variations of optical properties are dominated by the change of particle characteristics. The backscattering ratio can be used to discriminate particle types as a proxy of particulate bulk refractive index. In the YSBS, for waters with a backscattering ratio of less than 0.012, the variations of optical properties are dominated by the phytoplankton particles. For waters with a backscattering ratio greater than 0.012, the variations of optical properties are dominated by inorganic sediment particles. In addition, for the YSBS, the variations in optical properties of upper surface layer waters can be elucidated well by the vertical variations.

Predicting minimum uncertainties in the inversion of ocean color geophysical parameters based on Cramer-Rao bounds

We present an analytical approach based on Cramer-Rao Bounds (CRBs) to investigate the uncertainties in estimated ocean color parameters resulting from the propagation of uncertainties in the bio-optical reflectance modeling through the inversion process. Based on given bio-optical and noise probabilistic models, CRBs can be computed efficiently for any set of ocean color parameters and any sensor configuration, directly providing the minimum estimation variance that can be possibly attained by any unbiased estimator of any targeted parameter. Here, CRBs are explicitly developed using (1) two water reflectance models corresponding to deep and shallow waters, resp., and (2) four probabilistic models describing the environmental noises observed within four Sentinel-2 MSI, HICO, Sentinel-3 OLCI and MODIS images, resp. For both deep and shallow waters, CRBs are shown to be consistent with the experimental estimation variances obtained using two published remote-sensing methods, while not requiring one to perform any inversion. CRBs are also used to investigate to what extent perfect a priori knowledge on one or several geophysical parameters can improve the estimation of remaining unknown parameters. For example, using pre-existing knowledge of bathymetry (e.g., derived from LiDAR) within the inversion is shown to greatly improve the retrieval of bottom cover for shallow waters. Finally, CRBs are shown to provide valuable information on the best estimation performances that may be achieved with the MSI, HICO, OLCI and MODIS configurations for a variety of oceanic, coastal and inland waters. CRBs are thus demonstrated to be an informative and efficient tool to characterize minimum uncertainties in inverted ocean color geophysical parameters.

High colored dissolved organic matter (CDOM) absorption in surface waters of the central-eastern Arctic Ocean: Implications for biogeochemistry and ocean color algorithms

As consequences of global warming sea-ice shrinking, permafrost thawing and changes in fresh water and terrestrial material export have already been reported in the Arctic environment. These processes impact light penetration and primary production. To reach a better understanding of the current status and to provide accurate forecasts Arctic biogeochemical and physical parameters need to be extensively monitored. In this sense, bio-optical properties are useful to be measured due to the applicability of optical instrumentation to autonomous platforms, including satellites. This study characterizes the non-water absorbers and their coupling to hydrographic conditions in the poorly sampled surface waters of the central and eastern Arctic Ocean. Over the entire sampled area colored dissolved organic matter (CDOM) dominates the light absorption in surface waters. The distribution of CDOM, phytoplankton and non-algal particles absorption reproduces the hydrographic variability in this region of the Arctic Ocean which suggests a subdivision into five major bio-optical provinces: Laptev Sea Shelf, Laptev Sea, Central Arctic/Transpolar Drift, Beaufort Gyre and Eurasian/Nansen Basin. Evaluating ocean color algorithms commonly applied in the Arctic Ocean shows that global and regionally tuned empirical algorithms provide poor chlorophyll-a (Chl-a) estimates. The semi-analytical algorithms Generalized Inherent Optical Property model (GIOP) and Garver-Siegel-Maritorena (GSM), on the other hand, provide robust estimates of Chl-a and absorption of colored matter. Applying GSM with modifications proposed for the western Arctic Ocean produced reliable information on the absorption by colored matter, and specifically by CDOM. These findings highlight that only semi-analytical ocean color algorithms are able to identify with low uncertainty the distribution of the different optical water constituents in these high CDOM absorbing waters. In addition, a clustering of the Arctic Ocean into bio-optical provinces will help to develop and then select province-specific ocean color algorithms.

Spaceborne Lidar in the Study of Marine Systems

Satellite passive ocean color instruments have provided an unbroken ∼20-year record of global ocean plankton properties, but this measurement approach has inherent limitations in terms of spatial-temporal sampling and ability to resolve vertical structure within the water column. These limitations can be addressed by coupling ocean color data with measurements from a spaceborne lidar. Airborne lidars have been used for decades to study ocean subsurface properties, but recent breakthroughs have now demonstrated that plankton properties can be measured with a satellite lidar. The satellite lidar era in oceanography has arrived. Here, we present a review of the lidar technique, its applications in marine systems, a perspective on what can be accomplished in the near future with an ocean- and atmosphere-optimized satellite lidar, and a vision for a multiplatform virtual constellation of observational assets that would enable a three-dimensional reconstruction of global ocean ecosystems.

Remote Sensing of CDOM, CDOM Spectral Slope, and Dissolved Organic Carbon in the Global Ocean

A Global Ocean Carbon Algorithm Database (GOCAD) has been developed from over 500 oceanographic field campaigns conducted worldwide over the past 30 years including in situ reflectances and coincident satellite imagery, multi- and hyperspectral Chromophoric Dissolved Organic Matter (CDOM) absorption coefficients from 245–715 nm, CDOM spectral slopes in eight visible and ultraviolet wavebands, dissolved and particulate organic carbon (DOC and POC, respectively), and inherent optical, physical, and biogeochemical properties. From field optical and radiometric data and satellite measurements, several semi-analytical, empirical, and machine learning algorithms for retrieving global DOC, CDOM, and CDOM slope were developed, optimized for global retrieval, and validated. Global climatologies of satellite-retrieved CDOM absorption coefficient and spectral slope based on the most robust of these algorithms lag seasonal patterns of phytoplankton biomass belying Case 1 assumptions, and track terrestrial runoff on ocean basin scales. Variability in satellite retrievals of CDOM absorption and spectral slope anomalies are tightly coupled to changes in atmospheric and oceanographic conditions associated with El Niño Southern Oscillation (ENSO), strongly covary with the multivariate ENSO index in a large region of the tropical Pacific, and provide insights into the potential evolution and feedbacks related to sea surface dissolved carbon in a warming climate. Further validation of the DOC algorithm developed here is warranted to better characterize its limitations, particularly in mid-ocean gyres and the southern oceans.

An overview of approaches and challenges for retrieving marine inherent optical properties from ocean color remote sensing

Ocean color measured from satellites provides daily global, synoptic views of spectral waterleaving reflectances that can be used to generate estimates of marine inherent optical properties (IOPs). These reflectances, namely the ratio of spectral upwelled radiances to spectral downwelled irradiances, describe the light exiting a water mass that defines its color. IOPs are the spectral absorption and scattering characteristics of ocean water and its dissolved and particulate constituents. Because of their dependence on the concentration and composition of marine constituents, IOPs can be used to describe the contents of the upper ocean mixed layer. This information is critical to further our scientific understanding of biogeochemical oceanic processes, such as organic carbon production and export, phytoplankton dynamics, and responses to climatic disturbances. Given their importance, the international ocean color community has invested significant effort in improving the quality of satellite-derived IOP products, both regionally and globally. Recognizing the current influx of data products into the community and the need to improve current algorithms in anticipation of new satellite instruments (e.g., the global, hyperspectral spectroradiometer of the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission), we present a synopsis of the current state of the art in the retrieval of these core optical properties. Contemporary approaches for obtaining IOPs from satellite ocean color are reviewed and, for clarity, separated based their inversion methodology or the type of IOPs sought. Summaries of known uncertainties associated with each approach are provided, as well as common performance metrics used to evaluate them. We discuss current knowledge gaps and make recommendations for future investment for upcoming missions whose instrument characteristics diverge sufficiently from heritage and existing sensors to warrant reassessing current approaches.

Assessing uncertainties in scattering correction algorithms for reflective tube absorption measurements made with a WET Labs ac-9

In situ absorption measurements collected with a WET Labs ac-9 employing a reflective tube approach were scatter corrected using several possible methods and compared to reference measurements made by a PSICAM to assess performance. Overall, two correction methods performed best for the stations sampled: one using an empirical relationship between the ac-9 and PSICAM to derive the scattering error (ε) in the near-infrared (NIR), and one where ε was independently derived from concurrent measurements of the volume scattering function (VSF). Application of the VSF-based method may be more universally applicable, although difficult to routinely apply because of the lack of commercially available VSF instrumentation. The performance of the empirical approach is encouraging as it relies only on the ac meter measurement and may be readily applied to historical data, although there are inevitably some inherent assumptions about particle composition that hinder universal applicability. For even the best performing methods, residual errors of 20% or more were commonly observed for many water types. For clear ocean waters, a conventional baseline subtraction with the assumption of negligible near-IR absorption performed as well or better than the above methods because propagated uncertainties were lower than observed with the proportional method.

Simplified model of spectral absorption by non-algal particles and dissolved organic materials in aquatic environments

Absorption by non-algal particles (NAP, ad) and colored dissolved organic matter (CDOM, ag) are frequently modeled by exponential functions of wavelength, either separately or as a sum. We present a new representation of NAP-plus-CDOM absorption adg based on the stretched exponential function adg(λ) = A exp{-[s(λ - λo)]β}, whose parameter β can be considered a measure of optical heterogeneity. A double exponential representation of adg can be fit extremely well by a stretched exponential for all plausible parameter combinations, despite having one fewer free parameter than a double exponential. Fitting two published compilations of in situ adg data - one at low spectral resolution (n = 5, λ = 412-555 nm) and one at high spectral resolution (n = 201, λ = 300-700 nm) - the stretched exponential outperforms the single exponential, double exponential, and a power law. We thereby conclude that the stretched exponential is the preferred model for adg absorption in circumstances when NAP and CDOM cannot be separated, such as in remote sensing inversions.

Self-shading associated with a skylight-blocked approach system for the measurement of water-leaving radiance and its correction

Self-shading associated with a skylight-blocked approach (SBA) system for the measurement of water-leaving radiance (L_w) and its correction [Appl. Opt.52, 1693 (2013)APOPAI0003-693510.1364/AO.52.001693] is characterized by Monte Carlo simulations, and it is found that this error is in a range of ∼1%-20% under most water properties and solar positions. A model for estimating this shading error is further developed, and eventually a scheme to correct this error based on the shaded measurements is proposed and evaluated. It is found that the shade-corrected value in the visible domain is within 3% of the true value, which thus indicates that we can obtain not only high precision but also high accuracy L_w in the field with the SBA scheme.

From toes to top-of-atmosphere: Fowler's Sneaker Depth index of water clarity for the Chesapeake Bay

Fowler's Sneaker Depth (FSD), analogous to the well known Secchi disk depth (Zsd), is a visually discerned citizen scientist metric used to assess water clarity in the Patuxent River estuary. In this study, a simple remote sensing algorithm was developed to derive FSD from space-borne spectroradiometric imagery. An empirical model was formed that estimates FSD from red-end remote sensing reflectances at 645 nm, Rrs(645). The model is based on a hyperbolic function relating water clarity to Rrs(645) that was established using radiative transfer modeling and fine tuned using in-water FSD measurements and coincident Rrs(645) data observed by NASA's Moderate Resolution Imaging Spectroradiometer aboard the Aqua spacecraft (MODISA). The resultant FSD algorithm was applied to Landsat-8 Operational Land Imager data to derive a short time-series for the Patuxent River estuary from January 2015 to June 2016. Satellite-derived FSD had an inverse, statistically significant relationship (p<0.005) with total suspended sediment concentration (TSS). Further, a distinct negative relationship between FSD and chlorophyll concentration was discerned during periods of high biomass (> 4 μg L^-1). The complex nature of water quality in the mid-to-upper Chesapeake Bay was captured using a MODISA-based FSD time series (2002-2016). This study demonstrates how a citizen scientist-conceived observation can be coupled with remote sensing. With further refinement and validation, the FSD may be a useful tool for delivering scientifically relevant results and for informing and engaging local stakeholders and policy makers.

Developments in Earth observation for the assessment and monitoring of inland, transitional, coastal and shelf-sea waters

The Earth's surface waters are a fundamental resource and encompass a broad range of ecosystems that are core to global biogeochemical cycling and food and energy production. Despite this, the Earth's surface waters are impacted by multiple natural and anthropogenic pressures and drivers of environmental change. The complex interaction between physical, chemical and biological processes in surface waters poses significant challenges for in situ monitoring and assessment and often limits our ability to adequately capture the dynamics of aquatic systems and our understanding of their status, functioning and response to pressures. Here we explore the opportunities that Earth observation (EO) has to offer to basin-scale monitoring of water quality over the surface water continuum comprising inland, transition and coastal water bodies, with a particular focus on the Danube and Black Sea region. This review summarises the technological advances in EO and the opportunities that the next generation satellites offer for water quality monitoring. We provide an overview of algorithms for the retrieval of water quality parameters and demonstrate how such models have been used for the assessment and monitoring of inland, transitional, coastal and shelf-sea systems. Further, we argue that very few studies have investigated the connectivity between these systems especially in large river-sea systems such as the Danube-Black Sea. Subsequently, we describe current capability in operational processing of archive and near real-time satellite data. We conclude that while the operational use of satellites for the assessment and monitoring of surface waters is still developing for inland and coastal waters and more work is required on the development and validation of remote sensing algorithms for these optically complex waters, the potential that these data streams offer for developing an improved, potentially paradigm-shifting understanding of physical and biogeochemical processes across large scale river-sea systems including the Danube-Black Sea is considerable.

Absorption correction and phase function shape effects on the closure of apparent optical properties

We present a closure experiment between new inherent optical properties (IOPs: absorption a, scattering b, backscattering b_b) and apparent optical properties (AOPs: remote-sensing reflectance R_rs, irradiance reflectance R, and anisotropic factor at nadir Q_n) data of Ionian and Adriatic seawaters, from very clear to turbid waters, ranging across one order of magnitude in R_rs. The internal consistency of the IOP-AOP matchups was investigated though radiative transfer closure. Using the in situ IOPs, we predicted the AOPs with the commercial radiative transfer solver Hydrolight. Closure was limited by two unresolved issues, one regarding processing of in situ data and the other related to radiative transfer modeling. First, different correction methods of the absorption data measured by the Wetlabs ac-s produced high variations in simulated reflectances, reaching 40% for the highest reflectances in our dataset. Second, the lack of detailed volume scattering function measurements forces us to adopt analytical functions that are consistent with a given particle backscattering ratio. The analytical phase functions named Fournier-Forand and two-term Kopelevich presented reasonable angular shapes with respect to measurements at a few backward angles. Between these phase functions, induced changes were within 4% for R_rs, within 11% for R, and within 10% for Q_n. Additionally, closure of Q_n was generally not successful considering radiometric uncertainties. Simulated Q_n overestimated low values and underestimated high values, especially at 665 nm, where Hydrolight was unable to predict measured Q_n values greater than 6 sr. The physical nature of Q_n makes this mismatch almost independent of the measured IOPs, thus precluding Q_n tuning by varying the former. The non-closure of Q_n might be caused by an inaccurate phase function and, to a lesser extent, by the modeling of the incoming radiance. For the future, this remains the task of accurate absorption and phase function measurements, especially at red wavelengths.

Effects of inelastic radiative processes on the determination of water-leaving spectral radiance from extrapolation of underwater near-surface measurements

Extrapolation of near-surface underwater measurements is the most common method to estimate the water-leaving spectral radiance, L_w(λ) (where λ is the light wavelength in vacuum), and remote-sensing reflectance, R_rs(λ), for validation and vicarious calibration of satellite sensors, as well as for ocean color algorithm development. However, uncertainties in L_w(λ) arising from the extrapolation process have not been investigated in detail with regards to the potential influence of inelastic radiative processes, such as Raman scattering by water molecules and fluorescence by colored dissolved organic matter and chlorophyll-a. Using radiative transfer simulations, we examine high-depth resolution vertical profiles of the upwelling radiance, L_u(λ), and its diffuse attenuation coefficient, K_Lu (λ), within the top 10 m of the ocean surface layer and assess the uncertainties in extrapolated values of L_w(λ). The inelastic processes generally increase L_u and decrease K_Lu in the red and near-infrared (NIR) portion of the spectrum. Unlike K_Lu in the blue and green spectral bands, K_Lu in the red and NIR is strongly variable within the near-surface layer even in a perfectly homogeneous water column. The assumption of a constant K_Lu with depth that is typically employed in the extrapolation method can lead to significant errors in the estimate of L_w. These errors approach ∼100% at 900 nm, and the desired threshold of 5% accuracy or less cannot be achieved at wavelengths greater than 650 nm for underwater radiometric systems that typically take measurements at depths below 1 m. These errors can be reduced by measuring L_u within a much shallower surface layer of tens of centimeters thick or even less at near-infrared wavelengths longer than 800 nm, which suggests a requirement for developing appropriate radiometric instrumentation and deployment strategies.

Implementation of an analytical Raman scattering correction for satellite ocean-color processing

Raman scattering of photons by seawater molecules is an inelastic scattering process. This effect can contribute significantly to the water-leaving radiance signal observed by space-borne ocean-color spectroradiometers. If not accounted for during ocean-color processing, Raman scattering can cause biases in derived inherent optical properties (IOPs). Here we describe a Raman scattering correction (RSC) algorithm that has been integrated within NASA's standard ocean-color processing software. We tested the RSC with NASA's Generalized Inherent Optical Properties algorithm (GIOP). A comparison between derived IOPs and in situ data revealed that the magnitude of the derived backscattering coefficient and the phytoplankton absorption coefficient were reduced when the RSC was applied, whilst the absorption coefficient of colored dissolved and detrital matter remained unchanged. Importantly, our results show that the RSC did not degrade the retrieval skill of the GIOP. In addition, a time-series study of oligotrophic waters near Bermuda showed that the RSC did not introduce unwanted temporal trends or artifacts into derived IOPs.

Multi-method approach to quantify uncertainties in the measurements of light absorption by particles

Through technological and research advances, numerous methods and protocols have emerged to estimate spectral absorption of light by particles, a_p, in an aquatic medium. However, the level of agreement among measurements remains elusive. We employed a multi-method approach to estimate the measurement precision of measuring optical density of particles on a filter pad using two common spectrophotometric methods, and the determination precision, or uncertainty, of the computational techniques for estimating a_p for six ocean color wavelengths (412, 443, 490, 510, 555, 670 nm). The optical densities measured with the two methods exhibited a significant, positive correlation. Optical density measurement precision ranged from 0.061%-63% and exhibited a significant, positive correlation. Multi-method uncertainty ranged from 7.48%-119%. Values of a_p at 555 nm and 670 nm exhibited the highest values of uncertainty. Poor performance of modeled a_p compared to determined a_p suggest uncertainties are propagated into bio-optical algorithms.

An assessment of phytoplankton primary productivity in the Arctic Ocean from satellite ocean color/in situ chlorophyll-a based models

We investigated 32 net primary productivity (NPP) models by assessing skills to reproduce integrated NPP in the Arctic Ocean. The models were provided with two sources each of surface chlorophyll-a concentration (chlorophyll), photosynthetically available radiation (PAR), sea surface temperature (SST), and mixed-layer depth (MLD). The models were most sensitive to uncertainties in surface chlorophyll, generally performing better with in situ chlorophyll than with satellite-derived values. They were much less sensitive to uncertainties in PAR, SST, and MLD, possibly due to relatively narrow ranges of input data and/or relatively little difference between input data sources. Regardless of type or complexity, most of the models were not able to fully reproduce the variability of in situ NPP, whereas some of them exhibited almost no bias (i.e., reproduced the mean of in situ NPP). The models performed relatively well in low-productivity seasons as well as in sea ice-covered/deep-water regions. Depth-resolved models correlated more with in situ NPP than other model types, but had a greater tendency to overestimate mean NPP whereas absorption-based models exhibited the lowest bias associated with weaker correlation. The models performed better when a subsurface chlorophyll-a maximum (SCM) was absent. As a group, the models overestimated mean NPP, however this was partly offset by some models underestimating NPP when a SCM was present. Our study suggests that NPP models need to be carefully tuned for the Arctic Ocean because most of the models performing relatively well were those that used Arctic-relevant parameters.