Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cloud Layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA)

With their extensive coverage, marine low clouds greatly impact global climate. Presently, marine low clouds are poorly represented in global climate models, and the response of marine low clouds to changes in atmospheric greenhouse gases and aerosols remains the major source of uncertainty in climate simulations. The Eastern North Atlantic (ENA) is a region of persistent but diverse subtropical marine boundary layer clouds, whose albedo and precipitation are highly susceptible to perturbations in aerosol properties. In addition, the ENA is periodically impacted by continental aerosols, making it an excellent location to study the cloud condensation nuclei (CCN) budget in a remote marine region periodically perturbed by anthropogenic emissions, and to investigate the impacts of long-range transport of aerosols on remote marine clouds. The Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) campaign was motivated by the need of comprehensive in-situ measurements for improving the understanding of marine boundary layer CCN budget, cloud and drizzle microphysics, and the impact of aerosol on marine low cloud and precipitation. The airborne deployments took place from June 21 to July 20, 2017 and January 15 to February 18, 2018 in the Azores. The flights were designed to maximize the synergy between in-situ airborne measurements and ongoing long-term observations at a ground site. Here we present measurements, observation strategy, meteorological conditions during the campaign, and preliminary findings. Lastly, we discuss future analyses and modeling studies that improve the understanding and representation of marine boundary layer aerosols, clouds, precipitation, and the interactions among them.

54 ENVIRONMENTAL SCIENCES↗

Contrasting characteristics of open- and closed-cellular stratocumulus cloud in the eastern North Atlantic

Abstract. Extensive regions of marine boundary layer cloud impact the radiative balance through their significant shortwave albedo while having little impact on outgoing longwave radiation. Despite this importance, these cloud systems remain poorly represented in large-scale models due to difficulty in representing the processes that drive their life cycle and coverage. In particular, the mesoscale organization and cellular structure of marine boundary clouds have important implications for the subsequent cloud feedbacks. In this study, we use long-term (2013–2018) observations from the Atmospheric Radiation Measurement (ARM) Facility's Eastern North Atlantic (ENA) site on Graciosa Island, Azores, Portugal, to identify cloud cases with open- or closed-cellular organization. More than 500 h of each organization type are identified. The ARM observations are combined with reanalysis and satellite products to quantify the cloud, precipitation, aerosol, thermodynamic, and large-scale synoptic characteristics associated with these cloud types. Our analysis shows that both cloud organization populations occur during similar sea surface temperature conditions, but the open-cell cases are distinguished by stronger cold-air advection and large-scale subsidence compared to the closed-cell cases, consistent with their formation during cold-air outbreaks. We also find that the open-cell cases were associated with deeper boundary layers, stronger low-level winds, and higher rain rates compared to their closed-cell counterparts. Finally, raindrops with diameters larger than 1 mm were routinely recorded at the surface during both populations, with a higher number of large drops during the open-cellular cases. The similarities and differences noted herein provide important insights into the environmental and cloud characteristics during varying marine boundary layer cloud mesoscale organization and will be useful for the evaluation of model simulations for ENA marine clouds.

54 ENVIRONMENTAL SCIENCES↗

Identifying synoptic controls on boundary layer thermodynamic and cloud properties in a regional forecast model

Although most of our understanding of boundary layer cloudiness is based on idealized, subtropical, barotropic marine environments, boundary layer clouds exist across a range of conditions. In this study, we use the Naval Research Laboratory's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) and an automated cold-front-relative analysis framework to explore the boundary layer structure associated with low clouds across a transect through the cold front of a midlatitude synoptic cyclone. The model credibly captures boundary layer structure in line with conceptual models of midlatitude cyclones and ground-based observations at Graciosa Island in the Azores. The warm sector is conditionally unstable, with clouds that are too shallow and with too little liquid water, compared to cloud property retrievals from satellite and surface-based instruments. The cold-frontal region exhibits convection associated with weak stability and ascent. Northwest of the cold front, the boundary layer is well mixed, deeper, and capped by a strong inversion maintained by large-scale subsidence. Simulated clouds in frontal and post-frontal regions are mostly too thick, with too much liquid water and too little cloud base drizzle. The post-frontal clouds are associated with grid-scale updrafts, which appear to be the model's attempt to represent mesoscale organization of cellular convection typically observed in the cold sector of midlatitude cyclones. The deep, well-mixed post-frontal boundary layers and cloudiness are maintained by strong surface fluxes, as in cold air outbreaks. Our analysis framework serves as a unique approach to model verification, and our results offer insights into differences in boundary layer cloud behavior between subtropical and synoptic cold-sector regimes.

