Search

Open boundary conditions were developed for atmospheric large-eddy simulation (LES) models and implemented into the Dutch Atmospheric Large-Eddy Simulation model. The implementation was tested in a "Big Brother"-like setup, in which the simulation with open boundary conditions was forced by an identical control simulation with periodic boundary conditions. The results show that the open boundary implementation has minimal influence on the solution. Both the mean state and the turbulent structures are close to the control simulation, and disturbances at the in- and outflow boundaries are negligible. To emulate a setup in which the LES is coupled to a coarser model, the influence of coarse boundary input was tested by smoothing the output of the periodic control simulation both temporally and spatially before feeding it as input to the simulation with open boundary conditions. When smoothing is applied over larger spatial and longer temporal scales, disturbances start to form at the inflow boundary and an area exists where turbulence needs to develop. Adding synthetic turbulence to the smoothed input reduces the size of this area and the magnitude of the disturbances. [ABSTRACT FROM AUTHOR]

Cloud shadows can be detected in the radiance measurements of the TROPOMI instrument on board the Sentinel-5P satellite due to its high spatial resolution, and could possibly affect its air quality products. The cloud shadow induced signatures are, however, not always apparent and may depend on various cloud and scene parameters. Hence, the quantification of the cloud shadow impact requires the analysis of large data sets. Here we use the cloud shadow detection algorithm DARCLOS to detect cloud shadow pixels in the TROPOMI absorbing aerosol index (AAI) product over Europe during 8 months. For every shadow pixel, we automatically select cloud- and shadow-free neighbour pixels, in order to estimate the cloud shadow induced signature. In addition, we simulate the measured cloud shadow impact on the AAI with our newly developed 3D radiative transfer algorithm MONKI. Both the measurements and simulations show that the average cloud shadow impact on the AAI is close to zero (0.06 and 0.16, respectively). However, the top-of-atmosphere reflectance ratio between 340 and 380 nm, which is used to compute the AAI, is significantly increased in 95 % of the shadow pixels. So, cloud shadows are bluer than surrounding non-shadow pixels. Our simulations explain that the traditional AAI formula intrinsically already corrects for this cloud shadow effect, via the lower retrieved scene albedo. This cancellation of cloud shadow signatures is not always perfect, sometimes yielding second order low and high biases in the AAI which we also successfully reproduce with our simulations. We show that the magnitude of those second order cloud shadow effects depends on various cloud parameters which are difficult to determine for the shadows measured with TROPOMI. We conclude that a potential cloud shadow correction strategy for the TROPOMI AAI would therefore be complicated if not unnecessary. [ABSTRACT FROM AUTHOR]

Open boundary conditions were developed for atmospheric large eddy simulation (LES) models and implemented into the Dutch Atmospheric Large Eddy simulation model. The implementation was tested in a "Big Brother"-like setup, in which the simulation with open boundary conditions was forced by an identical control simulation with periodic boundary conditions. The results show that the open boundary implementation has minimal influence on the solution. Both the mean 5 state and the turbulent structures are close to the control simulation and disturbances at the in- and outflow boundaries are negligible. To emulate a setup in which the LES is coupled to a coarser model, the influence of coarse boundary input was tested by smoothing the output of the periodic control simulation both temporally and spatially before feeding it as input to the simulation with open boundary conditions. The results show that when the ratio between input and model resolution increases, disturbances start to form at the inflow boundary and an area exists where turbulence needs to develop. Adding synthetic 10 turbulence to the smoothed input reduces the size of this area and the magnitude of the disturbances. [ABSTRACT FROM AUTHOR]

Quantitative estimates of precipitation in a typical undisturbed trade wind region are derived from 2 months of radar reflectivity data and compared to the meteorological environment determined from soundings, surface flux, and airborne-lidar data. Shallow precipitation was ubiquitous, covering on average about 2% of the region and contributing to at least half of the total precipitation. Echo fractions on the scale of the radar domain range between 0% and 10% and vary greatly within a period from a few hours to a day. Variability in precipitation relates most strongly to variability in humidity and the zonal wind speed, although greater inversion heights and deeper clouds are also evident at times of more rain. The analysis herein suggests that subtle fluctuations in both the strength of the easterlies and in subsidence play a major role in regulating humidity and hence precipitation, even within a given meteorological regime (here, the undisturbed trades).

