Abstract Recent study indicates that the non-instantaneous interaction of convection and circulation is essential for evolution of large-scale convective systems. It is incorporated into cumulus parameterization (CP) by relating cloud-base mass flux of shallow convection to a composite of subcloud moisture convergence in the past 6 h. Three pairs of 19-yr simulations with original and modified CP schemes are conducted in a tropical channel model to verify their ability to reproduce the Madden–Julian oscillation (MJO). More coherent tropical precipitation and improved eastward propagation signal are observed in the simulations with the modified CP schemes based on the non-instantaneous interaction. It is found that enhanced feedback between shallow convection and low-level moisture convergence results in amplified shallow convective heating, and then generates reinforced moisture convergence, which transports more moisture upward. The improved simulations of eastward propagation of the MJO are largely attributed to higher specific humidity below 600 hPa in the free troposphere to the east of maximum rainfall center, which is related to stronger boundary layer moisture convergence forced by shallow convection. Large-scale horizontal advection causes asymmetric moisture tendencies relative to rainfall center (positive to the east and negative to the west) and also gives rise to eastward propagation. The zonal advection, especially the advection of anomalous specific humidity by mean zonal wind, is found to dominate the difference of horizontal advection between each pair of simulations. The results indicate the vital importance of non-instantaneous feedback between shallow convection and moisture convergence for convection organization and the eastward MJO propagation.
more »
« less
Computational Domain Size Effects on Large-Eddy Simulations of Precipitating Shallow Cumulus Convection
Idealized large-eddy simulations of shallow convection often utilize horizontally periodic computational domains. The development of precipitation in shallow cumulus convection changes the spatial structure of convection and creates large-scale organization. However, the limited periodic domain constrains the horizontal variability of the atmospheric boundary layer. Small computational domains cannot capture the mesoscale boundary layer organization and artificially constrain the horizontal convection structure. The effects of the horizontal domain size on large-eddy simulations of shallow precipitating cumulus convection are investigated using four computational domains, ranging from 40×40km2 to 320×320km2 and fine grid resolution (40 m). The horizontal variability of the boundary layer is captured in computational domains of 160×160km2. Small LES domains (≤40 km) cannot reproduce the mesoscale flow features, which are about 100km long, but the boundary layer mean profiles are similar to those of the larger domains. Turbulent fluxes, temperature and moisture variances, and horizontal length scales are converged with respect to domain size for domains equal to or larger than 160×160km2. Vertical velocity flow statistics, such as variance and spectra, are essentially identical in all domains and show minor dependence on domain size. Characteristic horizontal length scales (i.e., those relating to the mesoscale organization) of horizontal wind components, temperature and moisture reach an equilibrium after about hour 30.
more »
« less
- Award ID(s):
- 2143276
- NSF-PAR ID:
- 10489382
- Publisher / Repository:
- Atmosphere
- Date Published:
- Journal Name:
- Atmosphere
- Volume:
- 14
- Issue:
- 7
- ISSN:
- 2073-4433
- Page Range / eLocation ID:
- 1186
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract The goal of this study is to challenge a large eddy simulation model with a range of observations from a modern field campaign and to develop case studies useful to other modelers. The 2015 Cloud System Evolution in the Trades (CSET) field campaign provided a wealth of in situ and remote sensing observations of subtropical cloud transitions in the summertime Northeast Pacific. Two Lagrangian case studies based on these observations are used to validate the thermodynamic, radiative and microphysical properties of large eddy simulations (LES) of the stratocumulus to cumulus transition. The two cases contrast a relatively fast cloud transition in a clean, initially well-mixed boundary layer vs. a slower transition in an initially decoupled boundary layer with higher aerosol concentrations and stronger mean subsidence. For each case, simulations of two neighboring trajectories sample mesoscale variability and the coherence of the transition in adjacent air masses. In both cases, LES broadly reproduce satellite and aircraft observations of the transition. Simulations of the first case match observations more closely than for the second case, where simulations underestimate cloud cover early in the simulations and overestimate cloud top height later. For the first case, simulated cloud fraction and liquid water path increase if a larger cloud droplet number concentration is prescribed. In the second case, precipitation onset and inversion cloud breakup occurs earlier when the LES domain is chosen large enough to support strong mesoscale organization.more » « less
-
more » « less
Abstract A 3‐D cloud‐resolving model has been used to investigate the domain size dependence of simulations of convective self‐aggregation (CSA) in radiative‐convective equilibrium. We investigate how large a domain is needed to allow multiple convective clusters and also how the properties equilibrated CSA depend on domain size. We used doubly periodic square domains of widths 768, 1,536, 3,072, and 6,144 km, over 350 simulated days. In the 768‐, 1,536‐, and 3,072‐km domains, the simulations produced circular convective clusters surrounded by broader regions of dry, subsiding air. In the 6,144‐km domain, the convection ultimately forms two semiconnected bands. As the domain size increases, equilibrated CSA moistens in two ways. First, as the circulation widens, this leads to stronger boundary layer winds and a more humid boundary layer. Second, the stronger inflow into the convective region boundary layer is associated with a warmer convective region boundary layer, which leads to intensified deep convection, more melting and freezing near the freezing level, enhanced midlevel stability, increased congestus activity, and detrainment of moist air into the dry region. In the larger domains, the deep convection and congestus slowly oscillate out of phase with each other with a time period of about 25 to 30 days. We hypothesize that other important domain size sensitivities, including a decrease in net moist static energy export from the convective region, are fundamentally linked to the increasing relationship between domain size and boundary layer wind speed. Our results suggest that the statistics of CSA converge only for domains wider than about 3,000 km.
-
more » « less
Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement.
Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.
-
more » « less
A two‐column radiative–convective equilibrium (RCE) model is used to study the depth of convection that develops in the subsiding branch of a Walker‐like overturning circulation. The model numerically solves for two‐dimensional non‐rotating hydrostatic flow, which is damped by momentum diffusion in the boundary layer and model interior, and by convective momentum transport. Convection, clouds and radiative transfer are parametrized, and the convection scheme does not include explicit freezing or melting.
While integrating the model towards local RCE, the level of neutral buoyancy (LNB) fluctuates between mid‐ and high levels. Evaporation of detrained moisture at the LNB locally cools the environment, so that the final RCE state has a stable layer at mid‐levels (550 hPa ≈ 50–100 hPa below 0 °C), which is unrelated to melting of ice. Preferred detrainment at mid‐ and high levels leaves the middle‐to‐upper troposphere relatively dry.
A circulation is introduced by incrementally lowering the sea‐surface temperature in one column, which collapses convection: first to a congestus mode with tops near 550 hPa, below the dry layer created in RCE; then to congestus with tops near 650 hPa; and finally to shallow cumulus with tops near 850 hPa. Critical to stabilizing congestus near 650 hPa is large radiative cooling near moist cumulus tops under a dry upper atmosphere. This congestus mode is very sensitive, and only develops when horizontal temperature gradients created by evaporative and radiative cooling can persist against the work of gravity waves. This only happens in runs with ample momentum diffusion, which are those with convective momentum transport or large domains.
Compared to the shallow mode, the congestus mode produces a deep moist layer and more precipitation. This reduces radiative cooling in the cloud layer and enhances stability near the cloud base, which weakens the circulation, and leads to less precipitation over the warm ocean.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/atmos14071186
