Extending aquatic spectral information with the first radiometric IR-B field observations.

Planetary radiometric observations enable remote sensing of biogeochemical parameters to describe spatiotemporal variability in aquatic ecosystems. For approximately the last half century, the science of aquatic radiometry has established a knowledge base using primarily, but not exclusively, visible wavelengths. Scientific subdisciplines supporting aquatic radiometry have evolved hardware, software, and procedures to maximize competency for exploiting visible wavelength information. This perspective culminates with the science requirement that visible spectral resolution must be continually increased to extract more information. Other sources of information, meanwhile, remain underexploited, particularly information from nonvisible wavelengths. Herein, absolute radiometry is used to evaluate spectral limits for deriving and exploiting aquatic data products, specifically the normalized water-leaving radiance, Γ(λ), and its derivative products. Radiometric observations presented herein are quality assured for individual wavebands, and spectral verification is conducted by analyzing celestial radiometric results, comparing agreement of above- and in-water observations at applicable wavelengths, and evaluating consistency with bio-optical models and optical theory. The results presented include the first absolute radiometric field observations of Γ(λ) within the IR-B spectral domain (i.e. spanning 1400-3000 nm), which indicate that IR-B signals confer greater and more variable flux than formerly ascribed. Black-pixel processing, a routine correction in satellite and in situ aquatic radiometry wherein a spectrum is offset corrected relative to a nonvisible waveband (often IR-B or a shorter legacy waveband) set to a null value, is shown to degrade aquatic spectra and derived biogeochemical parameters.

Kelpwatch: A new visualization and analysis tool to explore kelp canopy dynamics reveals variable response to and recovery from marine heatwaves.

Giant kelp and bull kelp forests are increasingly at risk from marine heatwave events, herbivore outbreaks, and the loss or alterations in the behavior of key herbivore predators. The dynamic floating canopy of these kelps is well-suited to study via satellite imagery, which provides high temporal and spatial resolution data of floating kelp canopy across the western United States and Mexico. However, the size and complexity of the satellite image dataset has made ecological analysis difficult for scientists and managers. To increase accessibility of this rich dataset, we created Kelpwatch, a web-based visualization and analysis tool. This tool allows researchers and managers to quantify kelp forest change in response to disturbances, assess historical trends, and allow for effective and actionable kelp forest management. Here, we demonstrate how Kelpwatch can be used to analyze long-term trends in kelp canopy across regions, quantify spatial variability in the response to and recovery from the 2014 to 2016 marine heatwave events, and provide a local analysis of kelp canopy status around the Monterey Peninsula, California. We found that 18.6% of regional sites displayed a significant trend in kelp canopy area over the past 38 years and that there was a latitudinal response to heatwave events for each kelp species. The recovery from heatwave events was more variable across space, with some local areas like Bahía Tortugas in Baja California Sur showing high recovery while kelp canopies around the Monterey Peninsula continued a slow decline and patchy recovery compared to the rest of the Central California region. Kelpwatch provides near real time spatial data and analysis support and makes complex earth observation data actionable for scientists and managers, which can help identify areas for research, monitoring, and management efforts.

Spectrally simplified approach for leveraging legacy geostationary oceanic observations.

The use of multispectral geostationary satellites to study aquatic ecosystems improves the temporal frequency of observations and mitigates cloud obstruction, but no operational capability presently exists for the coastal and inland waters of the United States. The Advanced Baseline Imager (ABI) on the current iteration of the Geostationary Operational Environmental Satellites, termed the R Series (GOES-R), however, provides sub-hourly imagery and the opportunity to overcome this deficit and to leverage a large repository of existing GOES-R aquatic observations. The fulfillment of this opportunity is assessed herein using a spectrally simplified, two-channel aquatic algorithm consistent with ABI wave bands to estimate the diffuse attenuation coefficient for photosynthetically available radiation, Kd(PAR). First, an in situ ABI dataset was synthesized using a globally representative dataset of above- and in-water radiometric data products. Values of Kd(PAR) were estimated by fitting the ratio of the shortest and longest visible wave bands from the in situ ABI dataset to coincident, in situKd(PAR) data products. The algorithm was evaluated based on an iterative cross-validation analysis in which 80% of the dataset was randomly partitioned for fitting and the remaining 20% was used for validation. The iteration producing the median coefficient of determination (R2) value (0.88) resulted in a root mean square difference of 0.319m-1, or 8.5% of the range in the validation dataset. Second, coincident mid-day images of central and southern California from ABI and from the Moderate Resolution Imaging Spectroradiometer (MODIS) were compared using Google Earth Engine (GEE). GEE default ABI reflectance values were adjusted based on a near infrared signal. Matchups between the ABI and MODIS imagery indicated similar spatial variability (R2=0.60) between ABI adjusted blue-to-red reflectance ratio values and MODIS default diffuse attenuation coefficient for spectral downward irradiance at 490 nm, Kd(490), values. This work demonstrates that if an operational capability to provide ABI aquatic data products was realized, the spectral configuration of ABI would potentially support a sub-hourly, visible aquatic data product that is applicable to water-mass tracing and physical oceanography research.

Automated satellite remote sensing of giant kelp at the Falkland Islands (Islas Malvinas).

