Breakdown of hydrostatic balance at convective scales in the forecast errors in the Met Office Unified Model

For data assimilation in numerical weather prediction, the initial forecast-error covariance matrix Pf is required. For variational assimilation it is particularly important to prescribe an accurate initial matrix Pf, since Pf is either static (in the 3D-Var case) or constant at the beginning of each assimilation window (in the 4D-Var case). At large scales the atmospheric flow is well approximated by hydrostatic balance and this balance is strongly enforced in the initial matrix Pf used in operational variational assimilation systems such as that of the Met Office. However, at convective scales this balance does not necessarily hold any more. Here we examine the extent to which hydrostatic balance is valid in the vertical forecast-error covariances for high-resolution models in order to determine whether there is a need to relax this balance constraint in convective-scale data assimilation. We use the Met Office Global and Regional Ensemble Prediction System (MOGREPS) and a 1.5 km resolution version of the Unified Model for a case study characterized by the presence of convective activity. An ensemble of high-resolution forecasts valid up to three hours after the onset of convection is produced. We show that at 1.5 km resolution hydrostatic balance does not hold for forecast errors in regions of convection. This indicates that in the presence of convection hydrostatic balance should not be enforced in the covariance matrix used for variational data assimilation at this scale. The results show the need to investigate covariance models that may be better suited for convective-scale data assimilation. Finally, we give a measure of the balance present in the forecast perturbations as a function of the horizontal scale (from 3–90 km) using a set of diagnostics. Copyright © 2012 Royal Meteorological Society and British Crown Copyright, the Met Office

Departures from convective equilibrium with a rapidly-varying surface forcing

Convective equilibrium is a long-standing and useful concept for understanding many aspects of the behaviour of deep moist convection. For example, it is often invoked in developing parameterizations for large-scale models. However, the equilibrium assumption may begin to break down as models are increasingly used with shorter timesteps and finer resolutions. Here we perform idealized cloud-system resolving model simulations of deep convection with imposed time variations in the surface forcing. A range of rapid forcing timescales from 1 − 36hr are used, in order to induce systematic departures from equilibrium. For the longer forcing timescales, the equilibrium assumption remains valid, in at least the limited sense that cycle-integrated measures of convective activity are very similar from cycle to cycle. For shorter forcing timescales, cycle-integrated convection becomes more variable, with enhanced activity on one cycle being correlated with reduced activity on the next, suggesting a role for convective memory. Further investigation shows that the memory does not appear to be carried by the domain-mean thermodynamic fields but rather by structures on horizontal scales of 5 − 20km. Such structures are produced by the convective clouds and can persist beyond the lifetime of the cloud, even through to the next forcing cycle.

Mixing length controls on high resolution simulations of convective storms

We perform simulations of several convective events over the southern UK with the Met Office Unified Model (UM) at horizontal grid lengths ranging from 1.5 km to 200 m. Comparing the simulated storms on these days with the Met Office rainfall radar network allows us to apply a statistical approach to evaluate the properties and evolution of the simulated storms over a range of conditions. Here we present results comparing the storm morphology in the model and reality which show that the simulated storms become smaller as grid length decreases and that the grid length that fits the observations best changes with the size of the observed cells. We investigate the sensitivity of storm morphology in the model to the mixing length used in the subgrid turbulence scheme. As the subgrid mixing length is decreased, the number of small storms with high area-averaged rain rates increases. We show that by changing the mixing length we can produce a lower resolution simulation that produces similar morphologies to a higher resolution simulation.

The three-dimensional morphology of simulated and observed convective storms over southern England