Eissner, Jordan M. [Univ. of Kansas, Lawrence, KS ↗

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗

Cloud water adjustments to aerosol perturbations are buffered by solar heating in non-precipitating marine stratocumuli

Abstract. Marine low-level clouds are key to the Earth's energy budget due to their expansive coverage over global oceans and their high reflectance of incoming solar radiation. Their responses to anthropogenic aerosol perturbations remain the largest source of uncertainty in estimating the anthropogenic radiative forcing of climate. A major challenge is the quantification of the cloud water response to aerosol perturbations. In particular, the presence of feedbacks through microphysical, dynamical, and thermodynamical pathways at various spatial and temporal scales could augment or weaken the response. Central to this problem is the temporal evolution in cloud adjustment, governed by entangled feedback mechanisms. We apply an innovative conditional Monte Carlo subsampling approach to a large ensemble of diurnal large-eddy simulation of non-precipitating marine stratocumulus to study the role of solar heating in governing the evolution in the relationship between droplet number and cloud water. We find a persistent negative trend in this relationship at night, confirming that the role of microphysically enhanced cloud-top entrainment. After sunrise, the evolution in this relationship appears buffered and converges to ∼-0.2 in the late afternoon. This buffering effect is attributed to a strong dependence of cloud-layer shortwave absorption on cloud liquid water path. These diurnal cycle characteristics further demonstrate a tight connection between cloud brightening potential and the relationship between cloud water and droplet number at sunrise, which has implications for the impact of the timing of advertent aerosol perturbations.

Zhang, Jianhao (ORCID:0000000169882935)↗

Advancing the quantification of aerosol-cloud interactions with the CALIPSO-CloudSat-Aqua/MODIS record

Aerosol-cloud-precipitation interactions are assessed over the non-polar ocean using more than 11 years of combined Aqua-MODIS, CALIPSO-CALIOP, and CloudSat products. The analysis first shows the benefit of incorporating vertically resolved aerosol extinction coefficient (σext) in aerosol-cloud interactions (ACI) assessments, demonstrating that: σext vertically collocated with the cloud layer () correlates best with cloud droplet number concentration (Nd), column-integrated aerosol optical depth (AOD) cannot explain the Nd variability in the extratropics, and the S-shape of the AOD-Nd relationship reported in previous studies is not replicated when using instead of AOD, with a Nd- linearity more consistent with in-situ studies over the ocean. ACI metric, estimated as the log-scale regression between CALIOP and MODIS Nd reveals that the eastern Pacific is the region with the strongest ACI, followed by the Southern Ocean. The susceptibility of clouds to changes in their liquid water path (LWP) and frequency of precipitation followed a 2-step calculation by combining the Nd- regression (ACI) with the regression between these macrophysical variables and Nd. LWP susceptibility is negative (LWP decreases with aerosol loading) and statistically significant over the eastern Pacific, eastern Atlantic, and extratropics. In contrast, vast areas of the tropical and subtropical ocean feature negligible changes in LWP with aerosol. Precipitation frequency susceptibility is negative, but the values are only significant over the coastal eastern Pacific and Atlantic. The findings suggest that previous modeling assessments relying on AOD may need to be revisited by taking advantage of the synergy between passive and active sensors.

Li, Zhujun↗

The prevalence of precipitation from polar supercooled clouds

Abstract. Supercooled clouds substantially impact polar surface energy budgets, but large-scale models often underestimate their occurrence, which motivates accurately establishing metrics of basic processes. An analysis of long-term measurements at Utqiaġvik, Alaska, and McMurdo Station, Antarctica, combines lidar-validated use of soundings to identify supercooled cloud layers and colocated ground-based profiling radar measurements to quantify cloud base precipitation. We find that more than 85 % (75 %) of sampled supercooled layers are precipitating over the Arctic (Antarctic) site, with more than 75 % (50 %) precipitating continuously to the surface. Such high frequencies can be reconciled with substantially lesser spaceborne estimates by considering differences in radar hydrometeor detection sensitivity. While ice precipitation into supercooled clouds from aloft is common, we also find that the great majority of supercooled cloud layers without ice falling into them are themselves continuously generating precipitation. Such sustained primary ice formation is consistent with continuous activation of immersion-mode ice-nucleating particles (INPs), suggesting that supercooled cloud formation is a principal gateway to ice formation at temperatures greater than ∼-38 ∘C over polar regions. The prevalence of weak precipitation fluxes is also consistent with supercooled cloud longevity and with well-observed and widely simulated case studies. An analysis of colocated microwave radiometer retrievals suggests that weak precipitation fluxes can be nonetheless consequential to moisture budgets for supercooled clouds owing to small liquid water paths. The results here also demonstrate that the observed abundance of mixed-phase clouds can vary substantially with instrument sensitivity and methodology. Finally, we suggest that these ground-based precipitation rate statistics offer valuable guidance for improving the representation of polar cloud processes in large-scale models.

54 ENVIRONMENTAL SCIENCES↗

Distinct dynamical and structural properties of marine stratocumulus and shallow cumulus clouds in the Eastern North Atlantic

Low-level clouds in a marine environment are examined using multiple ground-based instruments. Here, we used Doppler radar and Doppler lidar measurements (1) to examine cloud properties including cloud base and top height and thickness and (2) to investigate vertical velocity (w) statistics throughout the subcloud and in-cloud layers. Marine stratocumulus (Sc) and shallow cumulus (Cu) clouds have distinct properties in their horizontal and vertical extent. We also find that the two cloud types have different thermodynamic (coupled and decoupled conditions) and dynamical characteristics for the first four moments of the w distributions (mean, variance, skewness, and kurtosis). We identify a strong dependence of subcloud mean w on near-surface horizontal wind direction, implying that topographic effects exert a substantial impact on boundary layer dynamics and cloud properties at Graciosa Island. Our findings strongly suggest that care must be taken when interpreting island-based observations as representative of a generic marine dynamical environment. This information also can be useful for guiding and constraining the future global model development using observational dataset.

54 ENVIRONMENTAL SCIENCES↗

Framework for idealized climate simulations with spatiotemporal stochastic clouds and planetary-scale circulations

In climate predictions, clouds are the leading source of uncertainty. This is partly because, to simulate the fluid dynamics of climate over the entire globe, a large grid spacing must be used, so clouds are a subgrid-scale parametrization rather than a resolved feature. Here, a framework is investigated with finer grid spacing of O(1) or O(10)km so that some clouds are not subgrid-scale; instead, clouds evolve on the numerical grid. This cloud evolution is achieved using stochastic modeling. Hence, the framework is idealized in the sense that the full fluid dynamics of cloud circulations is still not resolved, and simplified vertical structures are used. Nevertheless, the fluid dynamics model includes evolving clouds that interactively adjust in size, shape, lifetime, and regional coverage. In addition, different cloud types are included with different roles in the climate system, including deep convective clouds and also boundary-layer clouds such as shallow cumulus and stratocumulus clouds. Other basic aspects of the idealized climate system are planetary-scale circulations (e.g., Walker circulation) and radiation. With these ingredients (evolving clouds, planetary-scale circulations, and radiation), the framework has the potential for idealized investigations of climate change with interactive cloud–radiative feedback of individual clouds. Here, the formulation of the model equations is presented, and numerical simulations are shown to illustrate the model dynamics and climate change.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Impact of Slantwise Descending Dry Intrusions on the Marine Boundary Layer and Air-Sea Interface Over the ARM Eastern North Atlantic Site

Dry intrusions (DIs) are synoptic-scale slantwise descending airstreams from the midlatitude upper troposphere toward the boundary layer at lower latitudes. Typically occurring behind cold fronts, such intrusions of dry air often reach the boundary layer and cause its deepening, thereby affecting boundary-layer clouds. Although subsidence is generally an inherent feature of the subtropical marine boundary layer (MBL), it is unclear how the MBL reacts to the transient, dynamically distinct DI. In this study, reanalysis data were combined with observations from the Atmospheric Radiation Measurement Eastern North Atlantic (ENA) site (39.1 degrees N,28.0 degrees W) to characterize the impact of DIs on MBL characteristics and clouds. Specifically, an objective criterion is applied to the observations made during the winter months of 2016-2018 to identify the DI days from those before and following DIs, and reference periods without DIs. The analyses suggest substantial deepening of the well-mixed boundary layer accompanied by changes in the cloud, precipitation and thermodynamic properties during the DI events. During the DI, the lower troposphere cooled and dried substantially thereby inducing strong surface sensible and latent heat fluxes. All while a strong inversion builds up at the elevated MBL top affecting cloud occurrence. The results show DIs to affect the boundary layer and cloud structure at the ENA site similar to 21% of the time in winter months, with the response of the cloud fields to the DI-fronts substantially different than that to the non-DI fronts. Hence, the DI events should be considered while studying boundary layer and cloud processes in the region.

54 ENVIRONMENTAL SCIENCES↗

Mixed‐Phase Clouds Over the Southern Ocean as Observed From Satellite and Surface Based Lidar and Radar

Abstract This study investigates the occurrence of mixed‐phase clouds (MPC, i.e., cloud layers containing both liquid and ice water at sub‐freezing temperatures) over the Southern Ocean (SO) using space‐ and surface‐based lidar and radar observations. The occurrence of supercooled clouds is dominated by geometrically thin (<1 km) layers that rarely contain ice. We diagnose layers that are geometrically thicker than 1 km to contain ice ~65% and ~4% of the time from below by surface remote sensors and from above by orbiting remote sensors, respectively. We examine the discrepancy in MPC occurrence statistics as diagnosed from below and above the cloud layer. From above, we find that MPC occurrence has a gradient associated with the Antarctic Polar Front near 55°S with a rare occurrence of satellite‐derived MPC south of that latitude. In contrast, surface sensors find ice in 33% of supercooled liquid water layers. We infer using observing system simulation experiments and data analysis that space‐based lidar cannot identify the occurrence of MPC except when secondary ice‐forming processes operate in convection that is, sufficiently strong to loft ice crystals to cloud tops. We conclude that the CALIPSO phase statistics of MPC have a severe low bias in MPC occurrence. Based on surface‐based statistics in the SO, we present a parameterization of the frequency of MPC as a function of cloud top temperature that differs substantially from that used in recent climate model simulations.

54 ENVIRONMENTAL SCIENCES↗

Toward Understanding the Simulated Phase Partitioning of Arctic Single-Layer Mixed-Phase Clouds in E3SM

21 Mixed-phase clouds are widespread in the Arctic and have substantial impacts on the 22 surface energy budget, sea ice extent, and regional climate. Significant changes are found in the 23 modeled phase partitioning of Arctic mixed-phase clouds in the newly developed U.S. 24 Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model 25 version 1 (EAMv1) compared to its predecessor (CAM5/EAMv0). The mixed-phase clouds 26 simulated by EAMv1 is overly dominated by super-cool liquid which is often underestimated by 27 CAM5. Compared to CAM5, EAMv1 has made several important changes in its cloud physical 28 parameterizations, including using CLUBB to unify shallow convection, PBL turbulence, and 29 cloud macrophysics, as well as updating the MG two-moment cloud microphysics scheme and 30 using the CNT heterogeneous ice nucleation scheme to replace the Meyers scheme for mixed- 31 phase clouds. To understand which of the newly-implemented physical parameterizations is 32 responsible for the change of modeled mixed-phase cloud phase partitioning in EAMv1, 33 sensitivity tests are conducted by utilizing the short-term hindcast approach and results are 34 compared to the observations from the ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) 35 field campaign. 36 It is found that although each scheme contributes to the decrease of cloud ice water 37 content (IWC) and increase of cloud liquid water content (LWC), only the combined effect of 38 CLUBB and CNT can explain the significant underestimation of IWC. By examining the budgets 39 of cloud microphysical processes, CLUBB is found to play a dominant role due to the fact that 40 CLUBB does not diagnoses cloud ice phase water in its assumed joint PDF method and thus may 41 not provide sufficient cloud ice to initialize ice phase microphysical processes such as the growth 42 of ice crystals through vapor deposition. Although a more effective heterogeneous ice nucleation 43 may alleviate the underestimation of IWC, the ice phase process in CLUBB requires urgent 44 attention in the future model development.

