Journal of Advances in Modeling Earth Systems最新文献

Climate models have long-standing difficulties simulating the South Pacific Convergence Zone (SPCZ) and its variability. For example, the default Zhang-McFarlane (ZM) convection scheme in the Community Atmosphere Model version 5 (CAM5) produces too much light precipitation and too little heavy precipitation in the SPCZ, with this bias toward light precipitation even more pronounced in the SPCZ than in the tropics as a whole. Here, we show that implementing a recently developed convection scheme in the CAM5 yields significant improvements in the simulated SPCZ during austral summer and discuss the reasons behind these improvements. In addition to intensifying both mean rainfall and its variability in the SPCZ, the new scheme produces a larger heavy rainfall fraction that is more consistent with observations and state-of-the-art reanalyses. This shift toward heavier, more variable rainfall increases both the magnitude and altitude of diabatic heating associated with convective precipitation, intensifying lower tropospheric convergence and increasing the influence of convection on the upper-level circulation. Increased diabatic production of potential vorticity in the upper troposphere intensifies the distortion effect exerted by convection on transient Rossby waves that pass through the SPCZ. Weaker distortion effects in simulations using the ZM scheme allow waves to propagate continuously through the region rather than dissipating locally, further reducing updrafts and weakening convection in the SPCZ. Our results outline a dynamical framework for evaluating model representations of tropical–extratropical interactions within the SPCZ and clarify why convective parameterizations that produce “top-heavy” profiles of deep convective heating better represent the SPCZ and its variability.

The effects of small-scale topography-induced land surface heterogeneity are not well represented in current Earth System Models (ESMs). In this study, a new topography-based subgrid structure referred to as topographic units (TGU) designed to better capture subgrid topographic effects, and methods to downscale atmospheric forcing to the land TGUs have been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Effects of the subgrid scheme and downscaling methods on ELM simulated land surface processes are evaluated over the conterminous United States (CONUS). For this purpose, ELM simulations are performed using two configurations without (NoD ELM) and with (D ELM) downscaling, both using TGUs derived for the 0.5-degree grids and the same land surface parameters. Simulations using the two ELM configurations are compared over the CONUS domain, regional levels, and at observational sites (e.g., SNOTEL). The CONUS-level results suggest that D ELM simulates more snowfall and snow water equivalent (SWE), higher runoff, and less ET during spring and summer. Regional-level results suggest more pronounced impacts of downscaling over regions dominated by higher elevation TGUs and regions with maximum precipitation occurring during cool seasons. Results at the SNOTEL sites suggest that D ELM has superior capability of reproducing the observed SWE at 83% of the sites, with more pronounced performance over topographically heterogeneous TGUs with their maximum precipitation occurring during cool seasons. The results highlight the importance of improving representation of small-scale surface heterogeneity in ESMs and motivate future research to understand their effects on land-atmosphere interactions, streamflow, and water resources management over mountainous regions.

Machine learning (ML) is used to build a bulk microphysical parameterization including ice processes. Simulations of the Lagrangian super-particle model McSnow are used as training data. The ML performs a coarse-graining of the particle-resolved microphysics to multi-category two-moment bulk equations. Besides mass and number, prognostic particle properties (P3) like melt water, rime mass, and rime volume are predicted by the ML-based bulk model. The ML-based scheme is tested with simulations of increasing complexity. As a box model, the ML-based bulk scheme can reproduce the simulations of McSnow quite accurately. In 3d idealized squall line simulations, the ML-based P3-like scheme provides a more realistic extended stratiform region when compared to the standard two-moment bulk scheme in ICON. In a realistic case study, the ML-based scheme runs stably, but can not significantly improve the results. This shows that ML can be used to coarse-grain super-particle simulations to a bulk scheme of arbitrary complexity.