A set of high-resolution radar observations of convective storms has been collected to evaluate such storms in the UK Met Office Unified Model during the DYMECS project (Dynamical and Microphysical Evolution of Convective Storms). The 3-GHz Chilbolton Advanced Meteorological Radar was set up with a scan-scheduling algorithm to automatically track convective storms identified in real-time from the operational rainfall radar network. More than 1,000 storm observations gathered over fifteen days in 2011 and 2012 are used to evaluate the model under various synoptic conditions supporting convection. In terms of the detailed three-dimensional morphology, storms in the 1500-m grid-length simulations are shown to produce horizontal structures a factor 1.5–2 wider compared to radar observations. A set of nested model runs at grid lengths down to 100m show that the models converge in terms of storm width, but the storm structures in the simulations with the smallest grid lengths are too narrow and too intense compared to the radar observations. The modelled storms were surrounded by a region of drizzle without ice reflectivities above 0 dBZ aloft, which was related to the dominance of ice crystals and was improved by allowing only aggregates as an ice particle habit. Simulations with graupel outperformed the standard configuration for heavy-rain profiles, but the storm structures were a factor 2 too wide and the convective cores 2 km too deep.

The DYMECS project: a statistical approach for the evaluation of convective storms in high-resolution NWP models

A new frontier in weather forecasting is emerging by operational forecast models now being run at convection-permitting resolutions at many national weather services. However, this is not a panacea; significant systematic errors remain in the character of convective storms and rainfall distributions. The DYMECS project (Dynamical and Microphysical Evolution of Convective Storms) is taking a fundamentally new approach to evaluate and improve such models: rather than relying on a limited number of cases, which may not be representative, we have gathered a large database of 3D storm structures on 40 convective days using the Chilbolton radar in southern England. We have related these structures to storm life-cycles derived by tracking features in the rainfall from the UK radar network, and compared them statistically to storm structures in the Met Office model, which we ran at horizontal grid length between 1.5 km and 100 m, including simulations with different subgrid mixing length. We also evaluated the scale and intensity of convective updrafts using a new radar technique. We find that the horizontal size of simulated convective storms and the updrafts within them is much too large at 1.5-km resolution, such that the convective mass flux of individual updrafts can be too large by an order of magnitude. The scale of precipitation cores and updrafts decreases steadily with decreasing grid lengths, as does the typical storm lifetime. The 200-m grid-length simulation with standard mixing length performs best over all diagnostics, although a greater mixing length improves the representation of deep convective storms.

Convective updraught evaluation in high-resolution NWP simulations using single-Doppler radar measurements

This study presents an evaluation of the size and strength of convective updraughts in high-resolution simulations by the UK Met Office Unified Model (UM). Updraught velocities have been estimated from range–height indicator (RHI) Doppler velocity measurements using the Chilbolton advanced meteorological radar, as part of the Dynamical and Microphysical Evolution of Convective Storms (DYMECS) project. Based on mass continuity and the vertical integration of the observed radial convergence, vertical velocities tend to be underestimated for convective clouds due to the undetected cross-radial convergence. Velocity fields from the UM at a resolution corresponding to the radar observations are used to scale such estimates to mitigate the inherent biases. The analysis of more than 100 observed and simulated storms indicates that the horizontal scale of updraughts in simulations tend to decrease with grid length; the 200 m grid length agreed most closely with the observations. Typical updraught mass fluxes in the 500 m grid length simulations were up to an order of magnitude greater than observed, and greater still in the 1.5 km grid length simulations. The effect of increasing the mixing length in the sub-grid turbulence scheme depends on the grid length. For the 1.5 km simulations, updraughts were weakened though their horizontal scale remained largely unchanged. Progressively more so for the sub-kilometre grid lengths, updraughts were broadened and intensified; horizontal scale was now determined by the mixing length rather than the grid length. In general, simulated updraughts were found to weaken too quickly with height. The findings were supported by the analysis of the widths of reflectivity patterns in both the simulations and observations.

The sensitivity of convective aggregation to diabatic processes in idealized radiative-convective equilibrium simulations

