Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean sciences”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Accelerating the Lagrangian Particle Tracking in Hydrologic Modeling to Continental-Scale

Unprecedented climate change and anthropogenic activities have induced increasing ecohydrological problems, which have motivated the development of large-scale modeling for solutions. Water age/quality is as important as water quantity for understanding the water cycle. However, current scientific progress in tracking water parcels at large-scale with high spatiotemporal resolutions is far behind that in water balance/quantity owing to the lack of powerful tools. EcoSLIM is a particle tracking model that works with the hydrologic model ParFlow-CLM, which couples surface-subsurface hydrology with land surface processes. Here, we demonstrate a parallel framework to accelerate EcoSLIM to continental-scale on a distributed, multi-GPU platform with CUDA-Aware MPI. In tests from catchment-, to regional-, and then to continental-scale using 25-million to 1.6-billion particles, EcoSLIM shows significant speedup and excellent parallel performance. The parallel framework is portable to atmospheric and oceanic particle tracking models, where parallelization is inadequate and a standard parallel framework is absent. Parallelized EcoSLIM is a promising tool to accelerate our understanding of the terrestrial water cycle and the upscaling of subsurface hydrology to Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Machine learning prediction of connectivity, biodiversity and resilience in the Coral Triangle

Abstract Even optimistic climate scenarios predict catastrophic consequences for coral reef ecosystems by 2100. Understanding how reef connectivity, biodiversity and resilience are shaped by climate variability would improve chances to establish sustainable management practices. In this regard, ecoregionalization and connectivity are pivotal to designating effective marine protected areas. Here, machine learning algorithms and physical intuition are applied to sea surface temperature anomaly data over a twenty-four-year period to extract ecoregions and assess connectivity and bleaching recovery potential in the Coral Triangle and surrounding oceans. Furthermore, the impacts of the El Niño Southern Oscillation (ENSO) on biodiversity and resilience are quantified. We find that resilience is higher for reefs north of the Equator and that the extraordinary biodiversity of the Coral Triangle is dynamic in time and space, and benefits from ENSO. The large-scale exchange of genetic material is enhanced between the Indian Ocean and the Coral Triangle during La Niña years, and between the Coral Triangle and the central Pacific in neutral conditions. Through machine learning the outstanding biodiversity of the Coral Triangle, its evolution and the increase of species richness are contextualized through geological times, while offering new hope for monitoring its future.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Fire, dust, air and water: Improving aerosol biogeochemistry interactions in ACME. Final report

Aerosol-biogeochemistry interactions will have important consequences for human health, ecosystem services, and climate feedbacks over the next several decades. In this proposal we seek to link how human activities on land, including direct industrial emissions and landscape fires, impact the atmosphere and the biogeochemistry of the ocean and land. Our overall goal was to complement existing work at DOE labs to improve the capability and performance of the E3SM model in short-term, high-resolution projections of climate. Projections of ocean uptake of anthropogenic carbon are highly dependent on the details of the ocean physical and biogeochemistry parameterizations. Potentially equally important are the estimates of increases in iron and soluble nitrogen deposition in aerosols to the open ocean, but there are substantial uncertainties in these feedbacks. The E3SM land model simulates natural and anthropogenic fires, but improvement is needed to more accurately simulate fire responses to climate variability as well as fire contributions to terrestrial climate feedbacks. This proposal addressed four of the six priorities within the ACME/E3SM call by improving the representation of fire, industrial combustion sources, and desert dust and adding iron and soluble iron to the atmospheric chemistry (1). Our work focused on improving ocean biogeochemistry and the response of the oceans to anthropogenic atmospheric iron and nitrogen (3). We also improved the representation of land disturbance on fires and fire mediated carbon cycle feedbacks in the land model (4). Finally we had planned to couple the land, atmosphere and ocean biogeochemical systems (6) and explore how this coupling changes the feedbacks in the climate system. Our work under this proposal was separated into the following 5 tasks: i) improve representation of fire dynamics in both boreal and tropical ecosystems, ii) improve estimates of industrial emissions of iron and phosphorus, iii) improve representation of atmospheric iron, phosphorus and nitrogen transport, chemistry and deposition, iv) improve representation of iron, phosphorus and nitrogen deposition impacts on ocean biogeochemistry, and v) couple emissions into atmosphere, to deposition onto land and ocean biogeochemistry. For all 5 tasks, we conducted offline experiments of the historical period to enable comparison with available aerosol and land surface observations. We had planned to conduct fully coupled simulations, however, these will be conducted as part of the coupling in version 3 and 4 of the E3SM model, and thus have not been completed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Using ground-based lidar data to investigate the water–vapor budget in the daytime atmospheric boundary layer