Tropical forest diversity governs forest structures, compositions, and influences the ecosystem response to environmental changes. Better representation of forest diversity in ecosystem demography (ED) models within Earth system models is thus necessary to accurately capture and predict how tropical forests affect Earth system dynamics subject to climate changes. However, achieving forest coexistence in ED models is challenging due to their computational expense and limited understanding of the mechanisms governing forest functional diversity. This study applies the advanced Multi-Objective Population-based Parallel Local Surrogate-assisted search (MOPLS) optimization algorithm to simultaneously calibrate ecosystem fluxes and coexistence of two physiologically distinct tropical forest species in a size- and age-structured ED model with realistic representation of wood harvest. MOPLS exhibits satisfactory model performance, capturing hydrological and biogeochemical dynamics observed in Barro Colorado Island, Panama, and robustly achieving coexistence for the two representative forest species. This demonstrates its effectiveness in calibrating tropical forest coexistence. The optimal solution is applied to investigate the recovery trajectories of forest biomass after various intensities of clear-cut deforestation. We find that a 20% selective logging can take approximately 40 years for aboveground biomass to return to the initial level. This is due to the slow recovery rate of late successional trees, which only increases by 4% over the 40-year period. This study lays the foundation to calibrate coexistence in ED models. MOPLS can be an effective tool to help better represent tropical forest diversity in Earth system models and inform forest management practices.

This study proposes deterministic and stochastic energy-aware hybrid models that should enable simulations of idealized and primitive-equations Geophysical Fluid Dynamics (GFD) models at low resolutions without compromising on quality compared with high-resolution runs. Such hybrid models bridge the data-driven and physics-driven modeling paradigms by combining regional stability and classical GFD models at low resolution that cannot reproduce high-resolution reference flow features (large-scale flows and small-scale vortices) which are, however, resolved. Hybrid models use an energy-aware correction of advection velocity and extra forcing compensating for the drift of the low-resolution model away from the reference phase space. The main advantages of hybrid models are that they allow for physics-driven flow recombination within the reference energy band, reproduce resolved reference flow features, and produce more accurate ensemble forecasts than their classical GFD counterparts. Hybrid models offer appealing benefits and flexibility to the modeling and forecasting communities, as they are computationally cheap and can use both numerically-computed flows and observations from different sources. All these suggest that the hybrid approach has the potential to exploit low-resolution models for long-term weather forecasts and climate projections thus offering a new cost effective way of GFD modeling. The proposed hybrid approach has been tested on a three-layer quasi-geostrophic model for a beta-plane Gulf Stream flow configuration. The results show that the low-resolution hybrid model reproduces the reference flow features that are resolved on the coarse grid and also gives a more accurate ensemble forecast than the physics-driven model.

Data assimilation (DA) plays a pivotal role in diverse applications, ranging from weather forecasting to trajectory planning for autonomous vehicles. A prime example is the widely used ensemble Kalman filter (EnKF), which relies on the Kalman filter's linear update equation to correct each of the ensemble forecast member's state with incoming observations. Recent advancements have witnessed the emergence of deep learning approaches in this domain, primarily within a supervised learning framework. However, the adaptability of such models to untrained scenarios remains a challenge. In this study, we introduce a new DA strategy that utilizes reinforcement learning (RL) to apply state corrections using full or partial observations of the state variables. Our investigation focuses on demonstrating this approach to the chaotic Lorenz 63 and 96 systems, where the agent's objective is to maximize the geometric series with terms that are proportional to the negative root-mean-squared error (RMSE) between the observations and corresponding forecast states. Consequently, the agent develops a correction strategy, enhancing model forecasts based on available observations. Our strategy employs a stochastic action policy, enabling a Monte Carlo-based DA framework that relies on randomly sampling the policy to generate an ensemble of assimilated realizations. Numerical results demonstrate that the developed RL algorithm performs favorably when compared to the EnKF. Additionally, we illustrate the agent's capability to assimilate non-Gaussian observations, addressing one of the limitations of the EnKF.

Large-scale dynamical and thermodynamical processes are common environmental drivers of high-impact weather systems causing extreme weather events. However, such large-scale environmental conditions often display systematic biases in climate simulations, posing challenges to evaluating high-impact weather systems and extreme weather events. In this paper, a machine learning (ML) approach was employed to bias correct the large-scale wind, temperature, and humidity simulated by the atmospheric component of the Energy Exascale Earth System Model (E3SM) at ∼1° resolution. The usefulness of the ML approach for extreme weather analysis was demonstrated with a focus on three high-impact weather systems, including tropical cyclones (TCs), extratropical cyclones (ETCs), and atmospheric rivers (ARs). We show that the ML model can effectively reduce climate bias in large-scale wind, temperature, and humidity while preserving their responses to imposed climate change perturbations. The bias correction is found to directly improve water vapor transport associated with ARs, and representations of thermodynamical flows associated with ETCs. When the bias-corrected large-scale winds are used to drive a synthetic TC track forecast model over the Atlantic basin, the resulting TC track density agrees better with that of the TC track model driven by observed winds. In addition, the ML model insignificantly interferes with the mean climate change signals of large-scale storm environments as well as the occurrence and intensity of three weather systems. This study suggests that the proposed ML approach can be used to improve the downscaling of extreme weather events by providing more realistic large-scale storm environments simulated by low-resolution climate models.

