Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled models of the climate system”

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 343 records · Page 19

The Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICA)

To explore the various couplings across space, time and between ecosystems in a consistent manner, atmospheric modeling is moving away from the fractured limited-scale modeling strategy of the past towards a unification of the range of scales inherent in the Earth System. This paper describes the forward-looking MUlti-Scale Infrastructure for Chemistry and Aerosols (MUSICA), which is intended to become the next generation community infrastructure for research involving atmospheric chemistry and aerosols. MUSICA will be developed collaboratively by the National Center for Atmospheric Research (NCAR) and university and government researchers, with the goal of serving the international research and applications communities. The capability of unifying various spatio-temporal scales, coupling to other Earth System components and process-level modularization will allow advances in both fundamental and applied research in atmospheric composition, air quality and climate and is also envisioned to become a platform that addresses the needs of policy makers and stakeholders.

Pfister, Gabriele G.↗

High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) towards CMIP7

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. Consequently, the role of enhanced horizontal resolution in improved process representation in all components of the climate system continues to be of great interest. Recent simulations suggest the possibility of significant changes in both large-scale aspects of the ocean and atmospheric circulations and in the regional responses to climate change, as well as improvements in representations of small-scale processes and extremes, when resolution is enhanced. The first phase of the High-Resolution Model Intercomparison Project (HighResMIP1) was successful at producing a baseline multi-model assessment of global simulations with model grid spacings of 25–50 km in the atmosphere and 10–25 km in the ocean, a significant increase when compared to models with standard resolutions on the order of 1° that are typically used as part of the Coupled Model Intercomparison Project (CMIP) experiments. In addition to over 250 peer-reviewed manuscripts using the published HighResMIP1 datasets, the results were widely cited in the Intergovernmental Panel on Climate Change report and were the basis of a variety of derived datasets, including tracked cyclones (both tropical and extratropical), river discharge, storm surge, and impact studies. There were also suggestions from the few ocean eddy-rich coupled simulations that aspects of climate variability and change might be significantly influenced by improved process representation in such models. The compromises that HighResMIP1 made should now be revisited, given the recent major advances in modelling and computing resources. Aspects that will be reconsidered include experimental design and simulation length, complexity, and resolution. In addition, larger ensemble sizes and a wider range of future scenarios would enhance the applicability of HighResMIP. Therefore, we propose the High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) to improve and extend the previous work, to address new science questions, and to further advance our understanding of the role of horizontal resolution (and hence process representation) in state-of-the-art climate simulations. With further increases in high-performance computing resources and modelling advances, along with the ability to take full advantage of these computational resources, an enhanced investigation of the drivers and consequences of variability and change in both large- and synoptic-scale weather and climate is now possible. With the arrival of global cloud-resolving models (currently run for relatively short timescales), there is also an opportunity to improve links between such models and more traditional CMIP models, with HighResMIP providing a bridge to link understanding between these domains. HighResMIP also aims to link to other CMIP projects and international efforts such as the World Climate Research Program lighthouse activities and various digital twin initiatives. It also has the potential to be used as training and validation data for the fast-evolving machine learning climate models.

54 ENVIRONMENTAL SCIENCES↗

Understanding AMOC stability: the North Atlantic Hosing Model Intercomparison Project

Abstract. The Atlantic meridional overturning circulation (AMOC) is an important part of our climate system. The AMOC is predicted to weaken under climate change; however, theories suggest that it may have a tipping point beyond which recovery is difficult, hence showing quasi-irreversibility (hysteresis). Although hysteresis has been seen in simple models, it has been difficult to demonstrate in comprehensive global climate models. Here, we outline a set of experiments designed to explore AMOC hysteresis and sensitivity to additional freshwater input as part of the North Atlantic Hosing Model Intercomparison Project (NAHosMIP). These experiments include adding additional freshwater (hosing) for a fixed length of time to examine the rate and mechanisms of AMOC weakening and whether the AMOC subsequently recovers once hosing stops. Initial results are shown from eight climate models participating in the Sixth Coupled Model Intercomparison Project (CMIP6). The AMOC weakens in all models as a result of the freshening, but once the freshening ceases, the AMOC recovers in half of the models, and in the other half it stays in a weakened state. The difference in model behaviour cannot be explained by the ocean model resolution or type nor by details of subgrid-scale parameterisations. Likewise, it cannot be explained by previously proposed properties of the mean climate state such as the strength of the salinity advection feedback. Instead, the AMOC recovery is determined by the climate state reached when hosing stops, with those experiments where the AMOC is weakest not experiencing a recovery.

