Search

This study uses semi-idealized simulations to investigate multiscale processes induced by the heterogeneity of soil moisture observed during the 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign. The semi-idealized simulations have realistic land heterogeneity, but large-scale winds are removed. Analysis on isentropic coordinates enables the tracking of circulation that transports energy vertically and facilitates the identification of the primary convective processes induced by realistic land heterogeneity. The isentropes associated with upward motion are found to connect the ground characterized by high latent heat flux to cloud bases directly over the ground with high sensible heat flux, while isentropes associated with downward motion connect precipitation to the ground characterized by high sensible heat fluxes. The mixing of dry, warm parcels ascending from the ground with high sensible heat fluxes and moist parcels from high latent heat regions leads to cloud formation. This new mechanism explains how soil moisture heterogeneity provides the key ingredients such as buoyancy and moisture for shallow cloud formation. We also found that the submesoscale dominates upward energy transport in the boundary layer, while mesoscale circulations contribute to vertical energy transport above the boundary layer. Our novel method better illustrates and elucidates the nature of land atmospheric interactions under irregular and realistic soil moisture patterns. Significance Statement Models that resolve boundary layer turbulence and clouds have been used extensively to understand processes controlling land–atmosphere interactions, but many of their configurations and computational expense limit the use of variable land properties. This study aims to understand how heterogeneous land properties over multiple spatial scales affect energy redistribution by moist convection. Using a more realistic land representation and isentropic analyses, we found that high sensible heat flux regions are associated with relatively higher vertical velocity near the surface, and the high latent heat flux regions are associated with relatively higher moist energy. The mixing of parcels rising from these two regions results in the formation of shallow clouds.

Some of the climate research puzzles relate to a limited understanding of the critical factors governing the life cycle of cumulus clouds. These factors force the initiation and the various mixing processes during cloud life cycles. To shed some light into these processes, we tracked the life cycle of thousands of individual shallow cumulus clouds in a large-eddy simulation during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems field campaign in the U.S. southern Great Plains. Concurrent evolution of clouds is tracked and their respective neighboring clouds are examined. Results show that the clouds initially smaller than neighboring clouds can grow larger than the neighboring clouds by a factor of 2 within 20% of their lifetime. Two groups of the tracked clouds with growing and decaying neighboring clouds, respectively, show distinct characteristics in their life cycles. Clouds with growing neighboring clouds form above regions with larger surface heterogeneity, whereas clouds with decaying neighboring clouds are associated with less heterogeneous surfaces. Also, those with decaying neighboring clouds experience larger instability and a more humid boundary layer, indicating evaporation below the cloud base is likely occurring before those clouds are formed. Larger instability leads to higher vertical velocity and convergence within the cloud, which causes stronger surrounding downdrafts and water vapor removal in the surrounding area. The latter appears to be the reason for the decaying neighboring clouds. Understanding those processes provide insights into how cloud–cloud interactions modulate the evolution of cloud population and into how this evolution can be represented in future cumulus parameterizations.

For the Community Atmosphere Model version 6 (CAM6), an adjustment is needed to conserve dry air mass. This adjustment exposes an inconsistency in how CAM6’s energy budget incorporates water—in CAM6 water in the vapor phase has energy, but condensed phases of water do not. When water vapor condenses, only its latent energy is retained in the model, while its remaining internal, potential, and kinetic energy are lost. A global fixer is used in the default CAM6 model to maintain global energy conservation, but locally the energy tendency associated with water changing phase violates the divergence theorem. This error in energy tendency is intrinsically tied to the water vapor tendency, and reaches its highest values in regions of heavy rainfall, where the error can be as high as 40 W m−2 annually averaged. Several possible changes are outlined within this manuscript that would allow CAM6 to satisfy the divergence theorem locally. These fall into one of two categories: 1) modifying the surface flux to balance the local atmospheric energy tendency and 2) modifying the local atmospheric tendency to balance the surface plus top-of-atmosphere energy fluxes. To gauge which aspects of the simulated climate are most sensitive to this error, the simplest possible change—where condensed water still does not carry energy and a local energy fixer is used in place of the global one—is implemented within CAM6. Comparing this experiment with the default configuration of CAM6 reveals precipitation, particularly its variability, to be highly sensitive to the energy budget formulation. Significance Statement This study examines and explains spurious regional sources and sinks of energy in a widely used climate model. These energy errors result from not tracking energy associated with water after it transitions from the vapor phase to either liquid or ice. Instead, the model used a global fixer to offset the energy tendency related to the energy sources and sinks associated with condensed water species. We replace this global fixer with a local one to examine the model sensitivity to the regional energy error and find a large sensitivity in the simulated hydrologic cycle. This work suggests that the underlying thermodynamic assumptions in the model should be revisited to build confidence in the model-simulated regional-scale water and energy cycles.

Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