The moisture advection term in the water–vapor budget equation is investigated with a combination of a vertically-staring water–vapor lidar and Doppler lidar systems. These instruments make it possible to get the mean profile of moisture tendency and the latent heat flux (LHF) divergence. We use data of the Land–Atmosphere Feedback Experiment (LAFE) at the Atmospheric Radiation Measurement (ARM) Program’s Southern Great Plains (SGP) site, Oklahoma, USA, collected on 30 August 2017 between 15 and 24 UTC, which corresponds to 09 to 18 LT. The lidars provide turbulence resolving profiles of moisture and vertical wind fluctuations. The LHF profile is derived from the covariance of these moisture and vertical wind fluctuations. The mean boundary layer height z i is determined from the peak of the moisture variance. The results demonstrate that the combination of two remote sensing instruments can be applied for determining the dominant water–vapor budget terms, namely moisture tendency, latent heat flux divergence and moisture advection.

Advection↗
DOE OSTI2026

Time irreversibility as an indicator of approaching tipping points in Earth subsystems

With shifting environmental trends, many Earth system elements may be poised to undergo critical transitions or ‘tipping’. Reliable anticipation of these tipping elements is vital to inform policy decisions. Many of the current methods for tipping point detection are based on loss of resilience or ‘critical slowdown’ of the system as it approaches a tipping point. However, these methods are prone to false alarms; the detected slowdown may be an artifact of nonstationary noise unrelated to tipping behavior. Here, we explore the efficacy of early warning signs based on a nonequilibrium thermodynamics framework. The model-free detection method relies on the increased intrinsic time-irreversibility due to detailed balance breaking, preceding the onset of tipping or instabilities. We demonstrate that these EWSs are effective for tipping point detection and robust against false alarms due to nonstationary noise, using idealized models for two key elements of the Earth system that are prone to tipping: the Atlantic Meridional Overturning Circulation and Arctic sea-ice loss.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Examination of Global Midlatitude Atmospheric River Lifecycles Using an Object-Oriented Methodology

Tracking atmospheric rivers (ARs) across their lifecycles is a field of recent interest with a multitude of emerging methodologies. The CONNected-objECT (CONNECT) algorithm is adapted for the tracking of global midlatitude AR lifecycles and associated precipitation by implementing a seeded region growing segmentation algorithm, creating the AR-CONNECT algorithm. To facilitate the permissiveness of the methodology, AR-CONNECT is without hard-coded geometric criteria yet is still shown to extract synoptic-scale elongated objects >99.99% of the time. One of the consequences of the methodology is the ability to occasionally track atmospheric water vapor anomalies before evolving into AR geometries, effectively tracking AR genesis further back than other studies. With the aid of subdaily satellite-derived rain data, we investigate the climatology, trends, and patterns of AR lifecycles from 1983–2016 and compare with other AR tracking studies. We find that AR frequency, genesis, and terminus locations are in generally good agreement with other AR tracking methodologies, though with key differences, and that AR frequencies in each hemisphere are determined by the number of AR hotspots. Furthermore, we uncover evidence that certain AR characteristics, such as frequency, areal extent, and duration, show evidence of increasing trends. Midlatitude precipitation uncovered by AR-CONNECT shows contributions up to 50% over land and 65% over the ocean. Trend analysis of AR precipitation shows an increase in precipitation associated with ARs propagated by the Southern Jet Stream and ARs that traverse over the Sahara Desert, among others, but is determined not to be a driver of changes in global precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Impacts of a Weakened Atlantic Meridional Overturning Circulation on ENSO in a Warmer Climate

This study quantifies the impacts of the Atlantic Meridional Overturning Circulation (AMOC) on the El Niño–Southern Oscillation (ENSO) under anthropogenic warming by comparing climate change model simulations with declining and fixed strengths of the overturning. After the 1980s, a weakened AMOC is shown to reduce the strength of the annual cycle of sea surface temperature (SST) in the eastern equatorial Pacific and induce anomalous cross-equatorial northerly winds there, which strengthens ENSO variability by about 11%. An analysis of the Bjerknes stability index reveals that this intensification of ENSO results mainly from enhanced Ekman upwelling feedback due to amplified atmospheric wind response to SST anomalies and oceanic upwelling response to equatorial wind stress anomalies. The weakened AMOC also promotes the occurrence of Central Pacific El Niño events and reduces ENSO skewness. These AMOC impacts on ENSO magnitude and complexity throughout the twenty-first century are however smaller than ENSO variations expected from internal climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