Zhang, Meng↗

An Improved Algorithm for Estimating Surface Shortwave Radiation: Preliminary Evaluation With MODIS Products

Cloud parameters, as key inputs in radiative transfer algorithms, have a critical impact on surface shortwave radiation (SSR) computation. By introducing a parameterization of cloud transmittance and reflectance, based on radiative transfer simulations, this study improves the accuracy of an existing physically based model which severely underestimates SSR under thick cloud conditions. The cloud parameterization adopts the single-layer cloud model and simulates cloud transmittances and reflectances by varying cloud optical thickness, cloud particle size, and solar zenith angle. The revised model is applied to estimate instantaneous SSR using Moderate-resolution Imaging Spectroradiometer (MODIS) atmospheric and land products. The retrieved SSR is evaluated against observation data from 41 Baseline Surface Radiation Network (BSRN) stations and is also compared with the MODIS official SSR product. The root mean square error (RMSE) of the estimated instantaneous radiation is approximately 52 and 98 W m -2 under clear-sky and all-sky conditions, respectively. The accuracy of the improved parameterization is higher than that of the original model, and there is no obvious underestimation of SSR in the case of high cloud optical thickness. Therefore, the new algorithm improves the accuracy of SSR estimates in the presence of thick clouds. Retrievals with the improved model also achieve higher accuracy than the MODIS official SSR product (MCD18A1). To conclude, the reliable performance of the scheme at most BSRN stations illustrates that the improved model can be used to map SSR on a global scale.