The UK contribution to the Detection and Attribution Model Intercomparison Project (DAMIP), part of the sixth phase of the Climate Model Intercomparison Project (CMIP6), is described. The lower atmosphere and ocean resolution configuration of the latest Hadley Centre global environmental model, HadGEM3-GC3.1, is used to create simulations driven either with historical changes in anthropogenic well-mixed greenhouse gases, anthropogenic aerosols, or natural climate factors. Global mean near-surface air temperatures from the HadGEM3-GC31-LL simulations are consistent with CMIP6 model ensembles for the equivalent experiments. While the HadGEM3-GC31-LL simulations with anthropogenic and natural forcing factors capture the overall observed warming, the lack of marked simulated warming until the 1990s is diagnosed as due to aerosol cooling mostly offsetting the well-mixed greenhouse gas warming until then. The model has unusual temperature variability over the Southern Ocean related to occasional deep convection bringing heat to the surface. This is most prominent in the model's aerosol only simulations, which have the curious feature of warming in the high southern latitudes, while the rest of the globe cools, a behavior not seen in other CMIP6 models. This has implications for studies that assume model responses, from different climate drivers, can be linearly combined. While DAMIP was predominantly designed for detection and attribution studies, the experiments are also very valuable for understanding how different climate drivers influence a model, and thus for interpretating the responses of combined anthropogenic and natural driven simulations. We recommend institutions provide model simulations for the high priority DAMIP experiments.

Eddy viscosity is employed throughout the majority of numerical fluid dynamical models, and has been the subject of a vigorous body of research spanning a variety of disciplines. It has long been recognized that the proper description of eddy viscosity uses tensor mathematics, but in practice it is almost always employed as a scalar due to uncertainty about how to constrain the extra degrees of freedom and physical properties of its tensorial form. This manuscript borrows techniques from outside the realm of geophysical fluid dynamics to consider the eddy viscosity tensor using its eigenvalues and eigenvectors, establishing a new framework by which tensorial eddy viscosity can be tested. This is made possible by a careful analysis of an operation called tensor unrolling, which casts the eigenvalue problem for a fourth-order tensor into a more familiar matrix-vector form, whereby it becomes far easier to understand and manipulate. New constraints are established for the eddy viscosity coefficients that are guaranteed to result in energy dissipation, backscatter, or a combination of both. Finally, a testing protocol is developed by which tensorial eddy viscosity can be systematically evaluated across a wide range of fluid regimes.

Numerous Earth system models exhibit excessive aerosol effective forcing at the top of the atmosphere (TOA), including the Department of Energy's Energy Exascale Earth System Model (E3SM). Here, in the context of the E3SM version 3 effort, the predicted particle property (P3) stratiform cloud microphysics scheme and an enhanced deep convection parameterization suite (ZM_plus) are implemented into E3SM. The ZM_plus includes a convective cloud microphysics scheme, a multi-scale coherent structure parameterization for mesoscale convective systems, and a revised cloud base mass flux formulation considering impacts of the large-scale environment. The P3 scheme improved cloud and radiation particularly over the Northern Hemisphere and the frequency of heavy precipitation over the tropics, and the ZM_plus improved clouds in the tropics. P3 decreases aerosol effective forcing by 0.15 W m⁻², while the ZM_plus increases it by 0.27 W m⁻², resulting from excessive direct (0.31 W m⁻²) and indirect forcing (−1.79 W m⁻²). The excessive aerosol forcings are due to aerosol overestimation associated with insufficient aerosol wet removal. By improving the physical treatments in the aerosol wet removal, we effectively mitigate anthropogenic aerosol overestimation and thus attenuate direct (0.09 W m⁻²) and indirect aerosol forcing (−1.52 W m⁻²). Adjustment to primary organic matter hygroscopicity reduces direct and indirect forcing to more reasonable values: −0.13 W m⁻² and −1.31 W m⁻², respectively. On climatology, improved aerosol treatments mitigate overestimation of aerosol optical depth.