Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean dynamics”

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 289 records · Page 16

Effect of Nozzle Curvature on Supersonic Gas Jets Used in Laser-Plasma Acceleration

Supersonic gas jets produced by converging-diverging (C-D) nozzles are commonly used as targets for laser-plasma acceleration (LPA) experiments. A major point of interest for these targets is the gas density at the region of interaction where the laser ionizes the gas plume to create a plasma, providing the acceleration structure. Tuning the density profiles at this interaction region is crucial to LPA optimization. A "flat-top" density profile is desired at this line of interaction to control laser propagation and high energy electron acceleration, while a short high-density profile is often preferred for acceleration of lower-energy tightly-focused laser-plasma interactions. A particular design parameter of interest is the curvature of the nozzle's diverging section. We examine three nozzle designs with different curvatures: the concave "bell", straight conical and convex "trumpet" nozzles. We demonstrate that, at mm-scale distances from the nozzle exit, the trumpet and straight nozzles, if optimized, produce "flat-top" density profiles whereas the bell nozzle creates focused regions of gas with higher densities. An optimization procedure for the trumpet nozzle is derived and compared to the straight nozzle optimization process. We find that the trumpet nozzle, by providing an extra parameter of control through its curvature, is more versatile for creating flat-top profiles and its optimization procedure is more refined compared to the straight nozzle and the straight nozzle optimization process. Furthermore, we present results for different nozzle designs from computational fluid dynamics (CFD) simulations performed with the program ANSYS Fluent and verify them experimentally using neutral density interferometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Outcomes of the DOE Workshop on Atmospheric Challenges for the Wind Energy Industry

The U.S. Department of Energy-funded Mesoscale-to-Microscale Coupling (MMC) project team planned and conducted a virtual Workshop on Atmospheric Challenges for the Wind Energy Industry on October 19 and 20, 2020. The goal of the workshop was to forge a dialog with the community, including industry representatives, on how modeling tools are currently being used, the present active atmospheric modeling research in support of wind energy, and required advancements in capabilities and technology to continue to advance wind energy deployment. The workshop was planned in collaboration with an industry advisory panel that included representatives from wind power plant developers, turbine manufacturers, and companies that provide resource assessment and forecasting services. The format of the workshop included panels from government research sponsors, visionaries from industry, and mixed panels of researchers discussing research status and needs. A shared keynote presentation from the Technical University of Denmark experts anchored the second day of the workshop. An emphasis was placed on understanding the research needs in the offshore environment. In addition, breakout opportunities were provided each day. On the first day, the breakout discussions addressed predesigned questions configured to elicit participants’ thoughts on needed research directions. The second-day breakouts treated three important technical topics through a combination of presentations and group conversations. Each workshop participant chose their breakout preference from among downscaling details, modeling for turbines, and using artificial intelligence for atmospheric modeling. The discussions were robust and productive. The outcomes of the workshop include archiving a series of recommendations from industry and the research community on research directions required to further advance wind energy deployment. Discussions confirmed the need for high-fidelity modeling but that there are specific areas of applicability and other areas where the time and cost of computation is prohibitive. In those cases, the high-fidelity models can inform low-order models that are more practical for real-time or widely deployed applications. Industry must consider the financial cost of performing more expensive modeling approaches, but industry engineers and researchers are using these approaches where there appears to be a return on investment. An emerging type of low-order model is based on machine learning (ML). Participants confirmed that there are many atmospheric phenomena that need to be modeled better, including low-level jets, cold air outbreaks, land-sea induced circulations, diurnal variability, thin stable boundary layers, dynamic changes such as from frontal passage, interaction of wakes and blockage, and more. For the offshore environment, there is wide agreement that some level of ocean-wave-atmospheric coupling is necessary to capture variations in rotor-level winds needed to plan and operate offshore wind plants. Another recurring recommendation is that more observations are needed, particularly for the offshore environment. Those observations should consider the needs for model improvement, both for physically based models and for ML models. Observations must capture atmospheric profiles of variables that are important to understanding and modeling atmospheric and oceanic phenomena that impact boundary layer winds. Models must be validated with data and the uncertainty quantified, particularly those that are sensitive to initial and boundary conditions. Finally, a repeated request was to consider the holistic needs of hybrid plants of wind, solar, and storage resources because those types of plants are likely to be the wave of the future. In addition, industry wishes to understand impacts of the resource under a changing climate for long-term planning.