VPF-Class: taxonomic assignment and host prediction of uncultivated viruses based on viral protein families

Abstract Motivation Two key steps in the analysis of uncultured viruses recovered from metagenomes are the taxonomic classification of the viral sequences and the identification of putative host(s). Both steps rely mainly on the assignment of viral proteins to orthologs in cultivated viruses. Viral Protein Families (VPFs) can be used for the robust identification of new viral sequences in large metagenomics datasets. Despite the importance of VPF information for viral discovery, VPFs have not yet been explored for determining viral taxonomy and host targets. Results In this work, we classified the set of VPFs from the IMG/VR database and developed VPF-Class. VPF-Class is a tool that automates the taxonomic classification and host prediction of viral contigs based on the assignment of their proteins to a set of classified VPFs. Applying VPF-Class on 731K uncultivated virus contigs from the IMG/VR database, we were able to classify 363K contigs at the genus level and predict the host of over 461K contigs. In the RefSeq database, VPF-class reported an accuracy of nearly 100% to classify dsDNA, ssDNA and retroviruses, at the genus level, considering a membership ratio and a confidence score of 0.2. The accuracy in host prediction was 86.4%, also at the genus level, considering a membership ratio of 0.3 and a confidence score of 0.5. And, in the prophages dataset, the accuracy in host prediction was 86% considering a membership ratio of 0.6 and a confidence score of 0.8. Moreover, from the Global Ocean Virome dataset, over 817K viral contigs out of 1 million were classified. Availability and implementation The implementation of VPF-Class can be downloaded from https://github.com/biocom-uib/vpf-tools. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Global seamless tidal simulation using a 3D unstructured-grid model (SCHISM v5.10.0)

We present a new 3D unstructured-grid global ocean model to study both tidal and nontidal processes, with a focus on the total water elevation. Unlike existing global ocean models, the new model resolves estuaries and rivers down to ~8 m without the need for grid nesting. The model is validated with both satellite and in situ observations for elevation, temperature, and salinity. Tidal elevation solutions have a mean complex root-mean-square error (RMSE) of 4.2 cm for M2 and 5.4 cm for all five major constituents in the deep ocean. The RMSEs for the other four constituents, S2, N2, K1, and O1, are, respectively, 2.05, 0.93, 2.08, and 1.34 cm). The nontidal residual assessed by a tide gauge dataset (GESLA) has a mean RMSE of 7 cm. For the first time ever, we demonstrate the potential for seamless simulation on a single mesh from the global ocean into several estuaries along the US West Coast. The model is able to accurately capture the total elevation, even at some upstream stations. The model can therefore potentially serve as the backbone of a global tide surge and compound flooding forecasting framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Composition and Flux of Dissolved and Particulate Carbon and Nitrogen in the Lower Tocantins River

The Tocantins River contributes ~5 % of the total flux of water to the Amazon River plume in the Atlantic Ocean. Here, we evaluate monthly variability in the composition and abundance of carbon, nitrogen, and suspended sediment in the lower reaches of the Tocantins River from 2014-2016. Dissolved organic carbon (DOC) concentrations generally increased during periods of high discharge and are ~1.5 times lower than average concentrations at the mouth of the Amazon River. Dissolved inorganic carbon (DIC) similarly increased during periods of high discharge. Total dissolved nitrogen (TDN) and individual nitrogen species followed a similar temporal pattern, increasing during high water. $\text{NO}^-_3$ predominated the dissolved inorganic nitrogen (DIN) pool, followed by $\text{NH}_4^+$, and $\text{NO}_2^-$ , characteristic of environments with a relatively low anthropogenic impact. Dissolved fractions represented 92 % of the total carbon exported and 78 % of the total nitrogen exported. The suspended particulate sediment flux was 2.72 x 10 6 t yr -1 , with fine suspended sediment (FSS) dominating (71.3 %). Concentrations of carbon relative to nitrogen indicate a primarily terrigenous source of organic matter and CO 2 derived from in situ respiration of this material during the rainy season and a primarily algal/bacterial source of organic matter during the dry season. Considering past estimates of dissolved carbon and nitrogen fluxes from the Amazon River to the Atlantic Ocean, the Tocantins River contributes 3 and 3.7 % to total fluxes to the Amazon River plume region, respectively. While this contribution is relatively small, it may be influenced by future changes to the basin’s land use and hydrology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dynamics for El Niño-La Niña asymmetry constrain equatorial-Pacific warming pattern

