Search NASA⌕ Search

SEARCH · Search NASA

Results for “Air simulations”

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 73 records · Page 4

Comparing Cetane Number Measurement Methods

Cetane number is one of the most important fuel performance metrics for mixing controlled compression-ignition “diesel” engines, quantifying a fuel’s propensity for autoignition when injected into end-of-compression-type temperature and pressure conditions. The historical default and referee method on a Cooperative Fuel Research (CFR) engine configured with indirect fuel injection and variable compression ratio is cetane number (CN) rating. A subject fuel is evaluated against primary reference fuel blends, with heptamethylnonane defining a low-reactivity endpoint of CN = 15 and hexadecane defining a high-reactivity endpoint of CN = 100. While the CN scale covers the range from zero (0) to 100, typical testing is in the range of 30 to 65 CN. Alternatively, several constant-volume combustion chamber (CVCC)-based cetane rating devices have been developed to rate fuels with an equivalent derived cetane number (DCN) or indicated cetane number (ICN). These devices measure ignition delay for fuel injected into a fixed volume of high-temperature and high-pressure air to simulate end-of-compression-type conditions. In this study, a range of novel fuel compounds are evaluated across three CVCC methods: the Ignition Quality Tester (IQT), Fuel Ignition Tester (FIT), and Advanced Fuel Ignition Delay Analyzer (AFIDA). Resulting DCNs and ICNs are compared for fuels within the normal diesel fuel range of reactivity, as well as very high (~100) and very low DCNs/ICNs (~5). Distinct differences between results from various devices are discussed. This is important to consider because some new, high-efficiency advanced compression-ignition (CI) engine combustion strategies operate with more kinetically controlled distributed combustion as opposed to mixing controlled diffusion flames. These advanced combustion strategies may benefit from new fuel chemistries, but current rating methods of CN, DCN, and ICN may not fully describe their performance. In addition, recent evidence suggests ignition delay in modern on-road diesel engines with high-pressure common rail fuel injection systems may no longer directly correlate to traditional CN fuel ratings. Simulated end-of-compression conditions are compared for CN, DCN, and ICN and discussed in the context of modern diesel engines to provide additional insight. Results highlight the potential need for revised and/or multiple fuel test conditions to measure fuel performance for advanced CI strategies.

47 OTHER INSTRUMENTATION↗

Direct measurement of the muonic content of extensive air showers between 2 $\times$ 10 17 and 2 $\times$ 10 18 eV at the Pierre Auger Observatory

The hybrid design of the Pierre Auger Observatory allows for the measurement of the properties of extensive air showers initiated by ultra-high energy cosmic rays with unprecedented precision. By using an array of prototype underground muon detectors, we have performed the first direct measurement, by the Auger Collaboration, of the muon content of air showers between 2 $\times$ 10 17 and 2 $\times$ 10 18 eV. We have studied the energy evolution of the attenuation-corrected muon density, and compared it to predictions from air shower simulations. The observed densities are found to be larger than those predicted by models. We quantify this discrepancy by combining the measurements from the muon detector with those from the Auger fluorescence detector at 10 17.5 eV and 10 18 eV. We find that, for the models to explain the data, an increase in the muon density of 38% $±$4%(12%) $±^{21\%}_{18\%}$ for EPOS-LHC, and of 50%(53%) $±$4%(13%) $±^{23\%}_{20\%}$ for QGSJETII-04, is respectively needed.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating Total Methane Emissions from the Denver-Julesburg Basin Using Bottom-Up Approaches