17 WIND ENERGY↗

Estuarine hurricane wind can intensify surge-dominated extreme water level in shallow and converging coastal systems

Based on the projected increase in hurricane landfall frequency on the middle to lower US east coast, we examined the crucial role of the estuarine wind field in exacerbating coastal flooding. A regionally refined atmospheric and two high-resolution hydrology and ocean models are integrated to provide plausible and physically consistent ensembles of hurricane events and the associated flooding inside the Delaware Bay and River, a US mid-Atlantic estuary. Model results show that the hurricane propagation direction, estuarine geometry, remote surge from the open ocean, and direct nearshore upwind stress could magnify the flood magnitude. More specifically, inland-bound tracks that make landfall before reaching the mid-Atlantic coast produce a more significant surge within Delaware Bay than the shore-parallel tracks, where the estuarine wind direction plays the primary role in surge amplification. Ultimately, this study emphasized the need for integrated models to capture the nonlinear dynamics and interactions in flood hazard modeling.

54 ENVIRONMENTAL SCIENCES↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

Defensive symbiosis against giant viruses in amoebae

Protists are important regulators of microbial communities and key components in food webs with impact on nutrient cycling and ecosystem functioning. In turn, their activity is shaped by diverse intracellular parasites, including bacterial symbionts and viruses. Yet, bacteria–virus interactions within protists are poorly understood. Here, we studied the role of bacterial symbionts of free-living amoebae in the establishment of infections with nucleocytoplasmic large DNA viruses (Nucleocytoviricota). To investigate these interactions in a system that would also be relevant in nature, we first isolated and characterized a giant virus (Viennavirus, family Marseilleviridae) and a sympatric potential Acanthamoeba host infected with bacterial symbionts. Subsequently, coinfection experiments were carried out, using the fresh environmental isolates as well as additional amoeba laboratory strains. Employing fluorescence in situ hybridization and qPCR, we show that the bacterial symbiont, identified as Parachlamydia acanthamoebae, represses the replication of the sympatric Viennavirus in both recent environmental isolates as well as Acanthamoeba laboratory strains. In the presence of the symbiont, virions are still taken up, but viral factory maturation is inhibited, leading to survival of the amoeba host. The symbiont also suppressed the replication of the more complex Acanthamoeba polyphaga mimivirus and Tupanvirus deep ocean (Mimiviridae). Our work provides an example of an intracellular bacterial symbiont protecting a protist host against virus infections. The impact of virus–symbiont interactions on microbial population dynamics and eventually ecosystem processes requires further attention.

59 BASIC BIOLOGICAL SCIENCES↗

Soil microbial ecology and microbiome-metabolite linkages improve understanding of ecosystem states along terrestrial-aquatic interfaces

