Stochastic Convection Parameterizations

NASA Technical Reports Server (NTRS)

Teixeira, Joao; Reynolds, Carolyn; Suselj, Kay; Matheou, Georgios

2012-01-01

computational fluid dynamics, radiation, clouds, turbulence, convection, gravity waves, surface interaction, radiation interaction, cloud and aerosol microphysics, complexity (vegetation, biogeochemistry, radiation versus turbulence/convection stochastic approach, non-linearities, Monte Carlo, high resolutions, large-Eddy Simulations, cloud structure, plumes, saturation in tropics, forecasting, parameterizations, stochastic, radiation-clod interaction, hurricane forecasts

The Grell-Freitas Convection Parameterization: Recent Developments and Applications Within the NASA GEOS Global Model

NASA Technical Reports Server (NTRS)

Freitas, Saulo R.; Grell, Georg; Molod, Andrea; Thompson, Matthew A.

2017-01-01

We implemented and began to evaluate an alternative convection parameterization for the NASA Goddard Earth Observing System (GEOS) global model. The parameterization is based on the mass flux approach with several closures, for equilibrium and non-equilibrium convection, and includes scale and aerosol awareness functionalities. Recently, the scheme has been extended to a tri-modal spectral size approach to simulate the transition from shallow, mid, and deep convection regimes. In addition, the inclusion of a new closure for non-equilibrium convection resulted in a substantial gain of realism in model simulation of the diurnal cycle of convection over the land. Here, we briefly introduce the recent developments, implementation, and preliminary results of this parameterization in the NASA GEOS modeling system.

Improving microphysics in a convective parameterization: possibilities and limitations

NASA Astrophysics Data System (ADS)

Labbouz, Laurent; Heikenfeld, Max; Stier, Philip; Morrison, Hugh; Milbrandt, Jason; Protat, Alain; Kipling, Zak

2017-04-01

The convective cloud field model (CCFM) is a convective parameterization implemented in the climate model ECHAM6.1-HAM2.2. It represents a population of clouds within each ECHAM-HAM model column, simulating up to 10 different convective cloud types with individual radius, vertical velocities and microphysical properties. Comparisons between CCFM and radar data at Darwin, Australia, show that in order to reproduce both the convective cloud top height distribution and the vertical velocity profile, the effect of aerodynamic drag on the rising parcel has to be considered, along with a reduced entrainment parameter. A new double-moment microphysics (the Predicted Particle Properties scheme, P3) has been implemented in the latest version of CCFM and is compared to the standard single-moment microphysics and the radar retrievals at Darwin. The microphysical process rates (autoconversion, accretion, deposition, freezing, …) and their response to changes in CDNC are investigated and compared to high resolution CRM WRF simulations over the Amazon region. The results shed light on the possibilities and limitations of microphysics improvements in the framework of CCFM and in convective parameterizations in general.

Characteristics of different convective parameterization schemes on the simulation of intensity and track of severe extratropical cyclones over North Atlantic

NASA Astrophysics Data System (ADS)

Pradhan, P. K.; Liberato, Margarida L. R.; Ferreira, Juan A.; Dasamsetti, S.; Vijaya Bhaskara Rao, S.

2018-01-01

The role of the convective parameterization schemes (CPSs) in the ARW-WRF (WRF) mesoscale model is examined for extratropical cyclones (ETCs) over the North Atlantic Ocean. The simulation of very severe winter storms such as Xynthia (2010) and Gong (2013) are considered in this study. Most popular CPSs within WRF model, along with Yonsei University (YSU) planetary boundary layer (PBL) and WSM6 microphysical parameterization schemes are incorporated for the model experiments. For each storm, four numerical experiments were carried out using New Kain Fritsch (NKF), Betts-Miller-Janjic (BMJ), Grell 3D Ensemble (Gr3D) and no convection scheme (NCS) respectively. The prime objectives of these experiments were to recognize the best CPS that can forecast the intensity, track, and landfall over the Iberian Peninsula in advance of two days. The WRF model results such as central sea level pressure (CSLP), wind field, moisture flux convergence, geopotential height, jet stream, track and precipitation have shown sensitivity CPSs. The 48-hour lead simulations with BMJ schemes produce the best simulations both regarding ETCs intensity and track than Gr3D and NKF schemes. The average MAE and RMSE of intensities are least that (6.5 hPa in CSLP and 3.4 ms- 1 in the 10-m wind) found in BMJ scheme. The MAE and RMSE for and intensity and track error have revealed that NCS produces large errors than other CPSs experiments. However, for track simulation of these ETCs, at 72-, 48- and 24-hour means track errors were 440, 390 and 158 km respectively. In brevity, BMJ and Gr3D schemes can be used for short and medium range predictions of the ETCs over North Atlantic. For the evaluation of precipitation distributions using Gr3D scheme are good agreement with TRMM satellite than other CPSs.

Development and Testing of Coupled Land-surface, PBL and Shallow/Deep Convective Parameterizations within the MM5

NASA Technical Reports Server (NTRS)

Stauffer, David R.; Seaman, Nelson L.; Munoz, Ricardo C.

2000-01-01

The objective of this investigation was to study the role of shallow convection on the regional water cycle of the Mississippi and Little Washita Basins using a 3-D mesoscale model, the PSUINCAR MM5. The underlying premise of the project was that current modeling of regional-scale climate and moisture cycles over the continents is deficient without adequate treatment of shallow convection. It was hypothesized that an improved treatment of the regional water cycle can be achieved by using a 3-D mesoscale numerical model having a detailed land-surface parameterization, an advanced boundary-layer parameterization, and a more complete shallow convection parameterization than are available in most current models. The methodology was based on the application in the MM5 of new or recently improved parameterizations covering these three physical processes. Therefore, the work plan focused on integrating, improving, and testing these parameterizations in the MM5 and applying them to study water-cycle processes over the Southern Great Plains (SGP): (1) the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) described by Wetzel and Boone; (2) the 1.5-order turbulent kinetic energy (TKE)-predicting scheme of Shafran et al.; and (3) the hybrid-closure sub-grid shallow convection parameterization of Deng. Each of these schemes has been tested extensively through this study and the latter two have been improved significantly to extend their capabilities.

Radar and microphysical characteristics of convective storms simulated from a numerical model using a new microphysical parameterization

NASA Technical Reports Server (NTRS)

Ferrier, Brad S.; Tao, Wei-Kuo; Simpson, Joanne

1991-01-01

The basic features of a new and improved bulk-microphysical parameterization capable of simulating the hydrometeor structure of convective systems in all types of large-scale environments (with minimal adjustment of coefficients) are studied. Reflectivities simulated from the model are compared with radar observations of an intense midlatitude convective system. Simulated reflectivities using the novel four-class ice scheme with a microphysical parameterization rain distribution at 105 min are illustrated. Preliminary results indicate that this new ice scheme works efficiently in simulating midlatitude continental storms.

The Grell-Freitas Convective Parameterization: Recent Developments and Applications Within the NASA GEOS Global Model

NASA Astrophysics Data System (ADS)

Freitas, S.; Grell, G. A.; Molod, A.

2017-12-01

We implemented and began to evaluate an alternative convection parameterization for the NASA Goddard Earth Observing System (GEOS) global model. The parameterization (Grell and Freitas, 2014) is based on the mass flux approach with several closures, for equilibrium and non-equilibrium convection, and includes scale and aerosol awareness functionalities. Scale dependence for deep convection is implemented either through using the method described by Arakawa et al (2011), or through lateral spreading of the subsidence terms. Aerosol effects are included though the dependence of autoconversion and evaporation on the CCN number concentration.Recently, the scheme has been extended to a tri-modal spectral size approach to simulate the transition from shallow, congestus, and deep convection regimes. In addition, the inclusion of a new closure for non-equilibrium convection resulted in a substantial gain of realism in model simulation of the diurnal cycle of convection over the land. Also, a beta-pdf is employed now to represent the normalized mass flux profile. This opens up an additional venue to apply stochasticism in the scheme.

Sensitivity of U.S. summer precipitation to model resolution and convective parameterizations across gray zone resolutions

NASA Astrophysics Data System (ADS)

Gao, Yang; Leung, L. Ruby; Zhao, Chun; Hagos, Samson

2017-03-01

Simulating summer precipitation is a significant challenge for climate models that rely on cumulus parameterizations to represent moist convection processes. Motivated by recent advances in computing that support very high-resolution modeling, this study aims to systematically evaluate the effects of model resolution and convective parameterizations across the gray zone resolutions. Simulations using the Weather Research and Forecasting model were conducted at grid spacings of 36 km, 12 km, and 4 km for two summers over the conterminous U.S. The convection-permitting simulations at 4 km grid spacing are most skillful in reproducing the observed precipitation spatial distributions and diurnal variability. Notable differences are found between simulations with the traditional Kain-Fritsch (KF) and the scale-aware Grell-Freitas (GF) convection schemes, with the latter more skillful in capturing the nocturnal timing in the Great Plains and North American monsoon regions. The GF scheme also simulates a smoother transition from convective to large-scale precipitation as resolution increases, resulting in reduced sensitivity to model resolution compared to the KF scheme. Nonhydrostatic dynamics has a positive impact on precipitation over complex terrain even at 12 km and 36 km grid spacings. With nudging of the winds toward observations, we show that the conspicuous warm biases in the Southern Great Plains are related to precipitation biases induced by large-scale circulation biases, which are insensitive to model resolution. Overall, notable improvements in simulating summer rainfall and its diurnal variability through convection-permitting modeling and scale-aware parameterizations suggest promising venues for improving climate simulations of water cycle processes.

New Approaches to Parameterizing Convection

NASA Technical Reports Server (NTRS)

Randall, David A.; Lappen, Cara-Lyn

1999-01-01

Many general circulation models (GCMs) currently use separate schemes for planetary boundary layer (PBL) processes, shallow and deep cumulus (Cu) convection, and stratiform clouds. The conventional distinctions. among these processes are somewhat arbitrary. For example, in the stratocumulus-to-cumulus transition region, stratocumulus clouds break up into a combination of shallow cumulus and broken stratocumulus. Shallow cumulus clouds may be considered to reside completely within the PBL, or they may be regarded as starting in the PBL but terminating above it. Deeper cumulus clouds often originate within the PBL with also can originate aloft. To the extent that our models separately parameterize physical processes which interact strongly on small space and time scales, the currently fashionable practice of modularization may be doing more harm than good.

Simulation of heavy precipitation episode over eastern Peninsular Malaysia using MM5: sensitivity to cumulus parameterization schemes

NASA Astrophysics Data System (ADS)

Salimun, Ester; Tangang, Fredolin; Juneng, Liew

2010-06-01

A comparative study has been conducted to investigate the skill of four convection parameterization schemes, namely the Anthes-Kuo (AK), the Betts-Miller (BM), the Kain-Fritsch (KF), and the Grell (GR) schemes in the numerical simulation of an extreme precipitation episode over eastern Peninsular Malaysia using the Pennsylvania State University—National Center for Atmospheric Research Center (PSU-NCAR) Fifth Generation Mesoscale Model (MM5). The event is a commonly occurring westward propagating tropical depression weather system during a boreal winter resulting from an interaction between a cold surge and the quasi-stationary Borneo vortex. The model setup and other physical parameterizations are identical in all experiments and hence any difference in the simulation performance could be associated with the cumulus parameterization scheme used. From the predicted rainfall and structure of the storm, it is clear that the BM scheme has an edge over the other schemes. The rainfall intensity and spatial distribution were reasonably well simulated compared to observations. The BM scheme was also better in resolving the horizontal and vertical structures of the storm. Most of the rainfall simulated by the BM simulation was of the convective type. The failure of other schemes (AK, GR and KF) in simulating the event may be attributed to the trigger function, closure assumption, and precipitation scheme. On the other hand, the appropriateness of the BM scheme for this episode may not be generalized for other episodes or convective environments.

Evaluation of the Plant-Craig stochastic convection scheme in an ensemble forecasting system

NASA Astrophysics Data System (ADS)

Keane, R. J.; Plant, R. S.; Tennant, W. J.

2015-12-01

The Plant-Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic element only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant-Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant-Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.

Evaluating and Understanding Parameterized Convective Processes and their Role in the Development of Mesoscale Precipitation Systems

NASA Technical Reports Server (NTRS)

Fritsch, J. Michael; Kain, John S.

1997-01-01

Research efforts during the second year have centered on improving the manner in which convective stabilization is achieved in the Penn State/NCAR mesoscale model MM5. Ways of improving this stabilization have been investigated by (1) refining the partitioning between the Kain-Fritsch convective parameterization scheme and the grid scale by introducing a form of moist convective adjustment; (2) using radar data to define locations of subgrid-scale convection during a dynamic initialization period; and (3) parameterizing deep-convective feedbacks as subgrid-scale sources and sinks of mass. These investigations were conducted by simulating a long-lived convectively-generated mesoscale vortex that occurred during 14-18 Jul. 1982 and the 10-11 Jun. 1985 squall line that occurred over the Kansas-Oklahoma region during the PRE-STORM experiment. The long-lived vortex tracked across the central Plains states and was responsible for multiple convective outbreaks during its lifetime.

An Empirical Cumulus Parameterization Scheme for a Global Spectral Model

NASA Technical Reports Server (NTRS)

Rajendran, K.; Krishnamurti, T. N.; Misra, V.; Tao, W.-K.

2004-01-01

Realistic vertical heating and drying profiles in a cumulus scheme is important for obtaining accurate weather forecasts. A new empirical cumulus parameterization scheme based on a procedure to improve the vertical distribution of heating and moistening over the tropics is developed. The empirical cumulus parameterization scheme (ECPS) utilizes profiles of Tropical Rainfall Measuring Mission (TRMM) based heating and moistening derived from the European Centre for Medium- Range Weather Forecasts (ECMWF) analysis. A dimension reduction technique through rotated principal component analysis (RPCA) is performed on the vertical profiles of heating (Q1) and drying (Q2) over the convective regions of the tropics, to obtain the dominant modes of variability. Analysis suggests that most of the variance associated with the observed profiles can be explained by retaining the first three modes. The ECPS then applies a statistical approach in which Q1 and Q2 are expressed as a linear combination of the first three dominant principal components which distinctly explain variance in the troposphere as a function of the prevalent large-scale dynamics. The principal component (PC) score which quantifies the contribution of each PC to the corresponding loading profile is estimated through a multiple screening regression method which yields the PC score as a function of the large-scale variables. The profiles of Q1 and Q2 thus obtained are found to match well with the observed profiles. The impact of the ECPS is investigated in a series of short range (1-3 day) prediction experiments using the Florida State University global spectral model (FSUGSM, T126L14). Comparisons between short range ECPS forecasts and those with the modified Kuo scheme show a very marked improvement in the skill in ECPS forecasts. This improvement in the forecast skill with ECPS emphasizes the importance of incorporating realistic vertical distributions of heating and drying in the model cumulus scheme. This

Evaluating and Understanding Parameterized Convective Processes and Their Role in the Development of Mesoscale Precipitation Systems

NASA Technical Reports Server (NTRS)

Fritsch, J. Michael; Kain, John S.

1996-01-01

Research efforts focused on numerical simulations of two convective systems with the Penn State/NCAR mesoscale model. The first of these systems was tropical cyclone Irma, which occurred in 1987 in Australia's Gulf of Carpentaria during the AMEX field program. Comparison simulations of this system were done with two different convective parameterization schemes (CPS's), the Kain-Fritsch (KF) and the Betts-Miller (BM) schemes. The second system was the June 10-11, 1985 squall line simulation, which occurred over the Kansas-Oklahoma region during the PRE-STORM experiment. Simulations of this system using the KF scheme were examined in detail.

A Generalized Simple Formulation of Convective Adjustment Timescale for Cumulus Convection Parameterizations

EPA Science Inventory

Convective adjustment timescale (τ) for cumulus clouds is one of the most influential parameters controlling parameterized convective precipitation in climate and weather simulation models at global and regional scales. Due to the complex nature of deep convection, a pres...

Bias Reduction as Guidance for Developing Convection and Cloud Parameterization in GFDL AM4/CM4

NASA Astrophysics Data System (ADS)

Zhao, M.; Held, I.; Golaz, C.

2016-12-01

The representations of moist convection and clouds are challenging in global climate models and they are known to be important to climate simulations at all spatial and temporal scales. Many climate simulation biases can be traced to deficiencies in convection and cloud parameterizations. I will present some key biases that we are concerned about and the efforts that we have made to reduce the biases during the development of NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) new generation global climate model AM4/CM4. In particular, I will present a modified version of the moist convection scheme that is based on the University of Washington Shallow Cumulus scheme (UWShCu, Bretherton et. al 2004). The new scheme produces marked improvement in simulation of the Madden-Julian Oscillation (MJO) and the El Niño-Southern Oscillation (ENSO) compared to that used in AM3 and HIRAM. AM4/CM4 also produces high quality simulation of global distribution of cloud radiative effects and the precipitation with realistic mean climate state. This differs from models of improved MJO but with a much deteriorated mean state. The modifications to the UWShCu include an additional bulk plume for representing deep convection. The entrainment rate in the deep plume is parameterized to be a function of column-integrated relative humidity. The deep convective closure is based on relaxation of the convective available potential energy (CAPE) or cloud work function. The plumes' precipitation efficiency is optimized for better simulations of the cloud radiative effects. Precipitation re-evaporation is included in both shallow and deep plumes. In addition, a parameterization of convective gustiness is included with an energy source driven by cold pool derived from precipitation re-evaporation within the boundary layer and energy sink due to dissipation. I will present the motivations of these changes which are driven by reducing some aspects of the AM4/CM4 biases. Finally, I will also present

An improvement in mass flux convective parameterizations and its impact on seasonal simulations using a coupled model

NASA Astrophysics Data System (ADS)

Elsayed Yousef, Ahmed; Ehsan, M. Azhar; Almazroui, Mansour; Assiri, Mazen E.; Al-Khalaf, Abdulrahman K.

2017-02-01

A new closure and a modified detrainment for the simplified Arakawa-Schubert (SAS) cumulus parameterization scheme are proposed. In the modified convective scheme which is named as King Abdulaziz University (KAU) scheme, the closure depends on both the buoyancy force and the environment mean relative humidity. A lateral entrainment rate varying with environment relative humidity is proposed and tends to suppress convection in a dry atmosphere. The detrainment rate also varies with environment relative humidity. The KAU scheme has been tested in a single column model (SCM) and implemented in a coupled global climate model (CGCM). Increased coupling between environment and clouds in the KAU scheme results in improved sensitivity of the depth and strength of convection to environmental humidity compared to the original SAS scheme. The new scheme improves precipitation simulation with better representations of moisture and temperature especially during suppressed convection periods. The KAU scheme implemented in the Seoul National University (SNU) CGCM shows improved precipitation over the tropics. The simulated precipitation pattern over the Arabian Peninsula and Northeast African region is also improved.

Sensitivity of Tropical Cyclones to Parameterized Convection in the NASA GEOS5 Model

NASA Technical Reports Server (NTRS)

Lim, Young-Kwon; Schubert, Siegfried D.; Reale, Oreste; Lee, Myong-In; Molod, Andrea M.; Suarez, Max J.

2014-01-01

The sensitivity of tropical cyclones (TCs) to changes in parameterized convection is investigated to improve the simulation of TCs in the North Atlantic. Specifically, the impact of reducing the influence of the Relaxed Arakawa-Schubert (RAS) scheme-based parameterized convection is explored using the Goddard Earth Observing System version5 (GEOS5) model at 0.25 horizontal resolution. The years 2005 and 2006 characterized by very active and inactive hurricane seasons, respectively, are selected for simulation. A reduction in parameterized deep convection results in an increase in TC activity (e.g., TC number and longer life cycle) to more realistic levels compared to the baseline control configuration. The vertical and horizontal structure of the strongest simulated hurricane shows the maximum lower-level (850-950hPa) wind speed greater than 60 ms and the minimum sea level pressure reaching 940mb, corresponding to a category 4 hurricane - a category never achieved by the control configuration. The radius of the maximum wind of 50km, the location of the warm core exceeding 10 C, and the horizontal compactness of the hurricane center are all quite realistic without any negatively affecting the atmospheric mean state. This study reveals that an increase in the threshold of minimum entrainment suppresses parameterized deep convection by entraining more dry air into the typical plume. This leads to cooling and drying at the mid- to upper-troposphere, along with the positive latent heat flux and moistening in the lower-troposphere. The resulting increase in conditional instability provides an environment that is more conducive to TC vortex development and upward moisture flux convergence by dynamically resolved moist convection, thereby increasing TC activity.

Structural and parameteric uncertainty quantification in cloud microphysics parameterization schemes

NASA Astrophysics Data System (ADS)

van Lier-Walqui, M.; Morrison, H.; Kumjian, M. R.; Prat, O. P.; Martinkus, C.

2017-12-01

Atmospheric model parameterization schemes employ approximations to represent the effects of unresolved processes. These approximations are a source of error in forecasts, caused in part by considerable uncertainty about the optimal value of parameters within each scheme -- parameteric uncertainty. Furthermore, there is uncertainty regarding the best choice of the overarching structure of the parameterization scheme -- structrual uncertainty. Parameter estimation can constrain the first, but may struggle with the second because structural choices are typically discrete. We address this problem in the context of cloud microphysics parameterization schemes by creating a flexible framework wherein structural and parametric uncertainties can be simultaneously constrained. Our scheme makes no assuptions about drop size distribution shape or the functional form of parametrized process rate terms. Instead, these uncertainties are constrained by observations using a Markov Chain Monte Carlo sampler within a Bayesian inference framework. Our scheme, the Bayesian Observationally-constrained Statistical-physical Scheme (BOSS), has flexibility to predict various sets of prognostic drop size distribution moments as well as varying complexity of process rate formulations. We compare idealized probabilistic forecasts from versions of BOSS with varying levels of structural complexity. This work has applications in ensemble forecasts with model physics uncertainty, data assimilation, and cloud microphysics process studies.

Evaluating and Understanding Parameterized Convective Processes and Their Role in the Development of Mesoscale Precipitation Systems

NASA Technical Reports Server (NTRS)

Fritsch, J. Michael (Principal Investigator); Kain, John S.

1995-01-01