Methane is a powerful greenhouse gas with a 25 times higher 100-year warming potential than carbon dioxide and is a target for mitigation to achieve climate goals. To control and curb methane emissions, estimates are required from the sources and sectors which are typically generated using bottom-up methods. However, recent studies have shown that national and international bottom-up approaches can significantly underestimate emissions. In this study, we present three bottom-up approaches used to estimate methane emissions from all emission sectors in the Denver-Julesburg basin, CO, USA. Our data show emissions generated from all three methods are lower than historic measurements. A Tier 1/2 approach using IPCC emission factors estimated 2022 methane emissions of 358 Gg (0.8% of produced methane lost by the energy sector), while a Tier 3 EPA-based approach estimated emissions of 269 Gg (0.2%). Using emission factors informed by contemporary and region-specific measurement studies, emissions of 212 Gg (0.2%) were calculated. The largest difference in emissions estimates were a result of using the Mechanistic Air Emissions Simulator (MAES) for the production and transport of oil and gas in the DJ basin. The MAES accounts for changes to regulatory practice in the DJ basin, which include comprehensive requirements for compressors, pneumatics, equipment leaks, and fugitive emissions, which were implemented to reduce emissions starting in 2014. The measurement revealed that normalized gas loss is predicted to have been reduced by a factor of 20 when compared to 10-year-old normalization loss measurements and a factor of 10 less than a nearby oil and production area (Delaware basin, TX); however, we suggest that more measurements should be made to ensure that the long-tail emission distribution has been captured by the modeling. This study suggests that regulations implemented by the Colorado Department of Public Health and Environment could have reduced emissions by a factor of 20, but contemporary regional measurements should be made to ensure these bottom-up calculations are realistic.

03 NATURAL GAS↗

Simulation-Based Analysis of Impacts of Reduced Envelope and Duct Air Leakage on Indoor Air Pollutant Concentrations in Occupied Manufactured Homes

This report presents a simulation-based analysis that estimates the magnitude of indoor air pollutant concentration changes that could result from the proposed energy conservation standards. The analysis was conducted for a 1568 ft 2 , double-wide MH with variations in heating and cooling equipment (either furnace + air conditioner or heat pump) and various types of whole-house mechanical ventilation (continuous exhaust fan or central fan integrated supply), in three locations with varying climate conditions: Chicago IL, Fresno, CA; and Houston TX. The simulations tracked four air pollutants that can reach levels exceeding established safe target levels in homes: acrolein, formaldehyde, fine particulate matter (PM 2.5 ), and nitrogen dioxide (NO 2 ). The simulations considered acrolein and formaldehyde emitted from continuous indoor sources; acrolein, NO 2 and PM 2.5 from cooking; PM 2.5 from dispersed occupant activities; and NO 2 and PM 2.5 from outdoors, using historical data to identify typical levels. The impacts in homes operating or not operating whole-house mechanical ventilation equipment, kitchen and bath exhaust fans, and window opening as ventilation approaches were examined.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Refining seasonal performance metrics for room air-conditioning in emerging markets: Integrating building simulations with real-world equipment performance data

Buildings significantly impact worldwide energy consumption, emphasizing the need to reduce the cooling energy demand, especially in warm climates. Minimum energy performance standards (MEPS) and seasonal performance metrics such as the Cooling Seasonal Performance Factor (CSPF) are crucial for improving room air conditioning (RAC) efficiency. However, challenges remain, particularly in emerging markets like Brazil, where seasonal performance metrics have recently been introduced. This study assesses the factors influencing country-level seasonal efficiency metrics and proposes a framework to refine these calculations by considering local climates and expected RAC usage in real-world households via building simulations. Key considerations include outdoor air temperature binning for different climates, RAC usage patterns (i.e., daytime and nighttime usages), envelope thermal performance of households, and urban heat island (UHI) effects. The results reveal that CSPF values can vary significantly based on climate conditions, with observed CSPF ranging from 4.10 to 11.59 Wh/Wh across 577 Brazilian climates. The inclusion of UHI effects led to a reduction in CSPF values by up to 29% during nighttime operations in hot urban areas. Additionally, building envelope efficiency showed contrasting impacts on RAC performance, with CSPFs reaching up to 15.35 Wh/Wh under specific optimized conditions. These findings highlight the need for transparent policymaking in RAC performance databases, facilitating the application of approaches like those proposed in this study and supporting diverse stakeholders in decision-making.