Giant kelp populations that support productive and diverse coastal ecosystems at temperate and subpolar latitudes of both hemispheres are vulnerable to changing climate conditions as well as direct human impacts. Observations of giant kelp forests are spatially and temporally uneven, with disproportionate coverage in the northern hemisphere, despite the size and comparable density of southern hemisphere kelp forests. Satellite imagery enables the mapping of existing and historical giant kelp populations in understudied regions, but automating the detection of giant kelp using satellite imagery requires approaches that are robust to the optical complexity of the shallow, nearshore environment. We present and compare two approaches for automating the detection of giant kelp in satellite datasets: one based on crowd sourcing of satellite imagery classifications and another based on a decision tree paired with a spectral unmixing algorithm (automated using Google Earth Engine). Both approaches are applied to satellite imagery (Landsat) of the Falkland Islands or Islas Malvinas (FLK), an archipelago in the southern Atlantic Ocean that supports expansive giant kelp ecosystems. The performance of each method is evaluated by comparing the automated classifications with a subset of expert-annotated imagery (8 images spanning the majority of our continuous timeseries, cumulatively covering over 2,700 km of coastline, and including all relevant sensors). Using the remote sensing approaches evaluated herein, we present the first continuous timeseries of giant kelp observations in the FLK region using Landsat imagery spanning over three decades. We do not detect evidence of long-term change in the FLK region, although we observe a recent decline in total canopy area from 2017-2021. Using a nitrate model based on nearby ocean state measurements obtained from ships and incorporating satellite sea surface temperature products, we find that the area of giant kelp forests in the FLK region is positively correlated with the nitrate content observed during the prior year. Our results indicate that giant kelp classifications using citizen science are approximately consistent with classifications based on a state-of-the-art automated spectral approach. Despite differences in accuracy and sensitivity, both approaches find high interannual variability that impedes the detection of potential long-term changes in giant kelp canopy area, although recent canopy area declines are notable and should continue to be monitored carefully.

One- and Two-Band Sensors and Algorithms to Derive aCDOM(440) from Global Above- and In-Water Optical Observations.

The colored (or chromophoric, depending on the literature) dissolved organic matter (CDOM) spectral absorption coefficient, aCDOM(λ), is a variable of global interest that has broad application in the study of biogeochemical processes. Within the funding for scientific research, there is an overarching trend towards increasing the scale of observations both temporally and spatially, while simultaneously reducing the cost per sample, driving a systemic shift towards autonomous sensors and observations. Legacy aCDOM(λ) measurement techniques can be cost-prohibitive and do not lend themselves toward autonomous systems. Spectrally rich datasets carefully collected with advanced optical systems in diverse locations that span a global range of water bodies, in conjunction with appropriate quality assurance and processing, allow for the analysis of methods and algorithms to estimate aCDOM(440) from spectrally constrained one- and two-band subsets of the data. The resulting algorithms were evaluated with respect to established fit-for-purpose criteria as well as quality assured archival data. Existing and proposed optical sensors capable of exploiting the algorithms and intended for autonomous platforms are identified and discussed. One-band in-water algorithms and two-band above-water algorithms showed the most promise for practical use (accuracy of 3.0% and 6.5%, respectively), with the latter demonstrated for an airborne dataset.

Large-scale shift in the structure of a kelp forest ecosystem co-occurs with an epizootic and marine heatwave.

Climate change is responsible for increased frequency, intensity, and duration of extreme events, such as marine heatwaves (MHWs). Within eastern boundary current systems, MHWs have profound impacts on temperature-nutrient dynamics that drive primary productivity. Bull kelp (Nereocystis luetkeana) forests, a vital nearshore habitat, experienced unprecedented losses along 350 km of coastline in northern California beginning in 2014 and continuing through 2019. These losses have had devastating consequences to northern California communities, economies, and fisheries. Using a suite of in situ and satellite-derived data, we demonstrate that the abrupt ecosystem shift initiated by a multi-year MHW was preceded by declines in keystone predator population densities. We show strong evidence that northern California kelp forests, while temporally dynamic, were historically resilient to fluctuating environmental conditions, even in the absence of key top predators, but that a series of coupled environmental and biological shifts between 2014 and 2016 resulted in the formation of a persistent, altered ecosystem state with low primary productivity. Based on our findings, we recommend the implementation of ecosystem-based and adaptive management strategies, such as (1) monitoring the status of key ecosystem attributes: kelp distribution and abundance, and densities of sea urchins and their predators, (2) developing management responses to threshold levels of these attributes, and (3) creating quantitative restoration suitability indices for informing kelp restoration efforts.

Chlorophyll absorption and phytoplankton size information inferred from hyperspectral particulate beam attenuation.

Electromagnetic theory predicts spectral dependencies in extinction efficiency near a narrow absorption band for a particle with an index of refraction close to that of the medium in which it is immersed. These absorption band effects are anticipated in oceanographic beam-attenuation (beam-c) spectra, primarily due to the narrow red peak in absorption produced by the phytoplankton photopigment, chlorophyll a (Chl a). Here we present a method to obtain Chl a absorption and size information by analyzing an eigendecomposition of hyperspectral beam-c residuals measured in marine surface waters by an automatic underway system. We find that three principal modes capture more than 99% of the variance in beam-c residuals at wavelengths near the Chl a red absorption peak. The spectral shapes of the eigenvectors resemble extinction efficiency residuals attributed to the absorption band effects. Projection of the eigenvectors onto the beam-c residuals produces a time series of amplitude functions with absolute values that are strongly correlated to concurrent Chl a absorption line height (aLH) measurements (r values of 0.59 to 0.83) and hence provide a method to estimate Chl a absorption. Multiple linear regression of aLH on the amplitude functions enables an independent estimate of aLH, with RMSE of 3.19⋅10-3m-1 (3.3%) or log10-RMSE of 18.6%, and a raw-scale R2 value of 0.90 based on the Tara Oceans Expedition data. Relationships between the amplitude functions and the beam-c exponential slopes are in agreement with theory relating beam-c to the particle size distribution. Compared to multispectral analysis of beam-c slope, hyperspectral analysis of absorption band effects is anticipated to be relatively insensitive to the addition of nonpigmented particles and to monodispersion.