Research efforts during the first year focused on numerical simulations of two convective systems with the Penn State/NCAR mesoscale model. The first of these systems was tropical cyclone Irma, which occurred in 1987 in Australia's Gulf of Carpentaria during the AMEX field program. Comparison simulations of this system were done with two different convective parameterization schemes (CPS's), the Kain-Fritsch (1993 - KF) and the Betts-Miller (Betts 1986- BM) schemes. The second system was the June 10-11 1985 squall line simulation, which occurred over the Kansas-Oklahoma region during the PRE-STORM experiment. Simulations of this system using the KF scheme were examined in detail.

Spectral cumulus parameterization based on cloud-resolving model

NASA Astrophysics Data System (ADS)

Baba, Yuya

2018-02-01

We have developed a spectral cumulus parameterization using a cloud-resolving model. This includes a new parameterization of the entrainment rate which was derived from analysis of the cloud properties obtained from the cloud-resolving model simulation and was valid for both shallow and deep convection. The new scheme was examined in a single-column model experiment and compared with the existing parameterization of Gregory (2001, Q J R Meteorol Soc 127:53-72) (GR scheme). The results showed that the GR scheme simulated more shallow and diluted convection than the new scheme. To further validate the physical performance of the parameterizations, Atmospheric Model Intercomparison Project (AMIP) experiments were performed, and the results were compared with reanalysis data. The new scheme performed better than the GR scheme in terms of mean state and variability of atmospheric circulation, i.e., the new scheme improved positive bias of precipitation in western Pacific region, and improved positive bias of outgoing shortwave radiation over the ocean. The new scheme also simulated better features of convectively coupled equatorial waves and Madden-Julian oscillation. These improvements were found to be derived from the modification of parameterization for the entrainment rate, i.e., the proposed parameterization suppressed excessive increase of entrainment, thus suppressing excessive increase of low-level clouds.

Sensitivity of Pacific Cold Tongue and Double-ITCZ Bias to Convective Parameterization

NASA Astrophysics Data System (ADS)

Woelfle, M.; Bretherton, C. S.; Pritchard, M. S.; Yu, S.

2016-12-01

Many global climate models struggle to accurately simulate annual mean precipitation and sea surface temperature (SST) fields in the tropical Pacific basin. Precipitation biases are dominated by the double intertropical convergence zone (ITCZ) bias where models exhibit precipitation maxima straddling the equator while only a single Northern Hemispheric maximum exists in observations. The major SST bias is the enhancement of the equatorial cold tongue. A series of coupled model simulations are used to investigate the sensitivity of the bias development to convective parameterization. Model components are initialized independently prior to coupling to allow analysis of the transient response of the system directly following coupling. These experiments show precipitation and SST patterns to be highly sensitive to convective parameterization. Simulations in which the deep convective parameterization is disabled forcing all convection to be resolved by the shallow convection parameterization showed a degradation in both the cold tongue and double-ITCZ biases as precipitation becomes focused into off-equatorial regions of local SST maxima. Simulations using superparameterization in place of traditional cloud parameterizations showed a reduced cold tongue bias at the expense of additional precipitation biases. The equatorial SST responses to changes in convective parameterization are driven by changes in near equatorial zonal wind stress. The sensitivity of convection to SST is important in determining the precipitation and wind stress fields. However, differences in convective momentum transport also play a role. While no significant improvement is seen in these simulations of the double-ITCZ, the system's sensitivity to these changes reaffirm that improved convective parameterizations may provide an avenue for improving simulations of tropical Pacific precipitation and SST.

The Roles of Convection Parameterization in the Formation of Double ITCZ Syndrome in the NCAR CESM: I. Atmospheric Processes

NASA Astrophysics Data System (ADS)

Song, Xiaoliang; Zhang, Guang J.

2018-03-01

Several improvements are implemented in the Zhang-McFarlane (ZM) convection scheme to investigate the roles of convection parameterization in the formation of double intertropical convergence zone (ITCZ) bias in the NCAR CESM1.2.1. It is shown that the prominent double ITCZ biases of precipitation, sea surface temperature (SST), and wind stress in the standard CESM1.2.1 are largely eliminated in all seasons with the use of these improvements in convection scheme. This study for the first time demonstrates that the modifications of convection scheme can eliminate the double ITCZ biases in all seasons, including boreal winter and spring. Further analysis shows that the elimination of the double ITCZ bias is achieved not by improving other possible contributors, such as stratus cloud bias off the west coast of South America and cloud/radiation biases over the Southern Ocean, but by modifying the convection scheme itself. This study demonstrates that convection scheme is the primary contributor to the double ITCZ bias in the CESM1.2.1, and provides a possible solution to the long-standing double ITCZ problem. The atmospheric model simulations forced by observed SST show that the original ZM convection scheme tends to produce double ITCZ bias in high SST scenario, while the modified convection scheme does not. The impact of changes in each core component of convection scheme on the double ITCZ bias in atmospheric model is identified and further investigated.

A new approach to the convective parameterization of the regional atmospheric model BRAMS

NASA Astrophysics Data System (ADS)

Dos Santos, A. F.; Freitas, S. R.; de Campos Velho, H. F.; Luz, E. F.; Gan, M. A.; de Mattos, J. Z.; Grell, G. A.

2013-05-01

The summer characteristics of January 2010 was performed using the atmospheric model Brazilian developments on the Regional Atmospheric Modeling System (BRAMS). The convective parameterization scheme of Grell and Dévényi was used to represent clouds and their interaction with the large scale environment. As a result, the precipitation forecasts can be combined in several ways, generating a numerical representation of precipitation and atmospheric heating and moistening rates. The purpose of this study was to generate a set of weights to compute a best combination of the hypothesis of the convective scheme. It is an inverse problem of parameter estimation and the problem is solved as an optimization problem. To minimize the difference between observed data and forecasted precipitation, the objective function was computed with the quadratic difference between five simulated precipitation fields and observation. The precipitation field estimated by the Tropical Rainfall Measuring Mission satellite was used as observed data. Weights were obtained using the firefly algorithm and the mass fluxes of each closure of the convective scheme were weighted generating a new set of mass fluxes. The results indicated the better skill of the model with the new methodology compared with the old ensemble mean calculation.

A WRF-Chem Flash Rate Parameterization Scheme and LNOx Analysis of the 29-30 May 2012 Convective Event in Oklahoma During DC3

NASA Technical Reports Server (NTRS)

Cummings, Kristin A.; Pickering, Kenneth E.; Barth, M.; Weinheimer, A.; Bela, M.; Li, Y.; Allen, D.; Bruning, E.; MacGorman, D.; Rutledge, S.;

New Layer Thickness Parameterization of Diffusive Convection

NASA Astrophysics Data System (ADS)

Zhou, Sheng-Qi; Lu, Yuan-Zheng; Guo, Shuang-Xi; Song, Xue-Long; Qu, Ling; Cen, Xian-Rong; Fer, Ilker

2017-11-01

Double-diffusion convection is one of the most important non-mechanically driven mixing processes. Its importance has been particular recognized in oceanography, material science, geology, and planetary physics. Double-diffusion occurs in a fluid in which there are gradients of two (or more) properties with different molecular diffusivities and of opposing effects on the vertical density distribution. It has two primary modes: salt finger and diffusive convection. Recently, the importance of diffusive convection has aroused more interest due to its impact to the diapycnal mixing in the interior ocean and the ice and the ice-melting in the Arctic and Antarctic Oceans. In our recent work, we constructed a length scale of energy-containing eddy and proposed a new layer thickness parameterization of diffusive convection by using the laboratory experiment and in situ observations in the lakes and oceans. The new parameterization can well describe the laboratory convecting layer thicknesses (0.01 0.1 m) and those observed in oceans and lakes (0.1 1000 m). This work was supported by China NSF Grants (41476167,41406035 and 41176027), NSF of Guangdong Province, China (2016A030311042) and the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA11030302).

Improving Representation of Convective Transport for Scale-Aware Parameterization – Part I: Convection and Cloud Properties Simulated with Spectral Bin and Bulk Microphysics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fan, Jiwen; Liu, Yi-Chin; Xu, Kuan-Man

2015-04-27

The ultimate goal of this study is to improve representation of convective transport by cumulus parameterization for meso-scale and climate models. As Part I of the study, we perform extensive evaluations of cloud-resolving simulations of a squall line and mesoscale convective complexes in mid-latitude continent and tropical regions using the Weather Research and Forecasting (WRF) model with spectral-bin microphysics (SBM) and with two double-moment bulk microphysics schemes: a modified Morrison (MOR) and Milbrandt and Yau (MY2). Compared to observations, in general, SBM gives better simulations of precipitation, vertical velocity of convective cores, and the vertically decreasing trend of radar reflectivitymore » than MOR and MY2, and therefore will be used for analysis of scale-dependence of eddy transport in Part II. The common features of the simulations for all convective systems are (1) the model tends to overestimate convection intensity in the middle and upper troposphere, but SBM can alleviate much of the overestimation and reproduce the observed convection intensity well; (2) the model greatly overestimates radar reflectivity in convective cores (SBM predicts smaller radar reflectivity but does not remove the large overestimation); and (3) the model performs better for mid-latitude convective systems than tropical system. The modeled mass fluxes of the mid latitude systems are not sensitive to microphysics schemes, but are very sensitive for the tropical case indicating strong microphysics modification to convection. Cloud microphysical measurements of rain, snow and graupel in convective cores will be critically important to further elucidate issues within cloud microphysics schemes.« less

a Cumulus Parameterization Study with Special Attention to the Arakawa-Schubert Scheme

NASA Astrophysics Data System (ADS)

Kao, Chih-Yue Jim

Arakawa and Schubert (1974) developed a cumulus parameterization scheme in a framework that conceptually divides the mutual interaction of the cumulus convection and large-scale disturbance into the categories of large -scale budget requirements and the quasi-equilibrium assumption of cloud work function. We have applied the A-S scheme through a semi-prognostic approach to two different data sets: one is for an intense tropical cloud band event; the other is for tropical composite easterly wave disturbances. Both were observed in GATE. The cloud heating and drying effects predicted by the Arakawa-Schubert scheme are found to agree rather well with the observations. However, it is also found that the Arakawa-Schubert scheme underestimates both condensation and evaporation rates substantially when compared with the cumulus ensemble model results (Soong and Tao, 1980; Tao, 1983). An inclusion of the downdraft effects, as formulated by Johnson (1976), appears to alleviate this deficiency. In order to examine how the Arakawa-Schubert scheme works in a fully prognostic problem, a simulation of the evolution and structure of the tropical cloud band, mentioned above, under the influence of an imposed large-scale low -level forcing has been made, using a two-dimensional hydrostatic model with the inclusion of the Arakawa-Schubert scheme. Basically, the model result indicates that the meso-scale convective system is driven by the excess of the convective heating derived from the Arakawa-Schubert scheme over the adiabatic cooling due to the imposed large-scale lifting and induced meso-scale upward motion. However, as the convective system develops, the adiabatic warming due to the subsidence outside the cloud cluster gradually accumulates into a secondary temperature anomaly which subsequently reduces the original temperature contrast and inhibits the further development of the convective system. A 24 hour integration shows that the model is capable of simulating many

A Wrf-Chem Flash Rate Parameterization Scheme and LNO(x) Analysis of the 29-30 May 2012 Convective Event in Oklahoma During DC3

NASA Technical Reports Server (NTRS)

Cummings, Kristin A.; Pickering, Kenneth E.; Barth, M.; Weinheimer, A.; Bela, M.; Li, Y.; Allen, D.; Bruning, E.; MacGorman, D.; Rutledge, S.;

Boreal Winter MJO Teleconnection in the Community Atmosphere Model Version 5 with the Unified Convection Parameterization

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yoo, Changhyun; Park, Sungsu; Kim, Daehyun

2015-10-01

The Madden-Julian Oscillation (MJO), the dominant mode of tropical intraseasonal variability, influences weather and climate in the extratropics through atmospheric teleconnection. In this study, two simulations using the Community Atmosphere Model version 5 (CAM5) - one with the default shallow and deep convection schemes and the other with the Unified Convection scheme (UNICON) - are employed to examine the impacts of cumulus parameterizations on the simulation of the boreal wintertime MJO teleconnection in the Northern Hemisphere. We demonstrate that the UNICON substantially improves the MJO teleconnection. When the UNICON is employed, the simulated circulation anomalies associated with the MJO bettermore » resemble the observed counterpart, compared to the simulation with the default convection schemes. Quantitatively, the pattern correlation for the 300-hPa geopotential height anomalies between the simulations and observation increases from 0.07 for the default schemes to 0.54 for the UNICON. These circulation anomalies associated with the MJO further help to enhance the surface air temperature and precipitation anomalies over North America, although room for improvement is still evident. Initial value calculations suggest that the realistic MJO teleconnection with the UNICON is not attributed to the changes in the background wind, but primarily to the improved tropical convective heating associated with the MJO.« less

Assessing the Resolution Adaptability of the Zhang-McFarlane Cumulus Parameterization With Spatial and Temporal Averaging: RESOLUTION ADAPTABILITY OF ZM SCHEME

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yun, Yuxing; Fan, Jiwen; Xiao, Heng

Realistic modeling of cumulus convection at fine model resolutions (a few to a few tens of km) is problematic since it requires the cumulus scheme to adapt to higher resolution than they were originally designed for (~100 km). To solve this problem, we implement the spatial averaging method proposed in Xiao et al. (2015) and also propose a temporal averaging method for the large-scale convective available potential energy (CAPE) tendency in the Zhang-McFarlane (ZM) cumulus parameterization. The resolution adaptability of the original ZM scheme, the scheme with spatial averaging, and the scheme with both spatial and temporal averaging at 4-32more » km resolution is assessed using the Weather Research and Forecasting (WRF) model, by comparing with Cloud Resolving Model (CRM) results. We find that the original ZM scheme has very poor resolution adaptability, with sub-grid convective transport and precipitation increasing significantly as the resolution increases. The spatial averaging method improves the resolution adaptability of the ZM scheme and better conserves the total transport of moist static energy and total precipitation. With the temporal averaging method, the resolution adaptability of the scheme is further improved, with sub-grid convective precipitation becoming smaller than resolved precipitation for resolution higher than 8 km, which is consistent with the results from the CRM simulation. Both the spatial distribution and time series of precipitation are improved with the spatial and temporal averaging methods. The results may be helpful for developing resolution adaptability for other cumulus parameterizations that are based on quasi-equilibrium assumption.« less

An Extended Eddy-Diffusivity Mass-Flux Scheme for Unified Representation of Subgrid-Scale Turbulence and Convection

NASA Astrophysics Data System (ADS)

Tan, Zhihong; Kaul, Colleen M.; Pressel, Kyle G.; Cohen, Yair; Schneider, Tapio; Teixeira, João.

2018-03-01

Large-scale weather forecasting and climate models are beginning to reach horizontal resolutions of kilometers, at which common assumptions made in existing parameterization schemes of subgrid-scale turbulence and convection—such as that they adjust instantaneously to changes in resolved-scale dynamics—cease to be justifiable. Additionally, the common practice of representing boundary-layer turbulence, shallow convection, and deep convection by discontinuously different parameterizations schemes, each with its own set of parameters, has contributed to the proliferation of adjustable parameters in large-scale models. Here we lay the theoretical foundations for an extended eddy-diffusivity mass-flux (EDMF) scheme that has explicit time-dependence and memory of subgrid-scale variables and is designed to represent all subgrid-scale turbulence and convection, from boundary layer dynamics to deep convection, in a unified manner. Coherent up and downdrafts in the scheme are represented as prognostic plumes that interact with their environment and potentially with each other through entrainment and detrainment. The more isotropic turbulence in their environment is represented through diffusive fluxes, with diffusivities obtained from a turbulence kinetic energy budget that consistently partitions turbulence kinetic energy between plumes and environment. The cross-sectional area of up and downdrafts satisfies a prognostic continuity equation, which allows the plumes to cover variable and arbitrarily large fractions of a large-scale grid box and to have life cycles governed by their own internal dynamics. Relatively simple preliminary proposals for closure parameters are presented and are shown to lead to a successful simulation of shallow convection, including a time-dependent life cycle.

An Extended Eddy‐Diffusivity Mass‐Flux Scheme for Unified Representation of Subgrid‐Scale Turbulence and Convection

PubMed Central

Tan, Zhihong; Kaul, Colleen M.; Pressel, Kyle G.; Cohen, Yair; Teixeira, João

2018-01-01

Abstract Large‐scale weather forecasting and climate models are beginning to reach horizontal resolutions of kilometers, at which common assumptions made in existing parameterization schemes of subgrid‐scale turbulence and convection—such as that they adjust instantaneously to changes in resolved‐scale dynamics—cease to be justifiable. Additionally, the common practice of representing boundary‐layer turbulence, shallow convection, and deep convection by discontinuously different parameterizations schemes, each with its own set of parameters, has contributed to the proliferation of adjustable parameters in large‐scale models. Here we lay the theoretical foundations for an extended eddy‐diffusivity mass‐flux (EDMF) scheme that has explicit time‐dependence and memory of subgrid‐scale variables and is designed to represent all subgrid‐scale turbulence and convection, from boundary layer dynamics to deep convection, in a unified manner. Coherent up and downdrafts in the scheme are represented as prognostic plumes that interact with their environment and potentially with each other through entrainment and detrainment. The more isotropic turbulence in their environment is represented through diffusive fluxes, with diffusivities obtained from a turbulence kinetic energy budget that consistently partitions turbulence kinetic energy between plumes and environment. The cross‐sectional area of up and downdrafts satisfies a prognostic continuity equation, which allows the plumes to cover variable and arbitrarily large fractions of a large‐scale grid box and to have life cycles governed by their own internal dynamics. Relatively simple preliminary proposals for closure parameters are presented and are shown to lead to a successful simulation of shallow convection, including a time‐dependent life cycle. PMID:29780442

The Stochastic Parcel Model: A deterministic parameterization of stochastically entraining convection

DOE PAGES

Romps, David M.

2016-03-01

Convective entrainment is a process that is poorly represented in existing convective parameterizations. By many estimates, convective entrainment is the leading source of error in global climate models. As a potential remedy, an Eulerian implementation of the Stochastic Parcel Model (SPM) is presented here as a convective parameterization that treats entrainment in a physically realistic and computationally efficient way. Drawing on evidence that convecting clouds comprise air parcels subject to Poisson-process entrainment events, the SPM calculates the deterministic limit of an infinite number of such parcels. For computational efficiency, the SPM groups parcels at each height by their purity, whichmore » is a measure of their total entrainment up to that height. This reduces the calculation of convective fluxes to a sequence of matrix multiplications. The SPM is implemented in a single-column model and compared with a large-eddy simulation of deep convection.« less

ARM - Midlatitude Continental Convective Clouds

DOE Data Explorer

Jensen, Mike; Bartholomew, Mary Jane; Genio, Anthony Del; Giangrande, Scott; Kollias, Pavlos

2012-01-19

Convective processes play a critical role in the Earth's energy balance through the redistribution of heat and moisture in the atmosphere and their link to the hydrological cycle. Accurate representation of convective processes in numerical models is vital towards improving current and future simulations of Earths climate system. Despite improvements in computing power, current operational weather and global climate models are unable to resolve the natural temporal and spatial scales important to convective processes and therefore must turn to parameterization schemes to represent these processes. In turn, parameterization schemes in cloud-resolving models need to be evaluated for their generality and application to a variety of atmospheric conditions. Data from field campaigns with appropriate forcing descriptors have been traditionally used by modelers for evaluating and improving parameterization schemes.

Examining Chaotic Convection with Super-Parameterization Ensembles

NASA Astrophysics Data System (ADS)

Jones, Todd R.

This study investigates a variety of features present in a new configuration of the Community Atmosphere Model (CAM) variant, SP-CAM 2.0. The new configuration (multiple-parameterization-CAM, MP-CAM) changes the manner in which the super-parameterization (SP) concept represents physical tendency feedbacks to the large-scale by using the mean of 10 independent two-dimensional cloud-permitting model (CPM) curtains in each global model column instead of the conventional single CPM curtain. The climates of the SP and MP configurations are examined to investigate any significant differences caused by the application of convective physical tendencies that are more deterministic in nature, paying particular attention to extreme precipitation events and large-scale weather systems, such as the Madden-Julian Oscillation (MJO). A number of small but significant changes in the mean state climate are uncovered, and it is found that the new formulation degrades MJO performance. Despite these deficiencies, the ensemble of possible realizations of convective states in the MP configuration allows for analysis of uncertainty in the small-scale solution, lending to examination of those weather regimes and physical mechanisms associated with strong, chaotic convection. Methods of quantifying precipitation predictability are explored, and use of the most reliable of these leads to the conclusion that poor precipitation predictability is most directly related to the proximity of the global climate model column state to atmospheric critical points. Secondarily, the predictability is tied to the availability of potential convective energy, the presence of mesoscale convective organization on the CPM grid, and the directive power of the large-scale.

The sensitivity of Alpine summer convection to surrogate climate change: an intercomparison between convection-parameterizing and convection-resolving models

NASA Astrophysics Data System (ADS)

Keller, Michael; Kröner, Nico; Fuhrer, Oliver; Lüthi, Daniel; Schmidli, Juerg; Stengel, Martin; Stöckli, Reto; Schär, Christoph

2018-04-01

Climate models project an increase in heavy precipitation events in response to greenhouse gas forcing. Important elements of such events are rain showers and thunderstorms, which are poorly represented in models with parameterized convection. In this study, simulations with 12 km horizontal grid spacing (convection-parameterizing model, CPM) and 2 km grid spacing (convection-resolving model, CRM) are employed to investigate the change in the diurnal cycle of convection with warmer climate. For this purpose, simulations of 11 days in June 2007 with a pronounced diurnal cycle of convection are compared with surrogate simulations from the same period. The surrogate climate simulations mimic a future climate with increased temperatures but unchanged relative humidity and similar synoptic-scale circulation. Two temperature scenarios are compared: one with homogeneous warming (HW) using a vertically uniform warming and the other with vertically dependent warming (VW) that enables changes in lapse rate. The two sets of simulations with parameterized and explicit convection exhibit substantial differences, some of which are well known from the literature. These include differences in the timing and amplitude of the diurnal cycle of convection, and the frequency of precipitation with low intensities. The response to climate change is much less studied. We can show that stratification changes have a strong influence on the changes in convection. Precipitation is strongly increasing for HW but decreasing for the VW simulations. For cloud type frequencies, virtually no changes are found for HW, but a substantial reduction in high clouds is found for VW. Further, we can show that the climate change signal strongly depends upon the horizontal resolution. In particular, significant differences between CPM and CRM are found in terms of the radiative feedbacks, with CRM exhibiting a stronger negative feedback in the top-of-the-atmosphere energy budget.