Bavaresco, Mateus↗

Density scaling approximation for Monte-Carlo simulations of radioactive plumes

The release of radioactive gas into the atmosphere can diffuse into large volumes of air downwind from the point of release. The extent of radioactivity can cover thousands of cubic meters of air. For such large volumes, the weather models used to predict the down-wind distribution of the plume and the radiation transport models used to predict the radiation reaching ground-level from the plume can take tens of hours of computer time on multi-node institutional High-Performance Computing facilities. In this paper we focus on the radiation transport aspect of plume modeling. Here, we describe a phenomenological method for approximating the amounts of radiation that reach ground level from large volumes of a static radioactive plume that can be calculated on a stand-alone personal computer in much shorter computation times than those usually needed for such large volume evaluations. We refer to this method as the Density Scaling Approximation (DSA). Its ability to approximate ground-level count rates of large plumes comes from using a small-plume volume with a scaled-up value of air density to simulate the same number of scatterings that occur during transport in larger plume volumes at normal air density. We demonstrate the DSA by using a 100 m-diameter air-filled hemispherical dome geometry with a uniform volumetric activity of 135 Xe gas throughout the air-filled volume. The DSA for a larger dome diameter is obtained by evaluating the 100 m dome with an air density scaled up by the linear ratio of the larger diameter to the 100 m diameter. We find that this approximation works well for dome diameters up to 1200 m – the largest diameter studied and a size more than sufficient for accounting for all the radiation from 135 Xe. Moreover, most of our DSA results can be calculated over 500 times faster than corresponding full-sized geometry with normal air density. To help evaluate the accuracy of the DSA and gain insight into how well it can reproduce different regions of the spectra, we use three, easily understood regions of interest to compare the DSA results to the full-sized geometry at normal air density results. These regions are the full-energy peak, the region of single-Compton scattering, and the region of multiple-Compton scattering. We show how the dominance of the Compton scattering mechanism determines this division and thus provides insight into how Compton scattering is manifested in spectra from photon scattering through air in general, and how well the DSA approximation works.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ocean Energy Sandia TEAMER WEC Simulation Results

Computational fluid simulations for wave energy converters and supporting materials from Ocean Energy's WEC Buoy TEAMER project in collaboration with Sandia National Laboratories. Each file includes images and video of simulation along with the simulation files for the case. The videos and images include air and water simulations of the buoy in 2D and 3D views. Data spreadsheets of the simulation outputs are also included in the files. The data are represented as different cases in which the wave height(Hs) and peak period(Tp) were different.

16 TIDAL AND WAVE POWER↗

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗

Orientation and Conformation of Proteins at the Air–Water Interface Determined from Integrative Molecular Dynamics Simulations and Sum Frequency Generation Spectroscopy

Understanding the assembly of proteins at the air-water interface (AWI) informs the formation of protein films, emulsion properties, and protein aggregation. Determination of protein conformation and orientation at an interface is difficult to resolve with a single experimental or simulation technique alone. To date, the interfacial structure of even one of the most widely studied proteins, lysozyme, at the AWI remains unresolved. Here in this study, molecular dynamics (MD) simulations are used to determine if the protein adopts a side-on, head-on, or axial orientation at the AWI with two different forcefields, GROMOS-53a6 + SPC/E and a99SB-disp + TIP4P-D. Vibrational sum frequency generation (SFG) spectroscopy experiments and spectral SFG calculations validate consistency between the structure determined from MD and experiments. Overall, we show with strong agreement that lysozyme adopts an axial conformation at pH 7. Further, we provide molecular-level insight as to how pH influences the binding domains of lysozyme resulting in side-on adsorption near the isoelectric point of the lysozyme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Satellite and ARM Observations to Evaluate Cold Air Outbreak Cloud Transitions in E3SM Global Storm‐Resolving Simulations

