Carbon fluxes and primary magma CO2 contents along the global mid‐ocean ridge system
Explore the source record for details and available documents.
99 GENERAL AND MISCELLANEOUS↗
SEARCH · Search NASA
Results for “Global ocean”
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.
Evaluating the Benefits of Bayesian Hierarchical Methods for Analyzing Heterogeneous Environmental Datasets: A Case Study of Marine Organic Carbon Fluxes
Large compilations of heterogeneous environmental observations are increasingly available as public databases, allowing researchers to test hypotheses across datasets. Statistical complexities arise when analyzing compiled data due to unbalanced spatial sampling, variable environmental context, mixed measurement techniques, and other reasons. Hierarchical Bayesian modeling is increasingly used in environmental science to describe these complexities, however few studies explicitly compare the utility of hierarchical Bayesian models to simpler and more commonly applied methods. Here we demonstrate the utility of the hierarchical Bayesian approach with application to a large compiled environmental dataset consisting of 5,741 marine vertical organic carbon flux observations from 407 sampling locations spanning eight biomes across the global ocean. We fit a global scale Bayesian hierarchical model that describes the vertical profile of organic carbon flux with depth. Profile parameters within a particular biome are assumed to share a common deviation from the global mean profile. Individual station-level parameters are then modeled as deviations from the common biome-level profile. The hierarchical approach is shown to have several benefits over simpler and more common data aggregation methods. First, the hierarchical approach avoids statistical complexities introduced due to unbalanced sampling and allows for flexible incorporation of spatial heterogeneitites in model parameters. Second, the hierarchical approach uses the whole dataset simultaneously to fit the model parameters which shares information across datasets and reduces the uncertainty up to 95% in individual profiles. Third, the Bayesian approach incorporates prior scientific information about model parameters; for example, the non-negativity of chemical concentrations or mass-balance, which we apply here. We explicitly quantify each of these properties in turn. We emphasize the generality of the hierarchical Bayesian approach for diverse environmental applications and its increasing feasibility for large datasets due to recent developments in Markov Chain Monte Carlo algorithms and easy-to-use high-level software implementations.
54 ENVIRONMENTAL SCIENCES↗
Threat by marine heatwaves to adaptive large marine ecosystems in an eddy-resolving model
Marine heatwaves (MHWs), episodic periods of abnormally high sea surface temperature, severely affect marine ecosystems. Large marine ecosystems (LMEs) cover ~22% of the global ocean but account for 95% of global fisheries catches. Yet how climate change affects MHWs over LMEs remains unknown because such LMEs are confined to the coast where low-resolution climate models are known to have biases. Here, using a high-resolution Earth system model and applying a ‘future threshold’ that considers MHWs as anomalous warming above the long-term mean warming of sea surface temperatures, we find that future intensity and annual days of MHWs over the majority of the LMEs remain higher than in the present-day climate. Better resolution of ocean mesoscale eddies enables simulation of more realistic MHWs than low-resolution models. These increases in MHWs under global warming pose a serious threat to LMEs, even if resident organisms could adapt fully to the long-term mean warming.
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↗
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↗
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↗
Microbial Community Structure and Ecological Networks during Simulation of Diatom Sinking
Microbial-mediated utilization of particulate organic matter (POM) during its downward transport from the surface to the deep ocean constitutes a critical component of the global ocean carbon cycle. However, it remains unclear as to how high hydrostatic pressure (HHP) and low temperature (LT) with the sinking particles affects community structure and network interactions of the particle-attached microorganisms (PAM) and those free-living microorganisms (FLM) in the surrounding water. In this study, we investigated microbial succession and network interactions in experiments simulating POM sinking in the ocean. Diatom-derived 13C- and 12C-labeled POM were used to incubate surface water microbial communities from the East China Sea (ECS) under pressure (temperature) of 0.1 (25 °C), 20 (4 °C), and 40 (4 °C) MPa (megapascal). Our results show that the diversity and species richness of the PAM and FLM communities decreased significantly with HHP and LT. Microbial community analysis indicated an increase in the relative abundance of Bacteroidetes at high pressure (40 MPa), mostly at the expense of Gammaproteobacteria, Alphaproteobacteria, and Gracilibacteria at atmospheric pressure. Hydrostatic pressure and temperature affected lifestyle preferences between particle-attached (PA) and free-living (FL) microbes. Ecological network analysis showed that HHP and LT enhanced microbial network interactions and resulted in higher vulnerability to networks of the PAM communities and more resilience of those of the FLM communities. Most interestingly, the PAM communities occupied most of the module hubs of the networks, whereas the FLM communities mainly served as connectors of the modules, suggesting their different ecological roles of the two groups of microbes. These results provided novel insights into how HHP and LT affected microbial community dynamics, ecological networks during POM sinking, and the implications for carbon cycling in the ocean.
59 BASIC BIOLOGICAL SCIENCES↗
A Mineralogy–Based Anthropogenic Combustion–Iron Emission Inventory
Atmospheric supply of iron can modulate ocean biogeochemistry, due to its key role in global nitrogen and carbon cycles. Current estimates predict up to 20% of global ocean net primary productivity depends on an atmospheric iron source. Using a technology–based methodology, we revise total and soluble anthropogenic iron emissions and resolve iron into its mineral components, which allows modeling mineral–specific atmospheric reactions. We compare different methodologies for representing anthropogenic iron solubility: measured in mild and strong leaches and estimated using a mineralogy basis and identify the emissions that are most affected by such assumptions. The inclusion of metal smelting as an iron source increases iron emissions by up to 10 times higher in the fine aerosol fraction (smaller than 1 μm) than most previous inventories. Different solubility assumptions alter anthropogenic soluble iron emissions and deposition by a factor of 20 and 10, respectively. Using solubilities measured in mild leaches and calculated by mineralogy give 20–30 Gg/yr anthropogenic emissions and 40–50 Gg/yr deposition, while those measured in strong leaches give 80–440 Gg/yr emissions and 200–450 Gg/yr deposition. Here, this range of anthropogenic soluble iron deposition leads to global soluble iron deposition of 1,900–2,300 Gg/yr when dust, wildfires, and atmospheric processing are included, indicating such assumptions can affect global soluble iron supply by about 30%. In regions where marine primary productivity is iron limited, anthropogenic combustion–iron contributes up to half of the atmospheric soluble iron flux to the North Pacific Ocean but supplies less than 5% to the Southern Ocean.
54 ENVIRONMENTAL SCIENCES↗
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.
Great Bahama Bank carbonates↗
Impacts of Wildfire Aerosols on Global Energy Budget and Climate: The Role of Climate Feedbacks
Aerosols emitted from wildfires could significantly affect global climate through perturbing global radiation balance. In this study, the Community Earth System Model with prescribed daily fire aerosol emissions is used to investigate fire aerosols’ impacts on global climate with emphasis on the role of climate feedbacks. The total global fire aerosol radiative effect (RE) is estimated to be -0.78 ± 0.29 W m -2 , which is mostly from shortwave RE due to aerosol–cloud interactions (REaci; -0.70 ± 0.20 W m -2 ). The associated global annual-mean surface air temperature (SAT) change ΔT is -0.64 ± 0.16 K with the largest reduction in the Arctic regions where the shortwave REaci is strong. Associated with the cooling, the Arctic sea ice is increased, which acts to reamplify the Arctic cooling through a positive ice-albedo feedback. The fast response (irrelevant to ΔT) tends to decrease surface latent heat flux into atmosphere in the tropics to balance strong atmospheric fire black carbon absorption, which reduces the precipitation in the tropical land regions (southern Africa and South America). The climate feedback processes (associated with ΔT) lead to a significant surface latent heat flux reduction over global ocean areas, which could explain most (~80%) of the global precipitation reduction. The precipitation significantly decreases in deep tropical regions (5°N) but increases in the Southern Hemisphere tropical ocean, which is associated with the southward shift of the intertropical convergence zone and the weakening of Southern Hemisphere Hadley cell. Such changes could partly compensate the interhemispheric temperature asymmetry induced by boreal forest fire aerosol indirect effects, through intensifying the cross-equator atmospheric heat transport.
54 ENVIRONMENTAL SCIENCES↗
PMEL Passive Acoustics Research: Quantifying the Ocean Soundscape from Whales to Wave Energy
Passive acoustic monitoring of the global ocean has increased dramatically over the last decade, providing insights into seasonal sea ice and wind/wave variability, biodiversity, geophysical hazards, and anthropogenic noise impacts. All of these phenomena are sentinels of marine ecosystem health and ocean climate change. Recognizing the utility of underwater sound, the Pacific Marine Environmental Laboratory (PMEL) formed a passive acoustic research program with the goal of quantifying deep-ocean and coastal soundscapes in support of NOAA’s mission to conserve and manage marine ecosystems. PMEL Acoustics Program researchers have built a stable of novel ocean technologies, including autonomous stationary hydrophones, mobile platforms, and near-real-time surface buoys with satellite communication capability. These passive acoustic monitoring systems have been deployed in every major ocean basin on Earth, enabling significant advancements in understanding of natural and anthropogenic sounds. This progress includes evaluation of human-made sound levels across US waters, observations of ship noise fluctuations during the COVID-19 pandemic, and evaluation of noise levels from offshore wave-energy devices. Progress in natural sound research includes assessment of seasonal variability in the presence of endangered cetacean species due to population recovery and/or changing ocean temperatures as well as early detection of the collapse of an Antarctic ice shelf.
54 ENVIRONMENTAL SCIENCES↗
Comparing numerical accuracy and stability for different horizontal discretizations in MPAS-Ocean
This manuscript investigates the effectiveness of two possible horizontal discretizations for the global ocean model MPAS-Ocean, both applied to Spherical Centroidal Voronoi Tessellations (SCVTs). The first discretization is TRiSK, a C-grid, finite-volume method, that possesses many desirable mimetic properties, but has a low order accuracy. The second discretization was introduced for the first time by Peixoto (2016), and consists of modifications to the TRiSK scheme designed to achieve at least first-order accuracy in the L ∞ norm, with the loss of some of the mimetic properties. Tests on shallow-water and primitive-equation models show that the scheme due to Peixoto is indeed more accurate, but presents stability issues with respect to TRiSK. Here, TRiSK is indeed found to be often more stable in time and more robust with respect to errors in the geometric properties of the grid.
97 MATHEMATICS AND COMPUTING↗
Latitudinal gradient in the respiration quotient and the implications for ocean oxygen availability
Significance The loss of ocean oxygen caused by climate warming is a serious environmental issue and could lead to major declines in animal habitats. However, current Earth system models struggle in explaining observed trends in global ocean oxygen. Using a model, we demonstrate that marine oxygen concentrations depend sensitively on a rarely considered quantity known as the respiration quotient to prescribe the oxygen consumed per mole of organic carbon respired. Using a combination of direct chemical measurements across major biomes in the eastern Pacific Ocean plus changes in the global distribution of oxygen and carbon, we show how the respiration quotient increases systematically with temperature. The findings imply that we can experience more intense declines in ocean oxygen with warming.
54 ENVIRONMENTAL SCIENCES↗
LIGHT-bgcArgo-1.0: using synthetic float capabilities in E3SMv2 to assess spatiotemporal variability in ocean physics and biogeochemistry
Since their advent over 2 decades ago, autonomous Argo floats have revolutionized the field of oceanography, and, more recently, the addition of biogeochemical and biological sensors to these floats has greatly improved our understanding of carbon, nutrient, and oxygen cycling in the ocean. While Argo floats offer unprecedented horizontal, vertical, and temporal coverage of the global ocean, uncertainties remain about whether Argo sampling frequency and density capture the true spatiotemporal variability in physical, biogeochemical, and biological properties. As the true distributions of, e.g., temperature or oxygen are unknown, these uncertainties remain difficult to address with Argo floats alone. Numerical models with synthetic observing systems offer one potential avenue to address these uncertainties. Here, we implement synthetic biogeochemical Argo floats into the Energy Exascale Earth System Model version 2 (E3SMv2), which build on the Lagrangian In Situ Global High-Performance Particle Tracking (LIGHT) module in E3SMv2 (E3SMv2-LIGHT-bgcArgo-1.0). Since the synthetic floats sample the model fields at model run time, the end user defines the sampling protocol ahead of any model simulation, including the number and distribution of synthetic floats to be deployed, their sampling frequency, and the prognostic or diagnostic model fields to be sampled. Using a 6-year proof-of-concept simulation, we illustrate the utility of the synthetic floats in different case studies. In particular, we quantify the impact of (i) sampling density on the float-derived detection of deep-ocean change in temperature or oxygen and on float-derived estimates of phytoplankton phenology, (ii) sampling frequency and sea-ice cover on float trajectory lengths and hence float-derived estimates of current velocities, and (iii) short-term variability in ecosystem stressors on estimates of their seasonal variability.
54 ENVIRONMENTAL SCIENCES↗
A Gradient Based Subgrid-Scale Parameterization for Ocean Mesoscale Eddies
Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large-scale ocean flows. In this work, we introduce a new subgrid-scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high-resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse-resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non-dimensional parameters in our new parameterization through a non-ad hoc and tuning-free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy-permitting and non-eddying simulations.
54 ENVIRONMENTAL SCIENCES↗
Storm Surge Modeling as an Application of Local Time-Stepping in MPAS-Ocean
This paper presents the first practical application of local time-stepping (LTS) schemes in the Model for Prediction Across Scales-Ocean (MPAS-O). We use LTS schemes in a single-layer, global ocean model that predicts the storm surge around the eastern coast of the United States during Hurricane Sandy. The variable-resolution meshes used are of unprecedentedly high resolution in MPAS-O, containing cells as small as 125 m wide in Delaware Bay. It is shown that a particular, third-order LTS scheme (LTS3) produces sea-surface height solutions that are of comparable quality to solutions produced by the classical four-stage, fourth-order Runge-Kutta method (RK4) with a uniform time step on the same meshes. Furthermore, LTS3 is up to 35% faster in the best cases considered, where the number of cells using the coarse time-step relative to those using the fine time-step is as low as 1:1. This shows that LTS schemes are viable for use in MPAS-O with the added benefit of substantially less computational cost. The results of these performance experiments inform us of the requirements for efficient mesh design and configuration of LTS regions for LTS schemes. In particular, we see that for LTS to be efficient on a given mesh, it is important to have enough cells using the coarse time-step relative to those using the fine time-step, typically at least 1:5 to see an increase in performance.
54 ENVIRONMENTAL SCIENCES↗
Observationally constrained analysis of sulfur cycle in the marine atmosphere with NASA ATom measurements and AeroCom model simulations
The atmospheric sulfur cycle plays a key role in air quality, climate, and ecosystems, such as pollution, radiative forcing, new particle formation, and acid rain. In this study, we compare the spatially and temporally resolved measurements from the NASA Atmospheric Tomography (ATom) mission with simulations from five AeroCom III models for four sulfur species (dimethyl sulfide (DMS), sulfur dioxide (SO 2 ), particulate methanesulfonate (MSA), and particulate sulfate (SO 4 )). We focus on remote regions over the Pacific, Atlantic, and Southern oceans from near the surface to ~12 km altitude range covering all four seasons. In general, the differences among model results can be greater than 1 order of magnitude. Comparing with observations, model-simulated SO 2 is generally low, whereas SO 4 is generally high. Simulated DMS concentrations near the sea surface exceed observed levels by a factor of 5 in most cases, suggesting potential overestimation of DMS emissions in all models. With GEOS model simulations of tagging emission from anthropogenic, biomass burning, volcanic, and oceanic sources, we find that anthropogenic emissions are the dominant source of sulfate aerosol (40 %–60 % of the total amount) in the ATom measurements at almost all altitudes, followed by volcanic emissions (18 %–32 %) and oceanic sources (16 %–32 %). Similar source contributions can also be derived at broad ocean basins and on monthly scales, indicating the representativeness of ATom measurements for global ocean. Our work presents the first assessment of AeroCom sulfur study using ATom measurements, providing directions for improving sulfate simulations, which remain the largest uncertainty in radiative forcing estimates in aerosol climate models.
54 ENVIRONMENTAL SCIENCES↗
Underestimation of multi-decadal global O 2 loss due to an optimal interpolation method
Abstract. The global ocean's oxygen content has declined significantly over the past several decades and is expected to continue decreasing under global warming, with far-reaching impacts on marine ecosystems and biogeochemical cycling. Determining the oxygen trend, its spatial pattern, and uncertainties from observations is fundamental to our understanding of the changing ocean environment. This study uses a suite of CMIP6 Earth system models to evaluate the biases and uncertainties in oxygen distribution and trends due to sampling sparseness. Model outputs are sub-sampled according to the spatial and temporal distribution of the historical shipboard measurements, and the data gaps are filled by a simple optimal interpolation method using Gaussian covariance with a constant e-folding length scale. Sub-sampled results are compared to full model output, revealing the biases in global and basin-wise oxygen content trends. The simple optimal interpolation underestimates the modeled global deoxygenation trends, capturing approximately two-thirds of the full model trends. The North Atlantic and subpolar North Pacific are relatively well sampled, and the simple optimal interpolation is capable of reconstructing more than 80 % of the oxygen trend in the non-eddying CMIP models. In contrast, pronounced biases are found in the equatorial oceans and the Southern Ocean, where the sampling density is relatively low. The application of the simple optimal interpolation method to the historical dataset estimated the global oxygen loss to be 1.5 % over the past 50 years. However, the ratio of the global oxygen trend between the sub-sampled and full model output has increased the estimated loss rate in the range of 1.7 % to 3.1 % over the past 50 years, which partially overlaps with previous studies. The approach taken in this study can provide a framework for the intercomparison of different statistical gap-filling methods to estimate oxygen content trends and their uncertainties due to sampling sparseness.
Ito, Takamitsu (ORCID:000000019873099X)↗