Importance of convective parameterization in ENSO predictions

NASA Astrophysics Data System (ADS)

Zhu, Jieshun; Kumar, Arun; Wang, Wanqiu; Hu, Zeng-Zhen; Huang, Bohua; Balmaseda, Magdalena A.

2017-06-01

This letter explored the influence of atmospheric convection scheme on El Niño-Southern Oscillation (ENSO) predictions using a set of hindcast experiments. Specifically, a low-resolution version of the Climate Forecast System version 2 is used for 12 month hindcasts starting from each April during 1982-2011. The hindcast experiments are repeated with three atmospheric convection schemes. All three hindcasts apply the identical initialization with ocean initial conditions taken from the European Centre for Medium-Range Weather Forecasts and atmosphere/land initial states from the National Centers for Environmental Prediction. Assessments indicate a substantial sensitivity of the sea surface temperature prediction skill to the different convection schemes, particularly over the eastern tropical Pacific. For the Niño 3.4 index, the anomaly correlation skill can differ by 0.1-0.2 at lead times longer than 2 months. Long-term simulations are further conducted with the three convection schemes to understand the differences in prediction skill. By conducting heat budget analyses for the mixed-layer temperature anomalies, it is suggested that the convection scheme having the highest skill simulates stronger and more realistic coupled feedbacks related to ENSO. Particularly, the strength of the Ekman pumping feedback is better represented, which is traced to more realistic simulation of surface wind stress. Our results imply that improving the mean state simulations in coupled (ocean-atmosphere) general circulation model (e.g., ameliorating the Intertropical Convergence Zone simulation) might further improve our ENSO prediction capability.

Demonstration of Effects on Tropical Cyclone Forecasts with a High Resolution Global Model from Variation in Cumulus Convection Parameterization

NASA Technical Reports Server (NTRS)

Miller, Timothy L.; Robertson, Franklin R.; Cohen, Charles; Mackaro, Jessica

2009-01-01

The Goddard Earth Observing System Model, Version 5 (GEOS-5) is a system of models that have been developed at Goddard Space Flight Center to support NASA's earth science research in data analysis, observing system modeling and design, climate and weather prediction, and basic research. The work presented used GEOS-5 with 0.25o horizontal resolution and 72 vertical levels (up to 0.01 hP) resolving both the troposphere and stratosphere, with closer packing of the levels close to the surface. The model includes explicit (grid-scale) moist physics, as well as convective parameterization schemes. Results will be presented that will demonstrate strong dependence in the results of modeling of a strong hurricane on the type of convective parameterization scheme used. The previous standard (default) option in the model was the Relaxed Arakawa-Schubert (RAS) scheme, which uses a quasi-equilibrium closure. In the cases shown, this scheme does not permit the efficient development of a strong storm in comparison with observations. When this scheme is replaced by a modified version of the Kain-Fritsch scheme, which was originally developed for use on grids with intervals of order 25 km such as the present one, the storm is able to develop to a much greater extent, closer to that of reality. Details of the two cases will be shown in order to elucidate the differences in the two modeled storms.

ARM - Midlatitude Continental Convective Clouds (comstock-hvps)

DOE Data Explorer

Jensen, Mike; Comstock, Jennifer; Genio, Anthony Del; Giangrande, Scott; Kollias, Pavlos

2012-01-06

Convective processes play a critical role in the Earth's energy balance through the redistribution of heat and moisture in the atmosphere and their link to the hydrological cycle. Accurate representation of convective processes in numerical models is vital towards improving current and future simulations of Earths climate system. Despite improvements in computing power, current operational weather and global climate models are unable to resolve the natural temporal and spatial scales important to convective processes and therefore must turn to parameterization schemes to represent these processes. In turn, parameterization schemes in cloud-resolving models need to be evaluated for their generality and application to a variety of atmospheric conditions. Data from field campaigns with appropriate forcing descriptors have been traditionally used by modelers for evaluating and improving parameterization schemes.

Rapid Parameterization Schemes for Aircraft Shape Optimization

NASA Technical Reports Server (NTRS)

Li, Wu

2012-01-01

A rapid shape parameterization tool called PROTEUS is developed for aircraft shape optimization. This tool can be applied directly to any aircraft geometry that has been defined in PLOT3D format, with the restriction that each aircraft component must be defined by only one data block. PROTEUS has eight types of parameterization schemes: planform, wing surface, twist, body surface, body scaling, body camber line, shifting/scaling, and linear morphing. These parametric schemes can be applied to two types of components: wing-type surfaces (e.g., wing, canard, horizontal tail, vertical tail, and pylon) and body-type surfaces (e.g., fuselage, pod, and nacelle). These schemes permit the easy setup of commonly used shape modification methods, and each customized parametric scheme can be applied to the same type of component for any configuration. This paper explains the mathematics for these parametric schemes and uses two supersonic configurations to demonstrate the application of these schemes.

Stochastic Convection Parameterizations: The Eddy-Diffusivity/Mass-Flux (EDMF) Approach (Invited)

NASA Astrophysics Data System (ADS)

Teixeira, J.

2013-12-01

In this presentation it is argued that moist convection parameterizations need to be stochastic in order to be realistic - even in deterministic atmospheric prediction systems. A new unified convection and boundary layer parameterization (EDMF) that optimally combines the Eddy-Diffusivity (ED) approach for smaller-scale boundary layer mixing with the Mass-Flux (MF) approach for larger-scale plumes is discussed. It is argued that for realistic simulations stochastic methods have to be employed in this new unified EDMF. Positive results from the implementation of the EDMF approach in atmospheric models are presented.

Evaluation of the Plant-Craig stochastic convection scheme (v2.0) in the ensemble forecasting system MOGREPS-R (24 km) based on the Unified Model (v7.3)

NASA Astrophysics Data System (ADS)

Keane, Richard J.; Plant, Robert S.; Tennant, Warren J.

2016-05-01

The Plant-Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic scheme only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant-Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant-Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.

Exploring the potential of machine learning to break deadlock in convection parameterization

NASA Astrophysics Data System (ADS)

Pritchard, M. S.; Gentine, P.

2017-12-01

We explore the potential of modern machine learning tools (via TensorFlow) to replace parameterization of deep convection in climate models. Our strategy begins by generating a large ( 1 Tb) training dataset from time-step level (30-min) output harvested from a one-year integration of a zonally symmetric, uniform-SST aquaplanet integration of the SuperParameterized Community Atmosphere Model (SPCAM). We harvest the inputs and outputs connecting each of SPCAM's 8,192 embedded cloud-resolving model (CRM) arrays to its host climate model's arterial thermodynamic state variables to afford 143M independent training instances. We demonstrate that this dataset is sufficiently large to induce preliminary convergence for neural network prediction of desired outputs of SP, i.e. CRM-mean convective heating and moistening profiles. Sensitivity of the machine learning convergence to the nuances of the TensorFlow implementation are discussed, as well as results from pilot tests from the neural network operating inline within the SPCAM as a replacement to the (super)parameterization of convection.

Investigating the scale-adaptivity of a shallow cumulus parameterization scheme with LES

NASA Astrophysics Data System (ADS)

Brast, Maren; Schemann, Vera; Neggers, Roel

2017-04-01

In this study we investigate the scale-adaptivity of a new parameterization scheme for shallow cumulus clouds in the gray zone. The Eddy-Diffusivity Multiple Mass-Flux (or ED(MF)n ) scheme is a bin-macrophysics scheme, in which subgrid transport is formulated in terms of discretized size densities. While scale-adaptivity in the ED-component is achieved using a pragmatic blending approach, the MF-component is filtered such that only the transport by plumes smaller than the grid size is maintained. For testing, ED(MF)n is implemented in a large-eddy simulation (LES) model, replacing the original subgrid-scheme for turbulent transport. LES thus plays the role of a non-hydrostatic testing ground, which can be run at different resolutions to study the behavior of the parameterization scheme in the boundary-layer gray zone. In this range convective cumulus clouds are partially resolved. We find that at high resolutions the clouds and the turbulent transport are predominantly resolved by the LES, and the transport represented by ED(MF)n is small. This partitioning changes towards coarser resolutions, with the representation of shallow cumulus clouds becoming exclusively carried by the ED(MF)n. The way the partitioning changes with grid-spacing matches the results of previous LES studies, suggesting some scale-adaptivity is captured. Sensitivity studies show that a scale-inadaptive ED component stays too active at high resolutions, and that the results are fairly insensitive to the number of transporting updrafts in the ED(MF)n scheme. Other assumptions in the scheme, such as the distribution of updrafts across sizes and the value of the area fraction covered by updrafts, are found to affect the location of the gray zone.

Improving and Understanding Climate Models: Scale-Aware Parameterization of Cloud Water Inhomogeneity and Sensitivity of MJO Simulation to Physical Parameters in a Convection Scheme

NASA Astrophysics Data System (ADS)

Xie, Xin

Microphysics and convection parameterizations are two key components in a climate model to simulate realistic climatology and variability of cloud distribution and the cycles of energy and water. When a model has varying grid size or simulations have to be run with different resolutions, scale-aware parameterization is desirable so that we do not have to tune model parameters tailored to a particular grid size. The subgrid variability of cloud hydrometers is known to impact microphysics processes in climate models and is found to highly depend on spatial scale. A scale- aware liquid cloud subgrid variability parameterization is derived and implemented in the Community Earth System Model (CESM) in this study using long-term radar-based ground measurements from the Atmospheric Radiation Measurement (ARM) program. When used in the default CESM1 with the finite-volume dynamic core where a constant liquid inhomogeneity parameter was assumed, the newly developed parameterization reduces the cloud inhomogeneity in high latitudes and increases it in low latitudes. This is due to both the smaller grid size in high latitudes, and larger grid size in low latitudes in the longitude-latitude grid setting of CESM as well as the variation of the stability of the atmosphere. The single column model and general circulation model (GCM) sensitivity experiments show that the new parameterization increases the cloud liquid water path in polar regions and decreases it in low latitudes. Current CESM1 simulation suffers from the bias of both the pacific double ITCZ precipitation and weak Madden-Julian oscillation (MJO). Previous studies show that convective parameterization with multiple plumes may have the capability to alleviate such biases in a more uniform and physical way. A multiple-plume mass flux convective parameterization is used in Community Atmospheric Model (CAM) to investigate the sensitivity of MJO simulations. We show that MJO simulation is sensitive to entrainment rate

Towards a parameterization of convective wind gusts in Sahel

NASA Astrophysics Data System (ADS)

Largeron, Yann; Guichard, Françoise; Bouniol, Dominique; Couvreux, Fleur; Birch, Cathryn; Beucher, Florent