While the Madden–Julian Oscillation (MJO) has been shown to affect tropical cyclones (TCs) worldwide through its modulation of large-scale circulation in the atmosphere, little or no role for the ocean has been identified to date in this influence of MJO on TCs. Using observations and numerical model simulations, we demonstrate that MJO events substantially impact TCs over the Maritime Continent (MC) region through an oceanic pathway. While propagating across the MC region, MJO events cause significant sea surface cooling with an area-averaged value of about 0.35 ± 0.12 °C. Hence, TCs over the MC region immediately following the passage of MJO events encounter considerably cooler sea surface temperatures. Consequently, the enthalpy fluxes under the storms are reduced and the intensification rates decrease by more than 50% on average. These results highlight an important role played by the ocean in facilitating MJO-induced sub-seasonal variability in TC activity over the MC region.

The frequency of North Pacific atmospheric rivers (ARs) affects water supply and flood risk over western North America. Thus, understanding factors that affect the variability of landfalling AR frequency is of scientific and societal importance. This study aims at identifying the sources of the moisture for North Pacific ARs and assessing how different modes of variability modulate these sources. To this end, the sources and variability of the background divergent component of the integrated moisture flux (DIVT) in ARs are identified using MERRA reanalysis. It is shown that in the boreal winter, this background DIVT in ARs is related to the outflow from the subsidence over the subtropics which transports moisture northward, while in summer it is related to the Asian monsoon and it transports moisture northwestward. This leads to a seasonal northwest/southeast movement of the AR frequency climatology. At the intra-seasonal scale, propagation of the Madden-Julian Oscillation introduces an anti-clockwise rotation of the background DIVT, with northward transport in phases 1 and 2, westward in 3 and 4, southward in 5 and 6 and eastward in 7 and 8, making landfall over the west coast of North America most likely during the last two phases. Similarly, El Niño Southern Oscillation variability also affects the frequency of ARs through modulation of the westerly background DIVT, favoring landfall over the US west coast during strong El Niño phases. It is shown that in general the likelihood of AR landfall over the western US is correlated with the zonal background DIVT over northeastern Pacific.

It is well documented that over the tropical oceans, column-integrated precipitable water (pw) and precipitation (P) have a nonlinear relationship. In this study moisture budget analysis is used to examine this P–pw relationship in a normalized precipitable water framework. It is shown that the parameters of the nonlinear relationship depend on the vertical structure of moisture convergence. Specifically, the precipitable water values at which precipitation is balanced independently by evaporation versus by moisture convergence define a critical normalized precipitable water, pwnc. This is a measure of convective inhibition that separates tropical precipitation into two regimes: a local evaporation-controlled regime with widespread drizzle and a precipitable water–controlled regime. Most of the 17 CMIP6 historical simulations examined here have higher pwnc compared to ERA5, and more frequently they operate in the drizzle regime. When compared to observations, they overestimate precipitation over the high-evaporation oceanic regions off the equator, thereby producing a “double ITCZ” feature, while underestimating precipitation over the large tropical landmasses and over the climatologically moist oceanic regions near the equator. The responses to warming under the SSP585 scenario are also examined using the normalized precipitable water framework. It is shown that the critical normalized precipitable water value at which evaporation versus moisture convergence balance precipitation decreases as a result of the competing dynamic and thermodynamic responses to warming, resulting in an increase in drizzle and total precipitation. Statistically significant historical trends corresponding to the thermodynamic and dynamic changes are detected in ERA5 and in low-intensity drizzle precipitation in the PERSIANN precipitation dataset.