The propagation of optical and acoustical waves is affected by the atmospheric turbulence through the local instantaneous refractive index structure parameter in a volume of characteristic size r, denoted [C.sub.n,r.sup.2]. Often r is well within the inertial--convective range. In this study, a large-eddy simulation (LES) of spatial resolution [DELTA] is used to analyze the distribution of [C.sub.n,[DELTA].sup.2] in the convective boundary layer. The local formulation used to calculate [C.sub.n,[DELTA].sup.2] is described and is found to back up the subgrid parameterization of atmospheric LES. The mean vertical profile behaves according to the mixed layer similarity theory. The spatial organization relates to the presence of buoyant ascending plumes. This structure is associated with some bimodal probability density functions of the inertial-convective range variables, in agreement with experimental results. The standard model of jointly lognormal statistics is challenged. The intermittency of these variables is characterized, and its inertial-convective range component quantitatively agrees with experimental estimates. The present LES results are used to document the bias produced by averaging the dissipation rates of TKE and temperature variance to estimate the average structure parameters. Comparing with nonconvective turbulence measurements, the bias for [C.sub.T,r.sup.2] is found to depend on the stability. A physical explanation is offered that emphasizes the role of the large-scale turbulence under convective conditions.

A better conceptual understanding and more realistic parameterizations of convective boundary layers in climate and weather prediction models have been major challenges in meteorological research. In particular, parameterizations of the dry convective boundary layer, in spite of the absence of water phase-changes and its consequent simplicity as compared to moist convection, typically suffer from problems in attempting to represent realistically the boundary layer growth and what is often referred to as countergradient fluxes. The eddy-diffusivity (ED) approach has been relatively successful in representing some characteristics of neutral boundary layers and surface layers in general. The mass-flux (MF) approach, on the other hand, has been used for the parameterization of shallow and deep moist convection. In this paper, a new approach that relies on a combination of the ED and MF parameterizations (EDMF) is proposed for the dry convective boundary layer. It is shown that the EDMF approach follows naturally from a decomposition of the turbulent fluxes into 1) a part that includes strong organized updrafts, and 2) a remaining turbulent field. At the basis of the EDMF approach is the concept that nonlocal subgrid transport due to the strong updrafts is taken into account by the MF approach, while the remaining transport is taken into account by an ED closure. Large-eddy simulation (LES) results of the dry convective boundary layer are used to support the theoretical framework of this new approach and to determine the parameters of the EDMF model. The performance of the new formulation is evaluated against LES results, and it is shown that the EDMF closure is able to reproduce the main properties of dry convective boundary layers in a realistic manner. Furthermore, it will be shown that this approach has strong advantages over the more traditional countergradient approach, especially in the entrainment layer. As a result, this EDMF approach opens the way to parameterize the clear and cumulus-topped boundary layer in a simple and unified way.

This paper reports an intercomparison study on undisturbed trade wind cumulus convection under steady-state conditions as observed during the Barbados Oceanographic and Meteorological Experiment (BOMEX) with 10 large eddy simulation (LES) models. A main objective of this study is to obtain a quantitative assessment of the quality of the turbulent dynamics for this type of boundary layer clouds as produced by the different LES codes. A 6-h simulation shows excellent model-to-model agreement of the observed vertical thermodynamical structure, reasonable agreement of variances and turbulent fluxes, and good agreement of quantities conditionally sampled within the model clouds, such as cloud cover, liquid water, and cloud updraft strength. In the second part of this paper the LES dataset is used to evaluate simple models that are used in parameterizations of current general circulation models (GCMs). Finally, the relation of this work to subsequent LES studies of more complicated regimes is discussed, and guidance is given for the design of future observational studies of shallow cumulus boundary layers.

Cloud shadows are observed by the TROPOMI satellite instrument as a result of its high spatial resolution as compared to its predecessor instruments. These shadows contaminate TROPOMI's air quality measurements, because shadows are generally not taken into account in the models that are used for aerosol and trace gas retrievals. If the shadows are to be removed from the data, or if shadows are to be studied, an automatic detection of the shadow pixels is needed. We present the Detection AlgoRithm for CLOud Shadows (DARCLOS) for TROPOMI, which is the first cloud shadow detection algorithm for a spaceborne spectrometer. DARCLOS raises potential cloud shadow flags (PCSFs), and actual cloud shadow flags (ACSFs). The PCSFs indicate the TROPOMI ground pixels that are potentially affected by cloud shadows based on a geometric consideration with safety margins. The ACSFs are a refinement of the PCSFs using spectral reflectance information of the PCSF pixels, and identify the TROPOMI ground pixels that are confidently affected by cloud shadows. We validate DARCLOS with true color images made by the VIIRS instrument on board of Suomi NPP orbiting in close constellation with TROPOMI on board of Sentinel 5-P. We conclude that the PCSF can be used to exclude cloud shadow contamination from TROPOMI data, while the ACSF can be used to select pixels for the scientific analysis of cloud shadow effects. [ABSTRACT FROM AUTHOR]

