<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mar. Sci.</journal-id>
<journal-title>Frontiers in Marine Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mar. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-7745</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmars.2023.1267681</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Marine Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>What data are needed to detect wildfire effects on coastal ecosystems? A case study during the Thomas Fire</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Kramer</surname>
<given-names>Sasha J.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1973107"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Bisson</surname>
<given-names>Kelsey M.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/797724"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mitchell</surname>
<given-names>Catherine</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/317444"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Monterey Bay Aquarium Research Institute</institution>, <addr-line>Moss Landing, CA</addr-line>, <country>United States</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Botany &amp; Plant Pathology, Oregon State University</institution>, <addr-line>Corvallis, OR</addr-line>, <country>United States</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Bigelow Laboratory for Ocean Sciences</institution>, <addr-line>East Boothbay, ME</addr-line>, <country>United States</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Katherine Dafforn, Macquarie University, Australia</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Xuerong Sun, University of Exeter, United Kingdom; Joshua Picotte, United States Geological Survey (USGS), United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Sasha J. Kramer, <email xlink:href="mailto:skramer@mbari.org">skramer@mbari.org</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>10</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>10</volume>
<elocation-id>1267681</elocation-id>
<history>
<date date-type="received">
<day>26</day>
<month>07</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>03</day>
<month>10</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Kramer, Bisson and Mitchell</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Kramer, Bisson and Mitchell</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Wildfires are growing in frequency and severity worldwide, with anthropogenic climate change predicted to worsen the effects of wildfires in the future. While most wildfire impacts occur on land, coastal fires can also affect the ocean via smoke production and ash deposition. The impacts of wildfires on marine ecology and biogeochemistry have been studied infrequently, as it is difficult to conduct fieldwork rapidly and safely during unpredictable natural disasters. Increasingly, remote sensing measurements are used to study the impacts of wildfires on marine ecosystems through optical proxies. Given the optical impacts of smoke and in-water ash, these measurements may be limited in their scope and accuracy. Here, we evaluate the potential and limitations of remote sensing data collected from MODIS-Aqua to describe the effects of wildfires on optics and phytoplankton observations. Using samples collected in the Santa Barbara Channel (California, USA) during the Thomas Fire in December 2017, we found that MODIS-Aqua data were unsuited for interpreting ecosystem effects during a wildfire. Our results identified a persistent overestimation of chlorophyll-<italic>a</italic> concentration from MODIS-Aqua compared to <italic>in situ</italic> measurements. Optical models applied to <italic>in situ</italic> radiometry data overestimated the absorption by colored dissolved organic matter (CDOM) during the wildfire. Satellites will remain an important tool to measure the impacts of wildfires on marine ecosystems, but this analysis demonstrates the importance of <italic>in situ</italic> sampling to quantify the impacts of wildfires on ocean ecology and biogeochemistry due to the difficulty of interpreting remote sensing data during these events.</p>
</abstract>
<kwd-group>
<kwd>phytoplankton</kwd>
<kwd>wildfires</kwd>
<kwd>remote sensing</kwd>
<kwd>ocean color</kwd>
<kwd>phytoplankton absorption</kwd>
<kwd>CDOM absorption</kwd>
<kwd>Santa Barbara Channel</kwd>
</kwd-group>
<contract-num rid="cn001">NNX17AK04G, 80NSSC21K0015 , 80NSSC22K0049</contract-num>
<contract-num rid="cn002">Postdoctoral Fellowship in Marine Microbial Ecology</contract-num>
<contract-num rid="cn003">1821916</contract-num>
<contract-num rid="cn004">2101A</contract-num>
<contract-sponsor id="cn001">National Aeronautics and Space Administration<named-content content-type="fundref-id">10.13039/100000104</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Simons Foundation<named-content content-type="fundref-id">10.13039/100000893</named-content>
</contract-sponsor>
<contract-sponsor id="cn003">National Science Foundation<named-content content-type="fundref-id">10.13039/100000001</named-content>
</contract-sponsor>
<contract-sponsor id="cn004">North Pacific Research Board<named-content content-type="fundref-id">10.13039/100012635</named-content>
</contract-sponsor>
<contract-sponsor id="cn005">David and Lucile Packard Foundation<named-content content-type="fundref-id">10.13039/100000008</named-content>
</contract-sponsor>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="7"/>
<ref-count count="50"/>
<page-count count="12"/>
<word-count count="6760"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Ocean Observation</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Around the world, wildfires are increasing in magnitude and severity (<xref ref-type="bibr" rid="B15">Doerr and Sant&#xed;n, 2016</xref>; <xref ref-type="bibr" rid="B18">Hamilton et&#xa0;al., 2018</xref>). The conditions that lead to more intense and destructive wildfires are also increasing, in part due to the effects of anthropogenic climate change on prolonged drought conditions and frequent severe weather events (<xref ref-type="bibr" rid="B50">Westerling et&#xa0;al., 2006</xref>). These impacts are particularly true in the western United States, where wildfire impacts on inland waters and watersheds have been studied (e.g., <xref ref-type="bibr" rid="B12">Coombs and Melack, 2013</xref>; <xref ref-type="bibr" rid="B39">Urmy et&#xa0;al., 2016</xref>). However, while wildfires are frequently proximal to the ocean, there have been few studies examining the impact of wildfires on the ecology and biogeochemistry of marine systems until recently. One example hypothesized about the impact of wildfire smoke on coral reef productivity in the near-coastal zone (e.g., <xref ref-type="bibr" rid="B1">Abram et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B42">van Woesik, 2004</xref>). Additional inferences about the influence of wildfire smoke and ash on the ocean have also been drawn from the parallel example of volcanic ash and smoke. Volcanic impacts have measurable effects on the ocean, including fertilizing the nutrient-limited open ocean (<xref ref-type="bibr" rid="B16">Duggen et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B19">Hamme et&#xa0;al., 2010</xref>).</p>
<p>Like volcanic eruptions, wildfires are largely unpredictable events. While some regions may be more susceptible to wildfires in certain seasons, the timing, duration, and extent of a potential fire is largely unknown prior to that event. Consequently, coastal studies are easier to conduct than open ocean approaches and can take advantage of existing nearshore time series to compare wildfire impacts to historical norms. Given the planning required to conduct a successful oceanographic research cruise, virtually no <italic>in situ</italic> studies existed to compare wildfire effects on ocean ecosystems prior to the Thomas Fire, which occurred during a pre-existing cruise (<xref ref-type="bibr" rid="B6">Bisson et&#xa0;al., 2020</xref>) and allowed for the evaluation of phytoplankton community composition during the fire (<xref ref-type="bibr" rid="B23">Kramer et&#xa0;al., 2020</xref>). Since then, other studies in the coastal ocean have examined the impacts of wildfires on dissolved black carbon (<xref ref-type="bibr" rid="B43">Wagner et&#xa0;al., 2021</xref>), on the concentrations of fecal indicator bacteria (<xref ref-type="bibr" rid="B10">Cira et&#xa0;al., 2022</xref>), and on dissolved organic matter deposition (<xref ref-type="bibr" rid="B13">Coward et&#xa0;al., 2022</xref>). These studies relied on the availability of ships to access the sampling sites during wildfires and collect samples on timescales relevant to the disaster, which can provide a further logistical and funding challenge.</p>
<p>Increasingly, satellite remote sensing provides a (near) real-time tool to monitor the impacts of these sporadic wildfire events on the open ocean, as has been done with volcanic impacts on the ocean (<xref ref-type="bibr" rid="B49">Westberry et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B8">Bisson et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B17">Franz et al., 2023</xref>). Both ocean color remote sensing and autonomous profiling floats (e.g., Argo) have been used as tools to identify phytoplankton blooms following a wildfire, particularly in inaccessible, open ocean and/or high-latitude regions (<xref ref-type="bibr" rid="B3">Ardyna et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B44">Wang et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B45">Weis et&#xa0;al., 2022</xref>). All three of these studies examined the changes in chlorophyll-<italic>a</italic> concentration and phytoplankton carbon indirectly through optical proxies and determined that wildfire smoke and ash is potentially fertilizing the Southern and Arctic Oceans, similarly to the effects observed after volcanic eruptions. Furthermore, these studies used business-as-usual satellite products and assumed those aerosol models and atmospheric corrections were appropriate, which may not be true during a wildfire. Given the unpredictable nature of wildfires on land and the potentially rapid response by phytoplankton in the coastal and open ocean, this area of study urgently requires further inquiry. However, the remote sensing component must be carefully considered to avoid drawing improper conclusions from ocean color imagery impacted by insufficient atmospheric correction or unexpected organic matter absorption and scattering.</p>
<p>It is highly likely that the proximal wildfire smoke and <italic>in situ</italic> ash affected remote sensing measurements and/or optical proxies employed in these previous studies, because smoke and ash obscure the ocean signal and provide challenges to conventional atmospheric correction. The majority (&gt;90%) of the top-of-atmosphere radiance signal observed from passive satellites comes from atmospheric contributions. Even routine errors in atmospheric correction in the absence of smoke can create large downstream errors in ocean color products, and it has been demonstrated that errors in remote sensing reflectance (<inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) are up to an order of magnitude larger than was previously understood in some regions (<xref ref-type="bibr" rid="B7">Bisson et&#xa0;al., 2021</xref>). Additionally, since <inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is proportional to the ratio of backscattering (<inline-formula>
<mml:math display="inline" id="im3">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) to absorption (<inline-formula>
<mml:math display="inline" id="im4">
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>), any changes to the inherent optical properties of the ocean water due to ash deposited into the ocean will further impact the ocean color signal (Equations 1-3):</p>
<disp-formula>
<label>[1],</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x221d;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where</p> <disp-formula>
<label>[2]</label>
<mml:math display="block" id="im5">
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p><italic>and</italic></p>
<disp-formula><label>[3].</label>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>
<inline-formula>
<mml:math display="inline" id="im6">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> represents absorption by seawater and <inline-formula>
<mml:math display="inline" id="im7">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is backscattering by seawater. The particulate absorption component, <inline-formula>
<mml:math display="inline" id="im8">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, is composed of absorption by phytoplankton (<inline-formula>
<mml:math display="inline" id="im9">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) and absorption by non-algal particles (also referred to as depigmented particles) (<inline-formula>
<mml:math display="inline" id="im10">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>). Absorption by colored dissolved organic material (CDOM), <inline-formula>
<mml:math display="inline" id="im11">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, also contributes to the total absorption signal. Particulate backscattering, <inline-formula>
<mml:math display="inline" id="im12">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, includes backscattering by phytoplankton (<inline-formula>
<mml:math display="inline" id="im13">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) and depigmented particles (<inline-formula>
<mml:math display="inline" id="im14">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>). Deposited wildfire ash contributes non-algal particulates and CDOM with absorbing and scattering properties that will impact the water-leaving radiance signal in the same spectral range as absorption and scattering by phytoplankton (as has been seen with volcanic ash; e.g., <xref ref-type="bibr" rid="B9">Browning et&#xa0;al., 2015</xref>). Through these combined effects of atmospheric interference and surface ocean deposition, wildfire smoke and ash might impact the conclusions that can be drawn from ocean color data.</p>
<p>Here, we use the Thomas Fire as a case study to explore the data that are needed to detect optical changes during a coastal wildfire. When the Thomas Fire ignited in December 2017, our team was prepared to take advantage of a serendipitously planned research cruise in the Santa Barbara Channel that month. We were also able to contextualize the observations made during the fire through the historical optical data collected as part of the Plumes and Blooms timeseries (PnB), with high-resolution optical data collected near-monthly since 1996. Due to the aforementioned difficulties in collecting <italic>in situ</italic> data during a wildfire and placing that data in context, our study represents a rare opportunity to consider the most useful data for observing optical and biogeochemical changes during a coastal wildfire. We leverage these data with available satellite matchups from the MODIS-Aqua sensor to examine how wildfire ash affected the Santa Barbara Channel, and to consider the most relevant, accurate, and impactful data during this extreme event.</p>
<p>At the time of sampling in 2017, the Thomas Fire was the largest wildfire in California history&#x2014;it is now the eighth largest fire wildfire and the fifteenth most destructive (CAL Fire; <ext-link ext-link-type="uri" xlink:href="https://www.fire.ca.gov/our-impact/statistics">https://www.fire.ca.gov/our-impact/statistics</ext-link>). It is likely that coastal wildfires will continue to grow in intensity, frequency, and magnitude, and remote sensing measurements will remain an important, practical, and cost-effective tool to monitor the varied and poorly-constrained impacts of these extreme events on the ocean. However, satellite-based approaches will need to consider the potential impacts of the fire on the bio-optical regime of the water column, particularly through increasing CDOM and non-algal particles alongside potential increases in phytoplankton biomass. The results presented here suggest that satellite data do not reliably capture the complexity of ocean color changes during coastal wildfires, and that <italic>in situ</italic> data are required to describe both optical variability and ecological impacts of wildfire smoke and ash.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Sampling during the Thomas Fire</title>
<p>All <italic>in situ</italic> samples were collected in the Santa Barbara Channel, California, USA as part of the Across the Channel Investigating Diel Dynamics (ACIDD) cruise from December 16-23, 2017 (<xref ref-type="bibr" rid="B6">Bisson et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B23">Kramer et&#xa0;al., 2020</xref>). Samples were collected at local noon, plus or minus one hour, for correspondence with MODIS-Aqua overpasses. Due to the difficulty of obtaining coincident, unobstructed ocean color data (i.e., free from clouds or thick smoke) during the wildfire, match-up samples between <italic>in situ</italic> and remote sensing were only available from three days of sampling: December 19-21, 2017. The location of <italic>in situ</italic> samples overlaid with MODIS ocean color data is shown in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>MODIS chlorophyll-<italic>a</italic> concentration for December 19, December 20, and December 21, 2017. The location of <italic>in situ</italic> sampling for each day is shown as a black cross. The location of Plumes and Blooms Station 4 is shown in a white octagon. White patches on December 19 and 20 were masked out and flagged due to cloud cover, high polarization, or other product failures.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1267681-g001.tif"/>
</fig>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>
<italic>In situ</italic> measurements</title>
<sec id="s2_2_1">
<label>2.2.1</label>
<title>Compact-Optical Profiling System (C-OPS) reflectance</title>
<p>The Compact-Optical Profiling System (C-OPS, Biospherical Instruments, Inc.) was used to measure upwelling spectral radiance (<inline-formula>
<mml:math display="inline" id="im15">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) and downwelling spectral irradiance (<inline-formula>
<mml:math display="inline" id="im16">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) from 320 and 780 nm at 18 distinct wavelengths. Instrument deployment and data processing for the C-OPS in free-fall mode followed methods presented in <xref ref-type="bibr" rid="B38">Toole and Siegel (2001)</xref>; <xref ref-type="bibr" rid="B22">Kostadinov et&#xa0;al. (2007)</xref>; <xref ref-type="bibr" rid="B5">Barr&#xf3;n et&#xa0;al. (2014)</xref>, and <xref ref-type="bibr" rid="B20">Henderikx Freitas et&#xa0;al. (2017)</xref>. The instrument was deployed at noon from the stern of the ship to minimize the influence of ship shadow or solar zenith angle on the measurements. Remote sensing reflectance just below the surface, <inline-formula>
<mml:math display="inline" id="im17">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, was obtained following equation 4 (<xref ref-type="bibr" rid="B28">Mobley, 1994</xref>):</p>
<disp-formula>
<label>[4]</label>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac><mml:mo>.</mml:mo></mml:mrow>
</mml:math>
</disp-formula>
<p>
<inline-formula>
<mml:math display="inline" id="im18">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> was then converted to remote sensing reflectance just above the surface of the water, <inline-formula>
<mml:math display="inline" id="im19">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, following <xref ref-type="bibr" rid="B25">Lee et&#xa0;al. (2002)</xref>:</p>
<disp-formula>
<label>[5]</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>0.52</mml:mn>
<mml:mo>*</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1.7</mml:mn>
<mml:mo>*</mml:mo>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac><mml:mo>.</mml:mo></mml:mrow>
</mml:math>
</disp-formula>
<p>C-OPS samples were collected from December 17-21, 2017 (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;1A</bold>
</xref>). Replicate profiles were obtained to produce uncertainty estimates (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Remote sensing reflectance measured <bold>(A)</bold> <italic>in situ</italic> with the C-OPS and <bold>(B)</bold> from space with MODIS. Spectra were measured on December 19 (red), December 20 (green), and December 21 (blue), 2017. Error bars show +/- 1 standard deviation: for C-OPS, within the top 1m of water; for MODIS-Aqua, within the 5x5 pixel box surrounding the C-OPS cast site.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1267681-g002.tif"/>
</fig>
</sec>
<sec id="s2_2_2">
<label>2.2.2</label>
<title>Filter pad absorption by phytoplankton (<inline-formula>
<mml:math display="inline" id="im21">
<mml:mrow>
<mml:msub>
<mml:mi>aph</mml:mi>
<mml:mtext>&#xa0;</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>(<italic>&#x3bb;</italic>)) and non-algal particles (<inline-formula>
<mml:math display="inline" id="im103">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>(<italic>&#x3bb;</italic>))</title>
<p>Two-liter whole seawater samples (N = 24) were collected from the surface (&lt;5&#xa0;m) Niskin bottle from the CTD rosette and filtered through 0.7&#xb5;m 25&#xa0;mm GF/F filters (Whatman), which were then folded in half and stored in pre-labeled aluminum foil packets. Filters were stored immediately in liquid nitrogen for 1 week and at -80&#xb0;C for 5 months until analysis. Three blank samples were prepared at sea by filtering 2L MilliQ water through the same filter type at the same pump pressure and were stored in the same way as the discrete samples until the time of analysis. Just before analysis, one laboratory blank was prepared by filtering 2L MilliQ water. This sample was used as the reference blank in the spectrophotometer.</p>
<p>A Shimadzu UV-2401PB dual-beam spectrophotometer with a 60&#xa0;mm diameter ISR-2200 integrating sphere was used to measure optical densities of particulate material on the filter, <inline-formula>
<mml:math display="inline" id="im22">
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>. Samples were run in transmission-reflectance (T-R) mode. Frozen samples were thawed prior to analysis. Filters were kept hydrated and protected from light when not in the spectrophotometer. Samples were placed on the integrating sphere and run against the laboratory blank reference filter. All filters were scanned first for total particulate absorption (<inline-formula>
<mml:math display="inline" id="im23">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>). The filters were then extracted in methanol (following <xref ref-type="bibr" rid="B21">Kishino et&#xa0;al., 1985</xref>): 10ml hot methanol was poured carefully on each filter and allowed to sit for a few minutes before filtering through. The hot methanol rinse was repeated, and each filter was then rinsed with ~20ml boiling MilliQ water to remove phycobiliproteins, which do not extract in methanol (<xref ref-type="bibr" rid="B33">Roesler et&#xa0;al., 2018</xref>). Filters were then re-scanned for non-algal (depigmented) particle absorption (<inline-formula>
<mml:math display="inline" id="im24">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>). All samples and blanks were scanned 3 times (with a 90&#xb0; rotation between each scan) to obtain an average scan for all measured blanks and samples. The averaged blank filter absorbance was subtracted from the averaged absorbance of the corresponding sample filter. All spectra were zeroed at 850 nm as needed, where no absorbance is expected, following the beta-correction approach of <xref ref-type="bibr" rid="B31">Roesler (1998)</xref>.</p>
<p>Absorbance was converted to absorption and all spectra were corrected for pathlength amplification following the recommendations of the IOCCG protocols (<xref ref-type="bibr" rid="B33">Roesler et&#xa0;al., 2018</xref>) based on <xref ref-type="bibr" rid="B37">Tassan and Ferrari (1995)</xref> and <xref ref-type="bibr" rid="B36">Stramski et&#xa0;al. (2015)</xref> for a spectrophotometer in transmittance-reflectance (T-R) mode with an integrating sphere and externally-mounted samples, where the optical density of particulate material on the filter (<inline-formula>
<mml:math display="inline" id="im25">
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) is corrected to absorbance, <inline-formula>
<mml:math display="inline" id="im26">
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, following equation 6:</p>
<disp-formula>
<label>[6]</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mn>0.719</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mn>1.2287</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Absorbance (<inline-formula>
<mml:math display="inline" id="im27">
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) is then converted to absorption (<inline-formula>
<mml:math display="inline" id="im28">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) following equation 7:</p>
<disp-formula>
<label>[7]</label>
<mml:math display="block" id="M6">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2.303</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>*</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where the pathlength is the volume filtered (m<sup>3</sup>) divided by the area of the filter (m<sup>2</sup>), which was measured for each sample. This correction was applied to all particulate and depigmented particle absorption spectra. Finally, spectral absorption by phytoplankton (<inline-formula>
<mml:math display="inline" id="im29">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) is computed as the difference between total particulate absorption (<inline-formula>
<mml:math display="inline" id="im30">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) and depigmented particulate (<inline-formula>
<mml:math display="inline" id="im31">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) absorption coefficients.</p>
</sec>
<sec id="s2_2_3">
<label>2.2.3</label>
<title>Colored dissolved organic matter (CDOM) absorption (<inline-formula>
<mml:math display="inline" id="im32">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>(<italic>&#x3bb;</italic>))</title>
<p>CDOM absorption (<inline-formula>
<mml:math display="inline" id="im33">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) was measured following <xref ref-type="bibr" rid="B29">Nelson et&#xa0;al. (2007)</xref>. Whole seawater samples (N = 21) were collected from the surface Niskin bottle in 100ml amber glass vials. Samples were filtered immediately at sea using pre-extracted 0.2&#xb5;m polycarbonate filters (Nuclepore). The filtrate was collected in amber glass EPA vials and kept in the dark at 4&#xb0;C until analysis. Samples were equilibrated to room temperature before analysis. CDOM absorption (<inline-formula>
<mml:math display="inline" id="im34">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) was measured using an UltraPath single beam spectrophotometer (World Precision Instruments). Optical density (<inline-formula>
<mml:math display="inline" id="im35">
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) was measured and adjusted to <inline-formula>
<mml:math display="inline" id="im36">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> following equation 7, where the pathlength was 0.1m. Salinity corrections were applied as in <xref ref-type="bibr" rid="B29">Nelson et&#xa0;al. (2007)</xref>.</p>
</sec>
<sec id="s2_2_4">
<label>2.2.4</label>
<title>High Performance Liquid Chromatography (HPLC) pigments</title>
<p>Two-liter whole seawater samples were collected from the surface Niskin bottle. Samples were filtered through 0.7&#xb5;m 25&#xa0;mm Whatman GF/F filters and stored immediately in liquid nitrogen until analysis at NASA&#x2019;s Goddard Space Flight Center (GSFC) by the Ocean Biology Processing Group following <xref ref-type="bibr" rid="B40">Van Heukelem and Thomas (2001)</xref> and <xref ref-type="bibr" rid="B41">Van Heukelem and Hooker (2011)</xref>. Total chlorophyll-<italic>a</italic> (Tchla; the sum of monovinyl chlorophyll-<italic>a</italic>, divinyl chlorophyll a, chlorophyllide, and chlorophyll-<italic>a</italic> allomers and epimers) from HPLC was used here to compare to ocean color estimates (24 samples total were collected; 3 samples are shown here for comparison with exact matchups with ocean color data).</p>
</sec>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Generalized IOP algorithm (GIOP) applied to <italic>in situ</italic> <inline-formula>
<mml:math display="inline" id="im37">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>(<italic>&#x3bb;</italic>)</title>
<p>The Generalized IOP (GIOP) algorithm of <xref ref-type="bibr" rid="B47">Werdell et&#xa0;al. (2013)</xref> was used to invert <italic>in situ</italic> <inline-formula>
<mml:math display="inline" id="im38">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> measured with the C-OPS to retrieve spectral inherent optical properties. With the input of <inline-formula>
<mml:math display="inline" id="im39">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, GIOP uses the OC4 algorithm to estimate chlorophyll concentration (<ext-link ext-link-type="uri" xlink:href="https://oceancolor.gsfc.nasa.gov/atbd/chlor_a/">https://oceancolor.gsfc.nasa.gov/atbd/chlor_a/</ext-link>) and model particle backscattering (<inline-formula>
<mml:math display="inline" id="im40">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>), total absorption by CDOM and particles (<inline-formula>
<mml:math display="inline" id="im41">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>&#x261;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>), absorption by phytoplankton (<inline-formula>
<mml:math display="inline" id="im42">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>), and absorption by CDOM and depigmented particles (<inline-formula>
<mml:math display="inline" id="im43">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x261;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>). Since the GIOP in the default configuration was developed for the open ocean and the Santa Barbara Channel is a coastal site, the slope of <inline-formula>
<mml:math display="inline" id="im44">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> was adjusted from the default value of 0.018 to 0.012 to reflect the higher values typically found in coastal regions (<xref ref-type="bibr" rid="B32">Roesler et&#xa0;al., 1989</xref>) and compare the resulting model performance vs. <italic>in situ</italic> IOP data. All other parameters were used at their default values. While the GIOP is typically applied to satellite data, in our study, the GIOP was <italic>only</italic> applied to the C-OPS <inline-formula>
<mml:math display="inline" id="im45">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra. Satellite <inline-formula>
<mml:math display="inline" id="im46">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> data could not be inverted due to failures in the GIOP given the data quality.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>MODIS-Aqua Level 2 1-km products</title>
<p>MODIS-Aqua Level 2 1-km remote sensing reflectance products were downloaded from the NASA Ocean Biology Processing Group (<ext-link ext-link-type="uri" xlink:href="https://oceancolor.gsfc.nasa.gov">https://oceancolor.gsfc.nasa.gov</ext-link>). The most recent reprocessing (R2022) was used. While there were MODIS-Aqua images for each day of <italic>in situ</italic> sampling, there were only concurrent images of the location where C-OPS, inherent optical properties (IOPs), and HPLC samples were collected on December 19-21, 2017. Following the matchup criteria presented by <xref ref-type="bibr" rid="B4">Bailey and Werdell (2006)</xref>, <inline-formula>
<mml:math display="inline" id="im47">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> values were averaged for the 5x5 pixel box around the <italic>in situ</italic> observations. While many pixels in the match-up region were flagged due to high polarization or other product failures, no pixels contained any of the exclusion flags specified by <xref ref-type="bibr" rid="B4">Bailey and Werdell (2006)</xref>. It is of particular note that no spectra were flagged for a failed atmospheric correction algorithm (given by the ATMFAIL flag) and the product failures related to failed derived products (i.e., particulate inorganic carbon; PIC) rather than any other common flags (Sean Bailey, <italic>pers. comm.</italic>). However, the coefficient of variation was greater than 0.15 for multiple wavelengths of the December 19 and December 21 spectra, meaning that only the December 20 spectrum passed all <xref ref-type="bibr" rid="B4">Bailey and Werdell (2006)</xref> criteria. The spectra provide relevant context for the <italic>in situ</italic> samples and highlight the challenges of relying on ocean color when atmospheric interference is higher (e.g., during a coastal wildfire). In addition to <inline-formula>
<mml:math display="inline" id="im48">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra, we also retrieved the standard mass concentration of chlorophyll-<italic>a</italic> provided in the Level-2 data files (details at <ext-link ext-link-type="uri" xlink:href="https://oceancolor.gsfc.nasa.gov/atbd/chlor_a/">https://oceancolor.gsfc.nasa.gov/atbd/chlor_a/</ext-link>). Chlorophyll-<italic>a</italic> from MODIS Aqua compares favorably to <italic>in situ</italic> measurements in this region (<xref ref-type="bibr" rid="B34">Siegel et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B27">Li et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B2">Anderson et&#xa0;al., 2008</xref>).</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Plumes and Blooms sampling</title>
<p>Optical data from winter monthly cruises (defined as November to January) collected through the Plumes and Blooms sampling program (PnB; <ext-link ext-link-type="uri" xlink:href="http://www.oceancolor.ucsb.edu/plumes_and_blooms/">http://www.oceancolor.ucsb.edu/plumes_and_blooms/</ext-link>) were used to compare the historical context of the Santa Barbara Channel during years without wildfires. Plumes and Blooms has monitored the bio-optical variability of the Santa Barbara Channel along a consistent, 7-station transect approximately monthly since 1996. Samples were collected for <inline-formula>
<mml:math display="inline" id="im49">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im50">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im51">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> using the same methods detailed above. Only surface data from Station 4 (34&#x2019;15.01&#xb0;N, 119&#x2019;54.38&#xb0;W) were used for this comparison, as it was collocated with the sampling locations use throughout the Thomas Fire cruise. A total of 28 samples, collected in winter cruises between 2005-2016 and in 2018, were used here.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<p>The impact of wildfire smoke and ash on reflectance as measured by the C-OPS <italic>in situ</italic> was compared to reflectance measured by MODIS-Aqua via remote sensing. The shape and magnitude of C-OPS <inline-formula>
<mml:math display="inline" id="im52">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra were relatively consistent between the three days of sampling (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>). Relatively lower <inline-formula>
<mml:math display="inline" id="im53">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> values in the blue wavelengths (350-450nm) suggest high absorption by CDOM and NAP. A chlorophyll fluorescence peak is present at 670nm, consistent with the relatively high chlorophyll-<italic>a</italic> concentrations measured from HPLC on these days (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>). Alternately, the shape and magnitude of <inline-formula>
<mml:math display="inline" id="im54">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra measured from MODIS-Aqua are quite different across the three days of sampling and do not match the C-OPS spectra measured on those same days and times (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Chlorophyll-<italic>a</italic> concentration calculated using the OC3M algorithm (NASA GSFC; <uri xlink:href="https://oceancolor.gsfc.nasa.gov/atbd/chlor_a/">https://oceancolor.gsfc.nasa.gov/atbd/chlor_a/</uri>) for <bold>(A)</bold> C-OPS samples and <bold>(B)</bold> MODIS-Aqua samples, both overlaid on the NOMAD dataset for context (<xref ref-type="bibr" rid="B48">Werdell and Bailey, 2005</xref>). Chlorophyll-<italic>a</italic> concentrations from the three match-up days are compared in <bold>(C)</bold> to concentrations measured by HPLC.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1267681-g003.tif"/>
</fig>
<p>The discrepancies between satellite and <italic>in situ</italic> ocean color spectra necessarily result in different chlorophyll-<italic>a</italic> retrievals from the OC3M algorithm, which uses a ratio between the blue and green wavelengths to calculate chlorophyll-<italic>a</italic> concentration. The MODIS-Aqua spectra are theoretically unrealistic, in that these spectral shapes do not reflect realistic scenarios for radiative transfer in the ocean. When the derived chlorophyll-<italic>a</italic> concentrations from these satellite and <italic>in situ</italic> <inline-formula>
<mml:math display="inline" id="im55">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra are compared directly to each other and to chlorophyll-<italic>a</italic> measured from HPLC, chlorophyll-<italic>a</italic> from MODIS-Aqua exceeds chlorophyll-<italic>a</italic> from C-OPS and HPLC by a factor of 3 on December 19 and by a factor of 2 on December 21 (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>). On December 20, MODIS-Aqua chlorophyll-<italic>a</italic> is lower than HPLC and C-OPS chlorophyll-<italic>a</italic>. In all cases, the MODIS-Aqua chlorophyll-<italic>a</italic> estimates demonstrate the poorest relationship with HPLC chlorophyll-<italic>a</italic> or chlorophyll-<italic>a</italic> from C-OPS. The chlorophyll-<italic>a</italic> concentrations retrieved from C-OPS reflectances (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3A</bold>
</xref>) and MODIS-Aqua reflectances (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3B</bold>
</xref>) are within normal ranges, particularly when compared to the range of samples included in the NASA bio-Optical Marine Algorithm Dataset (NOMAD). However, the MODIS-Aqua values are on the outer limits of the relationship between observed between modeled and measured chlorophyll-<italic>a</italic> concentrations in the NOMAD dataset.</p>
<p>The bio-optical context in which these samples were collected is considered more holistically via the direct measurement of inherent optical properties (IOPs). Comparing the average <inline-formula>
<mml:math display="inline" id="im56">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im57">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im58">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> coefficients during the Thomas Fire (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4A</bold>
</xref>) to the mean values of these properties during a typical December (2005-2018) in the Santa Barbara Channel from PnB highlights the differences during the Thomas Fire (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4B</bold>
</xref>). Non-algal particulate absorption and CDOM absorption had similar magnitudes in December 2017 to historical Decembers, but phytoplankton absorption was relatively lower in December 2017. Accordingly, the <inline-formula>
<mml:math display="inline" id="im62">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> to <inline-formula>
<mml:math display="inline" id="im63">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> ratio and <inline-formula>
<mml:math display="inline" id="im64">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> to <inline-formula>
<mml:math display="inline" id="im65">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> ratio were much higher during the Thomas Fire than is expected in the Santa Barbara Channel in December (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4C</bold>
</xref>). The absorption budgets demonstrate that <inline-formula>
<mml:math display="inline" id="im66">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> dominated the overall absorption in the blue and green wavelengths in December 2017 compared to the relatively high phytoplankton absorption signal relative to other absorbing components in a typical December (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4D</bold>
</xref>).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Mean absorption by CDOM (<inline-formula>
<mml:math display="inline" id="im59">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, yellow), phytoplankton (<inline-formula>
<mml:math display="inline" id="im60">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, green), and non-algal particles (<inline-formula>
<mml:math display="inline" id="im61">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, brown) from <bold>(A)</bold> December 2017 and <bold>(B)</bold> other Decembers in the Plumes and Blooms dataset. Shading shows +/- one standard deviation of the mean value at each wavelength. The fractional contribution of each component to total absorption at each wavelength for these datasets is shown in <bold>(C)</bold> for December 2017 and <bold>(D)</bold> for other Decembers from Plumes and Blooms Station 4. Pure seawater absorption is also shown (<xref ref-type="bibr" rid="B30">Pope and Fry, 1997</xref>).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1267681-g004.tif"/>
</fig>
<p>Finally, the reflectance and IOP data were combined to consider the effects of wildfire smoke and ash on optical modeling. The measured IOP spectra were compared to the modeled IOP spectra from the inversion of C-OPS data measured in situ. GIOP retrieved the magnitude and relative shape of the mean <inline-formula>
<mml:math display="inline" id="im67">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra with high accuracy (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5A</bold>
</xref>), but underestimated the magnitude of CDOM+NAP absorption (<inline-formula>
<mml:math display="inline" id="im71">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x261;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>; <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5B</bold>
</xref>), resulting in an underestimation of the total summed particulate absorption term (<inline-formula>
<mml:math display="inline" id="im72">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>&#x261;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>; <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5C</bold>
</xref>). When the slope of <inline-formula>
<mml:math display="inline" id="im73">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> was adjusted to reflect the conditions expected in coastal regions (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>; <xref ref-type="bibr" rid="B32">Roesler et&#xa0;al., 1989</xref>; <xref ref-type="bibr" rid="B47">Werdell et&#xa0;al., 2013</xref>), the mean <inline-formula>
<mml:math display="inline" id="im74">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>&#x261;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra are the same magnitude as the mean measured <inline-formula>
<mml:math display="inline" id="im75">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra at most wavelengths (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2C</bold>
</xref>), but the magnitude of <inline-formula>
<mml:math display="inline" id="im76">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> was underestimated relative to measured values (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2A</bold>
</xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Mean measured (solid lines) and modeled (dashed lines) inherent optical properties. <bold>(A)</bold> Phytoplankton absorption (<inline-formula>
<mml:math display="inline" id="im68">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>), <bold>(B)</bold> summed CDOM and depigmented particulate absorption (<inline-formula>
<mml:math display="inline" id="im69">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>), and <bold>(C)</bold> total particulate and CDOM absorption (<inline-formula>
<mml:math display="inline" id="im70">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>). Shading shows +/- one standard deviation of the mean value at each wavelength.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1267681-g005.tif"/>
</fig>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>Our study leverages satellite data, <italic>in situ</italic> data, and time series data to explore the impacts of wildfires on ocean properties, including the potential effects of smoke, ash particle deposition, and increased dissolved organic matter from ash leachate. Of these three data sources, the satellite data were the easiest and cheapest data to collect. However, for our study region during this sampling period, these data were not of sufficient quality to evaluate ocean color or derived products. The variability in the MODIS-Aqua <inline-formula>
<mml:math display="inline" id="im77">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra, which were measured from the top of atmosphere to capture ocean color just above the surface of the water, are consistent with the variable smoke coverage over the sampling site on each day of sampling (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>; <xref ref-type="bibr" rid="B23">Kramer et&#xa0;al., 2020</xref>). These spectra result in excessively high (December 19, 21) or low (December 20) chlorophyll-<italic>a</italic> concentrations compared to HPLC (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>). The MODIS-Aqua <inline-formula>
<mml:math display="inline" id="im78">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra show three different but theoretically unrealistic spectral shapes compared to the consistent measurements of <inline-formula>
<mml:math display="inline" id="im79">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> from the C-OPS (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). Only the spectra measured on December 20 passed the <xref ref-type="bibr" rid="B4">Bailey and Werdell (2006)</xref> quality control assessment for validation-quality ocean color data, and none of these spectra were suitable for inversion using GIOP. Furthermore, these spectra were only flagged due to high polarization. More specifically, none of the flags masked in Level 3 data were triggered, meaning that these pixels would have been included in the generation of Level 3 data products. As Level 3 data are generally used for ocean color applications, these pixels could potentially be used to assess the impact of wildfires on marine ecosystems, despite clear data quality issues.</p>
<p>The difficulty of retrieving accurate chlorophyll concentrations from remote sensing methods has been observed in other optically-complex marine systems with impacts of dust, atmospheric aerosols, or high CDOM concentrations (<xref ref-type="bibr" rid="B11">Claustre et&#xa0;al., 2002</xref>; <xref ref-type="bibr" rid="B35">Stramska et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B26">Lewis and Arrigo, 2020</xref> and references therein). While these exact effects have not been observed <italic>in situ</italic> for coastal wildfire conditions prior to this study, our results are consistent with the optical properties of other ecosystems influenced by absorbing constituents in the atmosphere and/or ocean. In assessing the impact of the wildfire smoke and ash on the quality of the satellite <inline-formula>
<mml:math display="inline" id="im80">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> data, documented challenges with atmospheric correction of <inline-formula>
<mml:math display="inline" id="im81">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> data over the Santa Barbara Channel were also considered (e.g., <xref ref-type="bibr" rid="B20">Henderikx Freitas et&#xa0;al., 2017</xref>). When the spectra were assessed with the Quality Water Index Polynomial (QWIP; <xref ref-type="bibr" rid="B14">Dierssen et&#xa0;al., 2022</xref>), the remote sensing data collected over the Santa Barbara Channel in December 2017 during the Thomas Fire had a QWIP score that was comparable to the MODIS-Aqua timeseries over the same region in Decembers 2002-2022 (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;3</bold>
</xref>), suggesting that there was still the average amount of usable data collected in this region during the Thomas Fire as in wildfire-free years. Other aerosol sources may also contribute to poor atmospheric correction for remote sensing data. In this case study, the role of marine aerosol-producing phytoplankton (e.g., dimethylsulfide producers) can be ruled out following a high-resolution characterization of the phytoplankton community during this period (<xref ref-type="bibr" rid="B23">Kramer et&#xa0;al., 2020</xref>). Ultimately, given the quality of the satellite <inline-formula>
<mml:math display="inline" id="im82">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> data, these data were not useful for evaluating the color of the ocean or the potential ecological impact of the wildfire on the ocean.</p>
<p>In this study, we were fortunate to have a serendipitous research cruise already planned in the region where the fire broke out (<xref ref-type="bibr" rid="B6">Bisson et&#xa0;al., 2020</xref>). The <italic>in situ</italic> samples collected alongside <italic>in situ</italic> ocean color measurements in this study provide important points of comparison and more detailed information about the bio-optical regime of the Santa Barbara Channel during the Thomas Fire. The <italic>in situ</italic> reflectance data demonstrate the impact of the Thomas Fire on the coastal ocean. The C-OPS <inline-formula>
<mml:math display="inline" id="im83">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra, which were measured below the wildfire smoke plume and just below the surface of the water, are consistent in shape and magnitude, resulting in more accurate chlorophyll-<italic>a</italic> retrievals relative to HPLC (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>). These spectra were not impacted by the same quality control concerns as the satellite ocean color data, and were more suitable for both direct analysis (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>) and inversion approaches (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>). While the short length of this cruise resulted in relatively few direct matchups between in-water <inline-formula>
<mml:math display="inline" id="im84">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and satellite <inline-formula>
<mml:math display="inline" id="im85">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, we collected surface in-water optical samples in duplicate or triplicate, leading to a larger dataset against which to compare the satellite samples (21-24 samples, depending on the variable).</p>
<p>The <italic>in situ</italic> IOP data collected alongside the reflectance data help to further demonstrate some of the variability in different oceanic components during the Thomas Fire compared to a typical December in the Santa Barbara Channel. When the <italic>in situ</italic> absorption spectra are compared to the absorption by different optically significant components estimated by inversion modeling, the strengths and weaknesses of that type of ocean color model can be considered in a coastal ecosystem undergoing an extreme ecological event, such as during the Thomas Fire. The results shown here demonstrate that inversion models may capture some components of the reflectance signal (<inline-formula>
<mml:math display="inline" id="im86">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5A</bold>
</xref>) better than others (<inline-formula>
<mml:math display="inline" id="im87">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im88">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, <xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5B, C</bold>
</xref>). GIOP solves for <inline-formula>
<mml:math display="inline" id="im89">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> as a balance of absorption to backscattering (Equations 1-3). In this system, the unexpectedly high <italic>in situ</italic> <inline-formula>
<mml:math display="inline" id="im90">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> signal during the Thomas Fire relative to other absorbing constituents (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4C</bold>
</xref>) was not retrieved by inversion modeling using the GIOP in its default configuration. When a higher <inline-formula>
<mml:math display="inline" id="im91">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> slope was applied, the shape and magnitude of the <inline-formula>
<mml:math display="inline" id="im92">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x261;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra and <inline-formula>
<mml:math display="inline" id="im93">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>&#x261;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra were a better fit to the measured data (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2B, C</bold>
</xref>). However, the <inline-formula>
<mml:math display="inline" id="im94">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra underestimated the magnitude of the measured <inline-formula>
<mml:math display="inline" id="im95">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> data (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>), demonstrating the difficulty of capturing absorption by all optical components in a complex region with high concentrations of phytoplankton and other absorption constituents, such as the coastal ocean during a wildfire.</p>
<p>Finally, the PnB time series was leveraged to put the optical and pigment data in a broader historical and ecological context. In <xref ref-type="bibr" rid="B23">Kramer et&#xa0;al. (2020)</xref>, the PnB pigment data contextualized the abundance and magnitude of different phytoplankton groups during the Thomas Fire relative to a typical December. In this work, the <italic>in situ</italic> IOP data from PnB was compared to the data collected during the fire. The relatively low phytoplankton absorption coefficient in December 2017 (<xref ref-type="fig" rid="f4">
<bold>Figures&#xa0;4A, B</bold>
</xref>) was expected given the lower chlorophyll-<italic>a</italic> concentrations measured in December 2017 compared to the PnB average (<xref ref-type="bibr" rid="B23">Kramer et&#xa0;al., 2020</xref>), and the magnitude of the measured <inline-formula>
<mml:math display="inline" id="im96">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> in December 2017 matches the values retrieved from <italic>in situ</italic> ocean color inversion during this same period (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5A</bold>
</xref>). The relatively high <inline-formula>
<mml:math display="inline" id="im97">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>&#x261;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> coefficient during the Thomas Fire compared to other IOPs (<xref ref-type="fig" rid="f4">
<bold>Figures&#xa0;4C, D</bold>
</xref>) is consistent with recent studies showing increased deposition of black carbon and dissolved organic matter from wildfires to the coastal ocean (<xref ref-type="bibr" rid="B43">Wagner et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B13">Coward et&#xa0;al., 2022</xref>). Despite the relatively high chlorophyll-<italic>a</italic> concentrations during the Thomas Fire, the optical signal was not dominated by phytoplankton, but rather by CDOM absorption.</p>
<p>While satellite data are the cheapest and most readily-available form of data to leverage during an extreme event such as the Thomas Fire, for our study region and during this wildfire, these data did not have sufficient quality to allow evaluation of ocean color. Instead, <italic>in situ</italic> measurements of <inline-formula>
<mml:math display="inline" id="im98">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, IOPs, and chlorophyll-<italic>a</italic> allowed for a system-wide evaluation of the conditions in the SBC during the Thomas Fire. The <italic>in situ</italic> optical and pigment data from the PnB time series also contextualized these results within the broader context of the region and ecosystem. In many recent studies examining wildfire impacts on the ocean, particularly in more remote or open ocean environments, <italic>in situ</italic> samples are not available (<xref ref-type="bibr" rid="B3">Ardyna et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B44">Wang et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B45">Weis et&#xa0;al., 2022</xref>). Inversion modeling and derived products are essential in these cases: all three studies used chlorophyll-<italic>a</italic> from remote sensing to consider the possibility of a phytoplankton bloom occurring in the wake of a wildfire smoke plume. The relationship between <inline-formula>
<mml:math display="inline" id="im99">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and phytoplankton carbon (<inline-formula>
<mml:math display="inline" id="im100">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>y</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) was also used to further describe the impact of wildfires on primary productivity beyond the chlorophyll-<italic>a</italic> concentration (i.e., photoacclimation, physiological stress; <xref ref-type="bibr" rid="B44">Wang et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B45">Weis et&#xa0;al., 2022</xref>). The accuracy of these two parameters (chlorophyll, <inline-formula>
<mml:math display="inline" id="im101">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>y</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) is paramount to the interpretation of results and significance. The results shown here demonstrate that the derivation of these parameters is certainly also impacted by the wildfire, particularly from remote sensing. <italic>In situ</italic> sampling of IOPs and reflectance highlights an optical signature of the wildfire, but the accuracy of these measurements was only made possible by shipboard collection below the smoke plume, and the anomaly of these results was only made clear by comparison to time series data.</p>
<p>While ocean color has been and will remain an important, economical tool for considering the effects of wildfires on marine ecosystems, more high-quality measurements of IOPs and reflectance <italic>in situ</italic> during wildfires are urgently needed. As remote sensing approaches gain popularity for monitoring and evaluating the impacts of wildfires on the coastal and open ocean, it is essential to compare the performance of these approaches with <italic>in situ</italic> measurements in different ecosystems and regions. The impacts of wildfires on the ocean may often be nonlinear (e.g., wildfire smoke may block available light for photosynthesis, while leaching ash may provide limiting nutrients). It is important to have in-water measurements to evaluate these impacts in near real-time, particularly to compare to daily or weekly satellite measurements. Our results demonstrate that the impacts of wildfires on ocean ecology and optics may also vary based on the proximity to the wildfire. In the coastal ocean, increased CDOM from ash leaching may impact both ocean color (as seen in the shape and magnitude of the <italic>in situ</italic> C-OPS <inline-formula>
<mml:math display="inline" id="im102">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> spectra; <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>) and <italic>in situ</italic> IOPs (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5B, C</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figures&#xa0;2B, C</bold>
</xref>). <italic>In situ</italic> measurements will also allow for an evaluation of the impact of these other optically significant materials on the estimation of chlorophyll-<italic>a</italic> from remote sensing, particularly if further conclusions about ecosystem impacts will be drawn from an enhanced chlorophyll signal (e.g., <xref ref-type="bibr" rid="B3">Ardyna et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B44">Wang et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B45">Weis et&#xa0;al., 2022</xref>).</p>
<p>Finally, this study has demonstrated the ongoing challenge of &#x201c;known unknowns&#x201d; in work at the intersection of unpredictable extreme events (e.g., wildfires, volcanic eruptions) and ocean ecosystems (Bisson et&#xa0;al., <italic>in revision</italic>). While many elements of the system can be described or directly quantified (such as chlorophyll-<italic>a</italic> concentrations, absorption coefficients, reflectance spectra), there are many unknowns that cannot always be captured or investigated on the appropriate timescale&#x2014;the thickness of the smoke plume throughout the day, the amount of ash deposition and variability in ash composition, the rate of ash dissolution in natural seawater compared to laboratory studies, the retentiveness of the physical environment, etc. The combination of the dynamic ocean system and the unpredictable, ever-changing landscape of a wildfire may not be possible to capture in full detail, but the combination of remote sensing observations and <italic>in situ</italic> validation studies will ideally help to constrain the (near) real-time impacts of the fire on the ocean. As wildfires continue to increase in size and severity, particularly in coastal regions, future studies to quantify these impacts will be essential. Furthermore, the advent of satellite sensors with higher spectral resolution (e.g., NASA&#x2019;s Plankton Aerosol Cloud and ocean Ecosystem [PACE] sensor; <xref ref-type="bibr" rid="B46">Werdell et&#xa0;al., 2019</xref>) will allow for more information to be gathered from remote sensing during coastal wildfires, such as measurements of phytoplankton pigment concentrations that can link to community composition (<xref ref-type="bibr" rid="B24">Kramer et&#xa0;al., 2022</xref>). With these higher-resolution measurements and derived products, quality control and <italic>in situ</italic> samples for validation will remain essential, particularly during significant ecological disturbances.</p>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusions and recommendations</title>
<p>We close with three concrete recommendations for wildfire-ocean work going forward, based on the lessons learned and results shown here.</p>
<list list-type="simple">
<list-item>
<p>1) Coastal timeseries are essential for continuing observations and to provide historical context during a wildfire or other extreme event/ecological disturbance. In this study, the PnB dataset was compared to the higher-resolution sampling during the Thomas Fire to consider those results in the context of a typical December in the Santa Barbara Channel. As much as possible, we encourage the integration of optical measurements to these coastal time series, which are critical for remote sensing validation. Both reflectance and IOP measurements will help to confirm the patterns shown in ocean color observations and constrain the outcome of derived products.</p>
</list-item>
<list-item>
<p>2) Planning field work to collect <italic>in situ</italic> samples around an unpredictable event such as a wildfire may not be possible, and thus remote sensing will remain a vital tool for monitoring; however, whenever possible, we strongly recommend the collection of <italic>in situ</italic> samples to validate and confirm the results of these remote sensing measurements. Sample collection is more difficult in open ocean regions, but other autonomous resources beneath the smoke plume can be used in these cases when possible (e.g., floats, buoys, AERONET sensors, etc.).</p>
</list-item>
<list-item>
<p>3) We encourage funding agencies to develop creative opportunities for providing money and resources quickly when wildfires occur in coastal regions or when wildfire impacts in open ocean regions are observed from remote sensing (similar to the RAPID model, but ideally providing resources on timescales of days to weeks within the realistic course of most wildfires). Contingency plans should be in place for use of existing research vessels or smaller ships during wildfires or other extreme events, and where possible, funding should be set aside to mobilize for sample collection and analysis during these times. While these unpredictable events provide a logistical and financial challenge, more research is clearly needed in this area. The extremely limited datasets presented in existing <italic>in situ</italic> studies are likely only representative of those regions at that one moment in time, and a more holistic understanding of the impacts of wildfires on marine systems can only be reached with an investment of time, resources, and creative thinking.</p>
</list-item>
</list>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: <ext-link ext-link-type="uri" xlink:href="https://seabass.gsfc.nasa.gov/cruise/ACIDD_2017">https://seabass.gsfc.nasa.gov/cruise/ACIDD_2017</ext-link>, <ext-link ext-link-type="uri" xlink:href="https://seabass.gsfc.nasa.gov/experiment/Plumes_and_Blooms">https://seabass.gsfc.nasa.gov/experiment/Plumes_and_Blooms</ext-link>, <ext-link ext-link-type="uri" xlink:href="https://www.earthdata.nasa.gov/eosdis/daacs/obdaac">https://www.earthdata.nasa.gov/eosdis/daacs/obdaac</ext-link>.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>SK: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. KB: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Writing &#x2013; review &amp; editing. CM: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Visualization, Writing &#x2013; review &amp; editing.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The University of California&#x2019;s Ship Funds Program, the National Academies Keck Futures Initiative, UCSB&#x2019;s Coastal Fund, NASA grant #NNX17AK04G to Kyla Drushka, and NSF grant #1821916 to David Valentine provided funding for this project. The Plumes and Blooms time series is supported by NASA Ocean Biology and Biogeochemistry. SK is supported by NASA grant #80NSSC21K0015 to Colleen A. Durkin, the David and Lucille Packard Foundation, and a Simons Foundation Postdoctoral Fellowship in Marine Microbial Ecology. CM is supported by NPRB grant #2101A to CM. Funding for publication comes from NASA TWSC grant #80NSSC22K0049 to KB.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>This work could not have been accomplished without the enthusiastic support of the Captain and crew of R/V <italic>Sally Ride</italic> and all other members of the ACIDD science team. Thank you to Crystal Thomas at NASA GSFC for processing the ACIDD and PnB HPLC samples with high precision and accuracy. Thank you also to Paula Bontempi, Laura Lorenzoni, and Ivona Cetini&#x107; for support with HPLC sample processing. Collin Roesler and Aimee Neeley provided helpful comments on the phytoplankton absorption data. Jeremy Werdell, Amir Ibrahim, and Sean Bailey provided important context for the remote sensing data. Thank you to Nathalie Guillocheau and Dylan Catlett for PnB sample help. We are grateful for support from Dave Siegel. Comments from Robert Frouin and four reviewers helped to improve earlier versions of this manuscript.</p>
</ack>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmars.2023.1267681/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmars.2023.1267681/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet_1.pdf" id="SM1" mimetype="application/pdf"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abram</surname> <given-names>N. J.</given-names>
</name>
<name>
<surname>Gagan</surname> <given-names>M. K.</given-names>
</name>
<name>
<surname>McCulloch</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Chappell</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Hantoro</surname> <given-names>W. S.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Coral reef death during the 1997 Indian Ocean dipole linked to Indonesian wildfires</article-title>. <source>Science</source> <volume>301</volume> (<issue>5635</issue>), <fpage>952</fpage>&#x2013;<lpage>955</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1083841</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Anderson</surname> <given-names>C. R.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Kudela</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Brzezinski</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Empirical models of toxigenic Pseudo-nitzschia blooms: Potential use as a remote detection tool in the Santa Barbara Channel</article-title>. <source>Harmful. Algae.</source> <volume>8</volume> (<issue>3</issue>), <fpage>478</fpage>&#x2013;<lpage>492</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.hal.2008.10.005</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ardyna</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hamilton</surname> <given-names>D. S.</given-names>
</name>
<name>
<surname>Harmel</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Lacour</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Bernstein</surname> <given-names>D. N.</given-names>
</name>
<name>
<surname>Lalibert&#xe9;</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Wildfire aerosol deposition likely amplified a summertime Arctic phytoplankton bloom</article-title>. <source>Commun. Earth Environ.</source> <volume>3</volume> (<issue>1</issue>), <fpage>201</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s43247-022-00511-9</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bailey</surname> <given-names>S. W.</given-names>
</name>
<name>
<surname>Werdell</surname> <given-names>P. J.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>A multi-sensor approach for the on-orbit validation of ocean color satellite data products</article-title>. <source>Remote Sens. Environ.</source> <volume>102</volume>, <fpage>12</fpage>&#x2013;<lpage>23</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.rse.2006.01.015</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barr&#xf3;n</surname> <given-names>R. K.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Guillocheau</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Evaluating the importance of phytoplankton community structure to the optical properties of the Santa Barbara Channel, California</article-title>. <source>Limnol. Oceanogr.</source> <volume>59</volume> (<issue>3</issue>), <fpage>927</fpage>&#x2013;<lpage>946</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.2014.59.3.0927</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bisson</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Baetge</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Kramer</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Catlett</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Girling</surname> <given-names>G.</given-names>
</name>
<name>
<surname>McNair</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>California wildfire burns boundaries between science and art</article-title>. <source>Oceanography</source> <volume>33</volume> (<issue>1</issue>), <fpage>16</fpage>&#x2013;<lpage>19</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5670/oceanog.2020.110</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bisson</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Boss</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Werdell</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Ibrahim</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Frouin</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Behrenfeld</surname> <given-names>M. J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Seasonal bias in global ocean color observations</article-title>. <source>Appl. Optics.</source> <volume>60</volume> (<issue>23</issue>), <fpage>6978</fpage>&#x2013;<lpage>6988</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1364/AO.426137</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bisson</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Gass&#xf3;</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Mahowald</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Wagner</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Koffman</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Carn</surname> <given-names>S. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Observing ocean ecosystem responses to volcanic ash</article-title>. <source>Remote Sens. Environ.</source> <volume>296</volume>, <fpage>1</fpage>&#x2013;<lpage>13</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.rse.2023.113749</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Browning</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Stone</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Bouman</surname> <given-names>H. A.</given-names>
</name>
<name>
<surname>Mather</surname> <given-names>T. A.</given-names>
</name>
<name>
<surname>Pyle</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Moore</surname> <given-names>C. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Volcanic ash supply to the surface ocean&#x2014;remote sensing of biological responses and their wider biogeochemical significance</article-title>. <source>Front. Mar. Sci.</source> <volume>2</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmars.2015.00014</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cira</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Bafna</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Kong</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Holt</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Ginger</surname> <given-names>L.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Turbidity and fecal indicator bacteria in recreational marine waters increase following the 2018 Woolsey Fire</article-title>. <source>Sci. Rep.</source> <volume>12</volume> (<issue>1</issue>), <fpage>2428</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-022-05945-x</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Claustre</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Morel</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Hooker</surname> <given-names>S. B.</given-names>
</name>
<name>
<surname>Babin</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Antoine</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Oubelkheir</surname> <given-names>K.</given-names>
</name>
<etal/>
</person-group>. (<year>2002</year>). <article-title>Is desert dust making oligotrophic waters greener</article-title>? <source>Geophys. Res. Lett.</source> <volume>29</volume> (<issue>10</issue>), <fpage>107</fpage>&#x2013;<lpage>101</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2001GL014056</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Coombs</surname> <given-names>J. S.</given-names>
</name>
<name>
<surname>Melack</surname> <given-names>J. M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Initial impacts of a wildfire on hydrology and suspended sediment and nutrient export in California chaparral watersheds</article-title>. <source>Hydrol. Processes.</source> <volume>27</volume>, <fpage>3842</fpage>&#x2013;<lpage>3851</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/hyp.9508</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Coward</surname> <given-names>E. K.</given-names>
</name>
<name>
<surname>Seech</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Carter</surname> <given-names>M. L.</given-names>
</name>
<name>
<surname>Flick</surname> <given-names>R. E.</given-names>
</name>
<name>
<surname>Grassian</surname> <given-names>V. H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Of sea and smoke: evidence of marine dissolved organic matter deposition from 2020 western United States wildfires</article-title>. <source>Environ. Sci. Technol. Lett.</source> <volume>9</volume> (<issue>10</issue>), <fpage>869</fpage>&#x2013;<lpage>876</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1021/acs.estlett.2c00383</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dierssen</surname> <given-names>H. M.</given-names>
</name>
<name>
<surname>Vandermeulen</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Barnes</surname> <given-names>B. B.</given-names>
</name>
<name>
<surname>Castagna</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Knaeps</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Vanhellemont</surname> <given-names>Q.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>QWIP: A quantitative metric for quality control of aquatic reflectance spectral shape using the apparent visible wavelength</article-title>. <source>Front. Remote Sens.</source> <volume>3</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/frsen.2022.869611</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Doerr</surname> <given-names>S. H.</given-names>
</name>
<name>
<surname>Sant&#xed;n</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Global trends in wildfire and its impacts: perceptions versus realities in a changing world</article-title>. <source>Philos. Trans. R. Soc. B.</source> <volume>371</volume>, <fpage>1</fpage>&#x2013;<lpage>10</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1098/rstb.2015.0345</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Duggen</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Croot</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Schacht</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Hoffmann</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Subduction zone volcanic ash can fertilize the surface ocean and stimulate phytoplankton growth: Evidence from biogeochemical experiments and satellite data</article-title>. <source>Geophys. Res. Lett.</source> <volume>34</volume>, <fpage>1</fpage>&#x2013;<lpage>5</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2006GL027522</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Franz</surname> <given-names>B. A.</given-names>
</name>
<name>
<surname>Cetini&#x107;</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Ibrahim</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Sayer</surname> <given-names>A. M</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Anomalous trends in global ocean carbon concentrations following the 2022 eruptions of Hunga Tonga-Hunga Ha&#x2019;apai</article-title>. <fpage>1&#x2013;20</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.21203/rs.3.rs-3335677/v1</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hamilton</surname> <given-names>D. S.</given-names>
</name>
<name>
<surname>Hantson</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Scott</surname> <given-names>C. E.</given-names>
</name>
<name>
<surname>Kaplan</surname> <given-names>J. O.</given-names>
</name>
<name>
<surname>Pringle</surname> <given-names>K. J.</given-names>
</name>
<name>
<surname>Nieradzik</surname> <given-names>L. P.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Reassessment of pre-industrial fire emissions strongly affects anthropogenic aerosol forcing</article-title>. <source>Nat. Commun.</source> <volume>9</volume> (<issue>1</issue>), <fpage>3182</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41467-018-05592-9</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hamme</surname> <given-names>R. C.</given-names>
</name>
<name>
<surname>Webley</surname> <given-names>P. W.</given-names>
</name>
<name>
<surname>Crawford</surname> <given-names>W. R.</given-names>
</name>
<name>
<surname>Whitney</surname> <given-names>F. A.</given-names>
</name>
<name>
<surname>DeGrandpre</surname> <given-names>M. D.</given-names>
</name>
<name>
<surname>Emerson</surname> <given-names>S. R.</given-names>
</name>
<etal/>
</person-group>. (<year>2010</year>). <article-title>Volcanic ash fuels anomalous plankton bloom in subarctic northeast Pacific</article-title>. <source>Geophys. Res. Lett.</source> <volume>37</volume>, <fpage>1</fpage>&#x2013;<lpage>5</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2010GL044629</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Henderikx Freitas</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Maritorena</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Fields</surname> <given-names>E.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Satellite assessment of particulate matter and phytoplankton variations in the Santa Barbara Channel and its surrounding waters: Role of surface waves</article-title>. <source>J. Geophys. Res.: Oceans.</source> <volume>122</volume> (<issue>1</issue>), <fpage>355</fpage>&#x2013;<lpage>371</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2016JC012152</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kishino</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Takahashi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Okami</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Ichimura</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>Estimation of the spectral absorption coefficients of phytoplankton in the sea</article-title>. <source>Bull. Mar. Sci.</source> <volume>37</volume> (<issue>2</issue>), <fpage>634</fpage>&#x2013;<lpage>642</lpage>.</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kostadinov</surname> <given-names>T. S.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Maritorena</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Guillocheau</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Ocean color observations and modeling for an optically complex site: Santa Barbara Channel, California, USA</article-title>. <source>J. Geophys. Res.: Oceans.</source> <volume>112</volume> (<issue>C7</issue>), <fpage>1</fpage>&#x2013;<lpage>15</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2006JC003526</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kramer</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Bisson</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Fischer</surname> <given-names>A. D.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Observations of phytoplankton community composition in the Santa Barbara Channel during the Thomas Fire</article-title>. <source>J. Geophys. Res.: Oceans.</source> <volume>125</volume> (<issue>12</issue>), <fpage>1</fpage>&#x2013;<lpage>16</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2020JC016851</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kramer</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Maritorena</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Catlett</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Modeling surface ocean phytoplankton pigments from hyperspectral remote sensing reflectance on global scales</article-title>. <source>Remote Sens. Environ.</source> <volume>270</volume>, <elocation-id>112879</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.rse.2021.112879</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Carder</surname> <given-names>K. L.</given-names>
</name>
<name>
<surname>Arnone</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Deriving inherent optical properties from water color: a multiband quasi-analytical algorithm for optically deep waters</article-title>. <source>Appl. Optics.</source> <volume>41</volume> (<issue>27</issue>), <fpage>5755</fpage>&#x2013;<lpage>5772</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1364/AO.41.005755</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lewis</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Arrigo</surname> <given-names>K. R.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Ocean color algorithms for estimating chlorophyll a, CDOM absorption, and particle backscattering in the arctic ocean</article-title>. <source>J. Geophys. Res.: Oceans.</source> <volume>125</volume> (<issue>6</issue>), <elocation-id>e2019JC015706</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2019JC015706</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Stamnes</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Spurr</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Stamnes</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Simultaneous retrieval of aerosol and ocean properties by optimal estimation: SeaWiFS case studies for the Santa Barbara Channel</article-title>. <source>Int. J. Remote Sens.</source> <volume>29</volume> (<issue>19</issue>), <fpage>5689</fpage>&#x2013;<lpage>5698</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/01431160802007632</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Mobley</surname> <given-names>C. D.</given-names>
</name>
</person-group> (<year>1994</year>). <source>Light and Water: Radiative Transfer in Natural Waters</source> (<publisher-loc>New York</publisher-loc>: <publisher-name>Academic</publisher-name>).</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nelson</surname> <given-names>N. B.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Craig.</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Swan</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Smethie</surname> <given-names>W. M.</given-names>
</name>
<name>
<surname>Khatiwala</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Hydrography of chromophoric dissolved organic matter in the North Atlantic</article-title>. <source>Deep. Sea. Res. Part I.: Oceanogr. Res. Papers.</source> <volume>54</volume> (<issue>5</issue>), <fpage>710</fpage>&#x2013;<lpage>731</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.dsr.2007.02.006</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pope</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Fry</surname> <given-names>E. S.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Absorption spectrum (380&#x2013;700 nm) of pure water. II. Integrating cavity measurements</article-title>. <source>Appl. Optics.</source> <volume>36</volume> (<issue>33</issue>), <fpage>8710</fpage>&#x2013;<lpage>8723</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1364/AO.36.008710</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roesler</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Theoretical and experimental approaches to improve the accuracy of particulate absorption coefficients derived from the quantitative filter technique</article-title>. <source>Limnol. Oceanogr.</source> <volume>43</volume> (<issue>7</issue>), <fpage>1649</fpage>&#x2013;<lpage>1660</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.1998.43.7.1649</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roesler</surname> <given-names>C. S.</given-names>
</name>
<name>
<surname>Perry</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Carder</surname> <given-names>K. L.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Modeling in <italic>situ</italic> phytoplankton absorption from total absorption spectra in productive inland marine waters</article-title>. <source>Limnol. Oceanogr.</source> <volume>34</volume> (<issue>8</issue>), <fpage>1510</fpage>&#x2013;<lpage>1523</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.1989.34.8.1510</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Roesler</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Stramski</surname> <given-names>D.</given-names>
</name>
<name>
<surname>D&#x2019;Sa</surname> <given-names>E. J.</given-names>
</name>
<name>
<surname>R&#xf6;ttgers</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Reynolds</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>2018</year>). <source>Chapter 5: Spectrophotometric Measurements of Particulate Absorption Using Filter Pads (Ocean Optics &amp; Biogeochemistry Protocols for Satellite Ocean Colour Sensor Validation)</source> (<publisher-loc>Dartmouth, Canada</publisher-loc>: <publisher-name>International Ocean Colour Coordinating Group (IOCCG) in collaboration with National Aeronautics and Space Administration (NASA</publisher-name>), <fpage>35</fpage>.</citation>
</ref>
<ref id="B34">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Maritorena</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Nelson</surname> <given-names>N. B.</given-names>
</name>
</person-group> (<year>2003</year>). <source>Chapter 15. Plumes and Blooms: Modeling the Case II Waters of the Santa Barbara Channel (NASA Technical Memorandum No. 20040067991)</source> (<publisher-loc>Washington, D.C.</publisher-loc>: <publisher-name>University of California Santa Barbara</publisher-name>), <fpage>153</fpage>&#x2013;<lpage>159</lpage>. Available at: <uri xlink:href="https://ntrs.nasa.gov/api/citations/20040067991/downloads/20040067991.pdf">https://ntrs.nasa.gov/api/citations/20040067991/downloads/20040067991.pdf</uri>.</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stramska</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Stramski</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Cichocka</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Cieplak</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Wo&#x17a;niak</surname> <given-names>S. B.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Effects of atmospheric particles from Southern California on the optical properties of seawater</article-title>. <source>J. Geophys. Res.: Oceans.</source> <volume>113</volume> (<issue>C8</issue>), <fpage>1</fpage>&#x2013;<lpage>15</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2007JC004407</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stramski</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Reynolds</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Kaczmarek</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Uitz</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zheng</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Correction of pathlength amplification in the filter-pad technique for measurements of particulate absorption coefficient in the visible spectral region</article-title>. <source>Appl. Optics.</source> <volume>54</volume>, <fpage>6763</fpage>&#x2013;<lpage>6782</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1364/AO.54.006763</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tassan</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Ferrari</surname> <given-names>G. M.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>An alternative approach to absorption measurements of aquatic particles retained on filters</article-title>. <source>Limnol. Oceanogr.</source> <volume>40</volume>, <fpage>1358</fpage>&#x2013;<lpage>1368</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.1995.40.8.1358</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Toole</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>D. A.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Modes and mechanisms of ocean color variability in the Santa Barbara Channel</article-title>. <source>J. Geophys. Res.: Oceans.</source> <volume>106</volume> (<issue>C11</issue>), <fpage>26985</fpage>&#x2013;<lpage>27000</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2000JC000371</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Urmy</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Williamson</surname> <given-names>C. E.</given-names>
</name>
<name>
<surname>Leach</surname> <given-names>T. H.</given-names>
</name>
<name>
<surname>Schladow</surname> <given-names>S. G.</given-names>
</name>
<name>
<surname>Overholt</surname> <given-names>E. P.</given-names>
</name>
<name>
<surname>Warren</surname> <given-names>J. D.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Vertical redistribution of zooplankton in an oligotrophic lake associated with reduction in ultraviolet radiation by wildfire smoke</article-title>. <source>Geophys. Res. Lett.</source> <volume>43</volume>, <fpage>3746</fpage>&#x2013;<lpage>3753</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2016GL068533</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Van Heukelem</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Thomas</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Computer-assisted high-performance liquid chromatography method development with applications to the isolation and analysis of phytoplankton pigments</article-title>. <source>J. Chromatogr. A</source>, <fpage>910</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0378-4347(00)00603-4</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Van Heukelem</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Hooker</surname> <given-names>S. B.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>"The importance of a quality assurance plan for method validation and minimizing uncertainties in the HPLC analysis of phytoplankton pigments</article-title>. " in  <source>Phytoplankton Pigments: Characterization, Chemotaxonomy, and Applications in Oceanography</source>. Eds <person-group person-group-type="editor">
<name>
<surname>Roy</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Llewellyn</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Egeland</surname> <given-names>E. S.</given-names>
</name>
<name>
<surname>Johnsen</surname> <given-names>G.</given-names>
</name>
</person-group> (<publisher-loc>Cambridge, United Kingdom</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>), <fpage>195</fpage>&#x2013;<lpage>242</lpage>.</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>van Woesik</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Comment on &#x201c;Coral reef death during the 1997 Indian Ocean dipole linked to Indonesian wildfires</article-title>.&#x201d; <source>Science</source> <volume>303</volume> (<issue>5662</issue>), <fpage>1297</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1091983</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wagner</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Harvey</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Baetge</surname> <given-names>N.</given-names>
</name>
<name>
<surname>McNair</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Arrington</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Stubbins</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Investigating atmospheric inputs of dissolved black carbon to the santa barbara channel during the thomas fire (California, USA)</article-title>. <source>J. Geophys. Res.: Biogeosci.</source> <volume>126</volume> (<issue>8</issue>), <elocation-id>e2021JG006442</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2021JG006442</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.-H.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>R.</given-names>
</name>
<name>
<surname>He</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Xue</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Australian fire nourishes ocean phytoplankton bloom</article-title>. <source>Sci. Total. Environ.</source> <volume>807</volume>, <elocation-id>150775</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.scitotenv.2021.150775</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Weis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Schallenberg</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Chase</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Bowie</surname> <given-names>A. R.</given-names>
</name>
<name>
<surname>Wojtasiewicz</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Perron</surname> <given-names>M. M. G.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Southern ocean phytoplankton stimulated by wildfire emissions and sustained by iron recycling</article-title>. <source>Geophys. Res. Lett.</source> <volume>49</volume> (<issue>11</issue>), <elocation-id>e2021GL097538</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2021GL097538</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Werdell</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Behrenfeld</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Bontempi</surname> <given-names>P. S.</given-names>
</name>
<name>
<surname>Boss</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Cairns</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Davis</surname> <given-names>G. T.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>The Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission: Status, science, advances</article-title>. <source>Bull. Am. Meteorol. Soc.</source>, <fpage>1</fpage>&#x2013;<lpage>59</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1175/BAMS-D-18-0056.1</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Werdell</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Franz</surname> <given-names>B. A.</given-names>
</name>
<name>
<surname>Bailey</surname> <given-names>S. W.</given-names>
</name>
<name>
<surname>Feldman</surname> <given-names>G. C.</given-names>
</name>
<name>
<surname>Boss</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Brando</surname> <given-names>V. E.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>Generalized ocean color inversion model for retrieving marine inherent optical properties</article-title>. <source>Appl. Optics.</source> <volume>52</volume> (<issue>10</issue>), <fpage>2019</fpage>&#x2013;<lpage>2037</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1364/AO.52.002019</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Werdell</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Bailey</surname> <given-names>S. W.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>An improved in-situ bio-optical data set for ocean color algorithm development and satellite data product validation</article-title>. <source>Remote Sens. Environ.</source> <volume>98</volume>, <fpage>122</fpage>&#x2013;<lpage>140</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.rse.2005.07.001</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Westberry</surname> <given-names>T. K.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>Y. R.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Behrenfeld</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Remer</surname> <given-names>L. A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Satellite-detected ocean ecosystem response to volcanic eruptions in the subarctic Northeast Pacific Ocean</article-title>. <source>Geophys. Res. Lett.</source> <volume>46</volume>, <fpage>11270</fpage>&#x2013;<lpage>11280</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2019GL083977</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Westerling</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Hidalgo</surname> <given-names>H. G.</given-names>
</name>
<name>
<surname>Cayan</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Swetnam</surname> <given-names>T. W.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Warming and earlier spring increase Western U.S. forest wildfire activity</article-title>. <source>Science</source> <volume>313</volume> (<issue>5789</issue>), <fpage>940</fpage>&#x2013;<lpage>943</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1128834</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>