A record-breaking extreme rainfall event, the highest amount recorded since 1866, hit Indonesia's capital, Jakarta, in early January 2020. This torrential rainfall was mainly caused by the convergence of moisture-rich air due to an unusual blocking of cross-equatorial northerly surge over Northwest Java by the southerly winds induced by a cyclonic flow over the Indian Ocean. This condition caused a local increase in the amount of water vapor over Jakarta that eventually led to the formation of clouds and heavy precipitation. In addition, the concurrent occurrences of convectively active phases of equatorial waves (Kelvin, TD-type, and eastward propagating inertia-gravity waves) and Madden-Julian Oscillation during the event also partly contributed to the enhanced local moisture over the region by increasing low-level moisture flux convergence. Together, these large-scale dynamical drivers provided a convective environment that fostered the development of a massive rain-producing mesoscale convective system and, consequently, extreme rainfall over the region.

Realistic simulation of the Earth’s mean state climate remains a major challenge and yet it is crucial for predicting the climate system in transition. Deficiencies in models’ process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections as well as incorrect attribution of the physical mechanisms governing the past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and longstanding biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when running as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback as well as the fidelity of the simulated base climate state are important for constraining the climate in the past and future.

The present study examines the characteristics of the MJO events represented in the Energy Exascale Earth System Model version 1 (E3SMv1), DOE’s new Earth system model. The coupled E3SMv1 realistic...

Author(s): Balaguru, K; Patricola, CM; Hagos, SM; Leung, LR; Dong, L | Abstract: While El Nino–Southern Oscillation (ENSO) influences eastern North Pacific (ENP) tropical cyclones (TCs) through a variety of atmospheric processes when examined concurrently, ocean pathways dominate at longer lead times. The eastward displacement of the warm pool during an El Nino, which carries warm water into the ENP basin, is the primary oceanic mechanism. Despite this, the question of whether an accurate knowledge of preseason ENSO conditions enhances predictability of ENP TCs has not been addressed specifically. In this study, we show that relative to traditional indices of ENSO, the ENSO Longitude Index (ELI) captures changes in the location of deep convection and associated thermocline processes more accurately. Consequently, the ELI explains more variability in the upper-ocean heat content, and thus TC activity, at lead times of several months in the ENP basin. These results motivate the need to further explore the predictability of ENP TCs associated withnENSO.

In this study, we analyze the realism with which tropical cyclones (TCs) are simulated in the fully coupled low‐ and high‐resolution Energy Exascale Earth System Model (E3SM) version 1, with a focus on the latter. Compared to the low‐resolution (grid spacing of ∼1°), the representation of TCs improves considerably in the high‐resolution configuration (grid spacing of ∼0.25°). Significant improvements are found in the global TC frequency, TC lifetime maximum intensities, and the relative distribution of TCs among the different basins. However, at both resolutions, spurious TC activity is found in some basins, notably in the subtropical regions. Contrasting the simulated large‐scale TC environment with observations reveals that the model environment is unrealistically conducive for TC development in those regions. Further analysis indicates that these biases are likely related to those in thermodynamic potential intensity, caused by systematic SST biases, and vertical wind shear in the coupled model. TC‐ocean interaction is also examined in the high‐resolution configuration of the model. The salient features of the ocean's response to TC‐induced mixing and the ocean's impact on TC intensification are well‐reproduced. Finally, an evaluation of the influence of El Niño Southern Oscillation (ENSO) on TCs in the high‐resolution configuration of the model reveals that the ENSO‐TC relationship in the model has the right sign and is significant for the North Atlantic and Northwest Pacific, albeit weaker than in observations. In summary, the high‐resolution configuration of the E3SM model simulates TC activity reasonably and hence could be a useful tool for TC‐related research.