With increasing temperatures, it is likely that precipitation extremes increase as well. While, on larger spatial and longer temporal scales, the amplification of rainfall extremes often follows the Clausius‐Clapeyron relation, it has been shown that local short‐term convective precipitation extremes may well exceed the Clausius‐Clapeyron rate of around 6.5%/K. Most studies on this topic have focused exclusively on the intensity aspect, while only few have examined (with contradictory results) how warmer and moister conditions modulate the spatial characteristics of convective precipitation extremes and how these connect to increased intensities. Here we study this relation by using a large eddy simulation model. We simulate one diurnal cycle of heavy convective precipitation activity based on a realistic observation‐based strongly forced case setup. Systematically perturbed initial conditions of temperature and specific humidity enable an examination of the response of intensities and spatial characteristics of the precipitation field over an 8° dew point temperature range. We find that warmer and moister conditions result in an overall increase of both intensities and spatial extent of individual rain cells. Colder conditions favor the development of many but smaller rain cells. Under warmer conditions, we find a reduced number of individual cells, but their size significantly grows along with an increase of intensities over a large part of a rain cell. Combined, these factors lead to larger and more intense rain cells that can produce up to almost 20% more rain per degree warming and therefore have a large impact. Key Points: Warmer and moister conditions reshape the precipitation field to fewer but larger rain cellsLarger rain cells produce more intense rainRain rates aggregated over a rain cell can increase by up to almost 20% per degree warming [ABSTRACT FROM AUTHOR]

The fifth intercomparison of the Global Water and Energy Experiment Cloud System Studies Working Group 1 is used as a vehicle for better understanding the dynamics of trade wind cumuli capped by a strong inversion. The basis of the intercomparison is 10 simulations by 7 groups. These simulations are supplemented by many further sensitivity studies, including some with very refined grid meshes. The simulations help illustrate the turbulent dynamics of trade cumuli in such a regime. In many respects the dynamics are similar to those found in many previous simulations of trade cumuli capped by weaker inversions. The principal differences are the extent to which the cloud layer is quasi-steady in the current simulations, evidence of weak countergradient momentum transport within the cloud layer, and the development and influence of an incipient stratiform cloud layer at the top of the cloud layer. Although many elements of the turbulent structure (including the wind profiles, the evolution of cloud-base height, the statistics of the subcloud layer, and the nature of mixing in the lower and middle parts of the cloud layer) are robustly predicted, the representation of the stratiform cloud amount by the different simulations is remarkably sensitive to a number of factors. Chief among these are differences between numerical algorithms. These sensitivities persist even among simulations on relatively refined grid meshes. Part of this sensitivity is attributed to a physically realistic positive radiative feedback, whereby a propensity toward higher cloud fractions in any given simulation is amplified by longwave radiative cooling. The simulations also provide new insight into the dynamics of the transition layer at cloud base. In accord with observations, the simulations predict that this layer is most identifiable in terms of moisture variances and gradients. The simulations help illustrate the highly variable (in both height and thickness) nature of the transition layer, and we speculate that this variability helps regulate convection. Lastly the simulations are used to help evaluate simple models of trade wind boundary layers. In accord with previous studies, mass-flux models well represent the dynamics of the cloud layer, while mixing-length models well represent the subcloud layer. The development of the stratiform cloud layer is not, however, captured by the mass-flux models. The simulations indicate that future theoretical research needs to focus on interface rules, whereby the cloud layer is coupled to the subcloud layer below and the free atmosphere above. Future observational studies of this regime would be of most benefit if they could provide robust cloud statistics as a function of mean environmental conditions.

In many large-scale models mass-flux parameterizations are applied to prognose the effect of cumulus cloud convection on the large-scale environment. Key parameters in the mass-flux equations are the lateral entrainment and detrainment rates. The physical meaning of these parameters is that they quantify the mixing rate of mass across the thermal boundaries between the cloud and its environment. The prognostic equations for the updraft and downdraft value of a conserved variable are used to derive a prognostic variance equation in the mass-flux approach. The analogy between this equation and the Reynolds-averaged variance equation is discussed. It is demonstrated that the prognostic variance equation formulated in mass-flux variables contains a gradient-production, transport, and dissipative term. In the latter term, the sum of the lateral entrainment and detrainment rates represents an inverse timescale of the dissipation. Steady-state solutions of the variance equations are discussed. Expressions for the fractional entrainment and detrainment coefficients are derived. Also, solutions for the vertical flux of an arbitrary conserved variable are presented. For convection in which the updraft fraction equals the downdraft fraction, the vertical flux of the scalar flows down the local mean gradient. The turbulent mixing coefficient is given by the ratio of the vertical mass flux and the sum of the fractional entrainment and detrainment coefficients. For an arbitrary updraft fraction, however, flux correction terms are part of the solution. It is shown that for a convective boundary layer these correction terms can account for countergradient transport, which is illustrated from large eddy simulation results. In the cumulus convection limit the vertical flux flows down the 'cloud' gradient. It is concluded that in the mass-flux approach the turbulent mixing coefficients, and the correction terms that arise from the transport term, are very similar to closures applied to the Reynolds-averaged equations.