54 ENVIRONMENTAL SCIENCES↗

Changes in the Regional Water Cycle and Their Impact on Societies

ABSTRACT Changes in “blue water”, which is the total supply of fresh water available for human extraction over land, are quite closely related to changes in runoff or equivalently precipitation minus evaporation, . This article examines how climate change‐driven recent past and future changes in the regional water cycle relate to blue water availability and changes in human blue water demand. Although at the largest scales theoretical and numerical model predictions are in broad agreement with observations, at continental scales and below models predict large ranges of possible future and runoff especially at the scale of individual river catchments and for shorter timescale subseasonal floods and droughts. Nevertheless, it is expected that the occurrence and severity of floods will increase and that of droughts may increase, possibly compounded by human‐driven non‐climatic changes such as changes in land use, dam water impoundment, irrigation and extraction of groundwater. Contemporary assessments predict that increases in 21st century human water extraction in many highly‐populated regions are unlikely to be sustainable given projections of future . To reduce uncertainty in future predictions, there is an urgent need to improve modeling of atmospheric, land surface and human processes and how these components are coupled. This should be supported by maintaining the observing network and expanding it to improve measurements of land surface, oceanic and atmospheric variables. This includes the development of satellite observations stable over multiple decades and suitable for building reanalysis datasets appropriate for model evaluation.

54 ENVIRONMENTAL SCIENCES↗

Urban cells: Extending the energy hub concept to facilitate sector and spatial coupling

The rapid growth of urban areas and concerns over climate change make it vital to improve the energy sustainability of cities. Understanding the complex interactions within different sectors (sectoral) and localities (spatial) of cities plays a crucial role in improving efficiency and sustainability, which is extremely challenging due to the complex urban morphology. State-of-the-art energy concepts do not facilitate a detailed consideration of both sectoral and spatial coupling that energy infrastructure maintains at the urban scale. This has become a significant challenge when designing interconnected urban energy infrastructure. The Urban Cell concept is introduced to address this bottleneck. A novel computational model using a modular approach is introduced to create an interconnected urban infrastructure, including the energy, building, and transportation sectors. Optimal sizing of the distributed energy system (including renewables, energy storage, and dispatchable sources) and optimal urban morphology is determined within a modular unit. A game-theoretic approach is used to model the interactions between urban cells (modular units). The study revealed that the urban cell concept can reduce the net present value of the interconnected energy infrastructure by 37% while increasing the installed renewable energy capacity by 25%. This demonstrates the benefit potential of urban cells and the importance of considering interactions between different sectors and different parts within a city. The Urban Cell concept can be used to present the complex interactions maintained within a city.

Perera, ATD↗

ARM Aerial Instrument Workshop Report