The El Niño-Southern Oscillation (ENSO) results from the instability of and also modulates the strength of the tropical-Pacific cold tongue. While climate models reproduce observed ENSO amplitude relatively well, the majority still simulates its asymmetry between warm (El Niño) and cold (La Niña) phases very poorly. The causes of this major deficiency and consequences thereof are so far not well understood. Analysing both reanalyses and climate models, we here show that simulated ENSO asymmetry is largely proportional to subsurface nonlinear dynamical heating (NDH) along the equatorial Pacific thermocline. Most climate models suffer from too-weak NDH and too-weak linear dynamical ocean-atmosphere coupling. Nevertheless, a sizeable subset (about 1/3) having relatively realistic NDH shows that El Niño-likeness of the equatorial-Pacific warming pattern is linearly related to ENSO amplitude change in response to greenhouse warming. Therefore, better simulating the dynamics of ENSO asymmetry potentially reduces uncertainty in future projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Marginal Ice Zone as a dominant source region of atmospheric mercury during central Arctic summertime

Atmospheric gaseous elemental mercury (GEM) concentrations in the Arctic exhibit a clear summertime maximum, while the origin of this peak is still a matter of debate in the community. Based on summertime observations during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition and a modeling approach, we further investigate the sources of atmospheric Hg in the central Arctic. Simulations with a generalized additive model (GAM) show that long-range transport of anthropogenic and terrestrial Hg from lower latitudes is a minor contribution (~2%), and more than 50% of the explained GEM variability is caused by oceanic evasion. A potential source contribution function (PSCF) analysis further shows that oceanic evasion is not significant throughout the ice-covered central Arctic Ocean but mainly occurs in the Marginal Ice Zone (MIZ) due to the specific environmental conditions in that region. Our results suggest that this regional process could be the leading contributor to the observed summertime GEM maximum. In the context of rapid Arctic warming and the observed increase in width of the MIZ, oceanic Hg evasion may become more significant and strengthen the role of the central Arctic Ocean as a summertime source of atmospheric Hg.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton

ABSTRACT Metagenomics is a powerful method for interpreting the ecological roles and physiological capabilities of mixed microbial communities. Yet, many tools for processing metagenomic data are neither designed to consider eukaryotes nor are they built for an increasing amount of sequence data. EukHeist is an automated pipeline to retrieve eukaryotic and prokaryotic metagenome-assembled genomes (MAGs) from large-scale metagenomic sequence data sets. We developed the EukHeist workflow to specifically process large amounts of both metagenomic and/or metatranscriptomic sequence data in an automated and reproducible fashion. Here, we applied EukHeist to the large-size fraction data (0.8–2,000 µm) from Tara Oceans to recover both eukaryotic and prokaryotic MAGs, which we refer to as TOPAZ (Tara Oceans Particle-Associated MAGs). The TOPAZ MAGs consisted of >900 environmentally relevant eukaryotic MAGs and >4,000 bacterial and archaeal MAGs. The bacterial and archaeal TOPAZ MAGs expand upon the phylogenetic diversity of likely particle- and host-associated taxa. We use these MAGs to demonstrate an approach to infer the putative trophic mode of the recovered eukaryotic MAGs. We also identify ecological cohorts of co-occurring MAGs, which are driven by specific environmental factors and putative host-microbe associations. These data together add to a number of growing resources of environmentally relevant eukaryotic genomic information. Complementary and expanded databases of MAGs, such as those provided through scalable pipelines like EukHeist, stand to advance our understanding of eukaryotic diversity through increased coverage of genomic representatives across the tree of life. IMPORTANCE Single-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers’ efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The Impact of Parameterized Lateral Mixing on the Antarctic Circumpolar Current in a Coupled Climate Model