Idealized explicit convection simulations of the Met Office Unified Model exhibit spontaneous self-aggregation in radiative-convective equilibrium, as seen in other models in previous studies. This self-aggregation is linked to feedbacks between radiation, surface fluxes, and convection, and the organization is intimately related to the evolution of the column water vapor field. Analysis of the budget of the spatial variance of column-integrated frozen moist static energy (MSE), following Wing and Emanuel [2014], reveals that the direct radiative feedback (including significant cloud longwave effects) is dominant in both the initial development of self-aggregation and the maintenance of an aggregated state. A low-level circulation at intermediate stages of aggregation does appear to transport MSE from drier to moister regions, but this circulation is mostly balanced by other advective effects of opposite sign and is forced by horizontal anomalies of convective heating (not radiation). Sensitivity studies with either fixed prescribed radiative cooling, fixed prescribed surface fluxes, or both do not show full self-aggregation from homogeneous initial conditions, though fixed surface fluxes do not disaggregate an initialized aggregated state. A sensitivity study in which rain evaporation is turned off shows more rapid self-aggregation, while a run with this change plus fixed radiative cooling still shows strong self-aggregation, supporting a “moisture memory” effect found in Muller and Bony [2015]. Interestingly, self-aggregation occurs even in simulations with sea surface temperatures (SSTs) of 295 K and 290 K, with direct radiative feedbacks dominating the budget of MSE variance, in contrast to results in some previous studies.

Proper evaluation of the external convective heat transfer for the thermal analysis of cool roofs

Cool materials are characterized by high solar reflectance and high thermal emittance; when applied to the external surface of a roof, they make it possible to limit the amount of solar irradiance absorbed by the roof, and to increase the rate of heat flux emitted by irradiation to the environment, especially during nighttime. However, a roof also releases heat by convection on its external surface; this mechanism is not negligible, and an incorrect evaluation of its entity might introduce significant inaccuracy in the assessment of the thermal performance of a cool roof, in terms of surface temperature and rate of heat flux transferred to the indoors. This issue is particularly relevant in numerical simulations, which are essential in the design stage, therefore it deserves adequate attention. In the present paper, a review of the most common algorithms used for the calculation of the convective heat transfer coefficient due to wind on horizontal building surfaces is presented. Then, with reference to a case study in Italy, the simulated results are compared to the outcomes of a measurement campaign. Hence, the most appropriate algorithms for the convective coefficient are identified, and the errors deriving by an incorrect selection of this coefficient are discussed.

Climatic simulations of the eastern Andes low-level jet and its dependency on convective parameterizations

The South American low level jet (SALLJ) of the Eastern Andes is investigated with Regional Climate Model version 3 (RegCM3) simulations during the 2002-2003 austral summer using two convective parameterizations (Grell and Emanuel). The simulated SALLJ is compared with the special observations of SALLJEX (SALLJ Experiment). Both the Grell and Emanuel schemes adequately simulate the low level flow over South America. However, there are some intensity differences. Due to the larger (smaller) convective activity, the Emanuel (Grell) scheme simulates more intense (weaker) low level wind than analysis in the tropics and subtropics. The objectives criteria of Sugahara (SJ) and Bonner (BJ) were used for LLJ identification. When applied to the observations, both criteria suggest a larger frequency of the SALLJ in Santa Cruz, followed by Mariscal, Trinidad and Asuncin. In Mariscal and Asuncin, the diurnal cycle indicates that SJ occurs mainly at 12 UTCs (morning), while the BJ criterion presents the SALLJ as more homogenously distributed. The concentration into two of the four-times-a-day observations does not allow conclusions about the diurnal cycle in Santa Cruz and Trinidad. The simulated wind profiles result in a lower than observed frequency of SALLJ using both the SJ and BJ criteria, with fewer events obtained with the BJ. Due to the stronger simulated winds, the Emanuel scheme produces an equal or greater relative frequency of SALLJ than the Grell scheme. However, the Grell scheme using the SJ criterion simulates the SALLJ diurnal cycle closer to the observed one. Although some discrepancies between observed and simulated mean vertical profiles of the horizontal wind are noted, there is large agreement between the composites of the vertical structure of the SALLJ, especially when the SJ criterion is used with the Grell scheme. On an intraseasonal scale, a larger southward displacement of SALLJ in February and December when compared with January has been noted. The Grell and Emanuel schemes simulated this observed oscillation in the low-level flow. However, the spatial pattern and intensity of rainfall and circulation anomalies simulated by the Grell scheme are closer to the analyses than those obtained with the Emanuel scheme.