Abstract This study examines marine boundary layer cloud regime transition during a cold air outbreak (CAO) over the Norwegian Sea, simulated by a global storm‐resolving model (GSRM) known as the Simple Cloud‐Resolving Energy Exascale Earth System Model Atmosphere Model (SCREAM). By selecting observational references based on a combination of large‐scale conditions rather than strict time‐matched comparisons, this study finds that SCREAM qualitatively captures the CAO cloud transition, including boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also accurately locates the greatest ice and liquid in the mesoscale updrafts, however, underestimates supercooled liquid water in cumulus clouds. The model evaluation approach adopted by this study takes advantages of the existing computational‐expensive global simulations of GSRM and the available observations to understand model performance and can be applied to assessments of other cloud regimes in different regions. Such practice provides valuable guidance on the future effort to correct and improve biased model behaviors.

54 ENVIRONMENTAL SCIENCES↗

Low global warming potential (GWP) refrigerant supermarket refrigeration system modeling and its application

As an environmentally friendly low global warming potential (GWP) refrigerant, Carbon dioxide (CO 2 ) has continuously gained popularity and research attention as alternative refrigerant for supermarket refrigeration system. In this paper, to fulfill the increasing need of accurate Low-GWP supermarket refrigeration models for development of supervisory level control and optimization strategies, a high fidelity model is developed for CO 2 transcritical supermarket refrigeration system which includes compressor rack of low temperature (LT) compressors and medium temperature (MT) compressors, air-cooled gas cooler, evaporator, expansion valves and other auxiliary equipment. A resistance-capacity model structure is proposed to simulate the display cases. Semi-thermodynamic models are proposed to estimate reciprocating compressors volumetric efficiency and power consumption. Furthermore the zone modeling approach is used for evaporator simulation, and air-cooled gas cooler is modeled with distributed modeling method. The expansion valve simulation is based on orifice flow model. To calibrate these models, both manufacture performance data and experimental data are used. The experiments are conducted with a full instrumental CO 2 supermarket refrigeration system installed in Oak Ridge National Lab Flexible Research Platform (FRP). The simulation model can predict the system performance, including power consumption, cooling capacity, mass flow rate, temperature, and pressure, with high accuracy (within ±4%) compared to experimental data. In addition, this developed model has been used to create the system optimum high side pressure for high side expansion valve control, and to generate wide operating range simulation data for developing the batter-equivalent commercial refrigeration system model which can be used in grid interactive control development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Flexible and Generic Functional Mock-up Unit Based Threat Injection Framework for Grid-interactive Efficient Buildings: A Case Study in Modelica

Grid-interactive efficient buildings (GEBs) have been considered as an important asset to support the power grid reliability by utilizing the demand flexibility offered by GEBs. GEBs are enabled by advances in sensors and controls, and the communication between building equipment, whole buildings, and the grid. The integration of different building technologies and network-based communication system makes GEBs vulnerable to passive threats such as equipment failure and active threats such as cyber-attacks. Modeling and simulation is an effective way to evaluate the impact of threats on the system performance. This paper proposes a generic and flexible threat injection framework for commonly-used building energy simulators such as EnergyPlus and Modelica to support threat modeling and evaluation. This framework leverages functional mock-up unit (FMU) to develop a general modeling interface for threat injection and simulation. A numerical case study using Modelica as a building energy simulator is conducted to demonstrate the capability of the framework for supporting single/multiple-order threat modeling and simulation of a GEB. Four threats and their combinations are injected on a Modelica-based threat-free building energy and control system, including operating supply fan at its full speed, remotely cycling the chiller on and off, blocking the chiller from receiving the chilled water supply temperature setpoints, and hijacking the global zone air temperature setpoint. Simulation results show that the cyber-attack that leads to short-term signal blocking has small effects on the system operation due to the "self-healing" feature of the heating, ventilation, and air-conditioning (HVAC) interactive control system. The threat that takes control of resetting the global zone air temperature setpoints has the most adverse impact on the system energy use, peak power demand, thermal comfort and the provision of demand flexibility. The combination of four threats have aggregative effects on the system but the effects are less than the additive effects of the individual threat.

