This review paper explores the field of mesoscale to microscale modeling over complex terrain as it traverses multiple so-called gray zones. In an attempt to bridge the gap between previous large-scale and small-scale modeling efforts, atmospheric simulations are being run at an unprecedented range of resolutions. The gray zone is the range of grid resolutions where particular features are neither subgrid nor fully resolved, but rather are partially resolved. The definition of a gray zone depends strongly on the feature being represented and its relationship to the model resolution. This paper explores three gray zones relevant to simulations over complex terrain: turbulence, convection, and topography. Taken together, these may be referred to as the gray continuum. The focus is on horizontal grid resolutions from ∼10 km to ∼10 m. In each case, the challenges are presented together with recent progress in the literature. A common theme is to address cross-scale interaction and scale-awareness in parameterization schemes. How numerical models are designed to cross these gray zones is critical to complex terrain applications in numerical weather prediction, wind resource forecasting, and regional climate modeling, among others.

This study investigates how precipitation-driven cold pools aid the formation of wider clouds that are essential for a transition from shallow to deep convection. In connection with a temperature depression and a depletion of moisture inside developing cold pools, an accumulation of moisture in moist patches around the cold pools is observed. Convective clouds are formed on top of these moist patches. Larger moist patches form with time supporting more and larger clouds. Moreover, enhanced vertical lifting along the leading edges of the gravity current triggered by the cold pool is found. The interplay of moisture aggregation and lifting eventually promotes the formation of wider clouds that are less affected by entrainment and become deeper. These mechanisms are corroborated in a series of cloud-resolving model simulations representing different atmospheric environments. A positive feedback is observed in that, in an atmosphere in which cloud and rain formation is facilitated, stronger downdrafts will form. These stronger downdrafts lead to a stronger modification of the moisture field, which in turn favors further cloud development. This effect is not only observed in the transition phase but also active in prolonging the peak time of precipitation in the later stages of the diurnal cycle. These findings are used to propose a simple way for incorporating the effect of cold pools on cloud sizes and thereby entrainment rate into parameterization schemes for convection. Comparison of this parameterization to the cloud-resolving modeling output gives promising results.

The importance of soil moisture anomalies on airmass convection over semiarid regions has been recognized in several studies. The underlying mechanisms remain partly unclear. An open question is why wetter soils can result in either an increase or a decrease of precipitation (positive or negative soil moisture–precipitation feedback, respectively). Here an idealized cloud-resolving modeling framework is used to explore the local soil moisture–precipitation feedback. The approach is able to replicate both positive and negative feedback loops, depending on the environmental parameters.The mechanism relies on horizontal soil moisture variations, which may develop and intensify spontaneously. The positive expression of the feedback is associated with the initiation of convection over dry soil patches, but the convective cells then propagate over wet patches where they strengthen and preferentially precipitate. The negative feedback may occur when the wind profile is too weak to support the propagation of convective features from dry to wet areas. Precipitation is then generally weaker and falls preferentially over dry patches. The results highlight the role of the midtropospheric flow in determining the sign of the feedback. A key element of the positive feedback is the exploitation of both low convective inhibition (CIN) over dry patches (for the initiation of convection) and high CAPE over wet patches (for the generation of precipitation).

This study investigates the interplay between atmospheric moisture and deep convective clouds via cold‐pool dynamics in the absence of large‐scale forcing in a series of cloud‐resolving modelling studies. More specifically, the contribution of moisture advection, evaporation of rain and surface fluxes to the moisture budget over particular regions of the domain is investigated. This is done both for a continental case and an oceanic case although both cases show very similar behaviour. The accumulation of moisture in confined regions of the sub‐cloud layer that constitute preferred locations for future cloud development mainly results from the advection of moisture. The latter contributes ~86%, minor evaporation of precipitation contributes ~4%, whereas surface moisture fluxes yield ~11% in the continental case. In the oceanic case advection contributes ~125%, surface moisture fluxes ~−32% and evaporation of precipitation ~7%.To further identify the origin of the advected moisture, additional scalars marking moisture originating from the surface and from the evaporation of rain are introduced into the model. It is seen that the surface moisture and the evaporated rain water released within the last 2 h only make 55% of the moisture accumulated in the moist patches in the land case, the rest stemming from older moisture. In the ocean case, this share increases to 72%. The contribution of recently released moisture drops to 28% at cloud base in the continental case and to 56% in the ocean case. The contribution at cloud base is dominated by surface fluxes; the evaporation of rain is negligible.

On summertime fair-weather days, thermally driven wind systems play an important role in determining the initiation of convection and the occurrence of localized precipitation episodes over mountainous terrain. This study compares the mechanisms of convection initiation and precipitation development within a thermally driven flow over an idealized double-ridge system in large-eddy (LESs) and convection-resolving (CRM) simulations. First, LES at a horizontal grid spacing of 200 m is employed to analyze the developing circulations and associated clouds and precipitation. Second, CRM simulations at horizontal grid length of 1 km are conducted to evaluate the performance of a kilometer-scale model in reproducing the discussed mechanisms. Mass convergence and a weaker inhibition over the two ridges flanking the valley combine with water vapor advection by upslope winds to initiate deep convection. In the CRM simulations, the spatial distribution of clouds and precipitation is generally well captured. However, if the mountains are high enough to force the thermally driven flow into an elevated mixed layer, the transition to deep convection occurs faster, precipitation is generated earlier, and surface rainfall rates are higher compared to the LES. Vertical turbulent fluxes remain largely unresolved in the CRM simulations and are underestimated by the model, leading to stronger upslope winds and increased horizontal moisture advection toward the mountain summits. The choice of the turbulence scheme and the employment of a shallow convection parameterization in the CRM simulations change the strength of the upslope winds, thereby influencing the simulated timing and intensity of convective precipitation.

Extreme precipitation events along the Alpine south side (AS) are often forced by upper-level positive potential vorticity (PV) anomalies over western Europe. These so-called PV streamers go along with a dynamical forcing for upward motion, a reduction of the static stability in the troposphere (hence facilitating convection), and are associated with low-level winds that transport moisture toward the Alps.A case of heavy precipitation is examined using the 40-yr ECMWF Re-Analysis data. Piecewise PV inversion (PPVI) and the limited-area Climate High Resolution Model (CHRM) are used to assess the influences of mesoscale parts of the streamer on the precipitation event. The impacts on the vertical stability are quantified by the convective available potential energy (CAPE) and an index of static stability. Very sensitive areas in terms of the stability are located beneath the southern tip of the streamer; smaller changes in the stability are observed in the Alpine region.The moisture transport toward the Alps is sensitive to the amplitude of the streamer, which influences the amount of water that can be transported along its eastern flank.The impacts of the topography on the flow are assessed by calculating an average inverse Froude number. Whether or not the air parcels are blocked by or lifted over the barrier (going along with suppressed and enhanced precipitation, respectively) depends on the vertical stability and the impinging wind velocity, two parameters that are inherently linked to the PV streamer and its substructure.