42 ENGINEERING↗

Realism of Lagrangian Large Eddy Simulations Driven by Reanalysis Meteorology: Tracking a Pocket of Open Cells Under a Biomass Burning Aerosol Layer

An approach to drive Lagrangian large eddy simulation (LES) of boundary layer clouds with reanalysis data is presented and evaluated using satellite (Spinning Enhanced Visible and Infrared Imager, SEVIRI) and aircraft (Cloud-Aerosol-Radiation Interactions and Forcing, CLARIFY) measurements. The simulations follow trajectories of the boundary layer flow. They track the formation and evolution of a pocket of open cells (POC) underneath a biomass burning aerosol layer in the free troposphere. The simulations reproduce the evolution of observed stratocumulus cloud morphology, cloud optical depth, and cloud drop effective radius, and capture the timing of the cloud state transition from closed to open cells seen in the satellite imagery on the three considered trajectories. They reproduce a biomass burning aerosol layer identified by the in-situ aircraft measurements above the inversion of the POC. Entrainment of aerosol from the biomass burning layer into the POC is limited to the extent of having no impact on cloud- or boundary layer properties, in agreement with the CLARIFY observations. The two-moment bin microphysics scheme used in the simulations reproduces the in-situ cloud microphysical properties reasonably well. A two-moment bulk microphysics scheme reproduces the satellite observations in the non-precipitating closed-cell state, but overestimates liquid water path and cloud optical depth in the precipitating open-cell state due to insufficient surface precipitation. A boundary layer cold and dry bias occurring in LES can be counteracted by reducing the grid aspect ratio and by tightening the large scale wind speed nudging towards the surface.

54 ENVIRONMENTAL SCIENCES↗

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