Fu, Yanyang↗

Significant impact of forcing uncertainty in a large ensemble of climate model simulations

Forcing due to solar and volcanic variability, on the natural side, and greenhouse gas and aerosol emissions, on the anthropogenic side, are the main inputs to climate models. Reliable climate model simulations of past and future climate change depend crucially upon them. Here we analyze large ensembles of simulations using a comprehensive Earth System Model to quantify uncertainties in global climate change attributable to differences in prescribed forcings. The different forcings considered here are those used in the two most recent phases of the Coupled Model Intercomparison Project (CMIP), namely CMIP5 and CMIP6. We show significant differences in simulated global surface air temperature due to volcanic aerosol forcing in the second half of the 19th century and in the early 21st century. The latter arise from small-to-moderate eruptions incorporated in CMIP6 simulations but not in CMIP5 simulations. We also find significant differences in global surface air temperature and Arctic sea ice area due to anthropogenic aerosol forcing in the second half of the 20th century and early 21st century. These differences are as large as those obtained in different versions of an Earth System Model employing identical forcings. In simulations from 2015 to 2100, we find significant differences in the rates of projected global warming arising from CMIP5 and CMIP6 concentration pathways that differ slightly but are equivalent in terms of their nominal radiative forcing levels in 2100. Our results highlight the influence of assumptions about natural and anthropogenic aerosol loadings on carbon budgets, the likelihood of meeting Paris targets, and the equivalence of future forcing scenarios.

54 ENVIRONMENTAL SCIENCES↗

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗

Analysis into Asymptotic Convergence to Full Nonlinear Solutions and Exploration of the Implication of Numerical Operator Mutation of Differential Systems

A robust, sufficiently accurate and practical hydrodynamic simulation toolset is required as a key component of the modeling and simulation of air-gap electrostatic discharge events. This work was performed to complement these ongoing efforts. In particular, hydrodynamic simulations must be vetted to ensure they are robust and sufficiently accurate over relevant characteristic scales. Verification models were generated in order to cultivate the technical knowledge and expertise needed to properly create, implement and execute numerical simulations. Furthermore, this effort was utilized extensively to educate students on the mathematical and numerical principles underlying hydrodynamic simulations. This education opportunity, provided in a holistic and rigorous manner, has greatly benefited developing scientists and engineers with the necessary understandings and toolsets required to excel at accomplishing the task at hand, and, more generally, it has enabled them to generate key programmatic deliverables. This report articulates several subtilties; specifically, how perturbations, nonlinear behavior, and dissipative mechanisms influence numerical stability, how to properly structure mathematical and numerical solutions, and how to properly generate error estimation/assignment. A more rigorous discussion of the consequences of such topics can be found in the body of this report in Chapters 2 and 3 with qualitative findings discussed in Chapter 4.

97 MATHEMATICS AND COMPUTING↗

The Importance of Convective Stage on Out-of-Cloud Convectively Induced Turbulence from High-Resolution Simulations

Turbulence (clear-air, mountain wave, convectively induced) is an aviation hazard that is a challenge to forecast due to the coarse resolution ultilized in operational weather models. Turbulence indices are commonly used to aid pilots in avoiding turbulence, but these indices have been designed and calibrated for midlatitude clear-air turbulence prediction (e.g., the Ellrod index). A significant limitation with current convectively induced turbulence (CIT) prediction is the lack of storm stage dependency. In this study, six high-resolution simulations of tropical oceanic and midlatitude continental convection are performed to characterize the turbulent environment near various convective types during the developing and mature stages. Second-order structure functions, a diagnostic commonly used to identify turbulence in turbulence prediction systems, are used to characterize the probability of turbulence for various convective types. Turbulence likelihood was found to be independent of region (i.e., tropical vs midlatitude) but dependent on convective stage. The probability of turbulence increased near developing convection for the majority of cases. Additional analysis of static stability and vertical wind shear, indicators of turbulence potential, showed that the convective environment near developing convection was more favorable for turbulence production than mature convection. Near developing convection, static stability decreased and vertical wind shear increased. Vertical wind shear near mature and developing convection was found to be weakly correlated to turbulence intensity in both the tropics and the midlatitudes. Finally, this study emphasizes the need for turbulence avoidance guidelines for the aviation community that are dependent on convective stage.