This study examines the impact of changing the lateral diffusion coefficient ARedi on the transport of the Antarctic Circumpolar Current (ACC). The lateral diffusion coefficient ARedi is poorly constrained, with values ranging across an order of magnitude in climate models. The ACC is difficult to accurately simulate, and there is a large spread in eastward transport in the Southern Ocean (SO) in these models. This paper examines how much of that spread can be attributed to different eddy parameterization coefficients. A coarse-resolution, fully coupled model suite was run with A Redi = 400, 800, 1200, and 2400 m 2 s -1 . Additionally, two simulations were run with two-dimensional representations of the mixing coefficient based on satellite altimetry. Relative to the 400 m 2 s -1 case, the 2400 m 2 s -1 case exhibits 1) an 11% decrease in average wind stress from 50° to 65°S, 2) a 20% decrease in zonally averaged eastward transport in the SO, and 3) a 14% weaker transport through the Drake Passage. The decrease in transport is well explained by changes in the thermal current shear, largely due to increases in ocean density occurring on the northern side of the ACC. In intermediate waters these increases are associated with changes in the formation of intermediate waters in the North Pacific. We hypothesize that the deep increases are associated with changes in the wind stress curl allowing Antarctic Bottom Water to escape and flow northward.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Large Eddy Simulations of Turbulence below Antarctic Ice Shelves

Turbulence below ice shelves is of key importance for predicting ice-shelf melt rates and consequently the contribution of ice sheets to sea level rise. In this project we conducted Large-Eddy Simulations (LES) to improve our understanding of sub-ice-shelf ocean turbulence and the relationship between ice-shelf melt rates and ocean conditions. Over the course of the second and final year of this project (FY2020), we accomplished both major code developments for the PArallel Large-eddy simulation Model (PALM; Maronga et al., 2015) and conducted a suite of simulations which form the basis for a publication in preparation.

58 GEOSCIENCES↗
DOE OSTI2025

Advancing Understanding of Deep Convective Anvil Clouds

In this project we use ground-based DOE ARM observations, satellite observations, and Weather Research and Forecasting (WRF) model simulations to arrive at a more complete understanding (or, new theory) of how atmospheric heating profiles in thunderstorm cores drives how quickly air rises in the storms. The rate at which air motion increases or decreases with altitude directly relates to how much cloud mass is ejected laterally into the adjacent atmosphere, which then influences the development and maintenance of “stratiform raining cloud anvils”. Our work has led to new equations connecting the heating in these thunderstorms to these stratiform clouds, with the connection modulated by the stability structure of the atmosphere (i.e., if more stable, then a more rapid change in thunderstorm heating with altitude results in more stratiform cloud, independent of whether over land or ocean, or how organized the thunderstorms are). Improving knowledge of stratiform clouds is important, since they contribute to substantial accumulated rainfall (important to agriculture and local hydrology) as well as radiation, given their large areal extents and lifetime. The radiation impacts likely impact simulated cloud feedback in all Earth System models (ESMs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Propagation of the 2014–2016 Northeast Pacific Marine Heatwave Through the Salish Sea

Effects and impacts of the Northeast Pacific marine heatwave of 2014–2016 on the inner coastal estuarine waters of the Salish Sea were examined using a combination of monitoring data and an established three-dimensional hydrodynamic and biogeochemical model of the region. The anomalous high temperatures reached the U.S. Pacific Northwest continental shelf toward the end of 2014 and primarily entered the Salish Sea waters through an existing strong estuarine exchange. Elevated temperatures up to + 2.3°C were observed at the monitoring stations throughout 2015 and 2016 relative to 2013 before dissipating in 2017. The hydrodynamic and biogeochemical responses to this circulating high-temperature event were examined using the Salish Sea Model over a 5-year window from 2013 to 2017. Responses of conventional water-quality indicator variables, such as temperature and salinity, nutrients and phytoplankton, zooplankton, dissolved oxygen, and pH, were evaluated relative to a baseline without the marine heatwave forcing. The simulation results relative to 2014 show an increase in biological activity (+14%, and 6% Δ phytoplankton biomass, respectively) during the peak heatwave year 2015 and 2016 propagating toward higher zooplankton biomass (+14%, +18% Δ mesozooplankton biomass). However, sensitivity tests show that this increase was a direct result of higher freshwater and associated nutrient loads accompanied by stronger estuarine exchange with the Pacific Ocean rather than warming due to the heatwave. Strong vertical circulation and mixing provided mitigation with only ≈+0.6°C domain-wide annual average temperature increase within Salish Sea, and served as a physical buffer to keep waters cooler relative to the continental shelf during the marine heatwave.

54 ENVIRONMENTAL SCIENCES↗