2014-05-01

] who focused on the wet tropical Pacific region, and linked wind gusts to convective precipitation rates alone, here, we also analyse the subgrid wind distribution during convective events, and quantify the statistical moments (variance, skewness and kurtosis) in terms of mean wind speed and convective indexes such as DCAPE. Next step of the work will be to formulate a parameterization of the cold pool convective gust from those probability density functions and analytical formulaes obtained from basic energy budget models. References : [Carslaw et al., 2010] A review of natural aerosol interactions and feedbacks within the earth system. Atmospheric Chemistry and Physics, 10(4):1701{1737. [Engelstaedter et al., 2006] North african dust emissions and transport. Earth-Science Reviews, 79(1):73{100. [Knippertz and Todd, 2012] Mineral dust aerosols over the sahara: Meteorological controls on emission and transport and implications for modeling. Reviews of Geophysics, 50(1). [Marsham et al., 2011] The importance of the representation of deep convection for modeled dust-generating winds over west africa during summer.Geophysical Research Letters, 38(16). [Marticorena and Bergametti, 1995] Modeling the atmospheric dust cycle: 1. design of a soil-derived dust emission scheme. Journal of Geophysical Research, 100(D8):16415{16. [Menut, 2008] Sensitivity of hourly saharan dust emissions to ncep and ecmwf modeled wind speed. Journal of Geophysical Research: Atmospheres (1984{2012), 113(D16). [Pierre et al., 2012] Impact of vegetation and soil moisture seasonal dynamics on dust emissions over the sahel. Journal of Geophysical Research: Atmospheres (1984{2012), 117(D6). [Redelsperger et al., 2000] A parameterization of mesoscale enhancement of surface fluxes for large-scale models. Journal of climate, 13(2):402{421.

Resolution-dependent behavior of subgrid-scale vertical transport in the Zhang-McFarlane convection parameterization

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xiao, Heng; Gustafson, Jr., William I.; Hagos, Samson M.

2015-04-18

With this study, to better understand the behavior of quasi-equilibrium-based convection parameterizations at higher resolution, we use a diagnostic framework to examine the resolution-dependence of subgrid-scale vertical transport of moist static energy as parameterized by the Zhang-McFarlane convection parameterization (ZM). Grid-scale input to ZM is supplied by coarsening output from cloud-resolving model (CRM) simulations onto subdomains ranging in size from 8 × 8 to 256 × 256 km 2s.

Regional Climate Model sesitivity to different parameterizations schemes with WRF over Spain

NASA Astrophysics Data System (ADS)

García-Valdecasas Ojeda, Matilde; Raquel Gámiz-Fortis, Sonia; Hidalgo-Muñoz, Jose Manuel; Argüeso, Daniel; Castro-Díez, Yolanda; Jesús Esteban-Parra, María

2015-04-01

scheme is the most appropriate parameterization in term of temperatures because it better describes monthly minimum temperatures and seems to perform well for maximum temperatures. Regarding precipitation, ERA-Interim time series are slightly higher correlated with observations than WRF, but the bias and the RMSE are largely worse. These results also suggest that CAM V.5.1 2-moment 5-class microphysics schemes should not be used due to the computational cost with no apparent gain with respect to simpler schemes such as WRF single-moment 3-class. For the convection scheme, this study suggests that Betts-Miller-Janjic scheme is an appropriate choice due to its robustness and Kain-Fritsch cumulus scheme should not be used over this region. KEY WORDS: Regional climate modelling, physics schemes, parameterizations, WRF. ACKNOWLEDGEMENTS This work has been financed by the projects P11-RNM-7941 (Junta de Andalucía-Spain) and CGL2013-48539-R (MINECO-Spain, FEDER).

Redistribution of ice nuclei between cloud and rain droplets: Parameterization and application to deep convective clouds

NASA Astrophysics Data System (ADS)

Paukert, M.; Hoose, C.; Simmel, M.

2017-03-01

In model studies of aerosol-dependent immersion freezing in clouds, a common assumption is that each ice nucleating aerosol particle corresponds to exactly one cloud droplet. In contrast, the immersion freezing of larger drops—"rain"—is usually represented by a liquid volume-dependent approach, making the parameterizations of rain freezing independent of specific aerosol types and concentrations. This may lead to inconsistencies when aerosol effects on clouds and precipitation shall be investigated, since raindrops consist of the cloud droplets—and corresponding aerosol particles—that have been involved in drop-drop-collisions. Here we introduce an extension to a two-moment microphysical scheme in order to account explicitly for particle accumulation in raindrops by tracking the rates of selfcollection, autoconversion, and accretion. This provides a direct link between ice nuclei and the primary formation of large precipitating ice particles. A new parameterization scheme of drop freezing is presented to consider multiple ice nuclei within one drop and effective drop cooling rates. In our test cases of deep convective clouds, we find that at altitudes which are most relevant for immersion freezing, the majority of potential ice nuclei have been converted from cloud droplets into raindrops. Compared to the standard treatment of freezing in our model, the less efficient mineral dust-based freezing results in higher rainwater contents in the convective core, affecting both rain and hail precipitation. The aerosol-dependent treatment of rain freezing can reverse the signs of simulated precipitation sensitivities to ice nuclei perturbations.

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

NASA Astrophysics Data System (ADS)

Pergaud, Julien; Masson, Valéry; Malardel, Sylvie; Couvreux, Fleur

2009-07-01

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy DiffusivityMass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w *. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCSARM) and conserve a realistic evolution of stratocumulus (EUROCSFIRE).

Sensitivity of the summertime tropical Atlantic precipitation distribution to convective parameterization and model resolution in ECHAM6

DOE PAGES

Siongco, Angela Cheska; Hohenegger, Cathy; Stevens, Bjorn

2017-02-09

A realistic simulation of the tropical Atlantic precipitation distribution remains a challenge for atmospheric general circulation models, owing to their too coarse resolution that makes it necessary to parameterize convection. During boreal summer, models tend to underestimate the northward shift of the tropical Atlantic rain belt, leading to deficient precipitation over land and an anomalous precipitation maximum over the west Atlantic ocean. In this study, the model ECHAM6 is used to test the sensitivity of the precipitation biases to convective parameterization and horizontal resolution. Two sets of sensitivity experiments are performed. In the first set of experiments, modifications are appliedmore » to the convection scheme in order to investigate the relative roles of the trigger, entrainment, and closure formulations. In the second set, the model is run at high resolution with low-resolution boundary conditions in order to identify the relative contributions of a high-resolution atmosphere, orography, and surface. Results show that the dry bias over land in the model can be reduced by weakening the entrainment rate over land. Over ocean, it is found that the anomalous precipitation maximum occurs because of model choices that decrease the sensitivity of convection to the monsoon circulation in the east Atlantic. A reduction of the west Atlantic precipitation bias can be achieved by (i) using a moisture convergence closure, (ii) increasing the resolution of orography, or (iii) enhancing the production of deep convection in the east Atlantic. As a result, the biases over land and over ocean do not impact each other.« less

Uncertainty Quantification and Parameter Tuning: A Case Study of Convective Parameterization Scheme in the WRF Regional Climate Model

NASA Astrophysics Data System (ADS)

Qian, Y.; Yang, B.; Lin, G.; Leung, R.; Zhang, Y.

2012-04-01

The current tuning process of parameters in global climate models is often performed subjectively or treated as an optimization procedure to minimize model biases based on observations. The latter approach may provide more plausible values for a set of tunable parameters to approximate the observed climate, the system could be forced to an unrealistic physical state or improper balance of budgets through compensating errors over different regions of the globe. In this study, the Weather Research and Forecasting (WRF) model was used to provide a more flexible framework to investigate a number of issues related uncertainty quantification (UQ) and parameter tuning. The WRF model was constrained by reanalysis of data over the Southern Great Plains (SGP), where abundant observational data from various sources was available for calibration of the input parameters and validation of the model results. Focusing on five key input parameters in the new Kain-Fritsch (KF) convective parameterization scheme used in WRF as an example, the purpose of this study was to explore the utility of high-resolution observations for improving simulations of regional patterns and evaluate the transferability of UQ and parameter tuning across physical processes, spatial scales, and climatic regimes, which have important implications to UQ and parameter tuning in global and regional models. A stochastic important-sampling algorithm, Multiple Very Fast Simulated Annealing (MVFSA) was employed to efficiently sample the input parameters in the KF scheme based on a skill score so that the algorithm progressively moved toward regions of the parameter space that minimize model errors. The results based on the WRF simulations with 25-km grid spacing over the SGP showed that the precipitation bias in the model could be significantly reduced when five optimal parameters identified by the MVFSA algorithm were used. The model performance was found to be sensitive to downdraft- and entrainment

Prediction of heavy rainfall over Chennai Metropolitan City, Tamil Nadu, India: Impact of microphysical parameterization schemes

NASA Astrophysics Data System (ADS)

Singh, K. S.; Bonthu, Subbareddy; Purvaja, R.; Robin, R. S.; Kannan, B. A. M.; Ramesh, R.

2018-04-01

This study attempts to investigate the real-time prediction of a heavy rainfall event over the Chennai Metropolitan City, Tamil Nadu, India that occurred on 01 December 2015 using Advanced Research Weather Research and Forecasting (WRF-ARW) model. The study evaluates the impact of six microphysical (Lin, WSM6, Goddard, Thompson, Morrison and WDM6) parameterization schemes of the model on prediction of heavy rainfall event. In addition, model sensitivity has also been evaluated with six Planetary Boundary Layer (PBL) and two Land Surface Model (LSM) schemes. Model forecast was carried out using nested domain and the impact of model horizontal grid resolutions were assessed at 9 km, 6 km and 3 km. Analysis of the synoptic features using National Center for Environmental Prediction Global Forecast System (NCEP-GFS) analysis data revealed strong upper-level divergence and high moisture content at lower level were favorable for the occurrence of heavy rainfall event over the northeast coast of Tamil Nadu. The study signified that forecasted rainfall was more sensitive to the microphysics and PBL schemes compared to the LSM schemes. The model provided better forecast of the heavy rainfall event using the logical combination of Goddard microphysics, YSU PBL and Noah LSM schemes, and it was mostly attributed to timely initiation and development of the convective system. The forecast with different horizontal resolutions using cumulus parameterization indicated that the rainfall prediction was not well represented at 9 km and 6 km. The forecast with 3 km horizontal resolution provided better prediction in terms of timely initiation and development of the event. The study highlights that forecast of heavy rainfall events using a high-resolution mesoscale model with suitable representations of physical parameterization schemes are useful for disaster management and planning to minimize the potential loss of life and property.

Comparison of the Effects of RAS vs. Kain-Fritsch Convective Schemes on Katrina Forecasts with GEOS-5

NASA Technical Reports Server (NTRS)

Miller, Timothy L.; Cohen, Charles; Paxton, Jessica; Robertson, F. R. (Pete)

2009-01-01

Global forecasts were made with the 0.25-degree latitude version of GEOS-5, with the RAS scheme and with the Kain-Fritsch scheme. Examination was made of the Katrina (2005) hurricane simulation. Replacement of the RAS convective scheme with the K-F scheme results in a much more vigorous Katrina, closer to reality. Still, the result is not as vigorous as reality. In terms of wind maximum, the gap was closed by 50%. The result seems to be due to the RAS scheme drying out the boundary layer, thus hampering the grid-scale secondary circulation and attending cyclone development. The RAS case never developed a full warm core, whereas the K-F case did. Not shown here: The K-F scheme also resulted in a more vigorous storm than when GEOS-5 is run with no convective parameterization. Also not shown: An experiment in which the RAS firing level was moved up by 3 model levels resulted in a stronger, warm-core storm, though not as strong as the K-F case. Effects on storm track were noticed, but not studied.

Non-local Second Order Closure Scheme for Boundary Layer Turbulence and Convection

NASA Astrophysics Data System (ADS)

Meyer, Bettina; Schneider, Tapio

2017-04-01

There has been scientific consensus that the uncertainty in the cloud feedback remains the largest source of uncertainty in the prediction of climate parameters like climate sensitivity. To narrow down this uncertainty, not only a better physical understanding of cloud and boundary layer processes is required, but specifically the representation of boundary layer processes in models has to be improved. General climate models use separate parameterisation schemes to model the different boundary layer processes like small-scale turbulence, shallow and deep convection. Small scale turbulence is usually modelled by local diffusive parameterisation schemes, which truncate the hierarchy of moment equations at first order and use second-order equations only to estimate closure parameters. In contrast, the representation of convection requires higher order statistical moments to capture their more complex structure, such as narrow updrafts in a quasi-steady environment. Truncations of moment equations at second order may lead to more accurate parameterizations. At the same time, they offer an opportunity to take spatially correlated structures (e.g., plumes) into account, which are known to be important for convective dynamics. In this project, we study the potential and limits of local and non-local second order closure schemes. A truncation of the momentum equations at second order represents the same dynamics as a quasi-linear version of the equations of motion. We study the three-dimensional quasi-linear dynamics in dry and moist convection by implementing it in a LES model (PyCLES) and compare it to a fully non-linear LES. In the quasi-linear LES, interactions among turbulent eddies are suppressed but nonlinear eddy—mean flow interactions are retained, as they are in the second order closure. In physical terms, suppressing eddy—eddy interactions amounts to suppressing, e.g., interactions among convective plumes, while retaining interactions between plumes and the

Sensitivity of Coupled Tropical Pacific Model Biases to Convective Parameterization in CESM1

NASA Astrophysics Data System (ADS)

Woelfle, M. D.; Yu, S.; Bretherton, C. S.; Pritchard, M. S.

2018-01-01

Six month coupled hindcasts show the central equatorial Pacific cold tongue bias development in a GCM to be sensitive to the atmospheric convective parameterization employed. Simulations using the standard configuration of the Community Earth System Model version 1 (CESM1) develop a cold bias in equatorial Pacific sea surface temperatures (SSTs) within the first two months of integration due to anomalous ocean advection driven by overly strong easterly surface wind stress along the equator. Disabling the deep convection parameterization enhances the zonal pressure gradient leading to stronger zonal wind stress and a stronger equatorial SST bias, highlighting the role of pressure gradients in determining the strength of the cold bias. Superparameterized hindcasts show reduced SST bias in the cold tongue region due to a reduction in surface easterlies despite simulating an excessively strong low-level jet at 1-1.5 km elevation. This reflects inadequate vertical mixing of zonal momentum from the absence of convective momentum transport in the superparameterized model. Standard CESM1simulations modified to omit shallow convective momentum transport reproduce the superparameterized low-level wind bias and associated equatorial SST pattern. Further superparameterized simulations using a three-dimensional cloud resolving model capable of producing realistic momentum transport simulate a cold tongue similar to the default CESM1. These findings imply convective momentum fluxes may be an underappreciated mechanism for controlling the strength of the equatorial cold tongue. Despite the sensitivity of equatorial SST to these changes in convective parameterization, the east Pacific double-Intertropical Convergence Zone rainfall bias persists in all simulations presented in this study.

Coupling between lower‐tropospheric convective mixing and low‐level clouds: Physical mechanisms and dependence on convection scheme

PubMed Central

Bony, Sandrine; Dufresne, Jean‐Louis; Roehrig, Romain

2016-01-01

Abstract Several studies have pointed out the dependence of low‐cloud feedbacks on the strength of the lower‐tropospheric convective mixing. By analyzing a series of single‐column model experiments run by a climate model using two different convective parametrizations, this study elucidates the physical mechanisms through which marine boundary‐layer clouds depend on this mixing in the present‐day climate and under surface warming. An increased lower‐tropospheric convective mixing leads to a reduction of low‐cloud fraction. However, the rate of decrease strongly depends on how the surface latent heat flux couples to the convective mixing and to boundary‐layer cloud radiative effects: (i) on the one hand, the latent heat flux is enhanced by the lower‐tropospheric drying induced by the convective mixing, which damps the reduction of the low‐cloud fraction, (ii) on the other hand, the latent heat flux is reduced as the lower troposphere stabilizes under the effect of reduced low‐cloud radiative cooling, which enhances the reduction of the low‐cloud fraction. The relative importance of these two different processes depends on the closure of the convective parameterization. The convective scheme that favors the coupling between latent heat flux and low‐cloud radiative cooling exhibits a stronger sensitivity of low‐clouds to convective mixing in the present‐day climate, and a stronger low‐cloud feedback in response to surface warming. In this model, the low‐cloud feedback is stronger when the present‐day convective mixing is weaker and when present‐day clouds are shallower and more radiatively active. The implications of these insights for constraining the strength of low‐cloud feedbacks observationally is discussed. PMID:28239438

Coupling between lower-tropospheric convective mixing and low-level clouds: Physical mechanisms and dependence on convection scheme.

PubMed

Vial, Jessica; Bony, Sandrine; Dufresne, Jean-Louis; Roehrig, Romain

2016-12-01

Several studies have pointed out the dependence of low-cloud feedbacks on the strength of the lower-tropospheric convective mixing. By analyzing a series of single-column model experiments run by a climate model using two different convective parametrizations, this study elucidates the physical mechanisms through which marine boundary-layer clouds depend on this mixing in the present-day climate and under surface warming. An increased lower-tropospheric convective mixing leads to a reduction of low-cloud fraction. However, the rate of decrease strongly depends on how the surface latent heat flux couples to the convective mixing and to boundary-layer cloud radiative effects: (i) on the one hand, the latent heat flux is enhanced by the lower-tropospheric drying induced by the convective mixing, which damps the reduction of the low-cloud fraction, (ii) on the other hand, the latent heat flux is reduced as the lower troposphere stabilizes under the effect of reduced low-cloud radiative cooling, which enhances the reduction of the low-cloud fraction. The relative importance of these two different processes depends on the closure of the convective parameterization. The convective scheme that favors the coupling between latent heat flux and low-cloud radiative cooling exhibits a stronger sensitivity of low-clouds to convective mixing in the present-day climate, and a stronger low-cloud feedback in response to surface warming. In this model, the low-cloud feedback is stronger when the present-day convective mixing is weaker and when present-day clouds are shallower and more radiatively active. The implications of these insights for constraining the strength of low-cloud feedbacks observationally is discussed.

Assessment of the turbulence parameterization schemes for the Martian mesoscale simulations

NASA Astrophysics Data System (ADS)

Temel, Orkun; Karatekin, Ozgur; Van Beeck, Jeroen

2016-07-01

Turbulent transport within the Martian atmospheric boundary layer (ABL) is one of the most important physical processes in the Martian atmosphere due to the very thin structure of Martian atmosphere and super-adiabatic conditions during the diurnal cycle [1]. The realistic modeling of turbulent fluxes within the Martian ABL has a crucial effect on the many physical phenomena including dust devils [2], methane dispersion [3] and nocturnal jets [4]. Moreover, the surface heat and mass fluxes, which are related with the mass transport within the sub-surface of Mars, are being computed by the turbulence parameterization schemes. Therefore, in addition to the possible applications within the Martian boundary layer, parameterization of turbulence has an important effect on the biological research on Mars including the investigation of water cycle or sub-surface modeling. In terms of the turbulence modeling approaches being employed for the Martian ABL, the "planetary boundary layer (PBL) schemes" have been applied not only for the global circulation modeling but also for the mesoscale simulations [5]. The PBL schemes being used for Mars are the variants of the PBL schemes which had been developed for the Earth and these schemes are either based on the empirical determination of turbulent fluxes [6] or based on solving a one dimensional turbulent kinetic energy equation [7]. Even though, the Large Eddy Simulation techniques had also been applied with the regional models for Mars, it must be noted that these advanced models also use the features of these traditional PBL schemes for sub-grid modeling [8]. Therefore, assessment of these PBL schemes is vital for a better understanding the atmospheric processes of Mars. In this framework, this present study is devoted to the validation of different turbulence modeling approaches for the Martian ABL in comparison to Viking Lander [9] and MSL [10] datasets. The GCM/Mesoscale code being used is the PlanetWRF, the extended version

Role of a cumulus parameterization scheme in simulating atmospheric circulation and rainfall in the nine-layer Goddard Laboratory for Atmospheres General Circulation Model

NASA Technical Reports Server (NTRS)

Sud, Y. C.; Chao, Winston C.; Walker, G. K.

1992-01-01

The influence of a cumulus convection scheme on the simulated atmospheric circulation and hydrologic cycle is investigated by means of a coarse version of the GCM. Two sets of integrations, each containing an ensemble of three summer simulations, were produced. The ensemble sets of control and experiment simulations are compared and differentially analyzed to determine the influence of a cumulus convection scheme on the simulated circulation and hydrologic cycle. The results show that cumulus parameterization has a very significant influence on the simulation circulation and precipitation. The upper-level condensation heating over the ITCZ is much smaller for the experiment simulations as compared to the control simulations; correspondingly, the Hadley and Walker cells for the control simulations are also weaker and are accompanied by a weaker Ferrel cell in the Southern Hemisphere. Overall, the difference fields show that experiment simulations (without cumulus convection) produce a cooler and less energetic atmosphere.

Improving Convection and Cloud Parameterization Using ARM Observations and NCAR Community Atmosphere Model CAM5

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Guang J.

2016-11-07

The fundamental scientific objectives of our research are to use ARM observations and the NCAR CAM5 to understand the large-scale control on convection, and to develop improved convection and cloud parameterizations for use in GCMs.

An Eddy-Diffusivity Mass-flux (EDMF) closure for the unified representation of cloud and convective processes

NASA Astrophysics Data System (ADS)

Tan, Z.; Schneider, T.; Teixeira, J.; Lam, R.; Pressel, K. G.

2014-12-01

Sub-grid scale (SGS) closures in current climate models are usually decomposed into several largely independent parameterization schemes for different cloud and convective processes, such as boundary layer turbulence, shallow convection, and deep convection. These separate parameterizations usually do not converge as the resolution is increased or as physical limits are taken. This makes it difficult to represent the interactions and smooth transition among different cloud and convective regimes. Here we present an eddy-diffusivity mass-flux (EDMF) closure that represents all sub-grid scale turbulent, convective, and cloud processes in a unified parameterization scheme. The buoyant updrafts and precipitative downdrafts are parameterized with a prognostic multiple-plume mass-flux (MF) scheme. The prognostic term for the mass flux is kept so that the life cycles of convective plumes are better represented. The interaction between updrafts and downdrafts are parameterized with the buoyancy-sorting model. The turbulent mixing outside plumes is represented by eddy diffusion, in which eddy diffusivity (ED) is determined from a turbulent kinetic energy (TKE) calculated from a TKE balance that couples the environment with updrafts and downdrafts. Similarly, tracer variances are decomposed consistently between updrafts, downdrafts and the environment. The closure is internally coupled with a probabilistic cloud scheme and a simple precipitation scheme. We have also developed a relatively simple two-stream radiative scheme that includes the longwave (LW) and shortwave (SW) effects of clouds, and the LW effect of water vapor. We have tested this closure in a single-column model for various regimes spanning stratocumulus, shallow cumulus, and deep convection. The model is also run towards statistical equilibrium with climatologically relevant large-scale forcings. These model tests are validated against large-eddy simulation (LES) with the same forcings. The comparison of

Development of a nonlocal convective mixing scheme with varying upward mixing rates for use in air quality and chemical transport models.

PubMed

Mihailović, Dragutin T; Alapaty, Kiran; Sakradzija, Mirjana

2008-06-01

Asymmetrical convective non-local scheme (CON) with varying upward mixing rates is developed for simulation of vertical turbulent mixing in the convective boundary layer in air quality and chemical transport models. The upward mixing rate form the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. This scheme provides a less rapid mass transport out of surface layer into other layers than other asymmetrical convective mixing schemes. In this paper, we studied the performance of a nonlocal convective mixing scheme with varying upward mixing in the atmospheric boundary layer and its impact on the concentration of pollutants calculated with chemical and air-quality models. This scheme was additionally compared versus a local eddy-diffusivity scheme (KSC). Simulated concentrations of NO(2) and the nitrate wet deposition by the CON scheme are closer to the observations when compared to those obtained from using the KSC scheme. Concentrations calculated with the CON scheme are in general higher and closer to the observations than those obtained by the KSC scheme (of the order of 15-20%). Nitrate wet deposition calculated with the CON scheme are in general higher and closer to the observations than those obtained by the KSC scheme. To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO(2)) and nitrate wet deposition was compared for the year 2002. The comparison was made for the whole domain used in simulations performed by the chemical European Monitoring and Evaluation Programme Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.

Impact of Stochastic Parameterization Schemes on Coupled and Uncoupled Climate Simulations with the Community Earth System Model

NASA Astrophysics Data System (ADS)

Christensen, H. M.; Berner, J.; Coleman, D.; Palmer, T.

2015-12-01

Stochastic parameterizations have been used for more than a decade in atmospheric models to represent the variability of unresolved sub-grid processes. They have a beneficial effect on the spread and mean state of medium- and extended-range forecasts (Buizza et al. 1999, Palmer et al. 2009). There is also increasing evidence that stochastic parameterization of unresolved processes could be beneficial for the climate of an atmospheric model through noise enhanced variability, noise-induced drift (Berner et al. 2008), and by enabling the climate simulator to explore other flow regimes (Christensen et al. 2015; Dawson and Palmer 2015). We present results showing the impact of including the Stochastically Perturbed Parameterization Tendencies scheme (SPPT) in coupled runs of the National Center for Atmospheric Research (NCAR) Community Atmosphere Model, version 4 (CAM4) with historical forcing. The SPPT scheme accounts for uncertainty in the CAM physical parameterization schemes, including the convection scheme, by perturbing the parametrised temperature, moisture and wind tendencies with a multiplicative noise term. SPPT results in a large improvement in the variability of the CAM4 modeled climate. In particular, SPPT results in a significant improvement to the representation of the El Nino-Southern Oscillation in CAM4, improving the power spectrum, as well as both the inter- and intra-annual variability of tropical pacific sea surface temperatures. References: Berner, J., Doblas-Reyes, F. J., Palmer, T. N., Shutts, G. J., & Weisheimer, A., 2008. Phil. Trans. R. Soc A, 366, 2559-2577 Buizza, R., Miller, M. and Palmer, T. N., 1999. Q.J.R. Meteorol. Soc., 125, 2887-2908. Christensen, H. M., I. M. Moroz & T. N. Palmer, 2015. Clim. Dynam., doi: 10.1007/s00382-014-2239-9 Dawson, A. and T. N. Palmer, 2015. Clim. Dynam., doi: 10.1007/s00382-014-2238-x Palmer, T.N., R. Buizza, F. Doblas-Reyes, et al., 2009, ECMWF technical memorandum 598.

Best convective parameterization scheme within RegCM4 to downscale CMIP5 multi-model data for the CORDEX-MENA/Arab domain

NASA Astrophysics Data System (ADS)

Almazroui, Mansour; Islam, Md. Nazrul; Al-Khalaf, A. K.; Saeed, Fahad

2016-05-01

A suitable convective parameterization scheme within Regional Climate Model version 4.3.4 (RegCM4) developed by the Abdus Salam International Centre for Theoretical Physics, Trieste, Italy, is investigated through 12 sensitivity runs for the period 2000-2010. RegCM4 is driven with European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim 6-hourly boundary condition fields for the CORDEX-MENA/Arab domain. Besides ERA-Interim lateral boundary conditions data, the Climatic Research Unit (CRU) data is also used to assess the performance of RegCM4. Different statistical measures are taken into consideration in assessing model performance for 11 sub-domains throughout the analysis domain, out of which 7 (4) sub-domains give drier (wetter) conditions for the area of interest. There is no common best option for the simulation of both rainfall and temperature (with lowest bias); however, one option each for temperature and rainfall has been found to be superior among the 12 options investigated in this study. These best options for the two variables vary from region to region as well. Overall, RegCM4 simulates large pressure and water vapor values along with lower wind speeds compared to the driving fields, which are the key sources of bias in simulating rainfall and temperature. Based on the climatic characteristics of most of the Arab countries located within the study domain, the drier sub-domains are given priority in the selection of a suitable convective scheme, albeit with a compromise for both rainfall and temperature simulations. The most suitable option Grell over Land and Emanuel over Ocean in wet (GLEO wet) delivers a rainfall wet bias of 2.96 % and a temperature cold bias of 0.26 °C, compared to CRU data. An ensemble derived from all 12 runs provides unsatisfactory results for rainfall (28.92 %) and temperature (-0.54 °C) bias in the drier region because some options highly overestimate rainfall (reaching up to 200 %) and underestimate

Evaluation of predicted diurnal cycle of precipitation after tests with convection and microphysics schemes in the Eta Model

NASA Astrophysics Data System (ADS)

Gomes, J. L.; Chou, S. C.; Yaguchi, S. M.

2012-04-01

Physics parameterizations and the model vertical and horizontal resolutions, for example, can significantly contribute to the uncertainty in the numerical weather predictions, especially at regions with complex topography. The objective of this study is to assess the influences of model precipitation production schemes and horizontal resolution on the diurnal cycle of precipitation in the Eta Model . The model was run in hydrostatic mode at 3- and 5-km grid sizes, the vertical resolution was set to 50 layers, and the time steps to 6 and 10 s, respectively. The initial and boundary conditions were taken from ERA-Interim reanalysis. Over the sea the 0.25-deg sea surface temperature from NOAA was used. The model was setup to run for each resolution over Angra dos Reis, located in the Southeast region of Brazil, for the rainy period between 18 December 2009 and 01 de January 2010, the model simulation range was 48 hours. In one set of runs the cumulus parameterization was switched off, in this case the model precipitation was fully simulated by cloud microphysics scheme, and in the other set the model was run with weak cumulus convection. The results show that as the model horizontal resolution increases from 5 to 3 km, the spatial pattern of the precipitation hardly changed, although the maximum precipitation core increased in magnitude. Daily data from automatic station data was used to evaluate the runs and shows that the diurnal cycle of temperature and precipitation were better simulated for 3 km when compared against observations. The model configuration results without cumulus convection shows a small contraction in the precipitating area and an increase in the simulated maximum values. The diurnal cycle of precipitation was better simulated with some activity of the cumulus convection scheme. The skill scores for the period and for different forecast ranges are higher at weak and moderate precipitation rates.

Some issues in uncertainty quantification and parameter tuning: a case study of convective parameterization scheme in the WRF regional climate model

NASA Astrophysics Data System (ADS)

Yang, B.; Qian, Y.; Lin, G.; Leung, R.; Zhang, Y.

2011-12-01

The current tuning process of parameters in global climate models is often performed subjectively or treated as an optimization procedure to minimize model biases based on observations. While the latter approach may provide more plausible values for a set of tunable parameters to approximate the observed climate, the system could be forced to an unrealistic physical state or improper balance of budgets through compensating errors over different regions of the globe. In this study, the Weather Research and Forecasting (WRF) model was used to provide a more flexible framework to investigate a number of issues related uncertainty quantification (UQ) and parameter tuning. The WRF model was constrained by reanalysis of data over the Southern Great Plains (SGP), where abundant observational data from various sources was available for calibration of the input parameters and validation of the model results. Focusing on five key input parameters in the new Kain-Fritsch (KF) convective parameterization scheme used in WRF as an example, the purpose of this study was to explore the utility of high-resolution observations for improving simulations of regional patterns and evaluate the transferability of UQ and parameter tuning across physical processes, spatial scales, and climatic regimes, which have important implications to UQ and parameter tuning in global and regional models. A stochastic important-sampling algorithm, Multiple Very Fast Simulated Annealing (MVFSA) was employed to efficiently sample the input parameters in the KF scheme based on a skill score so that the algorithm progressively moved toward regions of the parameter space that minimize model errors. The results based on the WRF simulations with 25-km grid spacing over the SGP showed that the precipitation bias in the model could be significantly reduced when five optimal parameters identified by the MVFSA algorithm were used. The model performance was found to be sensitive to downdraft- and entrainment

Some issues in uncertainty quantification and parameter tuning: a case study of convective parameterization scheme in the WRF regional climate model

NASA Astrophysics Data System (ADS)

Yang, B.; Qian, Y.; Lin, G.; Leung, R.; Zhang, Y.

2012-03-01

The current tuning process of parameters in global climate models is often performed subjectively or treated as an optimization procedure to minimize model biases based on observations. While the latter approach may provide more plausible values for a set of tunable parameters to approximate the observed climate, the system could be forced to an unrealistic physical state or improper balance of budgets through compensating errors over different regions of the globe. In this study, the Weather Research and Forecasting (WRF) model was used to provide a more flexible framework to investigate a number of issues related uncertainty quantification (UQ) and parameter tuning. The WRF model was constrained by reanalysis of data over the Southern Great Plains (SGP), where abundant observational data from various sources was available for calibration of the input parameters and validation of the model results. Focusing on five key input parameters in the new Kain-Fritsch (KF) convective parameterization scheme used in WRF as an example, the purpose of this study was to explore the utility of high-resolution observations for improving simulations of regional patterns and evaluate the transferability of UQ and parameter tuning across physical processes, spatial scales, and climatic regimes, which have important implications to UQ and parameter tuning in global and regional models. A stochastic importance sampling algorithm, Multiple Very Fast Simulated Annealing (MVFSA) was employed to efficiently sample the input parameters in the KF scheme based on a skill score so that the algorithm progressively moved toward regions of the parameter space that minimize model errors. The results based on the WRF simulations with 25-km grid spacing over the SGP showed that the precipitation bias in the model could be significantly reduced when five optimal parameters identified by the MVFSA algorithm were used. The model performance was found to be sensitive to downdraft- and entrainment

Shortwave radiation parameterization scheme for subgrid topography

NASA Astrophysics Data System (ADS)

Helbig, N.; LöWe, H.

2012-02-01

Topography is well known to alter the shortwave radiation balance at the surface. A detailed radiation balance is therefore required in mountainous terrain. In order to maintain the computational performance of large-scale models while at the same time increasing grid resolutions, subgrid parameterizations are gaining more importance. A complete radiation parameterization scheme for subgrid topography accounting for shading, limited sky view, and terrain reflections is presented. Each radiative flux is parameterized individually as a function of sky view factor, slope and sun elevation angle, and albedo. We validated the parameterization with domain-averaged values computed from a distributed radiation model which includes a detailed shortwave radiation balance. Furthermore, we quantify the individual topographic impacts on the shortwave radiation balance. Rather than using a limited set of real topographies we used a large ensemble of simulated topographies with a wide range of typical terrain characteristics to study all topographic influences on the radiation balance. To this end slopes and partial derivatives of seven real topographies from Switzerland and the United States were analyzed and Gaussian statistics were found to best approximate real topographies. Parameterized direct beam radiation presented previously compared well with modeled values over the entire range of slope angles. The approximation of multiple, anisotropic terrain reflections with single, isotropic terrain reflections was confirmed as long as domain-averaged values are considered. The validation of all parameterized radiative fluxes showed that it is indeed not necessary to compute subgrid fluxes in order to account for all topographic influences in large grid sizes.

Intercomparison Project on Parameterizations of Large-Scale Dynamics for Simulations of Tropical Convection

NASA Astrophysics Data System (ADS)

Sobel, A. H.; Wang, S.; Bellon, G.; Sessions, S. L.; Woolnough, S.

2013-12-01

Parameterizations of large-scale dynamics have been developed in the past decade for studying the interaction between tropical convection and large-scale dynamics, based on our physical understanding of the tropical atmosphere. A principal advantage of these methods is that they offer a pathway to attack the key question of what controls large-scale variations of tropical deep convection. These methods have been used with both single column models (SCMs) and cloud-resolving models (CRMs) to study the interaction of deep convection with several kinds of environmental forcings. While much has been learned from these efforts, different groups' efforts are somewhat hard to compare. Different models, different versions of the large-scale parameterization methods, and experimental designs that differ in other ways are used. It is not obvious which choices are consequential to the scientific conclusions drawn and which are not. The methods have matured to the point that there is value in an intercomparison project. In this context, the Global Atmospheric Systems Study - Weak Temperature Gradient (GASS-WTG) project was proposed at the Pan-GASS meeting in September 2012. The weak temperature gradient approximation is one method to parameterize large-scale dynamics, and is used in the project name for historical reasons and simplicity, but another method, the damped gravity wave (DGW) method, will also be used in the project. The goal of the GASS-WTG project is to develop community understanding of the parameterization methods currently in use. Their strengths, weaknesses, and functionality in models with different physics and numerics will be explored in detail, and their utility to improve our understanding of tropical weather and climate phenomena will be further evaluated. This presentation will introduce the intercomparison project, including background, goals, and overview of the proposed experimental design. Interested groups will be invited to join (it will not be

Order of accuracy of QUICK and related convection-diffusion schemes

NASA Technical Reports Server (NTRS)

Leonard, B. P.

1993-01-01

This report attempts to correct some misunderstandings that have appeared in the literature concerning the order of accuracy of the QUICK scheme for steady-state convective modeling. Other related convection-diffusion schemes are also considered. The original one-dimensional QUICK scheme written in terms of nodal-point values of the convected variable (with a 1/8-factor multiplying the 'curvature' term) is indeed a third-order representation of the finite volume formulation of the convection operator average across the control volume, written naturally in flux-difference form. An alternative single-point upwind difference scheme (SPUDS) using node values (with a 1/6-factor) is a third-order representation of the finite difference single-point formulation; this can be written in a pseudo-flux difference form. These are both third-order convection schemes; however, the QUICK finite volume convection operator is 33 percent more accurate than the single-point implementation of SPUDS. Another finite volume scheme, writing convective fluxes in terms of cell-average values, requires a 1/6-factor for third-order accuracy. For completeness, one can also write a single-point formulation of the convective derivative in terms of cell averages, and then express this in pseudo-flux difference form; for third-order accuracy, this requires a curvature factor of 5/24. Diffusion operators are also considered in both single-point and finite volume formulations. Finite volume formulations are found to be significantly more accurate. For example, classical second-order central differencing for the second derivative is exactly twice as accurate in a finite volume formulation as it is in single-point.

Sensitivities of Summertime Mesoscale Circulations in the Coastal Carolinas to Modifications of the Kain-Fritsch Cumulus Parameterization.

PubMed

Sims, Aaron P; Alapaty, Kiran; Raman, Sethu

2017-01-01

Two mesoscale circulations, the Sandhills circulation and the sea breeze, influence the initiation of deep convection over the Sandhills and the coast in the Carolinas during the summer months. The interaction of these two circulations causes additional convection in this coastal region. Accurate representation of mesoscale convection is difficult as numerical models have problems with the prediction of the timing, amount, and location of precipitation. To address this issue, the authors have incorporated modifications to the Kain-Fritsch (KF) convective parameterization scheme and evaluated these mesoscale interactions using a high-resolution numerical model. The modifications include changes to the subgrid-scale cloud formulation, the convective turnover time scale, and the formulation of the updraft entrainment rates. The use of a grid-scaling adjustment parameter modulates the impact of the KF scheme as a function of the horizontal grid spacing used in a simulation. Results indicate that the impact of this modified cumulus parameterization scheme is more effective on domains with coarser grid sizes. Other results include a decrease in surface and near-surface temperatures in areas of deep convection (due to the inclusion of the effects of subgrid-scale clouds on the radiation), improvement in the timing of convection, and an increase in the strength of deep convection.

A ubiquitous ice size bias in simulations of tropical deep convection

NASA Astrophysics Data System (ADS)

Stanford, McKenna W.; Varble, Adam; Zipser, Ed; Strapp, J. Walter; Leroy, Delphine; Schwarzenboeck, Alfons; Potts, Rodney; Protat, Alain

2017-08-01

The High Altitude Ice Crystals - High Ice Water Content (HAIC-HIWC) joint field campaign produced aircraft retrievals of total condensed water content (TWC), hydrometeor particle size distributions (PSDs), and vertical velocity (w) in high ice water content regions of mature and decaying tropical mesoscale convective systems (MCSs). The resulting dataset is used here to explore causes of the commonly documented high bias in radar reflectivity within cloud-resolving simulations of deep convection. This bias has been linked to overly strong simulated convective updrafts lofting excessive condensate mass but is also modulated by parameterizations of hydrometeor size distributions, single particle properties, species separation, and microphysical processes. Observations are compared with three Weather Research and Forecasting model simulations of an observed MCS using different microphysics parameterizations while controlling for w, TWC, and temperature. Two popular bulk microphysics schemes (Thompson and Morrison) and one bin microphysics scheme (fast spectral bin microphysics) are compared. For temperatures between -10 and -40 °C and TWC > 1 g m-3, all microphysics schemes produce median mass diameters (MMDs) that are generally larger than observed, and the precipitating ice species that controls this size bias varies by scheme, temperature, and w. Despite a much greater number of samples, all simulations fail to reproduce observed high-TWC conditions ( > 2 g m-3) between -20 and -40 °C in which only a small fraction of condensate mass is found in relatively large particle sizes greater than 1 mm in diameter. Although more mass is distributed to large particle sizes relative to those observed across all schemes when controlling for temperature, w, and TWC, differences with observations are significantly variable between the schemes tested. As a result, this bias is hypothesized to partly result from errors in parameterized hydrometeor PSD and single particle

Numerical Study of the Role of Shallow Convection in Moisture Transport and Climate

NASA Technical Reports Server (NTRS)

Seaman, Nelson L.; Stauffer, David R.; Munoz, Ricardo C.

2001-01-01

The objective of this investigation was to study the role of shallow convection on the regional water cycle of the Mississippi and Little Washita Basins of the Southern Great Plains (SGP) using a 3-D mesoscale model, the PSU/NCAR MM5. The underlying premise of the project was that current modeling of regional-scale climate and moisture cycles over the continents is deficient without adequate treatment of shallow convection. At the beginning of the study, it was hypothesized that an improved treatment of the regional water cycle can be achieved by using a 3-D mesoscale numerical model having high-quality parameterizations for the key physical processes controlling the water cycle. These included a detailed land-surface parameterization (the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) sub-model of Wetzel and Boone), an advanced boundary-layer parameterization (the 1.5-order turbulent kinetic energy (TKE) predictive scheme of Shafran et al.), and a more complete shallow convection parameterization (the hybrid-closure scheme of Deng et al.) than are available in most current models. PLACE is a product of researchers working at NASA's Goddard Space Flight Center in Greenbelt, MD. The TKE and shallow-convection schemes are the result of model development at Penn State. The long-range goal is to develop an integrated suite of physical sub-models that can be used for regional and perhaps global climate studies of the water budget. Therefore, the work plan focused on integrating, improving, and testing these parameterizations in the MM5 and applying them to study water-cycle processes over the SGP. These schemes have been tested extensively through the course of this study and the latter two have been improved significantly as a consequence.

New, Improved Bulk-microphysical Schemes for Studying Precipitation Processes in WRF. Part 1; Comparisons with Other Schemes

NASA Technical Reports Server (NTRS)

Tao, W.-K.; Shi, J.; Chen, S. S> ; Lang, S.; Hong, S.-Y.; Thompson, G.; Peters-Lidard, C.; Hou, A.; Braun, S.;

Redistribution of ice nuclei between cloud and rain droplets: Parameterization and application to deep convective clouds: ICE NUCLEI IN RAIN DROPLETS

DOE PAGES

Paukert, M.; Hoose, C.; Simmel, M.

2017-02-21

In model studies of aerosol-dependent immersion freezing in clouds, a common assumption is that each ice nucleating aerosol particle corresponds to exactly one cloud droplet. Conversely, the immersion freezing of larger drops—“rain”—is usually represented by a liquid volume-dependent approach, making the parameterizations of rain freezing independent of specific aerosol types and concentrations. This may lead to inconsistencies when aerosol effects on clouds and precipitation shall be investigated, since raindrops consist of the cloud droplets—and corresponding aerosol particles—that have been involved in drop-drop-collisions. We introduce an extension to a two-moment microphysical scheme in order to account explicitly for particle accumulation inmore » raindrops by tracking the rates of selfcollection, autoconversion, and accretion. This also provides a direct link between ice nuclei and the primary formation of large precipitating ice particles. A new parameterization scheme of drop freezing is presented to consider multiple ice nuclei within one drop and effective drop cooling rates. In our test cases of deep convective clouds, we find that at altitudes which are most relevant for immersion freezing, the majority of potential ice nuclei have been converted from cloud droplets into raindrops. Compared to the standard treatment of freezing in our model, the less efficient mineral dust-based freezing results in higher rainwater contents in the convective core, affecting both rain and hail precipitation. The aerosol-dependent treatment of rain freezing can reverse the signs of simulated precipitation sensitivities to ice nuclei perturbations.« less

Redistribution of ice nuclei between cloud and rain droplets: Parameterization and application to deep convective clouds: ICE NUCLEI IN RAIN DROPLETS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Paukert, M.; Hoose, C.; Simmel, M.

In model studies of aerosol-dependent immersion freezing in clouds, a common assumption is that each ice nucleating aerosol particle corresponds to exactly one cloud droplet. Conversely, the immersion freezing of larger drops—“rain”—is usually represented by a liquid volume-dependent approach, making the parameterizations of rain freezing independent of specific aerosol types and concentrations. This may lead to inconsistencies when aerosol effects on clouds and precipitation shall be investigated, since raindrops consist of the cloud droplets—and corresponding aerosol particles—that have been involved in drop-drop-collisions. We introduce an extension to a two-moment microphysical scheme in order to account explicitly for particle accumulation inmore » raindrops by tracking the rates of selfcollection, autoconversion, and accretion. This also provides a direct link between ice nuclei and the primary formation of large precipitating ice particles. A new parameterization scheme of drop freezing is presented to consider multiple ice nuclei within one drop and effective drop cooling rates. In our test cases of deep convective clouds, we find that at altitudes which are most relevant for immersion freezing, the majority of potential ice nuclei have been converted from cloud droplets into raindrops. Compared to the standard treatment of freezing in our model, the less efficient mineral dust-based freezing results in higher rainwater contents in the convective core, affecting both rain and hail precipitation. The aerosol-dependent treatment of rain freezing can reverse the signs of simulated precipitation sensitivities to ice nuclei perturbations.« less

Why is the simulated climatology of tropical cyclones so sensitive to the choice of cumulus parameterization scheme in the WRF model?

NASA Astrophysics Data System (ADS)

Zhang, Chunxi; Wang, Yuqing

2018-01-01

The sensitivity of simulated tropical cyclones (TCs) to the choice of cumulus parameterization (CP) scheme in the advanced Weather Research and Forecasting Model (WRF-ARW) version 3.5 is analyzed based on ten seasonal simulations with 20-km horizontal grid spacing over the western North Pacific. Results show that the simulated frequency and intensity of TCs are very sensitive to the choice of the CP scheme. The sensitivity can be explained well by the difference in the low-level circulation in a height and sorted moisture space. By transporting moist static energy from dry to moist region, the low-level circulation is important to convective self-aggregation which is believed to be related to genesis of TC-like vortices (TCLVs) and TCs in idealized settings. The radiative and evaporative cooling associated with low-level clouds and shallow convection in dry regions is found to play a crucial role in driving the moisture-sorted low-level circulation. With shallow convection turned off in a CP scheme, relatively strong precipitation occurs frequently in dry regions. In this case, the diabatic cooling can still drive the low-level circulation but its strength is reduced and thus TCLV/TC genesis is suppressed. The inclusion of the cumulus momentum transport (CMT) in a CP scheme can considerably suppress genesis of TCLVs/TCs, while changes in the moisture-sorted low-level circulation and horizontal distribution of precipitation are trivial, indicating that the CMT modulates the TCLVs/TCs activities in the model by mechanisms other than the horizontal transport of moist static energy.

A Parameterization for the Triggering of Landscape Generated Moist Convection

NASA Technical Reports Server (NTRS)

Lynn, Barry H.; Tao, Wei-Kuo; Abramopoulos, Frank

1998-01-01

A set of relatively high resolution three-dimensional (3D) simulations were produced to investigate the triggering of moist convection by landscape generated mesoscale circulations. The local accumulated rainfall varied monotonically (linearly) with the size of individual landscape patches, demonstrating the need to develop a trigger function that is sensitive to the size of individual patches. A new triggering function that includes the effect of landscapes generated mesoscale circulations over patches of different sizes consists of a parcel's perturbation in vertical velocity (nu(sub 0)), temperature (theta(sub 0)), and moisture (q(sub 0)). Each variable in the triggering function was also sensitive to soil moisture gradients, atmospheric initial conditions, and moist processes. The parcel's vertical velocity, temperature, and moisture perturbation were partitioned into mesoscale and turbulent components. Budget equations were derived for theta(sub 0) and q(sub 0). Of the many terms in this set of budget equations, the turbulent, vertical flux of the mesoscale temperature and moisture contributed most to the triggering of moist convection through the impact of these fluxes on the parcel's temperature and moisture profile. These fluxes needed to be parameterized to obtain theta(sub 0) and q(sub 0). The mesoscale vertical velocity also affected the profile of nu(sub 0). We used similarity theory to parameterize these fluxes as well as the parcel's mesoscale vertical velocity.

Prediction of convective activity using a system of parasitic-nested numerical models

NASA Technical Reports Server (NTRS)

Perkey, D. J.

1976-01-01

A limited area, three dimensional, moist, primitive equation (PE) model is developed to test the sensitivity of quantitative precipitation forecasts to the initial relative humidity distribution. Special emphasis is placed on the squall-line region. To accomplish the desired goal, time dependent lateral boundaries and a general convective parameterization scheme suitable for mid-latitude systems were developed. The sequential plume convective parameterization scheme presented is designed to have the versatility necessary in mid-latitudes and to be applicable for short-range forecasts. The results indicate that the scheme is able to function in the frontally forced squallline region, in the gently rising altostratus region ahead of the approaching low center, and in the over-riding region ahead of the warm front. Three experiments are discussed.

Parameterized Radiative Convective Equilibrium Across a Range of Domains: A Unifying Tool for General Circulation Models and High Resolution Models

NASA Astrophysics Data System (ADS)

Silvers, L. G.; Stevens, B. B.; Mauritsen, T.; Marco, G. A.

2015-12-01

The characteristics of clouds in General Circulation Models (GCMs) need to be constrained in a consistent manner with theory, observations, and high resolution models (HRMs). One way forward is to base improvements of parameterizations on high resolution studies which resolve more of the important dynamical motions and allow for less parameterizations. This is difficult because of the numerous differences between GCMs and HRMs, both technical and theoretical. Century long simulations at resolutions of 20-250 km on a global domain are typical of GCMs while HRMs often simulate hours at resolutions of 0.1km-5km on domains the size of a single GCM grid cell. The recently developed mode ICON provides a flexible framework which allows many of these difficulties to be overcome. This study uses the ICON model to compute SST perturbation simulations on multiple domains in a state of Radiative Convective Equilibrium (RCE) with parameterized convection. The domains used range from roughly the size of Texas to nearly half of Earth's surface area. All simulations use a doubly periodic domain with an effective distance between cell centers of 13 km and are integrated to a state of statistical stationarity. The primary analysis examines the mean characteristics of the cloud related fields and the feedback parameter of the simulations. It is shown that the simulated atmosphere of a GCM in RCE is sufficiently similar across a range of domain sizes to justify the use of RCE to study both a GCM and a HRM on the same domain with the goal of improved constraints on the parameterized clouds. The simulated atmospheres are comparable to what could be expected at midday in a typical region of Earth's tropics under calm conditions. In particular, the differences between the domains are smaller than differences which result from choosing different physics schemes. Significant convective organization is present on all domain sizes with a relatively high subsidence fraction. Notwithstanding

Intercomparison of Martian Lower Atmosphere Simulated Using Different Planetary Boundary Layer Parameterization Schemes

NASA Technical Reports Server (NTRS)

Natarajan, Murali; Fairlie, T. Duncan; Dwyer Cianciolo, Alicia; Smith, Michael D.

2015-01-01

We use the mesoscale modeling capability of Mars Weather Research and Forecasting (MarsWRF) model to study the sensitivity of the simulated Martian lower atmosphere to differences in the parameterization of the planetary boundary layer (PBL). Characterization of the Martian atmosphere and realistic representation of processes such as mixing of tracers like dust depend on how well the model reproduces the evolution of the PBL structure. MarsWRF is based on the NCAR WRF model and it retains some of the PBL schemes available in the earth version. Published studies have examined the performance of different PBL schemes in NCAR WRF with the help of observations. Currently such assessments are not feasible for Martian atmospheric models due to lack of observations. It is of interest though to study the sensitivity of the model to PBL parameterization. Typically, for standard Martian atmospheric simulations, we have used the Medium Range Forecast (MRF) PBL scheme, which considers a correction term to the vertical gradients to incorporate nonlocal effects. For this study, we have also used two other parameterizations, a non-local closure scheme called Yonsei University (YSU) PBL scheme and a turbulent kinetic energy closure scheme called Mellor- Yamada-Janjic (MYJ) PBL scheme. We will present intercomparisons of the near surface temperature profiles, boundary layer heights, and wind obtained from the different simulations. We plan to use available temperature observations from Mini TES instrument onboard the rovers Spirit and Opportunity in evaluating the model results.

Sensitivity of CONUS Summer Rainfall to the Selection of Cumulus Parameterization Schemes in NU-WRF Seasonal Simulations

NASA Technical Reports Server (NTRS)

Iguchi, Takamichi; Tao, Wei-Kuo; Wu, Di; Peters-Lidard, Christa; Santanello, Joseph A.; Kemp, Eric; Tian, Yudong; Case, Jonathan; Wang, Weile; Ferraro, Robert;

Processing Doppler Lidar and Cloud Radar Observations for Analysis of Convective Mass Flux Parameterizations Using DYNAMO Direct Observations

DTIC Science & Technology

2014-09-30

for Analysis of Convective Mass Flux Parameterizations Using DYNAMO Direct Observations R. Michael Hardesty CIRES/University of Colorado/NOAA 325...the RV-Revell during legs 2 & 3 of the DYNAMO experiement to help characterize vertical transport through the boundary layer and to build statistics...obtained during DYNAMO , and to investigate whether cold pools that emanate from convection organize the interplay between humidity and convection and

Sensitivity of a Cumulus Parameterization Scheme to Precipitation Production Representation and Its Impact on a Heavy Rain Event over Korea

DOE Office of Scientific and Technical Information (OSTI.GOV)

Han, Ji-Young; Hong, Song-You; Sunny Lim, Kyo-Sun

The sensitivity of a cumulus parameterization scheme (CPS) to a representation of precipitation production is examined. To do this, the parameter that determines the fraction of cloud condensate converted to precipitation in the simplified Arakawa–Schubert (SAS) convection scheme is modified following the results from a cloud-resolving simulation. While the original conversion parameter is assumed to be constant, the revised parameter includes a temperature dependency above the freezing level, whichleadstolessproductionoffrozenprecipitating condensate with height. The revised CPS has been evaluated for a heavy rainfall event over Korea as well as medium-range forecasts using the Global/Regional Integrated Model system (GRIMs). The inefficient conversionmore » of cloud condensate to convective precipitation at colder temperatures generally leads to a decrease in pre-cipitation, especially in the category of heavy rainfall. The resultant increase of detrained moisture induces moistening and cooling at the top of clouds. A statistical evaluation of the medium-range forecasts with the revised precipitation conversion parameter shows an overall improvement of the forecast skill in precipitation and large-scale fields, indicating importance of more realistic representation of microphysical processes in CPSs.« less

Parameterization Interactions in Global Aquaplanet Simulations

NASA Astrophysics Data System (ADS)

Bhattacharya, Ritthik; Bordoni, Simona; Suselj, Kay; Teixeira, João.

2018-02-01

Global climate simulations rely on parameterizations of physical processes that have scales smaller than the resolved ones. In the atmosphere, these parameterizations represent moist convection, boundary layer turbulence and convection, cloud microphysics, longwave and shortwave radiation, and the interaction with the land and ocean surface. These parameterizations can generate different climates involving a wide range of interactions among parameterizations and between the parameterizations and the resolved dynamics. To gain a simplified understanding of a subset of these interactions, we perform aquaplanet simulations with the global version of the Weather Research and Forecasting (WRF) model employing a range (in terms of properties) of moist convection and boundary layer (BL) parameterizations. Significant differences are noted in the simulated precipitation amounts, its partitioning between convective and large-scale precipitation, as well as in the radiative impacts. These differences arise from the way the subcloud physics interacts with convection, both directly and through various pathways involving the large-scale dynamics and the boundary layer, convection, and clouds. A detailed analysis of the profiles of the different tendencies (from the different physical processes) for both potential temperature and water vapor is performed. While different combinations of convection and boundary layer parameterizations can lead to different climates, a key conclusion of this study is that similar climates can be simulated with model versions that are different in terms of the partitioning of the tendencies: the vertically distributed energy and water balances in the tropics can be obtained with significantly different profiles of large-scale, convection, and cloud microphysics tendencies.

Optimizing the Betts-Miller-Janjic cumulus parameterization with Intel Many Integrated Core (MIC) architecture

NASA Astrophysics Data System (ADS)

Huang, Melin; Huang, Bormin; Huang, Allen H.-L.

2015-10-01

The schemes of cumulus parameterization are responsible for the sub-grid-scale effects of convective and/or shallow clouds, and intended to represent vertical fluxes due to unresolved updrafts and downdrafts and compensating motion outside the clouds. Some schemes additionally provide cloud and precipitation field tendencies in the convective column, and momentum tendencies due to convective transport of momentum. The schemes all provide the convective component of surface rainfall. Betts-Miller-Janjic (BMJ) is one scheme to fulfill such purposes in the weather research and forecast (WRF) model. National Centers for Environmental Prediction (NCEP) has tried to optimize the BMJ scheme for operational application. As there are no interactions among horizontal grid points, this scheme is very suitable for parallel computation. With the advantage of Intel Xeon Phi Many Integrated Core (MIC) architecture, efficient parallelization and vectorization essentials, it allows us to optimize the BMJ scheme. If compared to the original code respectively running on one CPU socket (eight cores) and on one CPU core with Intel Xeon E5-2670, the MIC-based optimization of this scheme running on Xeon Phi coprocessor 7120P improves the performance by 2.4x and 17.0x, respectively.

ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow

NASA Technical Reports Server (NTRS)

Leonard, B. P.; Mokhtari, Simin

1990-01-01

For convection-dominated flows, classical second-order methods are notoriously oscillatory and often unstable. For this reason, many computational fluid dynamicists have adopted various forms of (inherently stable) first-order upwinding over the past few decades. Although it is now well known that first-order convection schemes suffer from serious inaccuracies attributable to artificial viscosity or numerical diffusion under high convection conditions, these methods continue to enjoy widespread popularity for numerical heat transfer calculations, apparently due to a perceived lack of viable high accuracy alternatives. But alternatives are available. For example, nonoscillatory methods used in gasdynamics, including currently popular TVD schemes, can be easily adapted to multidimensional incompressible flow and convective transport. This, in itself, would be a major advance for numerical convective heat transfer, for example. But, as is shown, second-order TVD schemes form only a small, overly restrictive, subclass of a much more universal, and extremely simple, nonoscillatory flux-limiting strategy which can be applied to convection schemes of arbitrarily high order accuracy, while requiring only a simple tridiagonal ADI line-solver, as used in the majority of general purpose iterative codes for incompressible flow and numerical heat transfer. The new universal limiter and associated solution procedures form the so-called ULTRA-SHARP alternative for high resolution nonoscillatory multidimensional steady state high speed convective modelling.

A Mass-Flux Scheme View of a High-Resolution Simulation of a Transition from Shallow to Deep Cumulus Convection.

NASA Astrophysics Data System (ADS)

Kuang, Zhiming; Bretherton, Christopher S.

2006-07-01

In this paper, an idealized, high-resolution simulation of a gradually forced transition from shallow, nonprecipitating to deep, precipitating cumulus convection is described; how the cloud and transport statistics evolve as the convection deepens is explored; and the collected statistics are used to evaluate assumptions in current cumulus schemes. The statistical analysis methodologies that are used do not require tracing the history of individual clouds or air parcels; instead they rely on probing the ensemble characteristics of cumulus convection in the large model dataset. They appear to be an attractive way for analyzing outputs from cloud-resolving numerical experiments. Throughout the simulation, it is found that 1) the initial thermodynamic properties of the updrafts at the cloud base have rather tight distributions; 2) contrary to the assumption made in many cumulus schemes, nearly undiluted air parcels are too infrequent to be relevant to any stage of the simulated convection; and 3) a simple model with a spectrum of entraining plumes appears to reproduce most features of the cloudy updrafts, but significantly overpredicts the mass flux as the updrafts approach their levels of zero buoyancy. A buoyancy-sorting model was suggested as a potential remedy. The organized circulations of cold pools seem to create clouds with larger-sized bases and may correspondingly contribute to their smaller lateral entrainment rates. Our results do not support a mass-flux closure based solely on convective available potential energy (CAPE), and are in general agreement with a convective inhibition (CIN)-based closure. The general similarity in the ensemble characteristics of shallow and deep convection and the continuous evolution of the thermodynamic structure during the transition provide justification for developing a single unified cumulus parameterization that encompasses both shallow and deep convection.