The role of convective parameterization in the simulation of a cyclone over the South Atlantic

Numerical simulations are carried out to examine the role of the Kuo and Kain-Fritsch (KF) cumulus parameterization schemes and dry dynamics on a cyclone development, in a weak baroclinic atmosphere, over subtropical South Atlantic Ocean. The initial phase of the cyclone development is investigated with a coarse horizontal mesh (75 km) and when the cyclone reaches the mature stage two different horizontal resolutions are used (75 and 25 km). The best performance simulation for the cyclone initial phase occurs when the Kuo convective scheme is applied, and this may be attributed to a greater diabatic warming in the troposphere. On the other hand, the dry simulation is not capable of simulating the correct location and intensity of the cyclone in its initial phase. During the mature phase, a cyclone over deepening occurs in the Kuo scheme experiment associated with larger latent heat release in a deep vertical column. The presence of downdraft currents in the KF scheme, which acts to cool and dry the lower levels, is essential to stabilize the atmosphere and to reproduce the nearest observation cyclone deepening rate. The largest cyclone deepening is found in the Kuo scheme high resolution experiment. This suggests that the KF convective scheme is less sensitive to the horizontal grid resolution. It was also revealed that the diabatic processes are crucial to simulate the observed features of this marine cyclone over subtropical region.

Saturated flow boiling heat transfer and critical heat flux in small horizontal flattened tubes

This paper presents experimental results for flow boiling heat transfer coefficient and critical heat flux (CHF) in small flattened tubes. The tested flattened tubes have the same equivalent internal diameter of 2.2 mm, but different aspect height/width ratios (H/W) of 1/4, 1/2, 2 and 4. The experimental data were compared against results for circular tubes using R134a and R245fa as working fluids at a nominal saturation temperature of 31 degrees C. For mass velocities higher than 200 kg/m(2)s, the flattened and circular tubes presented similar heat transfer coefficients. Such a behavior is related to the fact that stratification effects are negligible under conditions of higher mass velocities. Heat transfer correlations from the literature, usually developed using only circular-channel experimental data, predicted the flattened tube results for mass velocities higher than 200 kg/m(2)s with mean absolute error lower than 20% using the equivalent diameter to account for the geometry effect. Similarly, the critical heat flux results were found to be independent of the tube aspect ratio when the same equivalent length was kept. Equivalent length is a new parameter which takes into account the channel heat transfer area. The CHF correlations for round tubes predicted the flattened tube data relatively well when using the equivalent diameter and length. Furthermore, a new proposed CHF correlation predicted the present flattened tube data with a mean absolute error of 5%. (C) 2012 Elsevier Ltd. All rights reserved.

Two-phase flow patterns and pressure drop inside horizontal tubes containing twisted-tape inserts

This paper presents an experimental study on two-phase flow patterns and pressure drop of R134a inside a 15.9 mm ID tube containing twisted-tape inserts. Experimental results were obtained in a horizontal test section for twisted-tape ratios of 3, 4, 9 and 14, mass velocities ranging from 75 to 250 kg/m(2) s and saturation temperatures of 5 and 15 degrees C. An unprecedented discussion on two-phase flow patterns inside tubes containing twisted-tape inserts is presented and the flow pattern effects on the frictional pressure drop are carefully discussed. Additionally, a new method to predict the frictional pressure drop during two-phase flow inside tubes containing twisted-tape inserts is proposed. (C) 2012 Elsevier Ltd. All rights reserved.

Influences of deep convective cloud systems on tropospheric trace gases and photochemistry over the tropical West Pacific

