SEARCH · Engineering Papers
Results for “Ocean biogeochemical modeling”
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.
Amplified bottom water acidification rates on the Bering Sea shelf from 1970–2022
The Bering Sea shelf supports a highly productive marine ecosystem that is vulnerable to ocean acidification (OA) due to the cold, carbon-rich waters. Previous observational evidence suggests that bottom waters on the shelf are already seasonally undersaturated with respect to aragonite (i.e. Ω arag <1) and that OA will continue to increase the spatial extent, duration, and intensity of these conditions. Here, we use a regional ocean biogeochemical model to simulate changes in ocean carbon chemistry for the Bering Sea shelf from 1970–2022. Over this timeframe, model results suggest that surface Ω arag decreases by −0.043 per decade and surface pH by −0.014 per decade, comparable to observed global rates of OA. However, bottom water pH decreases at twice the rate of surface pH, while bottom [H + ] decreases at nearly 3 times the rate of surface [H + ]. This amplified bottom water acidification has emerged over the past 25 years and is likely driven by a combination of anthropogenic carbon accumulation and increasing primary productivity and subsurface respiration and remineralization. Due to this enhanced bottom water acidification, the spatial extent of bottom waters with Ω arag <1 has greatly expanded over the past 2 decades, along with pH conditions harmful to red king crab. Interannual variability in surface and bottom Ω arag , pH, and [H + ] has also increased over the past 2 decades, resulting in part from the increased physical climate variability. We also find that the Bering Sea shelf is a net annual carbon sink of 1.1–7.9 Tg C yr −1 , with the range resulting from the difference in the two different atmospheric forcing reanalysis products used. Seasonally, the shelf is a significant carbon sink from April–October but a somewhat weaker carbon source from November–March.
Plankton energy flows using a global size-structured and trait-based model
Plankton community models are critical tools for understanding the processes that shape marine plankton communities, how plankton communities impact biogeochemical cycles, and the feedbacks between community structure and function. Here, using the flexible Marine Biogeochemistry Library (MARBL), we present the Size-based Plankton ECological TRAits (MARBL-SPECTRA) model, which is designed to represent a diverse plankton community while remaining computationally tractable. MARBL-SPECTRA is composed of nine phytoplankton and six zooplankton size classes represented using allometric scaling relationships for physiological traits and interactions within multiple functional types. MARBL-SPECTRA is embedded within the global ocean component of the Community Earth System Model (CESM) and simulates large-scale, emergent patterns in phytoplankton growth limitation, plankton phenology, plankton generation time, and trophic transfer efficiency. The model qualitatively reproduces observed global patterns of surface nutrients, chlorophyll biomass, net primary production, and the biogeographies of a range of plankton size classes. In addition, the model simulates how predator:prey dynamics and trophic efficiency vary across gradients in total ecosystem productivity. Shorter food chains that export proportionally more carbon from the surface to the ocean interior occur in productive, eutrophic regions, whereas in oligotrophic regions, the food chains are relatively long and export less organic matter from the surface. The union of functional type modeling with size-resolved, trait-based modeling approaches allows MARBL-SPECTRA to capture both large-scale elemental cycles and the structure of planktonic food webs affecting trophic transfer efficiency.
Viroid-like “obelisk” agents are widespread in the ocean and exceed the abundance of RNA viruses in the prokaryotic fraction
Abstract “Obelisks” are recently discovered ribonucleic acid (RNA) viroid-like elements present in diverse environments with no phylogenetic similarity to any known biological agent. obelisks were first identified in the human gut and in a commensal bacterium acting as a replicative host. They have a circular ∼1 kb RNA genome, rod-like secondary structures, and the encoding of a protein superfamily called “Oblins”. We performed a large-scale search of obelisks in the ocean using the Pebblescout program and the transcriptomic Sequence Archive Read databases, revealing the biogeography and abundance of these viroid-like RNA elements. We detected 55 obelisk genomes resulting in 35 marine clusters at the species level. These obelisks were detected in the prokaryotic fraction and to a lesser extent in the eukaryotic fraction, and distributed across all the oceans from surface to mesopelagic including the Arctic, and even in the coldest seawater of Earth beneath the Antarctic Ross Ice Shelf. The obelisk hallmark protein Oblin-1 confirmed by 3D models was found in various marine samples. Some of the detected marine obelisks harbor hammerhead self-cleaving ribozymes in both polarities. In the prokaryotic, but not the eukaryotic, fraction of the Tara Ocean dataset, relative abundance of obelisks calculated by transcriptomic fragment recruitment indicated that they are abundant in marine samples, reaching or even exceeding the relative abundance of the previously discovered uncultured RNA viruses. In conclusion, obelisks are abundant and widespread viroid-like elements that should be included in ocean biogeochemical models.
Results from a multi-laboratory ocean metaproteomic intercomparison: effects of LC-MS acquisition and data analysis procedures
Metaproteomics is an increasingly popular methodology that provides information regarding the metabolic functions of specific microbial taxa and has potential for contributing to ocean ecology and biogeochemical studies. A blinded multi-laboratory intercomparison was conducted to assess comparability and reproducibility of taxonomic and functional results and their sensitivity to methodological variables. Euphotic zone samples from the Bermuda Atlantic Time-series Study (BATS) in the North Atlantic Ocean collected by in situ pumps and the autonomous underwater vehicle (AUV) Clio were distributed with a paired metagenome, and one-dimensional (1D) liquid chromatographic data-dependent acquisition mass spectrometry analysis was stipulated. Analysis of mass spectra from seven laboratories through a common bioinformatic pipeline identified a shared set of 1056 proteins from 1395 shared peptide constituents. Quantitative analyses showed good reproducibility: pairwise regressions of spectral counts between laboratories yielded R 2 values averaged 0.62±0.11, and a Sørensen similarity analysis of the top 1000 proteins revealed 70 %–80 % similarity between laboratory groups. Taxonomic and functional assignments showed good coherence between technical replicates and different laboratories. A bioinformatic intercomparison study, involving 10 laboratories using eight software packages, successfully identified thousands of peptides within the complex metaproteomic datasets, demonstrating the utility of these software tools for ocean metaproteomic research. Lessons learned and potential improvements in methods were described. Future efforts could examine reproducibility in deeper metaproteomes, examine accuracy in targeted absolute quantitation analyses, and develop standards for data output formats to improve data interoperability. Together, these results demonstrate the reproducibility of metaproteomic analyses and their suitability for microbial oceanography research, including integration into global-scale ocean surveys and ocean biogeochemical models.
Using neural network ensembles to separate ocean biogeochemical and physical drivers of phytoplankton biogeography in Earth system models
Abstract. Earth system models (ESMs) are useful tools for predicting and understanding past and future aspects of the climate system. However, the biological and physical parameters used in ESMs can have wide variations in their estimates. Even small changes in these parameters can yield unexpected results without a clear explanation of how a particular outcome was reached. The standard method for estimating ESM sensitivity is to compare spatiotemporal distributions of variables from different runs of a single ESM. However, a potential pitfall of this method is that ESM output could match observational patterns because of compensating errors. For example, if a model predicts overly weak upwelling and low nutrient concentrations, it might compensate for this by allowing phytoplankton to have a high sensitivity to nutrients. Recently, we demonstrated that neural network ensembles (NNEs) are capable of extracting relationships between predictor and target variables within ocean biogeochemical models. Being able to view the relationships between variables, along with spatiotemporal distributions, allows for a more mechanistically based examination of ESM outputs. Here, we investigated whether we could apply NNEs to help us determine why different ESMs produce different spatiotemporal distributions of phytoplankton biomass. We tested this using three cases. The first and second case used different runs of the same ESM, except that the physical circulations differed between them in the first case, while the biological equations differed between them in the second. Our results indicated that the NNEs were capable of extracting the relationships between variables for different runs of a single ESM, allowing us to distinguish between differences due to changes in circulation (which do not change relationships) from changes in biogeochemical formulation (which do change relationships). In the third case, we applied NNEs to two different ESMs. The results of the third case highlighted the capability of NNEs to contrast the apparent relationships of different ESMs and some of the challenges it presents. Although applied specifically to the ocean components of an ESM, our study demonstrates that Earth system modelers can use NNEs to separate the contributions of different components of ESMs. Specifically, this allows modelers to compare the apparent relationships across different ESMs and observational datasets.
Biogeochemical Equilibrium Responses to Maximal Productivity in High Nutrient Low Chlorophyll Regions
Whether ocean iron fertilization in high nutrient low chlorophyll (HNLC) regions is a practical geoengineering strategy to combat climate warming has been debated because the fate of absorbed carbon by the ocean remains unclear. Here, we use an optimized oceanic biogeochemical inverse model to explore the equilibrium responses to increased productivity in major HNLC regions. We find that an overall downward shift of nutrients and carbon from the surface and intermediate to the deep waters with maximal productivity in the HNLC Southern Ocean (SO), North Pacific (NP) and eastern equatorial Pacific (EP). The efficiency of soft tissue pump is increased by 37.8%, 2.7% and 5.4% in the SO, EP and NP simulation while atmospheric CO 2 concentration is drawn down by 45.1, 3.4 and 6.7 ppm, respectively. Meanwhile, global mean oxygen is decreased by 21.1%, 1.3% and 2.4% in the SO, EP and NP simulation, accompanied by an expansion of oxygen minimum zones by 400%, 26% and 63%, respectively. Our study indicates that iron fertilization has the potential to be an effective geoengineering approach to mitigate climate warming but also has a strong impact on the development of ocean hypoxia.
Ocean physical‐biogeochemical interactions in the CMIP6 and E3SM Earth System Models
The primary objectives of this project are to quantify the spatio‐temporal linkages between physical climate and ocean biogeochemical variables in the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM), in CMIP6 Earth System Models (ESMs), and in historical observations, to evaluate the E3SM’s ability to reproduce the observations, and to develop a solid diagnostic package to do so. It will revolve around four hypotheses and will focus on key physical and biogeochemical processes that regulate the upper ocean carbon and oxygen cycling: water mass distribution, ventilation, and biological production.
Implementation of a turbine farm model into the Energy Exascale Earth System Model for investigation and quantification of global climate impacts
Although there has been widespread deployment of wind farms in the United States, and plans to continue deployment into the future, the complete effects of wind farms on Earth systems are not well understood. The work performed here has incorporated wind farm models into the Energy Exascale Earth System Model (E3SM) capable of simulating the effects of extracting momentum from the atmospheric flow field using power generating wind farms. This new capability will allow scientists to quantify the impacts of wind farm induced changes on Earth systems by exploiting E3SM’s ability to couple atmospheric, oceanic, and biogeochemical (BGC) models on a global scale and monitor precipitation levels, extreme weather events, soil moisture content and jet stream location over decades-long time periods. This tool will be used to inform decision making on wind farm citing and will contribute to the Lab’s ability to assess energy technology impacts on the environment and evaluate the trade-offs between energy infrastructure investments and their impacts on natural systems.
Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models
Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O 2 ). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O 2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O 2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O 2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×10 14 [mol/decade] for OMIP1, -0.37 ± 0.15×10 14 [mol/decade] for OMIP2, and -1.06 ± 0.68×10 14 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is -0.98×10 14 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O 2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O 2 . Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O 2 changes.
Manganese Limitation of Phytoplankton Physiology and Productivity in the Southern Ocean
Although iron and light are understood to regulate the Southern Ocean biological carbon pump, observations have also indicated a possible role for manganese. Low concentrations in Southern Ocean surface waters suggest manganese limitation is possible, but its spatial extent remains poorly constrained and direct manganese limitation of the marine carbon cycle has been neglected by ocean models. Here, using available observations, we develop a new global biogeochemical model and find that phytoplankton in over half of the Southern Ocean cannot attain maximal growth rates because of manganese deficiency. Manganese limitation is most extensive in austral spring and depends on phytoplankton traits related to the size of photosynthetic antennae and the inhibition of manganese uptake by high zinc concentrations in Antarctic waters. Importantly, manganese limitation expands under the increased iron supply of past glacial periods, reducing the response of the biological carbon pump. Overall, these model experiments describe a mosaic of controls on Southern Ocean productivity that emerge from the interplay of light, iron, manganese and zinc, shaping the evolution of Antarctic phytoplankton since the opening of the Drake Passage.
Can Biogeochemical Tracer Observations Constrain Southern Ocean Diapycnal Mixing Rates?
Abstract Direct observations of background diapycnal mixing rates in the Southern Ocean (SO) are limited spatially and temporally, making the choice of an appropriate value to parameterize this mixing in Earth system models a challenge. However, the deployment of Argo floats throughout the SO has provided an extensive range of observations of both physical and biogeochemical parameters. We use an ocean state estimate run with various background diapycnal mixing coefficients to assess if biogeochemical tracer observations can be used to better constrain SO diapycnal mixing rates. We find that vertical tracer distributions in the SO are highly sensitive to the rate of background diapycnal mixing and can provide an upper limit on background mixing rates. This demonstrates the importance of biogeochemical tracer observations throughout the full depth of the water column to validate ocean models.
Evaluation of Ocean Biogeochemistry and Carbon Cycling in CMIP Earth System Models With the International Ocean Model Benchmarking (IOMB) Software System
Abstract The International Ocean Model Benchmarking (IOMB) software package is a new community resource that we use here to evaluate surface and upper ocean biogeochemical variables and integrated anthropogenic carbon uptake from earth system models (ESMs) contributing to the 5th and 6th phases of the Coupled Model Intercomparison Project (CMIP5 and CMIP6). IOMB generates graphics and tables for systematically comparing model predictions against multiple datasets. Our analysis reveals some improvement in the multi‐model mean from CMIP5 to CMIP6 for most of the variables we examined. Compared to data‐constrained estimates of ocean anthropogenic carbon uptake for the 1994–2007 period, negative biases exist for many models between 30 and 50°S. Global model estimates of anthropogenic carbon uptake for the same period do not change significantly from CMIP5 to CMIP6, with the combined ensemble mean estimate of 27.8 ± 0.5 Pg C lower than a data‐constrained estimate of 33.0 ± 4.0 Pg C. At the same time, the change in the natural carbon inventory from CMIP is estimated to be a source of 0.7 ± 0.3 Pg C, which is considerably smaller in magnitude than a data‐constrained estimate of 5.0 ± 3.0 Pg C. With chlorofluorocarbon (CFC) predictions available for several models, we demonstrate that negative anthropogenic dissolved inorganic carbon biases coincide with negative biases in CFC concentration, highlighting the importance of weak exchange between the surface and interior ocean in regulating rates of anthropogenic carbon uptake. To examine the robustness of this attribution across the CMIP models, we calculate the global vertical temperature gradient between 200 and 1,000 m as a metric for global stratification and exchange between the surface and deeper waters. We find a linear relationship between the bias of the vertical temperature gradients and the bias in global anthropogenic carbon uptake, consistent with the hypothesis that model biases in anthropogenic carbon uptake are related to biases in surface‐to‐interior exchange by physical processes.
Simulating physical and biogeochemical dynamics of the Arctic Ocean [Slides]
Ocean physics (currents, mixing, air-sea transfer of heat and gases) and biogeochemistry (carbon and oxygen cycles, ecosystem dynamics, nutrient distributions) are tightly coupled. We used the DOE-developed Energy Exascale Earth System Model (E3SM) to simulate physical and biogeochemical dynamics of the global ocean, with an emphasis on the western North Pacific. In particular, we used the grid refinement capability of E3SM’s ocean component (LANL-developed MPASOcean) to quantify the effect of model resolution on the carbon and oxygen cycles using a pair of century-long simulations that differ only in the western North Pacific grid. Regional grid refinement has a noticeable impact on biogeochemistry. Distribution of important nutrients (such as nitrate) changes the amount of photosynthesis by phytoplankton, which in turn affects air-sea transfer of CO2 and the local carbon cycle.
Assessing the fidelity of shallow-water carbonates as records of the Ni isotope composition of surface seawater
Nickel is a bioessential metal that is used in enzymes important to the C, N, and O cycles, and changes in its marine abundance and bioavailability may have affected the evolutionary trajectory of early life. Changes over time in the Ni isotope composition (δ 60 Ni) of surface seawater, which reflects biological demand for Ni, could allow for the reconstruction of the dynamics of Ni demand over Earth’s history, but this approach would require geologic records of surface seawater. Here, we investigate the fidelity of shallow-water carbonates as a record of the Ni isotope composition of surface seawater by determining how Ni is first partitioned into natural carbonates and then how post-depositional processes influence the Ni signal. Our samples come from the Great Bahama Bank, which is a well-studied, modern carbonate platform often used to study ancient platforms. We found that Ni is fractionated from seawater upon incorporation into carbonates capturing shallow (<18 cm), recent deposition (0.1 ‰–0.4 ‰ lighter than seawater). Variation among these [Ni] and δ 60 Ni values may be controlled by variation in mineral proportions. Meteoric diagenesis shifts δ 60 Ni to lower values, which we attribute to isotopically light meteoric fluids. In contrast, carbonates that experienced sediment-buffered marine diagenesis with respect to Ca isotopes and Sr/Ca ratios do not appear to differ in δ 60 Ni values from sediments generally representative of their initial deposition. The sensitivity of δ 60 Ni to diagenetic reset in these samples appears comparable to the sensitivities of Ca isotopes and Sr/Ca ratios, to first order. Thus, in general, carbonates that experienced sediment-buffered marine diagenesis with respect to these elements may hold the most promise as a record of the δ 60 Ni of coeval surface seawater. Additionally, we use our results to infer that the fraction of Ni removed from seawater into carbonates is less than 10 % of the total Ni output from the global oceans and incorporation of this Ni sink into global biogeochemical models will only have a minor impact on the modeled modern Ni budget.
Trace metal contents of autotrophic flagellates from contrasting open‐ocean ecosystems
Abstract Nano‐ and picophytoplankton are a major component of open‐ocean ecosystems and one of the main plankton functional types in biogeochemical models, yet little is known about their trace metal contents. In cultures of the picoeukaryote Ostreococcus lucimarinus , iron limitation reduced iron quotas by 68%, a fraction of the plasticity known in diatoms. In contrast, a commonly co‐occurring cyanobacterium, Prochlorococcus , showed variable iron contents with iron availability in culture. Synchrotron X‐ray fluorescence was used to measure single‐cell metal (Mn, Fe, Co, Ni, Zn) quotas of autotrophic flagellates (1.4–16.8‐ μ m diameter) collected from four ocean regions. Iron quotas were tightly constrained and showed little response to iron availability, similar to cultured Ostreococcus . Zinc quotas also did not vary with zinc availability but appeared to vary with phosphorus availability. These results suggest that macronutrient and metal availability may be equally important for controlling metal contents of small eukaryotic open‐ocean phytoplankton.
Modelling the Sensitivity of Yukon River Biogeochemical Dynamics to Environmental and Chemical Drivers: Implications for Dissolved Organic Carbon
Riverine dissolved organic carbon (DOC) is a critical biogeochemical component that transmits information from Arctic soils to the Arctic Ocean, significantly influencing carbon dynamics in this unique ecosystem. As DOC travels downstream, it undergoes transformations that alter its composition and fate. The Yukon River serves as an effective testbed for modelling these dynamics, offering sufficient scale to capture key biogeochemical processes while having a simpler hydrology than other major Arctic rivers, as well as long-term DOC observational data for model validation. To investigate DOC transformations during transit in the Yukon River, we adapted our Arctic Riverine Organic Macromolecular Model by applying regional-specific parameterisations. Our model simulates the transport and transformation of 15 organic macromolecules, including CDOM (coloured dissolved organic matter), proteins, polysaccharides, lipids, lignin phenols, and humic substances. Initial DOC concentrations were derived from observed soil organic carbon stocks in the surrounding watershed, while chemical transformations and hydrological dynamics were modelled along the river's course. Sensitivity and uncertainty analyses were conducted using a Monte Carlo approach under two experimental setups. Results revealed that variability in DOC and CDOM concentrations at the river mouth were predominantly driven by initial DOC concentration (~70% of variability explained) and dilution at confluence points (~10%). The refractory fraction of DOC explained 21%–88% of the variability in 14 macromolecular concentrations and ranked in the top five sensitive parameters for all outputs when a uniform parameter distribution was assumed. However, when a more likely variability was applied to this parameter, its influence on DOC and CDOM decreased. Given that refractory DOC accounts for ~80% of total DOC in Arctic Rivers, this suggests that most DOC resists degradation and retains its chemical composition during transport to the coastal environment. River velocity, which determines residence time, explained 8%–47% of the variability in protein, polysaccharide, lipid, pigments, and lignin phenols at the river mouth. In contrast, chemical turnover times contributed only 1%–5% to output variability. Our findings underscore the need for improved land-specific headwater observations, including seasonal soil moisture and lateral transport dynamics that control the initial tributary-specific DOC inputs. With accelerated permafrost thaw and increasing river discharge, extending our model to other Arctic River systems and seasons will enhance understanding of Arctic riverine carbon fluxes and their contributions to the Arctic Ocean.
A New Framework for the Attribution of Air‐Sea CO 2 Exchange
Abstract The air‐sea transfer of carbon dioxide can be viewed as a dynamic system through which atmospheric and oceanic processes push surface waters away from thermodynamic equilibrium, while diffusive gas transfer pulls them back toward local equilibrium. These push/pull processes drive significant sub‐seasonal, seasonal, and interannual variability in air‐sea carbon fluxes, the quantification of which is critical both for diagnosing the ocean response to fossil fuel emissions and for attempts to mitigate anthropogenic climate disruption through intentional modification of surface ocean biogeochemistry. In this study, we present a new approach for attributing air‐sea carbon fluxes to specific mechanisms. The new framework is first applied to a two‐box ocean nutrient and carbon cycle model as an illustrative example. Next, outputs from a regional eddy‐resolving model of the Southern Ocean are analyzed. The roles of multiple physical and biogeochemical processes are identified. The decomposition of the seasonal air‐sea carbon flux shows the dominant role of biological carbon pumps that are partially compensated by the transport convergence. Finally, the framework is used to diagnose the response to mesoscale iron and alkalinity release, explicitly quantifying transport feedback and eventual impacts on net air‐sea carbon flux. Ocean carbon transport has divergent influences between iron and alkalinity release, due to opposing near‐surface gradients of dissolved inorganic carbon. We suggest that our attribution framework may be a useful analytical technique for monitoring natural ocean carbon fluxes and quantifying the impacts of human intervention on the ocean carbon cycle.