Challenges in Parameterizing the Lifecycle of Cumulus Convection in Global Climate Models

NASA Astrophysics Data System (ADS)

Del Genio, A. D.

2012-12-01

Moist convection exerts a strong influence on Earth's general circulation, energy cycle, and water cycle and has long been considered among the most difficult processes to represent in global climate models. Historically, convection has been portrayed in models as a collection of individual cells, and most of the attention has focused on deep precipitating convection that adjusts quickly to large-scale processes that destabilize the atmosphere. Only in the past decade has the need to represent the full convective lifecycle been recognized by the global climate modeling community, although many of the relevant features have been observed in field experiments for decades. Progress has accelerated in recent years with the aid of insights gained from cloud-resolving models and new satellite and surface remote sensing datasets. There has also been a welcome trend away from emphasis on the mean state and toward understanding of major modes of convective variability such as the Madden-Julian Oscillation and the continental diurnal cycle. On one end of the lifecycle, the need to capture the gradual transition from shallow to congestus to deep convection has renewed interest in understanding the process of entrainment and the previously underappreciated sensitivity of convection to the humidity of the free troposphere. On the other end, the tendency for convection to organize on the mesoscale in favorable humidity and shear conditions is only now beginning to receive attention in the parameterization community. Approaches to representing downdraft cold pools, which stimulate further convection and trigger organization, are now being implemented in GCMs. The subsequent evolution from convective cells to organized clusters with stratiform precipitation, which shifts the heating profile upward, extends the lifetime of convective systems, and can change the sign of convective momentum transport, remains a challenge, especially as model resolution increases.