Traditional parameterizations of the interaction between convection and the environment have relied on an assumption that the slowly-varying large-scale environment is in statistical equilibrium with a large number of small and short-lived convective clouds. They fail to capture non-equilibrium transitions such as the diurnal cycle and theformation of meso-scale convective systems as well asobserved precipitation statisticsand extremes. Informed by analysis of radar observations, cloud-permitting model simulation, theory and machine learning, this work presents a new stochastic cloud population dynamics model for characterizing the interactions between convective and stratiform clouds,with the ultimate goal of informing the representation ofthese interactions in global climate models. 15 wet seasons of precipitating cloud observations by a C-band radar at Darwin, Australia are fed into a machine learning algorithm to obtain transition functions that close a set of coupled equation relating large-scale forcing, mass flux, the convective cell size distribution and the stratiform area. Under realistic large-scale forcing, the derived transition functions show that, on the one hand, interactions with stratiform clouds act to dampen the variability in the size and number of convective cells and therefore in the convective mass flux. On the other hand, for a given convective area fraction, a larger number of smaller cells is more favorable for the growth of stratiform area than a smaller number of larger cells. The combination of these two factors gives rise to solutions with a number of convective cells embedded in a large stratiform area, reminiscent of mesoscale convective systems.

Regional climate simulations over the continental United States were conducted for the 2011 warm season using the Weather Research and Forecasting model at convection‐permitting resolution (4 km) with two commonly used microphysics parameterizations (Thompson and Morrison). Sensitivities of the simulated mesoscale convective system (MCS) properties and feedbacks to large‐scale environments are systematically examined against high‐resolution geostationary satellite and 3‐D mosaic radar observations. MCS precipitation including precipitation amount, diurnal cycle, and distribution of hourly precipitation intensity are reasonably captured by the two simulations despite significant differences in their simulated MCS properties. In general, the Thompson simulation produces better agreement with observations for MCS upper level cloud shield and precipitation area, convective feature horizontal and vertical extents, and partitioning between convective and stratiform precipitation. More importantly, Thompson simulates more stratiform rainfall, which agrees better with observations and results in top‐heavier heating profiles from robust MCSs compared to Morrison. A stronger dynamical feedback to the large‐scale environment is therefore seen in Thompson, wherein an enhanced mesoscale vortex behind the MCS strengthens the synoptic‐scale trough and promotes advection of cool and dry air into the rear of the MCS region. The latter prolongs the MCS lifetimes in the Thompson relative to the Morrison simulations. Hence, different treatment of cloud microphysics not only alters MCS convective‐scale dynamics but also has significant impacts on their macrophysical properties such as lifetime and precipitation. As long‐lived MCSs produced 2–3 times the amount of rainfall compared to short‐lived ones, cloud microphysics parameterizations have profound impact in simulating extreme precipitation and the hydrologic cycle.

The characteristics of Madden–Julian oscillation (MJO) events that strengthen and weaken over the Maritime Continent (MC) are examined. The real-time multivariate MJO (RMM) index is used to assess changes in global MJO amplitude over the MC. The MJO weakens at least twice as often as it strengthens over the MC, with weakening MJOs being twice as likely during El Niño compared to La Niña years and the reverse for strengthening events. MJO weakening shows a pronounced seasonal cycle that has not been previously documented. During the Northern Hemisphere (NH) summer and fall the RMM index can strengthen over the MC. MJOs that approach the MC during the NH winter typically weaken according to the RMM index. This seasonal cycle corresponds to whether the MJO crosses the MC primarily north or south of the equator. Because of the seasonal cycle, weakening MJOs are characterized by positive sea surface temperature and moist-static energy anomalies in the Southern Hemisphere (SH) of the MC compared to strengthening events. Analysis of the outgoing longwave radiation (OLR) MJO index (OMI) shows that MJO precipitation weakens when it crosses the MC along the equator. A possible explanation of this based on previous results is that the MJO encounters more landmasses and taller mountains when crossing along the equator or in the SH. The new finding of a seasonal cycle in MJO weakening over the MC highlights the importance of sampling MJOs throughout the year in future field campaigns designed to study MJO–MC interactions.