54 ENVIRONMENTAL SCIENCES↗

Improved Simulation of East Asian Summer Monsoon in the High-resolution CESM1 and Its Causes

Based on the high- and low-resolution Community Earth System Model, version 1 (CESM1), and corresponding simulations from phase 6 of the Coupled Model Intercomparison Project (CMIP6), we compare the interannual variability of the East Asian summer monsoon (EASM). The EASM interannual variability is characterized by the anomalous western North Pacific anticyclone (WNPAC) circulation and the dipole rainfall pattern with a negative southern lobe over the western North Pacific and a positive northern lobe along the Meiyu–Baiu region, which is better reproduced by the high-resolution models. Here, the reason for the improvement in the high-resolution models has been attributed to the better simulation of the warm temperature advection from the wind anomalies on the climatological temperature gradient. Positive sea surface temperature (SST) anomalies over the tropical Indian Ocean are the key to the improved wind anomalies featuring a WNPAC in the high-resolution models. The warm SST anomalies over the tropical Indian Ocean strengthen the WNPAC by triggering a Kelvin-wave response to the enhanced heat release induced by the increased precipitation. Based on the mixed-layer heat budget analysis, the warm SST anomalies over the western Indian Ocean in the high-resolution CESM1 are tied to the anomalous easterly wind along the equator, which reduces surface evaporation and upwelling. Therefore, the better simulations of air-sea feedback and the oceanic mesoscale eddy over the western Indian Ocean are the key for the improved simulation of the EASM interannual variations in the high-resolution CESM1.

54 ENVIRONMENTAL SCIENCES↗

Simulations With the Marine Biogeochemistry Library (MARBL)

Abstract The Marine Biogeochemistry Library (MARBL) is a prognostic ocean biogeochemistry model that simulates marine ecosystem dynamics and the coupled cycles of carbon, nitrogen, phosphorus, iron, silicon, and oxygen. MARBL is a component of the Community Earth System Model (CESM); it supports flexible ecosystem configuration of multiple phytoplankton and zooplankton functional types; it is also portable, designed to interface with multiple ocean circulation models. Here, we present scientific documentation of MARBL, describe its configuration in CESM2 experiments included in the Coupled Model Intercomparison Project version 6 (CMIP6), and evaluate its performance against a number of observational data sets. The model simulates present‐day air‐sea CO 2 flux and many aspects of the carbon cycle in good agreement with observations. However, the simulated integrated uptake of anthropogenic CO 2 is weak, which we link to poor thermocline ventilation, a feature evident in simulated chlorofluorocarbon distributions. This also contributes to larger‐than‐observed oxygen minimum zones. Moreover, radiocarbon distributions show that the simulated circulation in the deep North Pacific is extremely sluggish, yielding extensive oxygen depletion and nutrient trapping at depth. Surface macronutrient biases are generally positive at low latitudes and negative at high latitudes. CESM2 simulates globally integrated net primary production (NPP) of 48 Pg C yr −1 and particulate export flux at 100 m of 7.1 Pg C yr −1 . The impacts of climate change include an increase in globally integrated NPP, but substantial declines in the North Atlantic. Particulate export is projected to decline globally, attributable to decreasing export efficiency associated with changes in phytoplankton community composition.

Meteorology & Atmospheric Sciences↗