Evaluating Cloud Initialization in a Convection-permit NWP Model

NASA Astrophysics Data System (ADS)

Li, Jia; Chen, Baode

2015-04-01

In general, to avoid "double counting precipitation" problem, in convection permit NWP models, it was a common practice to turn off convective parameterization. However, if there were not any cloud information in the initial conditions, the occurrence of precipitation could be delayed due to spin-up of cloud field or microphysical variables. In this study, we utilized the complex cloud analysis package from the Advanced Regional Prediction System (ARPS) to adjust the initial states of the model on water substance, such as cloud water, cloud ice, rain water, et al., that is, to initialize the microphysical variables (i.e., hydrometers), mainly based on radar reflectivity observations. Using the Advanced Research WRF (ARW) model, numerical experiments with/without cloud initialization and convective parameterization were carried out at grey-zone resolutions (i.e. 1, 3, and 9 km). The results from the experiments without convective parameterization indicate that model ignition with radar reflectivity can significantly reduce spin-up time and accurately simulate precipitation at the initial time. In addition, it helps to improve location and intensity of predicted precipitation. With grey-zone resolutions (i.e. 1, 3, and 9 km), using the cumulus convective parameterization scheme (without radar data) cannot produce realistic precipitation at the early time. The issues related to microphysical parametrization associated with cloud initialization were also discussed.

Assessing impacts of PBL and surface layer schemes in simulating the surface–atmosphere interactions and precipitation over the tropical ocean using observations from AMIE/DYNAMO

DOE PAGES

Qian, Yun; Yan, Huiping; Berg, Larry K.; ...

2016-10-28

Accuracy of turbulence parameterization in representing Planetary Boundary Layer (PBL) processes in climate models is critical for predicting the initiation and development of clouds, air quality issues, and underlying surface-atmosphere-cloud interactions. In this study, we 1) evaluate WRF model-simulated spatial patterns of precipitation and surface fluxes, as well as vertical profiles of potential temperature, humidity, moist static energy and moisture tendency terms as simulated by WRF at various spatial resolutions and with PBL, surface layer and shallow convection schemes against measurements, 2) identify model biases by examining the moisture tendency terms contributed by PBL and convection processes through nudging experiments,more » and 3) evaluate the dependence of modeled surface latent heat (LH) fluxes onPBL and surface layer schemes over the tropical ocean. The results show that PBL and surface parameterizations have surprisingly large impacts on precipitation, convection initiation and surface moisture fluxes over tropical oceans. All of the parameterizations tested tend to overpredict moisture in PBL and free atmosphere, and consequently result in larger moist static energy and precipitation. Moisture nudging tends to suppress the initiation of convection and reduces the excess precipitation. The reduction in precipitation bias in turn reduces the surface wind and LH flux biases, which suggests that the model drifts at least partly because of a positive feedback between precipitation and surface fluxes. The updated shallow convection scheme KF-CuP tends to suppress the initiation and development of deep convection, consequently decreasing precipitation. The Eta surface layer scheme predicts more reasonable LH fluxes and the LH-Wind Speed relationship than the MM5 scheme, especially when coupled with the MYJ scheme. By examining various parameterization schemes in WRF, we identify sources of biases and weaknesses of current PBL, surface layer and

Tropical Cumulus Convection and Upward Propagating Waves in Middle Atmospheric GCMs

NASA Technical Reports Server (NTRS)

Horinouchi, T.; Pawson, S.; Shibata, K.; Langematz, U.; Manzini, E.; Giorgetta, M. A.; Sassi, F.; Wilson, R. J.; Hamilton, K. P.; deGranpre, J.;

Comparison of Gravity Wave Temperature Variances from Ray-Based Spectral Parameterization of Convective Gravity Wave Drag with AIRS Observations

NASA Technical Reports Server (NTRS)

Choi, Hyun-Joo; Chun, Hye-Yeong; Gong, Jie; Wu, Dong L.

2012-01-01

The realism of ray-based spectral parameterization of convective gravity wave drag, which considers the updated moving speed of the convective source and multiple wave propagation directions, is tested against the Atmospheric Infrared Sounder (AIRS) onboard the Aqua satellite. Offline parameterization calculations are performed using the global reanalysis data for January and July 2005, and gravity wave temperature variances (GWTVs) are calculated at z = 2.5 hPa (unfiltered GWTV). AIRS-filtered GWTV, which is directly compared with AIRS, is calculated by applying the AIRS visibility function to the unfiltered GWTV. A comparison between the parameterization calculations and AIRS observations shows that the spatial distribution of the AIRS-filtered GWTV agrees well with that of the AIRS GWTV. However, the magnitude of the AIRS-filtered GWTV is smaller than that of the AIRS GWTV. When an additional cloud top gravity wave momentum flux spectrum with longer horizontal wavelength components that were obtained from the mesoscale simulations is included in the parameterization, both the magnitude and spatial distribution of the AIRS-filtered GWTVs from the parameterization are in good agreement with those of the AIRS GWTVs. The AIRS GWTV can be reproduced reasonably well by the parameterization not only with multiple wave propagation directions but also with two wave propagation directions of 45 degrees (northeast-southwest) and 135 degrees (northwest-southeast), which are optimally chosen for computational efficiency.

Atmospheric numerical modeling resource enhancement and model convective parameterization/scale interaction studies

NASA Technical Reports Server (NTRS)

Cushman, Paula P.

1993-01-01

Research will be undertaken in this contract in the area of Modeling Resource and Facilities Enhancement to include computer, technical and educational support to NASA investigators to facilitate model implementation, execution and analysis of output; to provide facilities linking USRA and the NASA/EADS Computer System as well as resident work stations in ESAD; and to provide a centralized location for documentation, archival and dissemination of modeling information pertaining to NASA's program. Additional research will be undertaken in the area of Numerical Model Scale Interaction/Convective Parameterization Studies to include implementation of the comparison of cloud and rain systems and convective-scale processes between the model simulations and what was observed; and to incorporate the findings of these and related research findings in at least two refereed journal articles.

The Super Tuesday Outbreak: Forecast Sensitivities to Single-Moment Microphysics Schemes

NASA Technical Reports Server (NTRS)

Molthan, Andrew L.; Case, Jonathan L.; Dembek, Scott R.; Jedlovec, Gary J.; Lapenta, William M.

2008-01-01

Forecast precipitation and radar characteristics are used by operational centers to guide the issuance of advisory products. As operational numerical weather prediction is performed at increasingly finer spatial resolution, convective precipitation traditionally represented by sub-grid scale parameterization schemes is now being determined explicitly through single- or multi-moment bulk water microphysics routines. Gains in forecasting skill are expected through improved simulation of clouds and their microphysical processes. High resolution model grids and advanced parameterizations are now available through steady increases in computer resources. As with any parameterization, their reliability must be measured through performance metrics, with errors noted and targeted for improvement. Furthermore, the use of these schemes within an operational framework requires an understanding of limitations and an estimate of biases so that forecasters and model development teams can be aware of potential errors. The National Severe Storms Laboratory (NSSL) Spring Experiments have produced daily, high resolution forecasts used to evaluate forecast skill among an ensemble with varied physical parameterizations and data assimilation techniques. In this research, high resolution forecasts of the 5-6 February 2008 Super Tuesday Outbreak are replicated using the NSSL configuration in order to evaluate two components of simulated convection on a large domain: sensitivities of quantitative precipitation forecasts to assumptions within a single-moment bulk water microphysics scheme, and to determine if these schemes accurately depict the reflectivity characteristics of well-simulated, organized, cold frontal convection. As radar returns are sensitive to the amount of hydrometeor mass and the distribution of mass among variably sized targets, radar comparisons may guide potential improvements to a single-moment scheme. In addition, object-based verification metrics are evaluated for

The sensitivity of WRF daily summertime simulations over West Africa to alternative parameterizations. Part 2: Precipitation.

PubMed

Noble, Erik; Druyan, Leonard M; Fulakeza, Matthew

2016-01-01

This paper evaluates the performance of the Weather and Research Forecasting (WRF) model as a regional-atmospheric model over West Africa. It tests WRF sensitivity to 64 configurations of alternative parameterizations in a series of 104 twelve-day September simulations during eleven consecutive years, 2000-2010. The 64 configurations combine WRF parameterizations of cumulus convection, radiation, surface-hydrology, and PBL. Simulated daily and total precipitation results are validated against Global Precipitation Climatology Project (GPCP) and Tropical Rainfall Measuring Mission (TRMM) data. Particular attention is given to westward-propagating precipitation maxima associated with African Easterly Waves (AEWs). A wide range of daily precipitation validation scores demonstrates the influence of alternative parameterizations. The best WRF performers achieve time-longitude correlations (against GPCP) of between 0.35-0.42 and spatiotemporal variability amplitudes only slightly higher than observed estimates. A parallel simulation by the benchmark Regional Model-v.3 achieves a higher correlation (0.52) and realistic spatiotemporal variability amplitudes. The largest favorable impact on WRF precipitation validation is achieved by selecting the Grell-Devenyi convection scheme, resulting in higher correlations against observations than using the Kain-Fritch convection scheme. Other parameterizations have less obvious impact. Validation statistics for optimized WRF configurations simulating the parallel period during 2000-2010 are more favorable for 2005, 2006, and 2008 than for other years. The selection of some of the same WRF configurations as high scorers in both circulation and precipitation validations supports the notion that simulations of West African daily precipitation benefit from skillful simulations of associated AEW vorticity centers and that simulations of AEWs would benefit from skillful simulations of convective precipitation.

The sensitivity of WRF daily summertime simulations over West Africa to alternative parameterizations. Part 2: Precipitation

PubMed Central

Noble, Erik; Druyan, Leonard M.; Fulakeza, Matthew

2018-01-01

This paper evaluates the performance of the Weather and Research Forecasting (WRF) model as a regional-atmospheric model over West Africa. It tests WRF sensitivity to 64 configurations of alternative parameterizations in a series of 104 twelve-day September simulations during eleven consecutive years, 2000–2010. The 64 configurations combine WRF parameterizations of cumulus convection, radiation, surface-hydrology, and PBL. Simulated daily and total precipitation results are validated against Global Precipitation Climatology Project (GPCP) and Tropical Rainfall Measuring Mission (TRMM) data. Particular attention is given to westward-propagating precipitation maxima associated with African Easterly Waves (AEWs). A wide range of daily precipitation validation scores demonstrates the influence of alternative parameterizations. The best WRF performers achieve time-longitude correlations (against GPCP) of between 0.35–0.42 and spatiotemporal variability amplitudes only slightly higher than observed estimates. A parallel simulation by the benchmark Regional Model-v.3 achieves a higher correlation (0.52) and realistic spatiotemporal variability amplitudes. The largest favorable impact on WRF precipitation validation is achieved by selecting the Grell-Devenyi convection scheme, resulting in higher correlations against observations than using the Kain-Fritch convection scheme. Other parameterizations have less obvious impact. Validation statistics for optimized WRF configurations simulating the parallel period during 2000–2010 are more favorable for 2005, 2006, and 2008 than for other years. The selection of some of the same WRF configurations as high scorers in both circulation and precipitation validations supports the notion that simulations of West African daily precipitation benefit from skillful simulations of associated AEW vorticity centers and that simulations of AEWs would benefit from skillful simulations of convective precipitation. PMID:29563651

A unified parameterization of clouds and turbulence using CLUBB and subcolumns in the Community Atmosphere Model

DOE PAGES

Thayer-Calder, K.; Gettelman, A.; Craig, C.; ...

2015-06-30

Most global climate models parameterize separate cloud types using separate parameterizations. This approach has several disadvantages, including obscure interactions between parameterizations and inaccurate triggering of cumulus parameterizations. Alternatively, a unified cloud parameterization uses one equation set to represent all cloud types. Such cloud types include stratiform liquid and ice cloud, shallow convective cloud, and deep convective cloud. Vital to the success of a unified parameterization is a general interface between clouds and microphysics. One such interface involves drawing Monte Carlo samples of subgrid variability of temperature, water vapor, cloud liquid, and cloud ice, and feeding the sample points into amore » microphysics scheme.This study evaluates a unified cloud parameterization and a Monte Carlo microphysics interface that has been implemented in the Community Atmosphere Model (CAM) version 5.3. Results describing the mean climate and tropical variability from global simulations are presented. The new model shows a degradation in precipitation skill but improvements in short-wave cloud forcing, liquid water path, long-wave cloud forcing, precipitable water, and tropical wave simulation. Also presented are estimations of computational expense and investigation of sensitivity to number of subcolumns.« less

A unified parameterization of clouds and turbulence using CLUBB and subcolumns in the Community Atmosphere Model

DOE PAGES

Thayer-Calder, Katherine; Gettelman, A.; Craig, Cheryl; ...

2015-12-01

Most global climate models parameterize separate cloud types using separate parameterizations.This approach has several disadvantages, including obscure interactions between parameterizations and inaccurate triggering of cumulus parameterizations. Alternatively, a unified cloud parameterization uses one equation set to represent all cloud types. Such cloud types include stratiform liquid and ice cloud, shallow convective cloud, and deep convective cloud. Vital to the success of a unified parameterization is a general interface between clouds and microphysics. One such interface involves drawing Monte Carlo samples of subgrid variability of temperature, water vapor, cloud liquid, and cloud ice, and feeding the sample points into a microphysicsmore » scheme. This study evaluates a unified cloud parameterization and a Monte Carlo microphysics interface that has been implemented in the Community Atmosphere Model (CAM) version 5.3. Results describing the mean climate and tropical variability from global simulations are presented. In conclusion, the new model shows a degradation in precipitation skill but improvements in short-wave cloud forcing, liquid water path, long-wave cloud forcing, perceptible water, and tropical wave simulation. Also presented are estimations of computational expense and investigation of sensitivity to number of subcolumns.« less

Modelling heterogeneous ice nucleation on mineral dust and soot with parameterizations based on laboratory experiments

NASA Astrophysics Data System (ADS)

Hoose, C.; Hande, L. B.; Mohler, O.; Niemand, M.; Paukert, M.; Reichardt, I.; Ullrich, R.

2016-12-01

Between 0 and -37°C, ice formation in clouds is triggered by aerosol particles acting as heterogeneous ice nuclei. At lower temperatures, heterogeneous ice nucleation on aerosols can occur at lower supersaturations than homogeneous freezing of solutes. In laboratory experiments, the ability of different aerosol species (e.g. desert dusts, soot, biological particles) has been studied in detail and quantified via various theoretical or empirical parameterization approaches. For experiments in the AIDA cloud chamber, we have quantified the ice nucleation efficiency via a temperature- and supersaturation dependent ice nucleation active site density. Here we present a new empirical parameterization scheme for immersion and deposition ice nucleation on desert dust and soot based on these experimental data. The application of this parameterization to the simulation of cirrus clouds, deep convective clouds and orographic clouds will be shown, including the extension of the scheme to the treatment of freezing of rain drops. The results are compared to other heterogeneous ice nucleation schemes. Furthermore, an aerosol-dependent parameterization of contact ice nucleation is presented.

Improvement of Systematic Bias of mean state and the intraseasonal variability of CFSv2 through superparameterization and revised cloud-convection-radiation parameterization