CMIP5 models exhibit a mean dry bias and a large inter-model spread in simulating South Asian monsoon precipitation but the origins of the bias and spread are not well understood. Using moisture and energy budget analysis that exploits the weak temperature gradients in the tropics, we derived a non-linear relationship between the normalized precipitation and normalized precipitable water that is similar to the non-linear relationship between precipitation and precipitable water found in previous observational studies. About half of the 21 models analyzed fall in the steep gradient of the non-linear relationship where small differences in the normalized precipitable water in the equatorial Indian Ocean (EIO) manifest in large differences in normalized precipitation in the region. Models with larger normalized precipitable water in the EIO during spring contribute disproportionately to the large inter-model spread and multi-model mean dry bias in monsoon precipitation through perturbations of the large-scale winds. Thus the intermodel spread in precipitable water over EIO leads to the dry bias in the multi-model mean South Asian monsoon precipitation. The models with high normalized precipitable water over EIO also project larger response to warming and dominate the inter-model spread in the multi-model projections of monsoon rainfall. Conversely, models on the flat side of the relationship between normalized precipitation and precipitable water are in better agreement with each other and with observations. On average these models project a smaller increase in the projected monsoon precipitation than that from multi-model mean. This study identified the normalized precipitable water over EIO, which is determined by the relationship between the profiles of convergence and moisture and therefore is an essential outcome of the treatment of convection, as a key metric for understanding model biases and differentiating model skill in simulating South Asian monsoon precipitation.

Sea surface temperature (SST) is one of the most important parameters for tropical cyclone (TC) intensification. Here, it is shown that the relationship between SST and TC intensification varies considerably from basin to basin, with SST explaining less than 4% of the variance in TC intensification rates in the Atlantic, 12% in the western North Pacific, and 23% in the eastern Pacific. Several factors are shown to be responsible for these interbasin differences. First, variability of SST along TCs’ tracks is lower in the Atlantic. This is due to smaller horizontal SST gradients in the Atlantic, compared to the Pacific, and stronger damping of prestorm SST’s contribution to TC intensification by the storm-induced cold SST wake in the Atlantic. The damping occurs because SST tends to vary in phase with TC-induced SST cooling: in the Gulf of Mexico and northwestern Atlantic, where SSTs are highest, TCs tend to be strongest and their translations slowest, resulting in the strongest storm-induced cooling. The tendency for TCs to be more intense over the warmest SST in the Atlantic also limits the usefulness of SST as a predictor since stronger storms are less likely to experience intensification. Finally, SST tends to vary out of phase with vertical wind shear and outflow temperature in the western Pacific. This strengthens the relationship between SST and TC intensification more in the western Pacific than in the eastern Pacific or Atlantic. Combined, these factors explain why prestorm SST is such a poor predictor of TC intensification in the Atlantic, compared to the eastern and western North Pacific.

Sea surface temperature (SST) and tropical cyclone heat potential (TCHP) are metrics used to incorporate the ocean’s influence on hurricane intensification into the National Hurricane Center’s Statistical Hurricane Intensity Prediction Scheme (SHIPS). While both SST and TCHP serve as useful measures of the upper-ocean heat content, they do not accurately represent ocean stratification effects. Here, it is shown that replacing SST within the SHIPS framework with a dynamic temperature Tdy, which accounts for the oceanic negative feedback to the hurricane’s intensity arising from storm-induced vertical mixing and sea surface cooling, improves the model performance. While the model with SST and TCHP explains about 41% of the variance in 36-h intensity changes, replacing SST with Tdy increases the variance explained to nearly 44%. These results suggest that representation of the oceanic feedback, even through relatively simple formulations such as Tdy, may improve the performance of statistical hurricane intensity prediction models such as SHIPS.