These data are from Bandopadhyay et al., "Soil microbial ecology and microbiome-metabolite linkages improve understanding of ecosystem states along terrestrial-aquatic interfaces". This study aims to understand the soil microbial ecology along terrestrial-aquatic interfaces of a freshwater and estuarine region and how it relates to organic matter. We analyzed soil microbial (16S rRNA gene) and organic matter (Fourier-transform ion cyclotron resonance mass spectrometry, FTICR-MS) composition from upland (forested), transition (stressed forest), and wetland positions at three sites in each of the Lake Erie (freshwater) and Chesapeake Bay (estuarine) regions. This dataset includes 16S rRNA gene amplicon data (only processed file types included here) and organic matter composition from FTICR-MS data (raw and processed files included here) from upland (forested), transition (stressed forest), and wetland positions at three sites in each of the Lake Erie and Chesapeake Bay regions. These sites are part of the COMPASS-FME project (https://compass.pnnl.gov/FME/COMPASSFME). File formats and software needed to access files: 16S rRNA gene amplicon data: These files follow the format reported here https://ess-dive.gitbook.io/amplicon-sequencing-reporting-format#updates-in-v1.0.1. As per this format, there are four file types reported: 1. Taxon tables (also called sequence-by-sample or OTU (operational taxonomic unit)/ESV (exact sequence variant) tables) : available in a .txt file format and accessible using TextEdit or MS Excel. 2. Representative sequences (also called consensus sequences) : available in a .fasta format and accessible using TextEdit. 3. Sequencing metadata : available in a MS Excel workbook file format and CSV file format 4. Bioinformatic metadata : available in a MS Excel workbook file format and CSV file format FTICR-MS data: 1. Raw data converted to a processed file with intensities of the peaks in the given samples : available in a MS Excel CSV file format 2. Processed file used in analyses and visualizations (appended as icr_long_) : available in a MS Excel CSV file format 3. Metadata file for ICR features (appended as icr_meta) : available in a MS Excel CSV file format

54 ENVIRONMENTAL SCIENCES↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

The black carbon cycle and its role in the Earth system

Black carbon (BC) is produced by incomplete combustion of biomass by wildfires and burning of fossil fuels. BC is environmentally persistent over centuries to millennia, sequestering carbon in marine and terrestrial environments. However, its production, storage and dynamics, and therefore its role in the broader carbon cycling during global change, are poorly understood. Here, in this Review, we discuss BC cycling across the land-to-ocean continuum. Wildfires are the main source of BC, producing 128 ± 84 teragrams per year. Negative climate–BC feedbacks could arise as wildfire increases with anthropogenic warming, producing more BC, which in turn will sequester carbon, but the magnitude of these effects are unknown. Most BC is stored in terrestrial systems with some transported to the ocean via rivers and the atmosphere. However, the oceanic BC budget is not balanced, with known BC removal fluxes exceeding BC inputs. We demonstrate these observed inconsistencies using a simple ocean box model, which highlights key areas of future research. Measurements of BC mineralization and export rates along the land-to-ocean continuum and quantification of previously unexplored sources of oceanic BC are needed to close the global BC budget.

aquatic↗

Pronounced Impact of Salinity on Rapidly Intensifying Tropical Cyclones

We show the importance of salinity for rapidly intensifying Atlantic tropical cyclones and demonstrate the potential for improved prediction of rapid intensification through the inclusion of salinity. Tropical Cyclone (TC) rapid intensification (RI) is difficult to predict and poses a formidable threat to coastal populations. A warm upper ocean is well-known to favor RI, but the role of ocean salinity is less clear. This study shows a strong inverse relationship between salinity and TC RI in the eastern Caribbean and western tropical Atlantic due to near-surface freshening from the Amazon-Orinoco River system. In this region, rapidly intensifying TCs induce a much stronger surface enthalpy flux compared to more weakly intensifying storms, in part due to a reduction in SST cooling caused by salinity stratification. This reduction has a noticeable positive impact on TCs undergoing RI, but the impact of salinity on more weakly intensifying storms is insignificant. These statistical results are confirmed through experiments with an ocean mixed layer model, which show that the salinity-induced reduction in SST cold wakes increases significantly as the storm’s intensification rate increases. Currently, operational statistical-dynamical RI models do not use salinity as a predictor. Through experiments with a statistical RI prediction scheme, it is found that the inclusion of surface salinity significantly improves the RI detection skill, offering promise for improved operational RI prediction. Satellite surface salinity may be valuable for this purpose, given its global coverage and availability in near real-time.

Balaguru, Karthik↗

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics↗

Autonomous thermal tracking reveals spatiotemporal patterns of seabird activity relevant to interactions with floating offshore wind facilities

Planning is underway for placement of infrastructure needed to begin offshore wind (OSW) energy generation along the West Coast of the United States and elsewhere in the Pacific Ocean. In contrast to the primarily nearshore windfarms currently in the North Atlantic, the seabird communities inhabiting Pacific Wind Energy Areas (WEAs) include significant populations of species that fly by dynamic soaring, a behavior dependent on wind and in which flight height increases steeply with wind speed. Therefore, a more precise and detailed assessment of their 3D airspace use is needed to better understand the potential collision risks that OSW turbines may present to these seabirds. Toward this end, a novel technology called the ThermalTracker-3D (TT3D), which uses thermal imaging and stereo vision, was developed to render high-resolution (on average within ±5 m) flight tracks and related behavior of seabirds. The technology was developed and deployed on a wind-profiling LiDAR buoy in the Humboldt WEA, located 34 to 57 km off California’s coast. During the at-sea deployment between 24 May and 13 August 2021, the TT3D successfully tracked birds moving between 10 and 500 m from the device, around the clock, and in all weather conditions; a total of 1407 detections and their corresponding 3D flight trajectories were recorded. Mean altitudes of detections ranged 6-295 m above sea level (asl). Considering the degree of overlap with anticipated rotor swept zones (RSZ), which extend 25-260 m asl, 79% of detected birds (per m 3 of airspace) moved below the RSZ, 21% moved at heights overlapping the RSZ, and another 0.04% occurred at heights exceeding the RSZ. The high-resolution tracks provided valuable insight into seabird space use, especially at heights that make them vulnerable to collision during various environmental conditions (e.g., darkness, strong winds). Observations made by the TT3D will be useful in filling critical knowledge gaps related to estimating collision and avoidance between seabirds and OSW facilities in the Pacific and elsewhere. Future research will focus on enhancing the TT3D’s identification capabilities to the lowest taxon through validation studies and artificial intelligence, further contributing to seabird conservation efforts associated with OSW.

17 WIND ENERGY↗

Outcomes and Insights From Simplified Analytic Trajectory Optimization for a Tethered Underwater Kite

This letter formulates and solves a periodic trajectory optimization problem for a tethered underwater kite. The goal is to maximize the average mechanical power harvested by the kite. The type of kite considered in this letter extracts electricity from ocean currents by moving cross-current as its reel away from its base station, and consumes electricity to reel back. The problem of optimizing this kite’s trajectory is challenging due to the high dimensionality and nonlinearity of its dynamics. To tackle this challenge, the literature often separates the problem into two subproblems focusing on optimizing the cross-current and the reel-in/reel-out components of the trajectory, respectively, which may be sub-optimal. In contrast, this letter solves for the combined cross-current and reel-in/reel-out trajectory by linearizing the dynamics of the kite around a zero-power reference equilibrium trajectory in spherical coordinates. This allows the trajectory optimization problem to be solved analytically for simple sinusoidal input perturbations from equilibrium. Here we use linear quadratic regulation to enable the nonlinear kite model to track the optimized trajectory. The result is a computationally efficient approach that achieves an attractive Loyd factor of 19.9%, while providing important insights into the nature of the optimal trajectory.

42 ENGINEERING↗

Predictability and empirical dynamics of fisheries time series in the North Pacific

Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.

60 APPLIED LIFE SCIENCES↗

The role of the tropical Atlantic in tropical Pacific climate variability

Abstract Interactions between Atlantic and Pacific Oceans can affect tropical Pacific variability and its global impacts at both interannual and decadal timescales. Thus, a deepened understanding of the coupled Atlantic-Pacific interactions is needed. While possible dynamical mechanisms by which the Atlantic can influence the Pacific have been identified, the effectiveness of those mechanisms is difficult to establish using climate model simulations where Atlantic sea surface temperatures (SSTs) are prescribed and Pacific feedbacks cannot be realistically included. As an alternative approach, here we use a Linear Inverse Model (LIM) trained on observations and capable of correctly reproducing the observed statistics, to assess the relative role of the Atlantic-to-Pacific and Pacific-to-Atlantic influences on tropical Pacific variability. Our results indicate that Atlantic internal variability can enhance interannual SST anomalies in the eastern equatorial Pacific, and decadal SST anomalies in the central equatorial Pacific, while Pacific influences on the Atlantic significantly damp tropical Pacific decadal variability. This methodological framework could also be used to assess climate model fidelity in representing tropical basin interactions, helping to reconcile existing differences among models’ results.

Meteorology & Atmospheric Sciences↗

Morning boundary layer conditions for shallow to deep convective cloud evolution during the dry season in the central Amazon

Observations of the boundary layer (BL) processes are analyzed statistically for dry seasons of 2 years and in detail, as case studies, for 4 shallow convective days (ShCu) and 4 shallow-to-deep convective days (ShDeep) using a suite of ground-based measurements from the Observation and Modeling of the Green Ocean Amazon (GoAmazon 2014/5) Experiment. The BL stages in ShDeep days, from the nighttime to the cloudy mixing layer stage, are then described in comparison with ShCu days. Atmospheric thermodynamics and dynamics, environmental profiles, and surface turbulent fluxes were employed to compare these two distinct situations for each stage of the BL evolution. Particular attention is given to the morning transition stage, in which the BL changes from stable to unstable conditions in the early morning hours. Results show that the decrease in time duration of the morning transition on ShDeep days is associated with high humidity and well-established vertical wind shear patterns. Higher humidity since nighttime not only contributes to lowering the cloud base during the rapid growth of the BL but also contributes to the balance between radiative cooling and turbulent mixing during nighttime, resulting in higher sensible heat flux in the early morning. The sensible heat flux promotes rapid growth of the well-mixed layer, thus favoring the deeper BL starting from around 08:00 LST (UTC-4 h). Under these conditions, the time duration of morning transition is used to promote convection, having an important effect on the convective BL strength and leading to the formation of shallow cumulus clouds and their subsequent evolution into deep convective clouds. Statistical analysis was used to validate the conceptual model obtained from the case studies. Despite the case-to-case variability, the statistical analyses of the processes in the BL show that the described processes are very representative of cloud evolution during the dry season.

54 ENVIRONMENTAL SCIENCES↗

The Energy Exascale Earth System Model Version 3: 2. Overview of the Coupled System

The Energy Exascale Earth System Model version 3 (E3SMv3) represents the latest advancement in Earth system modeling developed by the U.S. Department of Energy (DOE). Building upon previous versions, E3SMv3 introduces significant updates across its coupled components to enhance capability and improve fidelity. The atmosphere component incorporates advancements in chemistry, aerosol-cloud interactions, convection, and microphysics. The ocean features a new time-stepping scheme and a higher-resolution unstructured mesh with sub-ice-shelf cavities, while the sea ice model integrates advanced snow and ice physics for more realistic cryospheric simulations. The land model introduces prognostic vegetation dynamics and a new sub-grid topographic treatment of solar radiation. A new tri-grid configuration harmonizes the horizontal grids of the land and river components for improved process coupling. It is enabled by a new non-linear remapping between the atmosphere and land. E3SMv3 underwent extensive testing through a comprehensive simulation campaign, including pre-industrial control, idealized CO 2 experiments, and historical simulations spanning 1850–2024. The model demonstrates significant improvements in simulating the evolution of the historical surface temperature, particularly addressing the “pothole cooling” bias in earlier versions. Reduced aerosol-related forcing contributes to more realistic radiative forcing and better alignment with the observational record. Ocean heat content (OHC) and sea ice trends are also improved as a result.

54 ENVIRONMENTAL SCIENCES↗

An Assessment of Nonhydrostatic and Hydrostatic Dynamical Cores at Seasonal Time Scales in the Energy Exascale Earth System Model (E3SM)

Abstract In global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central‐southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

58 GEOSCIENCES↗

Biophysical model of eelgrass and water quality in Coos Bay, OR shows greater mitigation potential for ocean acidification than hypoxia

Seagrass beds provide important ecosystem services and are valued, in part, for their potential to mediate stressors such as ocean acidification and hypoxia (OAH) for sensitive species. However, the susceptibility of seagrasses to anthropogenic impacts and recent declines motivate the need to better understand the drivers of seagrass and the water quality consequences that occur with variation in seagrass abundance. To meet this need, we leveraged existing monitoring data (water quality and seagrass), hydrodynamic circulation model, and biogeochemical model framework with seagrass submodel, to produce a biophysical model of Coos Bay estuary, Oregon, U.S. The model includes biogeochemical processes involving water quality, plankton, seagrass, and sediment-water interactions. Ecosystem models like this are useful for evaluating complex estuarine systems because they allow us to extend our understanding of system dynamics beyond existing observations and perform experiments to identify the processes driving observed patterns. We used the biophysical model of Coos Bay to evaluate the dynamics of water quality and native eelgrass (Zostera marina) under three eelgrass abundance scenarios (zero eelgrass, current extent, and maximum observed extent) to elucidate the relationship between eelgrass and OAH. Including eelgrass in the Coos Bay model produced results that more closely resembled water quality observations - dissolved oxygen (DO) and pH were more dynamic in simulations with eelgrass, often having both higher highs and lower lows. While there were some areas of the estuary where DO improved with the addition of eelgrass to the model there was overall a small net increase in harmful DO conditions (based on a salmon physiological threshold). In contrast, ocean acidification conditions, pH and calcium carbonate saturation state for aragonite (Ω), were improved (based on oyster requirements) with the addition of eelgrass - although the magnitude of improvement differed seasonally and spatially. Our new model represents a useful tool - one which accounts for and controls the relevant physical and biogeochemical processes - to evaluate conditions that confer resilience or enhance vulnerability to OAH in an important Pacific Northwest coastal estuary and results can inform the OAH-related dynamics occurring in other eastern boundary current estuaries.

FVCOM-ICM↗