NASA Astrophysics Data System (ADS)

Mukhopadhyay, P.; Phani Murali Krishna, R.; Goswami, Bidyut B.; Abhik, S.; Ganai, Malay; Mahakur, M.; Khairoutdinov, Marat; Dudhia, Jimmy

2016-05-01

Inspite of significant improvement in numerical model physics, resolution and numerics, the general circulation models (GCMs) find it difficult to simulate realistic seasonal and intraseasonal variabilities over global tropics and particularly over Indian summer monsoon (ISM) region. The bias is mainly attributed to the improper representation of physical processes. Among all the processes, the cloud and convective processes appear to play a major role in modulating model bias. In recent times, NCEP CFSv2 model is being adopted under Monsoon Mission for dynamical monsoon forecast over Indian region. The analyses of climate free run of CFSv2 in two resolutions namely at T126 and T382, show largely similar bias in simulating seasonal rainfall, in capturing the intraseasonal variability at different scales over the global tropics and also in capturing tropical waves. Thus, the biases of CFSv2 indicate a deficiency in model's parameterization of cloud and convective processes. Keeping this in background and also for the need to improve the model fidelity, two approaches have been adopted. Firstly, in the superparameterization, 32 cloud resolving models each with a horizontal resolution of 4 km are embedded in each GCM (CFSv2) grid and the conventional sub-grid scale convective parameterization is deactivated. This is done to demonstrate the role of resolving cloud processes which otherwise remain unresolved. The superparameterized CFSv2 (SP-CFS) is developed on a coarser version T62. The model is integrated for six and half years in climate free run mode being initialised from 16 May 2008. The analyses reveal that SP-CFS simulates a significantly improved mean state as compared to default CFS. The systematic bias of lesser rainfall over Indian land mass, colder troposphere has substantially been improved. Most importantly the convectively coupled equatorial waves and the eastward propagating MJO has been found to be simulated with more fidelity in SP-CFS. The reason of

Chemical transport models: the combined non-local diffusion and mixing schemes, and calculation of in-canopy resistance for dry deposition fluxes.

PubMed

Mihailovic, Dragutin T; Alapaty, Kiran; Podrascanin, Zorica

2009-03-01

Improving the parameterization of processes in the atmospheric boundary layer (ABL) and surface layer, in air quality and chemical transport models. To do so, an asymmetrical, convective, non-local scheme, with varying upward mixing rates is combined with the non-local, turbulent, kinetic energy scheme for vertical diffusion (COM). For designing it, a function depending on the dimensionless height to the power four in the ABL is suggested, which is empirically derived. Also, we suggested a new method for calculating the in-canopy resistance for dry deposition over a vegetated surface. The upward mixing rate forming the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. The vertical eddy diffusivity is parameterized using the mean turbulent velocity scale that is obtained by the vertical integration within the ABL. In-canopy resistance is calculated by integration of inverse turbulent transfer coefficient inside the canopy from the effective ground roughness length to the canopy source height and, further, from its the canopy height. This combination of schemes provides a less rapid mass transport out of surface layer into other layers, during convective and non-convective periods, than other local and non-local schemes parameterizing mixing processes in the ABL. The suggested method for calculating the in-canopy resistance for calculating the dry deposition over a vegetated surface differs remarkably from the commonly used one, particularly over forest vegetation. In this paper, we studied the performance of a non-local, turbulent, kinetic energy scheme for vertical diffusion combined with a non-local, convective mixing scheme with varying upward mixing in the atmospheric boundary layer (COM) and its impact on the concentration of pollutants

A Double-Moment Multiple-Phase Four-Class Bulk Ice Scheme. Part II: Simulations of Convective Storms in Different Large-Scale Environments and Comparisons with other Bulk Parameterizations.

NASA Astrophysics Data System (ADS)

Schoenberg Ferrier, Brad; Tao, Wei-Kuo; Simpson, Joanne

1995-04-01

Part I of this study described a detailed four-class bulk ice scheme (4ICE) developed to simulate the hydro-meteor profiles of convective and stratiform precipitation associated with mesoscale convective systems. In Part II, the 4ICE scheme is incorporated into the Goddard Cumulus Ensemble (GCE) model and applied without any `tuning' to two squall lines occurring in widely different environments, namely, one over the `Pica) ocean in the Global Atmospheric Research Program's (GARP) Atlantic Tropical Experiment (GATE) and the other over a midlatitude continent in the Cooperative Huntsville Meteorological Experiment (COHMEX). Comparisons were made both with earlier three-class ice formulations and with observations. In both cases, the 4ICE scheme interacted with the dynamics so as to resemble the observations much more closely than did the model runs with either of the three-class ice parameterizations. The following features were well simulated in the COHMEX case: a lack of stratiform rain at the surface ahead of the storm, reflectivity maxima near 60 dBZ in the vicinity of the melting level, and intense radar echoes up to near the tropopause. These features were in strong contrast with the GATE simulation, which showed extensive trailing stratiform precipitation containing a horizontally oriented radar bright band. Peak reflectivities were below the melting level, rarely exceeding 50 dBz, with a steady decrease in reflectivity with height above. With the other bulk formulations, the large stratiform rain areas were not reproduced in the GATE conditions.The microphysical structure of the model clouds in both environments were more realistic than that of earlier modeling efforts. Number concentrations of ice of O(100 L1) occurred above 6 km in the GATE model clouds as a result of ice enhancement and rime splintering in the 4ICE runs. These processes were more effective in the GATE simulation, because near the freezing level the weaker updrafts were comparable in

Sensitivity of Numerical Simulations of a Mesoscale Convective System to Ice Hydrometeors in Bulk Microphysical Parameterization

NASA Astrophysics Data System (ADS)

Pu, Zhaoxia; Lin, Chao; Dong, Xiquan; Krueger, Steven K.

2018-01-01

Mesoscale convective systems (MCSs) and their associated cloud properties are the important factors that influence the aviation activities, yet they present a forecasting challenge in numerical weather prediction. In this study, the sensitivity of numerical simulations of an MCS over the US Southern Great Plains to ice hydrometeors in bulk microphysics (MP) schemes has been investigated using the Weather Research and Forecasting (WRF) model. It is found that the simulated structure, life cycle, cloud coverage, and precipitation of the convective system as well as its associated cold pools are sensitive to three selected MP schemes, namely, the WRF single-moment 6-class (WSM6), WRF double-moment 6-class (WDM6, with the double-moment treatment of warm-rain only), and Morrison double-moment (MORR, with the double-moment representation of both warm-rain and ice) schemes. Compared with observations, the WRF simulation with WSM6 only produces a less organized convection structure with a short lifetime, while WDM6 can produce the structure and length of the MCS very well. Both simulations heavily underestimate the precipitation amount, the height of the radar echo top, and stratiform cloud fractions. With MORR, the model performs well in predicting the lifetime, cloud coverage, echo top, and precipitation amount of the convection. Overall results demonstrate the importance of including double-moment representation of ice hydrometeors along with warm-rain. Additional experiments are performed to further examine the role of ice hydrometeors in numerical simulations of the MCS. Results indicate that replacing graupel with hail in the MORR scheme improves the prediction of the convective structure, especially in the convective core region.

Experimental validation of convection-diffusion discretisation scheme employed for computational modelling of biological mass transport

PubMed Central

2010-01-01

Background The finite volume solver Fluent (Lebanon, NH, USA) is a computational fluid dynamics software employed to analyse biological mass-transport in the vasculature. A principal consideration for computational modelling of blood-side mass-transport is convection-diffusion discretisation scheme selection. Due to numerous discretisation schemes available when developing a mass-transport numerical model, the results obtained should either be validated against benchmark theoretical solutions or experimentally obtained results. Methods An idealised aneurysm model was selected for the experimental and computational mass-transport analysis of species concentration due to its well-defined recirculation region within the aneurysmal sac, allowing species concentration to vary slowly with time. The experimental results were obtained from fluid samples extracted from a glass aneurysm model, using the direct spectrophometric concentration measurement technique. The computational analysis was conducted using the four convection-diffusion discretisation schemes available to the Fluent user, including the First-Order Upwind, the Power Law, the Second-Order Upwind and the Quadratic Upstream Interpolation for Convective Kinetics (QUICK) schemes. The fluid has a diffusivity of 3.125 × 10-10 m2/s in water, resulting in a Peclet number of 2,560,000, indicating strongly convection-dominated flow. Results The discretisation scheme applied to the solution of the convection-diffusion equation, for blood-side mass-transport within the vasculature, has a significant influence on the resultant species concentration field. The First-Order Upwind and the Power Law schemes produce similar results. The Second-Order Upwind and QUICK schemes also correlate well but differ considerably from the concentration contour plots of the First-Order Upwind and Power Law schemes. The computational results were then compared to the experimental findings. An average error of 140% and 116% was demonstrated

An evaluation of the schemes of ocean surface albedo parameterization in shortwave radiation estimation

NASA Astrophysics Data System (ADS)

Niu, Hailin; Zhang, Xiaotong; Liu, Qiang; Feng, Youbin; Li, Xiuhong; Zhang, Jialin; Cai, Erli

2015-12-01

The ocean surface albedo (OSA) is a deciding factor on ocean net surface shortwave radiation (ONSSR) estimation. Several OSA schemes have been proposed successively, but there is not a conclusion for the best OSA scheme of estimating the ONSSR. On the base of analyzing currently existing OSA parameterization, including Briegleb et al.(B), Taylor et al.(T), Hansen et al.(H), Jin et al.(J), Preisendorfer and Mobley(PM86), Feng's scheme(F), this study discusses the difference of OSA's impact on ONSSR estimation in condition of actual downward shortwave radiation(DSR). Then we discussed the necessity and applicability for the climate models to integrate the more complicated OSA scheme. It is concluded that the SZA and the wind speed are the two most significant effect factor to broadband OSA, thus the different OSA parameterizations varies violently in the regions of both high latitudes and strong winds. The OSA schemes can lead the ONSSR results difference of the order of 20 w m-2. The Taylor's scheme shows the best estimate, and Feng's result just following Taylor's. However, the accuracy of the estimated instantaneous OSA changes at different local time. Jin's scheme has the best performance generally at noon and in the afternoon, and PM86's is the best of all in the morning, which indicate that the more complicated OSA schemes reflect the temporal variation of OWA better than the simple ones.

Collaborative Research: Reducing tropical precipitation biases in CESM — Tests of unified parameterizations with ARM observations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Larson, Vincent; Gettelman, Andrew; Morrison, Hugh

In state-of-the-art climate models, each cloud type is treated using its own separate cloud parameterization and its own separate microphysics parameterization. This use of separate schemes for separate cloud regimes is undesirable because it is theoretically unfounded, it hampers interpretation of results, and it leads to the temptation to overtune parameters. In this grant, we are creating a climate model that contains a unified cloud parameterization and a unified microphysics parameterization. This model will be used to address the problems of excessive frequency of drizzle in climate models and excessively early onset of deep convection in the Tropics over land.more » The resulting model will be compared with ARM observations.« less

Performance of multi-physics ensembles in convective precipitation events over northeastern Spain

NASA Astrophysics Data System (ADS)

García-Ortega, E.; Lorenzana, J.; Merino, A.; Fernández-González, S.; López, L.; Sánchez, J. L.

2017-07-01

Convective precipitation with hail greatly affects southwestern Europe, causing major economic losses. The local character of this meteorological phenomenon is a serious obstacle to forecasting. Therefore, the development of reliable short-term forecasts constitutes an essential challenge to minimizing and managing risks. However, deterministic outcomes are affected by different uncertainty sources, such as physics parameterizations. This study examines the performance of different combinations of physics schemes of the Weather Research and Forecasting model to describe the spatial distribution of precipitation in convective environments with hail falls. Two 30-member multi-physics ensembles, with two and three domains of maximum resolution 9 and 3km each, were designed using various combinations of cumulus, microphysics and radiation schemes. The experiment was evaluated for 10 convective precipitation days with hail over 2005-2010 in northeastern Spain. Different indexes were used to evaluate the ability of each ensemble member to capture the precipitation patterns, which were compared with observations of a rain-gauge network. A standardized metric was constructed to identify optimal performers. Results show interesting differences between the two ensembles. In two domain simulations, the selection of cumulus parameterizations was crucial, with the Betts-Miller-Janjic scheme the best. In contrast, the Kain-Fristch cumulus scheme gave the poorest results, suggesting that it should not be used in the study area. Nevertheless, in three domain simulations, the cumulus schemes used in coarser domains were not critical and the best results depended mainly on microphysics schemes. The best performance was shown by Morrison, New Thomson and Goddard microphysics.

CONVECTION IN CONDENSIBLE-RICH ATMOSPHERES

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ding, F.; Pierrehumbert, R. T., E-mail: fding@uchicago.edu

2016-05-01

Condensible substances are nearly ubiquitous in planetary atmospheres. For the most familiar case—water vapor in Earth’s present climate—the condensible gas is dilute, in the sense that its concentration is everywhere small relative to the noncondensible background gases. A wide variety of important planetary climate problems involve nondilute condensible substances. These include planets near or undergoing a water vapor runaway and planets near the outer edge of the conventional habitable zone, for which CO{sub 2} is the condensible. Standard representations of convection in climate models rely on several approximations appropriate only to the dilute limit, while nondilute convection differs in fundamentalmore » ways from dilute convection. In this paper, a simple parameterization of convection valid in the nondilute as well as dilute limits is derived and used to discuss the basic character of nondilute convection. The energy conservation properties of the scheme are discussed in detail and are verified in radiative-convective simulations. As a further illustration of the behavior of the scheme, results for a runaway greenhouse atmosphere for both steady instellation and seasonally varying instellation corresponding to a highly eccentric orbit are presented. The latter case illustrates that the high thermal inertia associated with latent heat in nondilute atmospheres can damp out the effects of even extreme seasonal forcing.« less

Comparison of 2-3D convection models with parameterized thermal evolution models: Application to Mars

NASA Astrophysics Data System (ADS)

Thiriet, M.; Plesa, A. C.; Breuer, D.; Michaut, C.

2017-12-01

To model the thermal evolution of terrestrial planets, 1D parametrized models are often used as 2 or 3D mantle convection codes are very time-consuming. In these parameterized models, scaling laws that describe the convective heat transfer rate as a function of the convective parameters are derived from 2-3D steady state convection models. However, so far there has been no comprehensive comparison whether they can be applied to model the thermal evolution of a cooling planet. Here we compare 2D and 3D thermal evolution models in the stagnant lid regime with 1D parametrized models and use parameters representing the cooling of the Martian mantle. For the 1D parameterized models, we use the approach of Grasset and Parmentier (1998) and treat the stagnant lid and the convecting layer separately. In the convecting layer, the scaling law for a fluid with constant viscosity is valid with Nu (Ra/Rac) ?, with Rac the critical Rayleigh number at which the thermal boundary layers (TBL) - top or bottom - destabilize. ? varies between 1/3 and 1/4 depending on the heating mode and previous studies have proposed intermediate values of b 0.28-0.32 according to their model set-up. The base of the stagnant lid is defined by the temperature at which the mantle viscosity has increased by a factor of 10; it thus depends on the rate of viscosity change with temperature multiplied by a factor? , whose value appears to vary depending on the geometry and convection conditions. In applying Monte Carlo simulations, we search for the best fit to temperature profiles and heat flux using three free parameters, i.e. ? of the upper TBL, ? and the Rac of the lower TBL. We find that depending on the definition of the stagnant lid thickness in the 2-3D models several combinations of ? and ? for the upper TBL can retrieve suitable fits. E.g. combinations of ? = 0.329 and ? = 2.19 but also ? = 0.295 and ? = 2.97 are possible; Rac of the lower TBL is 10 for all best fits. The results show that

ED(MF)n: Humidity-Convection Feedbacks in a Mass Flux Scheme Based on Resolved Size Densities

NASA Astrophysics Data System (ADS)

Neggers, R.

2014-12-01

Cumulus cloud populations remain at least partially unresolved in present-day numerical simulations of global weather and climate, and accordingly their impact on the larger-scale flow has to be represented through parameterization. Various methods have been developed over the years, ranging in complexity from the early bulk models relying on a single plume to more recent approaches that attempt to reconstruct the underlying probability density functions, such as statistical schemes and multiple plume approaches. Most of these "classic" methods capture key aspects of cumulus cloud populations, and have been successfully implemented in operational weather and climate models. However, the ever finer discretizations of operational circulation models, driven by advances in the computational efficiency of supercomputers, is creating new problems for existing sub-grid schemes. Ideally, a sub-grid scheme should automatically adapt its impact on the resolved scales to the dimension of the grid-box within which it is supposed to act. It can be argued that this is only possible when i) the scheme is aware of the range of scales of the processes it represents, and ii) it can distinguish between contributions as a function of size. How to conceptually represent this knowledge of scale in existing parameterization schemes remains an open question that is actively researched. This study considers a relatively new class of models for sub-grid transport in which ideas from the field of population dynamics are merged with the concept of multi plume modelling. More precisely, a multiple mass flux framework for moist convective transport is formulated in which the ensemble of plumes is created in "size-space". It is argued that thus resolving the underlying size-densities creates opportunities for introducing scale-awareness and scale-adaptivity in the scheme. The behavior of an implementation of this framework in the Eddy Diffusivity Mass Flux (EDMF) model, named ED(MF)n, is

Evaluation of simulated tropical convective updraft hydrometeor properties using aircraft observations

NASA Astrophysics Data System (ADS)

Stanford, McKenna W.

The High Altitude Ice Crystals - High Ice Water Content (HAIC-HIWC) field campaign produced aircraft retrievals of total condensed water content (TWC), hydrometeor particle size distributions, and vertical velocity (w) in high ice water content regions of tropical mesoscale convective systems (MCSs). These observations are used to evaluate deep convective updraft properties in high-resolution nested Weather Research and Forecasting (WRF) simulations of observed MCSs. Because simulated hydrometeor properties are highly sensitive to the parameterization of microphysics, three commonly used microphysical parameterizations are tested, including two bulk schemes (Thompson and Morrison) and one bin scheme (Fast Spectral Bin Microphysics). A commonly documented bias in cloud-resolving simulations is the exaggeration of simulated radar reflectivities aloft in tropical MCSs. This may result from overly strong convective updrafts that loft excessive condensate mass and from simplified approximations of hydrometeor size distributions, properties, species separation, and microphysical processes. The degree to which the reflectivity bias is a separate function of convective dynamics, condensate mass, and hydrometeor size has yet to be addressed. This research untangles these components by comparing simulated and observed relationships between w, TWC, and hydrometer size as a function of temperature. All microphysics schemes produce median mass diameters that are generally larger than observed for temperatures between -10 °C and -40 °C and TWC > 1 g m-3. Observations produce a prominent mode in the composite mass size distribution around 300 microm, but under most conditions, all schemes shift the distribution mode to larger sizes. Despite a much greater number of samples, all simulations fail to reproduce observed high TWC or high w conditions between -20 °C and -40 °C in which only a small fraction of condensate mass is found in relatively large particle sizes. Increasing

Impacts of Microphysical Scheme on Convective and Stratiform Characteristics in Two High Precipitation Squall Line Events

NASA Technical Reports Server (NTRS)

Wu, Di; Dong, Xiquan; Xi, Baike; Feng, Zhe; Kennedy, Aaron; Mullendore, Gretchen; Gilmore, Matthew; Tao, Wei-Kuo

2013-01-01

This study investigates the impact of snow, graupel, and hail processes on simulated squall lines over the Southern Great Plains in the United States. The Weather Research and Forecasting (WRF) model is used to simulate two squall line events in Oklahoma during May 2007, and the simulations are validated against radar and surface observations. Several microphysics schemes are tested in this study, including the WRF 5-Class Microphysics (WSM5), WRF 6-Class Microphysics (WSM6), Goddard Cumulus Ensemble (GCE) Three Ice (3-ice) with graupel, Goddard Two Ice (2-ice), and Goddard 3-ice hail schemes. Simulated surface precipitation is sensitive to the microphysics scheme when the graupel or hail categories are included. All of the 3-ice schemes overestimate the total precipitation with WSM6 having the largest bias. The 2-ice schemes, without a graupel/hail category, produce less total precipitation than the 3-ice schemes. By applying a radar-based convective/stratiform partitioning algorithm, we find that including graupel/hail processes increases the convective areal coverage, precipitation intensity, updraft, and downdraft intensities, and reduces the stratiform areal coverage and precipitation intensity. For vertical structures, simulations have higher reflectivity values distributed aloft than the observed values in both the convective and stratiform regions. Three-ice schemes produce more high reflectivity values in convective regions, while 2-ice schemes produce more high reflectivity values in stratiform regions. In addition, this study has demonstrated that the radar-based convective/stratiform partitioning algorithm can reasonably identify WRF-simulated precipitation, wind, and microphysical fields in both convective and stratiform regions.

A parameterization scheme for the x-ray linear attenuation coefficient and energy absorption coefficient.

PubMed

Midgley, S M

2004-01-21

A novel parameterization of x-ray interaction cross-sections is developed, and employed to describe the x-ray linear attenuation coefficient and mass energy absorption coefficient for both elements and mixtures. The new parameterization scheme addresses the Z-dependence of elemental cross-sections (per electron) using a simple function of atomic number, Z. This obviates the need for a complicated mathematical formalism. Energy dependent coefficients describe the Z-direction curvature of the cross-sections. The composition dependent quantities are the electron density and statistical moments describing the elemental distribution. We show that it is possible to describe elemental cross-sections for the entire periodic table and at energies above the K-edge (from 6 keV to 125 MeV), with an accuracy of better than 2% using a parameterization containing not more than five coefficients. For the biologically important elements 1 < or = Z < or = 20, and the energy range 30-150 keV, the parameterization utilizes four coefficients. At higher energies, the parameterization uses fewer coefficients with only two coefficients needed at megavoltage energies.

Prototype Mcs Parameterization for Global Climate Models

NASA Astrophysics Data System (ADS)

Moncrieff, M. W.

2017-12-01