A numerical model for studying the influences of deep convective cloud systems on photochemistry was developed based on a non-hydrostatic meteorological model and chemistry from a global chemistry transport model. The transport of trace gases, the scavenging of soluble trace gases, and the influences of lightning produced nitrogen oxides (NOx=NO+NO2) on the local ozone-related photochemistry were investigated in a multi-day case study for an oceanic region located in the tropical western Pacific. Model runs considering influences of large scale flows, previously neglected in multi-day cloud resolving and single column model studies of tracer transport, yielded that the influence of the mesoscale subsidence (between clouds) on trace gas transport was considerably overestimated in these studies. The simulated vertical transport and scavenging of highly soluble tracers were found to depend on the initial profiles, reconciling contrasting results from two previous studies. Influences of the modeled uptake of trace gases by hydrometeors in the liquid and the ice phase were studied in some detail for a small number of atmospheric trace gases and novel aspects concerning the role of the retention coefficient (i.e. the fraction of a dissolved trace gas that is retained in the ice phase upon freezing) on the vertical transport of highly soluble gases were illuminated. Including lightning NOx production inside a 500 km 2-D model domain was found to be important for the NOx budget and caused small to moderate changes in the domain averaged ozone concentrations. A number of sensitivity studies yielded that the fraction of lightning associated NOx which was lost through photochemical reactions in the vicinity of the lightning source was considerable, but strongly depended on assumptions about the magnitude and the altitude of the lightning NOx source. In contrast to a suggestion from an earlier study, it was argued that the near zero upper tropospheric ozone mixing ratios which were observed close to the study region were most probably not caused by the formation of NO associated with lightning. Instead, it was argued in agreement with suggestions from other studies that the deep convective transport of ozone-poor air masses from the relatively unpolluted marine boundary layer, which have most likely been advected horizontally over relatively large distances (both before and after encountering deep convection) probably played a role. In particular, it was suggested that the ozone profiles observed during CEPEX (Central Equatorial Pacific Experiment) were strongly influenced by the deep convection and the larger scale flow which are associated with the intra-seasonal oscillation.

Radar based improvements and diagnosis of the convective capabilities of the cosmo-2 model

In the last few years the resolution of numerical weather prediction (nwp) became higher and higher with the progresses of technology and knowledge. As a consequence, a great number of initial data became fundamental for a correct initialization of the models. The potential of radar observations has long been recognized for improving the initial conditions of high-resolution nwp models, while operational application becomes more frequent. The fact that many nwp centres have recently taken into operations convection-permitting forecast models, many of which assimilate radar data, emphasizes the need for an approach to providing quality information which is needed in order to avoid that radar errors degrade the model's initial conditions and, therefore, its forecasts. Environmental risks can can be related with various causes: meteorological, seismical, hydrological/hydraulic. Flash floods have horizontal dimension of 1-20 Km and can be inserted in mesoscale gamma subscale, this scale can be modeled only with nwp model with the highest resolution as the COSMO-2 model. One of the problems of modeling extreme convective events is related with the atmospheric initial conditions, in fact the scale dimension for the assimilation of atmospheric condition in an high resolution model is about 10 Km, a value too high for a correct representation of convection initial conditions. Assimilation of radar data with his resolution of about of Km every 5 or 10 minutes can be a solution for this problem. In this contribution a pragmatic and empirical approach to deriving a radar data quality description is proposed to be used in radar data assimilation and more specifically for the latent heat nudging (lhn) scheme. Later the the nvective capabilities of the cosmo-2 model are investigated through some case studies. Finally, this work shows some preliminary experiments of coupling of a high resolution meteorological model with an Hydrological one.

Hydrothermal waves and corotating rolls in laterally heated convection in simple liquids.

The stability of a liquid layer with an undeformable interface open to the atmo- sphere, subjected to a horizontal temperature gradient, is theoretically analysed. Buoyancy and surface tension forces give rise to a basic flow for any temperature dif- ference applied on the system. Depending on the liquid depth, this basic flow is desta- bilised either by an oscillatory instability, giving rise to the so-called hydrothermal waves, or by a stationary instability leading to corotating rolls. Oscillatory perturba- tions are driven by the basic flow and therefore one must distinguish between convec- tive and absolute thresholds. The instability mechanisms as well as the di¿erent re- gimes observed in experiments are discussed. The calculations are performed for a fluid used in recent experiments, namely silicone oil of 0.65 cSt ðPr 1?4 10Þ. In partic- ular, it is shown that two branches of absolute instability exist, which may be related to the two types of hydrothermal waves observed experimentally