An official website of the United States government
Convective parameterization is the long-lasting bottleneck of global climate modelling and one of the most difficult problems in atmospheric sciences. Uncertainty in convective parameterization is the leading cause of the widespread climate sensitivity in IPCC global warming projections. This paper reviews the observations and parameterizations of atmospheric convection with emphasis on the cloud structure, bulk effects, and closure assumption. The representative state-of-the-art convection schemes are presented, including the ECMWF convection scheme, the Grell scheme used in NCEP model and WRF model, the Zhang-MacFarlane scheme used in NCAR and DOE models, and parameterizations of shallow moist convection. The observed convection has self-suppression mechanisms caused by entrainment in convective updrafts, surface cold pool generated by unsaturated convective downdrafts, and warm and dry lower troposphere created by mesoscale downdrafts. The post-convection environment is often characterized by “diamond sounding” suggesting an over-stabilization rather than barely returning to neutral state. Then the pre-convection environment is characterized by slow moistening of lower troposphere triggered by surface moisture convergence and other mechanisms. The over-stabilization and slow moistening make the convection events episodic and decouple the middle/upper troposphere from the boundary layer, making the state-type quasi-equilibrium hypothesis invalid. Right now, unsaturated convective downdrafts and especially mesoscale downdrafts are missing in most convection schemes, while some schemes are using undiluted convective updrafts, all of which favour easily turned-on convection linked to double-ITCZ (inter-tropical convergence zone), overly weak MJO (Madden-Julian Oscillation) and precocious diurnal precipitation maximum. We propose a new strategy for convection scheme development using reanalysis-driven model experiments such as the assimilation runs in weather prediction centres and the decadal prediction runs in climate modelling centres, aided by satellite simulators evaluating key characteristics such as the lifecycle of convective cloud-top distribution and stratiform precipitation fraction.
more »
« less
This content will become publicly available on August 1, 2026
Edge-Averaged Virtual Element Methods for Convection-Diffusion and Convection-Dominated Problems
More Like this
-
-
null (Ed.)Abstract ‘Horizontal convection’ (HC) is the generic name for the flow resulting from a buoyancy variation imposed along a horizontal boundary of a fluid. We study the effects of rotation on three-dimensional HC numerically in two stages: first, when baroclinic instability is suppressed and, second, when it ensues and baroclinic eddies are formed. We concentrate on changes to the thickness of the near-surface boundary layer, the stratification at depth, the overturning circulation and the flow energetics during each of these stages. Our results show that, for moderate flux Rayleigh numbers ( $$O(1{0}^{11} )$$ ), rapid rotation greatly alters the steady-state solution of HC. When the flow is constrained to be uniform in the transverse direction, rapidly rotating solutions do not support a boundary layer, exhibit weaker overturning circulation and greater stratification at all depths. In this case, diffusion is the dominant mechanism for lateral buoyancy flux and the consequent buildup of available potential energy leads to baroclinically unstable solutions. When these rapidly rotating flows are perturbed, baroclinic instability develops and baroclinic eddies dominate both the lateral and vertical buoyancy fluxes. The resulting statistically steady solution supports a boundary layer, larger values of deep stratification and multiple overturning cells compared with non-rotating HC. A transformed Eulerian-mean approach shows that the residual circulation is dominated by the quasi-geostrophic eddy streamfunction and that the eddy buoyancy flux has a non-negligible interior diabatic component. The kinetic and available potential energies are greater than in the non-rotating case and the mixing efficiency drops from $${\sim }0. 7$$ to $${\sim }0. 17$$ . The eddies play an important role in the formation of the thermal boundary layer and, together with the negatively buoyant plume, help establish deep stratification. These baroclinically active solutions have characteristics of geostrophic turbulence.more » « less
-
This work employs single-mode equations to study convection and double-diffusive convection in a porous medium where the Darcy law provides large-scale damping. We first consider thermal convection with salinity as a passive scalar. The single-mode solutions resembling steady convection rolls reproduce the qualitative behavior of root-mean-square and mean temperature profiles of time-dependent states at high Rayleigh numbers from direct numerical simulations (DNS). We also show that the single-mode solutions are consistent with the heat-exchanger model that describes well the mean temperature gradient in the interior. The Nusselt number predicted from the single-mode solutions exhibits a scaling law with Rayleigh number close to that followed by exact 2D steady convection rolls, although large aspect ratio DNS results indicate a faster increase. However, the single-mode solutions at a high wavenumber predict Nusselt numbers close to the DNS results in narrow domains. We also employ the single-mode equations to analyze the influence of active salinity, introducing a salinity contribution to the buoyancy, but with a smaller diffusivity than the temperature. The single-mode solutions are able to capture the stabilizing effect of an imposed salinity gradient and describe the standing and traveling wave behaviors observed in DNS. The Sherwood numbers obtained from single-mode solutions show a scaling law with the Lewis number that is close to the DNS computations with passive or active salinity. This work demonstrates that single-mode solutions can be successfully applied to this system whenever periodic or no-flux boundary conditions apply in the horizontal.more » « less
-
null (Ed.)We use well resolved numerical simulations with the lattice Boltzmann method to study Rayleigh–Bénard convection in cells with a fractal boundary in two dimensions for $Pr = 1$ and $$Ra \in \left [10^7, 10^{10}\right ]$$ , where Pr and Ra are the Prandtl and Rayleigh numbers. The fractal boundaries are functions characterized by power spectral densities $S(k)$ that decay with wavenumber, $$k$$ , as $$S(k) \sim k^{p}$$ ( $p < 0$ ). The degree of roughness is quantified by the exponent $$p$$ with $p < -3$ for smooth (differentiable) surfaces and $$-3 \le p < -1$$ for rough surfaces with Hausdorff dimension $$D_f=\frac {1}{2}(p+5)$$ . By computing the exponent $$\beta$$ using power law fits of $$Nu \sim Ra^{\beta }$$ , where $Nu$ is the Nusselt number, we find that the heat transport scaling increases with roughness through the top two decades of $$Ra \in \left [10^8, 10^{10}\right ]$$ . For $$p$$ $= -3.0$ , $-2.0$ and $-1.5$ we find $$\beta = 0.288 \pm 0.005, 0.329 \pm 0.006$$ and $$0.352 \pm 0.011$$ , respectively. We also find that the Reynolds number, $Re$ , scales as $$Re \sim Ra^{\xi }$$ , where $$\xi \approx 0.57$$ over $$Ra \in \left [10^7, 10^{10}\right ]$$ , for all $$p$$ used in the study. For a given value of $$p$$ , the averaged $Nu$ and $Re$ are insensitive to the specific realization of the roughness.more » « less