Excellent progress has been made with observational, numerical and theoretical studies of MCS processes but the parameterization of those processes remain in a dire state and are missing from GCMs. The perceived complexity of the distribution, type, and intensity of organized precipitation systems has arguably daunted attention and stifled the development of adequate parameterizations. TRMM observations imply links between convective organization and large-scale meteorological features in the tropics and subtropics that are inadequately treated by GCMs. This calls for improved physical-dynamical treatment of organized convection to enable the next-generation of GCMs to reliably address a slew of challenges. The multiscale coherent structure parameterization (MCSP) paradigm is based on the fluid-dynamical concept of coherent structures in turbulent environments. The effects of vertical shear on MCS dynamics implemented as 2nd baroclinic convective heating and convective momentum transport is based on Lagrangian conservation principles, nonlinear dynamical models, and self-similarity. The prototype MCS parameterization, a minimalist proof-of-concept, is applied in the NCAR Community Climate Model, Version 5.5 (CAM 5.5). The MCSP generates convectively coupled tropical waves and large-scale precipitation features notably in the Indo-Pacific warm-pool and Maritime Continent region, a center-of-action for weather and climate variability around the globe.

Impact of Vegetation Cover Fraction Parameterization schemes on Land Surface Temperature Simulation in the Tibetan Plateau

NASA Astrophysics Data System (ADS)

Lv, M.; Li, C.; Lu, H.; Yang, K.; Chen, Y.

2017-12-01

The parameterization of vegetation cover fraction (VCF) is an important component of land surface models. This paper investigates the impacts of three VCF parameterization schemes on land surface temperature (LST) simulation by the Common Land Model (CoLM) in the Tibetan Plateau (TP). The first scheme is a simple land cover (LC) based method; the second one is based on remote sensing observation (hereafter named as RNVCF) , in which multi-year climatology VCFs is derived from Moderate-resolution Imaging Spectroradiometer (MODIS) NDVI (Normalized Difference Vegetation Index); the third VCF parameterization scheme derives VCF from the LAI simulated by LSM and clump index at every model time step (hereafter named as SMVCF). Simulated land surface temperature(LST) and soil temperature by CoLM with three VCF parameterization schemes were evaluated by using satellite LST observation and in situ soil temperature observation, respectively, during the period of 2010 to 2013. The comparison against MODIS Aqua LST indicates that (1) CTL produces large biases for both four seasons in early afternoon (about 13:30, local solar time), while the mean bias in spring reach to 12.14K; (2) RNVCF and SMVCF reduce the mean bias significantly, especially in spring as such reduce is about 6.5K. Surface soil temperature observed at 5 cm depth from three soil moisture and temperature monitoring networks is also employed to assess the skill of three VCF schemes. The three networks, crossing TP from West to East, have different climate and vegetation conditions. In the Ngari network, located in the Western TP with an arid climate, there are not obvious differences among three schemes. In Naqu network, located in central TP with a semi-arid climate condition, CTL shows a severe overestimates (12.1 K), but such overestimations can be reduced by 79% by RNVCF and 87% by SMVCF. In the third humid network (Maqu in eastern TP), CoLM performs similar to Naqu. However, at both Naqu and Maqu networks

Convective and microphysics parameterization impact on simulating heavy rainfall in Semarang (case study on February 12th, 2015)

NASA Astrophysics Data System (ADS)

Faridatussafura, Nurzaka; Wandala, Agie

2018-05-01

The meteorological model WRF-ARW version 3.8.1 is used for simulating the heavy rainfall in Semarang that occurred on February 12th, 2015. Two different convective schemes and two different microphysics scheme in a nested configuration were chosen. The sensitivity of those schemes in capturing the extreme weather event has been tested. GFS data were used for the initial and boundary condition. Verification on the twenty-four hours accumulated rainfall using GSMaPsatellite data shows that Kain-Fritsch convective scheme and Lin microphysics scheme is the best combination scheme among the others. The combination also gives the highest success ratio value in placing high intensity rainfall area. Based on the ROC diagram, KF-Lin shows the best performance in detecting high intensity rainfall. However, the combination still has high bias value.

The Sensitivity of WRF Daily Summertime Simulations over West Africa to Alternative Parameterizations. Part 1: African Wave Circulation

NASA Technical Reports Server (NTRS)

Noble, Erik; Druyan, Leonard M.; Fulakeza, Matthew

2014-01-01

The performance of the NCAR Weather Research and Forecasting Model (WRF) as a West African regional-atmospheric model is evaluated. The study tests the sensitivity of WRF-simulated vorticity maxima associated with African easterly waves to 64 combinations of alternative parameterizations in a series of simulations in September. In all, 104 simulations of 12-day duration during 11 consecutive years are examined. The 64 combinations combine WRF parameterizations of cumulus convection, radiation transfer, surface hydrology, and PBL physics. Simulated daily and mean circulation results are validated against NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA) and NCEP/Department of Energy Global Reanalysis 2. Precipitation is considered in a second part of this two-part paper. A wide range of 700-hPa vorticity validation scores demonstrates the influence of alternative parameterizations. The best WRF performers achieve correlations against reanalysis of 0.40-0.60 and realistic amplitudes of spatiotemporal variability for the 2006 focus year while a parallel-benchmark simulation by the NASA Regional Model-3 (RM3) achieves higher correlations, but less realistic spatiotemporal variability. The largest favorable impact on WRF-vorticity validation is achieved by selecting the Grell-Devenyi cumulus convection scheme, resulting in higher correlations against reanalysis than simulations using the Kain-Fritch convection. Other parameterizations have less-obvious impact, although WRF configurations incorporating one surface model and PBL scheme consistently performed poorly. A comparison of reanalysis circulation against two NASA radiosonde stations confirms that both reanalyses represent observations well enough to validate the WRF results. Validation statistics for optimized WRF configurations simulating the parallel period during 10 additional years are less favorable than for 2006.

Effects of Planetary Boundary Layer Parameterizations on CWRF Regional Climate Simulation

NASA Astrophysics Data System (ADS)

Liu, S.; Liang, X.

2011-12-01

Planetary Boundary Layer (PBL) parameterizations incorporated in CWRF (Climate extension of the Weather Research and Forecasting model) are first evaluated by comparing simulated PBL heights with observations. Among the 10 evaluated PBL schemes, 2 (CAM, UW) are new in CWRF while the other 8 are original WRF schemes. MYJ, QNSE and UW determine the PBL heights based on turbulent kinetic energy (TKE) profiles, while others (YSU, ACM, GFS, CAM, TEMF) are from bulk Richardson criteria. All TKE-based schemes (MYJ, MYNN, QNSE, UW, Boulac) substantially underestimate convective or residual PBL heights from noon toward evening, while others (ACM, CAM, YSU) well capture the observed diurnal cycle except for the GFS with systematic overestimation. These differences among the schemes are representative over most areas of the simulation domain, suggesting systematic behaviors of the parameterizations. Lower PBL heights simulated by the QNSE and MYJ are consistent with their smaller Bowen ratios and heavier rainfalls, while higher PBL tops by the GFS correspond to warmer surface temperatures. Effects of PBL parameterizations on CWRF regional climate simulation are then compared. The QNSE PBL scheme yields systematically heavier rainfall almost everywhere and throughout the year; this is identified with a much greater surface Bowen ratio (smaller sensible versus larger latent heating) and wetter soil moisture than other PBL schemes. Its predecessor MYJ scheme shares the same deficiency to a lesser degree. For temperature, the performance of the QNSE and MYJ schemes remains poor, having substantially larger rms errors in all seasons. GFS PBL scheme also produces large warm biases. Pronounced sensitivities are also found to the PBL schemes in winter and spring over most areas except the southern U.S. (Southeast, Gulf States, NAM); excluding the outliers (QNSE, MYJ, GFS) that cause extreme biases of -6 to +3°C, the differences among the schemes are still visible (±2°C), where the

Atmospheric parameterization schemes for satellite cloud property retrieval during FIRE IFO 2

NASA Technical Reports Server (NTRS)

Titlow, James; Baum, Bryan A.

1993-01-01

Satellite cloud retrieval algorithms generally require atmospheric temperature and humidity profiles to determine such cloud properties as pressure and height. For instance, the CO2 slicing technique called the ratio method requires the calculation of theoretical upwelling radiances both at the surface and a prescribed number (40) of atmospheric levels. This technique has been applied to data from, for example, the High Resolution Infrared Radiometer Sounder (HIRS/2, henceforth HIRS) flown aboard the NOAA series of polar orbiting satellites and the High Resolution Interferometer Sounder (HIS). In this particular study, four NOAA-11 HIRS channels in the 15-micron region are used. The ratio method may be applied to various channel combinations to estimate cloud top heights using channels in the 15-mu m region. Presently, the multispectral, multiresolution (MSMR) scheme uses 4 HIRS channel combination estimates for mid- to high-level cloud pressure retrieval and Advanced Very High Resolution Radiometer (AVHRR) data for low-level (is greater than 700 mb) cloud level retrieval. In order to determine theoretical upwelling radiances, atmospheric temperature and water vapor profiles must be provided as well as profiles of other radiatively important gas absorber constituents such as CO2, O3, and CH4. The assumed temperature and humidity profiles have a large effect on transmittance and radiance profiles, which in turn are used with HIRS data to calculate cloud pressure, and thus cloud height and temperature. For large spatial scale satellite data analysis, atmospheric parameterization schemes for cloud retrieval algorithms are usually based on a gridded product such as that provided by the European Center for Medium Range Weather Forecasting (ECMWF) or the National Meteorological Center (NMC). These global, gridded products prescribe temperature and humidity profiles for a limited number of pressure levels (up to 14) in a vertical atmospheric column. The FIRE IFO 2

A factorial assessment of the sensitivity of the BATS land-surface parameterization scheme. [BATS (Biosphere-Atmosphere Transfer Scheme)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Henderson-Sellers, A.

Land-surface schemes developed for incorporation into global climate models include parameterizations that are not yet fully validated and depend upon the specification of a large (20-50) number of ecological and soil parameters, the values of which are not yet well known. There are two methods of investigating the sensitivity of a land-surface scheme to prescribed values: simple one-at-a-time changes or factorial experiments. Factorial experiments offer information about interactions between parameters and are thus a more powerful tool. Here the results of a suite of factorial experiments are reported. These are designed (i) to illustrate the usefulness of this methodology andmore » (ii) to identify factors important to the performance of complex land-surface schemes. The Biosphere-Atmosphere Transfer Scheme (BATS) is used and its sensitivity is considered (a) to prescribed ecological and soil parameters and (b) to atmospheric forcing used in the off-line tests undertaken. Results indicate that the most important atmospheric forcings are mean monthly temperature and the interaction between mean monthly temperature and total monthly precipitation, although fractional cloudiness and other parameters are also important. The most important ecological parameters are vegetation roughness length, soil porosity, and a factor describing the sensitivity of the stomatal resistance of vegetation to the amount of photosynthetically active solar radiation and, to a lesser extent, soil and vegetation albedos. Two-factor interactions including vegetation roughness length are more important than many of the 23 specified single factors. The results of factorial sensitivity experiments such as these could form the basis for intercomparison of land-surface parameterization schemes and for field experiments and satellite-based observation programs aimed at improving evaluation of important parameters.« less

Characterizing the degree of convective clustering using radar reflectivity and its application to evaluating model simulations

NASA Astrophysics Data System (ADS)

Cheng, W. Y.; Kim, D.; Rowe, A.; Park, S.

2017-12-01

Despite the impact of mesoscale convective organization on the properties of convection (e.g., mixing between updrafts and environment), parameterizing the degree of convective organization has only recently been attempted in cumulus parameterization schemes (e.g., Unified Convection Scheme UNICON). Additionally, challenges remain in determining the degree of convective organization from observations and in comparing directly with the organization metrics in model simulations. This study addresses the need to objectively quantify the degree of mesoscale convective organization using high quality S-PolKa radar data from the DYNAMO field campaign. One of the most noticeable aspects of mesoscale convective organization in radar data is the degree of convective clustering, which can be characterized by the number and size distribution of convective echoes and the distance between them. We propose a method of defining contiguous convective echoes (CCEs) using precipitating convective echoes identified by a rain type classification algorithm. Two classification algorithms, Steiner et al. (1995) and Powell et al. (2016), are tested and evaluated against high-resolution WRF simulations to determine which method better represents the degree of convective clustering. Our results suggest that the CCEs based on Powell et al.'s algorithm better represent the dynamical properties of the convective updrafts and thus provide the basis of a metric for convective organization. Furthermore, through a comparison with the observational data, the WRF simulations driven by the DYNAMO large-scale forcing, similarly applied to UNICON Single Column Model simulations, will allow us to evaluate the ability of both WRF and UNICON to simulate convective clustering. This evaluation is based on the physical processes that are explicitly represented in WRF and UNICON, including the mechanisms leading to convective clustering, and the feedback to the convective properties.

Influence of the vertical mixing parameterization on the modeling results of the Arctic Ocean hydrology

NASA Astrophysics Data System (ADS)

Iakshina, D. F.; Golubeva, E. N.

2017-11-01

The vertical distribution of the hydrological characteristics in the upper ocean layer is mostly formed under the influence of turbulent and convective mixing, which are not resolved in the system of equations for large-scale ocean. Therefore it is necessary to include additional parameterizations of these processes into the numerical models. In this paper we carry out a comparative analysis of the different vertical mixing parameterizations in simulations of climatic variability of the Arctic water and sea ice circulation. The 3D regional numerical model for the Arctic and North Atlantic developed in the ICMMG SB RAS (Institute of Computational Mathematics and Mathematical Geophysics of the Siberian Branch of the Russian Academy of Science) and package GOTM (General Ocean Turbulence Model1,2, http://www.gotm.net/) were used as the numerical instruments . NCEP/NCAR reanalysis data were used for determination of the surface fluxes related to ice and ocean. The next turbulence closure schemes were used for the vertical mixing parameterizations: 1) Integration scheme based on the Richardson criteria (RI); 2) Second-order scheme TKE with coefficients Canuto-A3 (CANUTO); 3) First-order scheme TKE with coefficients Schumann and Gerz4 (TKE-1); 4) Scheme KPP5 (KPP). In addition we investigated some important characteristics of the Arctic Ocean state including the intensity of Atlantic water inflow, ice cover state and fresh water content in Beaufort Sea.

Cloud Simulations in Response to Turbulence Parameterizations in the GISS Model E GCM

NASA Technical Reports Server (NTRS)

Yao, Mao-Sung; Cheng, Ye

2013-01-01

The response of cloud simulations to turbulence parameterizations is studied systematically using the GISS general circulation model (GCM) E2 employed in the Intergovernmental Panel on Climate Change's (IPCC) Fifth Assessment Report (AR5).Without the turbulence parameterization, the relative humidity (RH) and the low cloud cover peak unrealistically close to the surface; with the dry convection or with only the local turbulence parameterization, these two quantities improve their vertical structures, but the vertical transport of water vapor is still weak in the planetary boundary layers (PBLs); with both local and nonlocal turbulence parameterizations, the RH and low cloud cover have better vertical structures in all latitudes due to more significant vertical transport of water vapor in the PBL. The study also compares the cloud and radiation climatologies obtained from an experiment using a newer version of turbulence parameterization being developed at GISS with those obtained from the AR5 version. This newer scheme differs from the AR5 version in computing nonlocal transports, turbulent length scale, and PBL height and shows significant improvements in cloud and radiation simulations, especially over the subtropical eastern oceans and the southern oceans. The diagnosed PBL heights appear to correlate well with the low cloud distribution over oceans. This suggests that a cloud-producing scheme needs to be constructed in a framework that also takes the turbulence into consideration.

The Double ITCZ Syndrome in GCMs: A Coupled Problem among Convection, Atmospheric and Ocean Circulations

NASA Astrophysics Data System (ADS)

Zhang, G. J.; Song, X.

2017-12-01

The double ITCZ bias has been a long-standing problem in coupled atmosphere-ocean models. A previous study indicates that uncertainty in the projection of global warming due to doubling of CO2 is closely related to the double ITCZ biases in global climate models. Thus, reducing the double ITCZ biases is not only important to getting the current climate features right, but also important to narrowing the uncertainty in future climate projection. In this work, we will first review the possible factors contributing to the ITCZ problem. Then, we will focus on atmospheric convection, presenting recent progress in alleviating the double ITCZ problem and its sensitivity to details of convective parameterization, including trigger conditions for convection onset, convective memory, entrainment rate, updraft model and closure in the NCAR CESM1. These changes together can result in dramatic improvements in the simulation of ITCZ. Results based on both atmospheric only and coupled simulations with incremental changes of convection scheme will be shown to demonstrate the roles of convection parameterization and coupled interaction between convection, atmospheric circulation and ocean circulation in the simulation of ITCZ.

Analysis and evaluation of WRF microphysical schemes for deep moist convection over south-eastern South America (SESA) using microwave satellite observations and radiative transfer simulations

NASA Astrophysics Data System (ADS)

Sol Galligani, Victoria; Wang, Die; Alvarez Imaz, Milagros; Salio, Paola; Prigent, Catherine

2017-10-01

In the present study, three meteorological events of extreme deep moist convection, characteristic of south-eastern South America, are considered to conduct a systematic evaluation of the microphysical parameterizations available in the Weather Research and Forecasting (WRF) model by undertaking a direct comparison between satellite-based simulated and observed microwave radiances. A research radiative transfer model, the Atmospheric Radiative Transfer Simulator (ARTS), is coupled with the WRF model under three different microphysical parameterizations (WSM6, WDM6 and Thompson schemes). Microwave radiometry has shown a promising ability in the characterization of frozen hydrometeors. At high microwave frequencies, however, frozen hydrometeors significantly scatter radiation, and the relationship between radiation and hydrometeor populations becomes very complex. The main difficulty in microwave remote sensing of frozen hydrometeor characterization is correctly characterizing this scattering signal due to the complex and variable nature of the size, composition and shape of frozen hydrometeors. The present study further aims at improving the understanding of frozen hydrometeor optical properties characteristic of deep moist convection events in south-eastern South America. In the present study, bulk optical properties are computed by integrating the single-scattering properties of the Liu(2008) discrete dipole approximation (DDA) single-scattering database across the particle size distributions parameterized by the different WRF schemes in a consistent manner, introducing the equal mass approach. The equal mass approach consists of describing the optical properties of the WRF snow and graupel hydrometeors with the optical properties of habits in the DDA database whose dimensions might be different (Dmax') but whose mass is conserved. The performance of the radiative transfer simulations is evaluated by comparing the simulations with the available coincident

A controlled variation scheme for convection treatment in pressure-based algorithm

NASA Technical Reports Server (NTRS)

Shyy, Wei; Thakur, Siddharth; Tucker, Kevin

1993-01-01

Convection effect and source terms are two primary sources of difficulties in computing turbulent reacting flows typically encountered in propulsion devices. The present work intends to elucidate the individual as well as the collective roles of convection and source terms in the fluid flow equations, and to devise appropriate treatments and implementations to improve our current capability of predicting such flows. A controlled variation scheme (CVS) has been under development in the context of a pressure-based algorithm, which has the characteristics of adaptively regulating the amount of numerical diffusivity, relative to central difference scheme, according to the variation in local flow field. Both the basic concepts and a pragmatic assessment will be presented to highlight the status of this work.

Analysis of Cloud-resolving Simulations of a Tropical Mesoscale Convective System Observed during TWP-ICE: Vertical Fluxes and Draft Properties in Convective and Stratiform Regions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mrowiec, Agnieszka A.; Rio, Catherine; Fridlind, Ann

2012-10-02

We analyze three cloud-resolving model simulations of a strong convective event observed during the TWP-ICE campaign, differing in dynamical core, microphysical scheme or both. Based on simulated and observed radar reflectivity, simulations roughly reproduce observed convective and stratiform precipitating areas. To identify the characteristics of convective and stratiform drafts that are difficult to observe but relevant to climate model parameterization, independent vertical wind speed thresholds are calculated to capture 90% of total convective and stratiform updraft and downdraft mass fluxes. Convective updrafts are fairly consistent across simulations (likely owing to fixed large-scale forcings and surface conditions), except that hydrometeor loadingsmore » differ substantially. Convective downdraft and stratiform updraft and downdraft mass fluxes vary notably below the melting level, but share similar vertically uniform draft velocities despite differing hydrometeor loadings. All identified convective and stratiform downdrafts contain precipitation below ~10 km and nearly all updrafts are cloudy above the melting level. Cold pool properties diverge substantially in a manner that is consistent with convective downdraft mass flux differences below the melting level. Despite differences in hydrometeor loadings and cold pool properties, convective updraft and downdraft mass fluxes are linearly correlated with convective area, the ratio of ice in downdrafts to that in updrafts is ~0.5 independent of species, and the ratio of downdraft to updraft mass flux is ~0.5-0.6, which may represent a minimum evaporation efficiency under moist conditions. Hydrometeor loading in stratiform regions is found to be a fraction of hydrometeor loading in convective regions that ranges from ~10% (graupel) to ~90% (cloud ice). These findings may lead to improved convection parameterizations.« less

Ensemble superparameterization versus stochastic parameterization: A comparison of model uncertainty representation in tropical weather prediction

NASA Astrophysics Data System (ADS)

Subramanian, Aneesh C.; Palmer, Tim N.

2017-06-01

Stochastic schemes to represent model uncertainty in the European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble prediction system has helped improve its probabilistic forecast skill over the past decade by both improving its reliability and reducing the ensemble mean error. The largest uncertainties in the model arise from the model physics parameterizations. In the tropics, the parameterization of moist convection presents a major challenge for the accurate prediction of weather and climate. Superparameterization is a promising alternative strategy for including the effects of moist convection through explicit turbulent fluxes calculated from a cloud-resolving model (CRM) embedded within a global climate model (GCM). In this paper, we compare the impact of initial random perturbations in embedded CRMs, within the ECMWF ensemble prediction system, with stochastically perturbed physical tendency (SPPT) scheme as a way to represent model uncertainty in medium-range tropical weather forecasts. We especially focus on forecasts of tropical convection and dynamics during MJO events in October-November 2011. These are well-studied events for MJO dynamics as they were also heavily observed during the DYNAMO field campaign. We show that a multiscale ensemble modeling approach helps improve forecasts of certain aspects of tropical convection during the MJO events, while it also tends to deteriorate certain large-scale dynamic fields with respect to stochastically perturbed physical tendencies approach that is used operationally at ECMWF.