The mission of the U.S. Department of Energy’s (DOE) Biological and Environmental Research (BER) program is to “support transformative science and scientific user facilities to achieve a predictive understanding of complex biological, earth, and environmental systems for energy and infrastructure security, independence, and prosperity.” (https://science.osti.gov/ber) Aligned with the BER central mission, the Earth and Environmental Systems Sciences Division (EESSD) plays a vital role in supporting the fundamental research to understand and predict Earth’s climate and environmental systems, and is also in a unique position to inform the development of sustainable solutions to the nation’s energy and environmental challenges. Specifically, EESSD manages two scientific user facilities: the Atmospheric Radiation Measurement (ARM) user facility and the Environmental Molecular Sciences Laboratory (EMSL). These facilities provide the broader scientific community with scientific expertise, technical capabilities, and unique data sets to facilitate science in areas of importance to DOE. As a multi-platform scientific user facility, ARM aims to fulfill the needs predominantly within the EESSD Atmospheric System Research (ASR) and the Earth and Environmental System Modeling (EESM) mission areas, and provide the critical measurements required to improve understanding of aerosol and cloud life cycles and their interactions, and their coupling with the Earth’s surface. Over the years, ARM has carried out piloted and unmanned aircraft campaigns under different organizational and operational paradigms (Schmid et al. 2014, 2016). Building on its success, the ARM Aerial Facility (AAF) continues to complement the ground-based observations with airborne in situ cloud, aerosol, and trace gas observations as well as measurements of atmospheric state and atmospheric radiation. During the past three years, ARM has managed field campaigns using unmanned aerial systems (UAS) and tethered balloon systems (TBS) at Oliktok Point in Alaska to improve understanding of atmospheric processes in the Arctic. In 2019, following a careful evaluation of scientific community needs, ARM acquired a Bombardier Challenger 850 regional jet to replace the vintage Grumman Gulfstream-159 turboprop aircraft previously used by AAF. With this new “laboratory in the sky”, AAF is evaluating its current and future aerial observation capabilities to continue satisfying the needs of the research community.

54 ENVIRONMENTAL SCIENCES↗

A Positive Low Cloud–Sea Surface Temperature Feedback in the East Asian Marginal Seas during El Niño Mature Winters and Their Following Spring

During the winter and subsequent spring of an El Niño year, the East Asian marginal sea (EAMS) exhibits positive sea surface temperature anomalies (SSTAs) and fewer low clouds, while the western North Pacific experiences negative SSTAs. In this study, we suggest that the positive SSTAs in EAMS are maintained by a positive low cloud–SST feedback. In neutral winters and springs, the EAMS is covered by low clouds, which have a cooling effect on surface temperatures. During an El Niño year, a western North Pacific anomalous anticyclone is established, and along its northwestern flank, there are favorable conditions for convergence of moisture and weaker surface latent heat flux over the EAMS. Once a positive SSTA has been established, a further reduction of turbulent mixing results in less low cloud and enhanced solar heating of the ocean mixed layer; this reinforces and maintains both the positive SSTA and the lack of low cloud via a positive feedback mechanism. Further, the concurrent increase of low cloud–SST feedback and anticyclone circulation strengths is evident in the coupled-model simulations from phase 6 of the Coupled Model Intercomparison Project. Furthermore, sensitivity experiments, performed with the atmospheric components of Community Earth System Model (CESM2), reveal that a positive SSTA helps to maintain the western North Pacific anomalous anticyclone. Four pacemaker-coupled experiments by CESM2, with sea surface temperature in the equatorial Pacific restored to the observational anomalies plus the model climatology and altered low cloud feedback over EAMS, suggest that the low cloud–SST feedback results in more than the maintenance of a positive SSTA over the EAMS: the positive feedback is also a previously overlooked mechanism for the maintenance of the western North Pacific anomalous anticyclone.

54 ENVIRONMENTAL SCIENCES↗

Application of Mobile Energy Storage for Enhancing Power Grid Resilience: A Review

Natural disasters can lead to large-scale power outages, affecting critical infrastructure and causing social and economic damages. These events are exacerbated by climate change, which increases their frequency and magnitude. Improving power grid resilience can help mitigate the damages caused by these events. Mobile energy storage systems, classified as truck-mounted or towable battery storage systems, have recently been considered to enhance distribution grid resilience by providing localized support to critical loads during an outage. Compared to stationary batteries and other energy storage systems, their mobility provides operational flexibility to support geographically dispersed loads across an outage area. This paper provides a comprehensive and critical review of academic literature on mobile energy storage for power system resilience enhancement. As mobile energy storage is often coupled with mobile emergency generators or electric buses, those technologies are also considered in the review. Allocation of these resources for power grid resilience enhancement requires modeling of both the transportation system constraints and the power grid operational constraints. These aspects are discussed, along with a discussion on the cost–benefit analysis of mobile energy resources. The paper concludes by presenting research gaps, associated challenges, and potential future directions to address these challenges.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimal Gas-Electric Energy System Decarbonization Planning

As energy utilities implement climate change mitigation policies, system planners require strategies for achieving affordable emissions reductions. Coordinated planning of electric power and natural gas systems will allow synergistic investments to address cross-sector operational constraints, competing uses for net-zero emissions fuels, and shifts in energy demands across energy carriers. In this study, we develop a novel optimization program that finds the cost-minimizing mix of infrastructure expansion or reduction across gas and electric systems to satisfy sector-specific emissions constraints. Alongside energy supply resources, our framework allows for central-planning of end-use equipment stocks to allow switching between gas and electric appliances upon failure or premature replacement. The proposed model is used to simulate case study scenarios for a benchmark 24-pipe gas network coupled to a 24-node power system test network. We find that electrification of greater than 80% of core gas demands is a component of the least-cost solution for modeled energy systems. Despite this substitution, the gas system is maintained to service difficult-to-electrify customers and to deliver net-zero emissions gas to electricity generators in times of peak electricity demand. Restricting electrification of gas appliances increases reliance on power- to-gas technologies and increases annual costs by 15% in 2040. Neglecting constraints on pipeline blending of hydrogen can produce a misleading result that relies on hydrogen blend fractions of greater than 50%. In all cases, we find the average costs of delivered gas increase nearly 5-fold across the decarbonization transition, highlighting the importance of future work investigating cost-allocation strategies for ensuring an equitable energy transition.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Using Convolutional Neural Network to Emulate Seasonal Tropical Cyclone Activity

Abstract It has been widely recognized that tropical cyclone (TC) genesis requires favorable large‐scale environmental conditions. Based on these linkages, numerous efforts have been made to establish an empirical relationship between seasonal TC activities and large‐scale environmental favorability in a quantitative way, which lead to conceptual functions such as the TC genesis index. However, due to the limited amount of reliable TC observations and complexity of the climate system, a simple analytic function may not be an accurate portrait of the empirical relationship between TCs and their ambiences. In this research, we use convolution neural networks (CNNs) to disentangle this complex relationship. To circumvent the limited amount of seasonal TC observation records, we implement transfer‐learning technique to train ensemble of CNNs first on suites of high‐resolution climate model simulations with realistic seasonal TC activities and large‐scale environmental conditions, and then on a state‐of‐the‐art reanalysis from 1950 to 2019. The trained CNNs can well reproduce the historical TC records and yields significant seasonal prediction skills when the large‐scale environmental inputs are provided by operational climate forecasts. Furthermore, by inputting the ensemble CNNs with 20th century reanalysis products and Phase 6 of the Coupled Model Intercomparison Project (CMIP6) simulations, we investigated TC variability and its changes in the past and future climates. Specifically, our ensemble CNNs project a decreasing trend of global mean TC activity in the future warming scenario, which is consistent with our future projections using high‐resolution climate model.

Meteorology & Atmospheric Sciences↗

Observational benchmarks inform representation of soil organic carbon dynamics in land surface models

Abstract. Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon–climate feedbacks. Machine learning models can help identify dominant environmental controllers and establish their functional relationships with SOC stocks. The resulting knowledge can be integrated into ESMs to reduce uncertainty and improve predictions of SOC dynamics over space and time. In this study, we used a large number of SOC field observations (n=54 000), geospatial datasets of environmental factors (n=46), and two machine learning approaches (namely random forest, RF, and generalized additive modeling, GAM) to (1) identify dominant environmental controllers of global and biome-specific SOC stocks, (2) derive functional relationships between environmental controllers and SOC stocks, and (3) compare the identified environmental controllers and predictive relationships with those in models used in Phase 6 of the Coupled Model Intercomparison Project (CMIP6). Our results showed that the diurnal temperature, drought index, cation exchange capacity, and precipitation were important observed environmental predictors of global SOC stocks. While the RF model identified 14 environmental factors that describe climatic, vegetation, and edaphic conditions as important predictors of global SOC stocks (R2=0.61, RMSE = 0.46 kg m−2), current ESMs oversimplify the relationships between environmental factors and SOC, with precipitation, temperature, and net primary productivity explaining > 96 % of the variability in ESM-modeled SOC stocks. Further, our study revealed notable disparities among the functional relationships between environmental factors and SOC stocks simulated by ESMs compared with observed relationships. To improve SOC representations in ESMs, it is imperative to incorporate additional environmental controls, such as the cation exchange capacity, and refine the functional relationships to align more closely with observations.

54 ENVIRONMENTAL SCIENCES↗

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↗

Multi-decadal variation of ENSO forecast skill since the late 1800s

Diagnosing El Niño-Southern Oscillation (ENSO) predictability within operational forecast models is hindered by computational expense and the need for initialization with three-dimensional fields generated by global data assimilation. We instead examine multi-year ENSO predictability since the late 1800s using the model-analog technique, which has neither limitation. We first draw global coupled model states from pre-industrial control simulations, from the Coupled Model Intercomparison Project Phase 6, that are chosen to initially match observed monthly sea surface temperature and height anomalies in the Tropics. Their subsequent 36-month model evolution are the hindcasts, whose 20th century ENSO skill is comparable to twice-yearly hindcasts generated by a state-of-the-art European operational forecasting system. Despite the so-called spring predictability barrier, present throughout the record, there is substantial second-year ENSO skill, especially after 1960. Overall, ENSO exhibited notably high values of both amplitude and skill towards the end of the 19th century, and again in recent decades.

54 ENVIRONMENTAL SCIENCES↗

Aerosol Impacts on Mesoscale Convective Systems Forming Under Different Vertical Wind Shear Conditions

Following our previous study of wind shear effect on mesoscale convective system (MCS) organization under a clean atmospheric condition using the Weather Research and Forecasting model coupled with Spectral-Bin Microphysics, we conduct sensitivity simulations by increasing cloud condensation nuclei (CCN) concentration and investigate aerosol impacts on MCSs forming under different wind shear conditions. We find that, aerosols induce stronger updrafts and downdrafts in all MCSs. The stronger updrafts and enlarged convective core area contribute to larger vertical mass fluxes and enhance precipitation, indicating convective invigoration. Increased updraft speed below 8-km altitude mainly results from enhanced condensational heating, and increased water loading effect is not large enough to overcome the increase in thermal buoyancy. Interestingly, above 8-km altitude we see reduced updraft speed by the increased aerosols due to reduced vertical pressure perturbation gradient force, and weakened depositional growth due to competition of water vapor by enhanced condensation. The accumulated rainfall and mean rain rate are increased, with an increased occurrence frequency of heavy rain. The increase of rain rate is seen in both convective and stratiform regions, with a more significant increase in convective rain rate. In general, we see a higher frequency of deep convective cloud in the polluted condition because of invigorated convection, and more stratiform/anvil clouds, but a lower frequency of shallow warm clouds. The unanimously invigorated MCSs by aerosols under various wind shear conditions revealed by this study has an important implication in weather and climate at warm and humid regions that are influenced by pollution.

54 ENVIRONMENTAL SCIENCES↗

CMIP7 Data Request: atmosphere priorities and opportunities

This paper presents a comprehensive overview of the Coupled Model Intercomparison Project Phase 7 (CMIP7) request for data unlocking key research avenues in atmospheric science and provides justification for the resources needed to produce this data. Topics within the CMIP7 Atmosphere Theme centre around processes and feedbacks in atmospheric science such as clouds, aerosols and atmospheric chemistry, atmospheric circulation, temperature variability and extremes, radiative forcings, and Earth system model evaluation. These topics are summarised in this paper as scientific “opportunities” which will be realised through CMIP7 experiments and Earth system model outputs. These opportunities were submitted by a thematic group of atmospheric science community representatives combined with an extended consultation process. The production of these variables will close key gaps and uncertainties identified during previous rounds of CMIP, and will be broadly used by scientific, policy, governmental, industry, and other communities that rely on climate model projections for research and decision making, including supporting the 7th Intergovernmental Panel on Climate Change Assessment Report (AR7). As an author group, we also reflect on the process used to collate this data request and make recommendations to future CMIP governance on implementing a consultation on this scale in the future.

58 GEOSCIENCES↗

A Causal Inference Model Based on Random Forests to Identify the Effect of Soil Moisture on Precipitation

Soil moisture influences precipitation mainly through its impact on land–atmosphere interactions. Understanding and correctly modeling soil moisture–precipitation (SM–P) coupling is crucial for improving weather forecasting and subseasonal to seasonal climate predictions, especially when predicting the persistence and magnitude of drought. However, the sign and spatial structure of SM–P feedback are still being debated in the climate research community, mainly due to the difficulty in establishing causal relationships and the high degree of nonlinearity in land–atmosphere processes. To this end, we developed a causal inference model based on the Granger causality analysis and a nonlinear machine learning model. This model includes three steps: nonlinear anomaly decomposition, nonlinear Granger causality analysis, and evaluation of the quality of SM–P feedback, which eliminates the nonlinear response of interannual and seasonal variability and the memory effects of climatic factors and isolates the causal relationship of local SM–P feedback. We applied this model by using National Climate Assessment–Land Data Assimilation System (NCA-LDAS) datasets over the United States. Here, the results highlight the importance of nonlinear atmosphere responses in land–atmosphere interactions. In addition, the strong feedback over the southwestern United States and the Great Plains both highlight the impacts of topographic factors rather than only the sensitivity of evapotranspiration to soil moisture. Furthermore, the SM–P index defined by our framework is used to benchmark Earth system models (ESMs), which provides a new metric for efficiently identifying potential model biases in modeling local land–atmosphere interactions and may help the development of ESMs in improving simulations of water cycle variability and extremes.

54 ENVIRONMENTAL SCIENCES↗

Abundant Nitrate and Nitric Acid Aerosol in the Upper Troposphere and Lower Stratosphere

The tropospheric and stratospheric nitrate aerosol is simulated by a sectional aerosol model coupled to the Community Earth System Model. The simulated nitrate mass fractional contribution to aerosols is significantly higher in the upper troposphere and lower stratosphere (UTLS) than that at the surface. Both in situ measurements and simulations show that nitrate aerosol accounts for about 30%–40% of the aerosol mass at the tropopause of the Asian summer monsoon (ASM) region. Furthermore, simulated condensed nitric acid particles account for ~20% of the annual mean aerosol mass at the tropical tropopause, and over 95% in the UTLS at the South Pole in June-July-August. Here, our study suggests that the extremely cold ambient conditions in the UTLS of the tropics, ASM and polar regions thermodynamically favor the condensation of ammonia and nitric acid. The widely distributed nitrate aerosol in the global UTLS may be overlooked by climate models.

54 ENVIRONMENTAL SCIENCES↗

Supporting hierarchical soil biogeochemical modeling: version 2 of the Biogeochemical Transport and Reaction model (BeTR-v2)

Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in Earth system models (ESMs). The the Biogeochemical Transport and Reaction model version 1 (BeTR-v1) code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multiphase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria, and fungi) dynamics. Here, we describe the new version, Biogeochemical Transport and Reaction model version 2 (BeTR-v2), which adopts more robust numerical solvers for multiphase diffusion and advection and coupling between biogeochemical reactions and improves code modularization over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run for a single topsoil layer or a vertically resolved soil column, and it allows the model to be fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the capability of BeTR-v2 with benchmark cases and example soil biogeochemical (BGC) implementations. By taking advantage of BeTR-v2's generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant–soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants – capabilities that the open-source BeTR-v2 provides. We contend that Earth system models should strive to minimize this uncertainty by applying better numerical solvers.

58 GEOSCIENCES↗