Abstract
A new system that resolves the stratosphere was implemented for operational medium-range weather forecasts at the Canadian Meteorological Centre. The model lid was raised from 10 to 0.1 hPa, ...parameterization schemes for nonorographic gravity wave tendencies and methane oxidation were introduced, and a new radiation scheme was implemented. Because of the higher lid height of 0.1 hPa, new measurements between 10 and 0.1 hPa were also added. This new high-top system resulted not only in dramatically improved forecasts of the stratosphere, but also in large improvements in medium-range tropospheric forecast skill. Pairs of assimilation experiments reveal that most of the stratospheric and tropospheric forecast improvement is obtained without the extra observations in the upper stratosphere. However, these observations further improve forecasts in the winter hemisphere but not in the summer hemisphere. Pairs of forecast experiments were run in which initial conditions were the same for each experiment but the forecast model differed. The large improvements in stratospheric forecast skill are found to be due to the higher lid height of the new model. The new radiation scheme helps to improve tropospheric forecasts. However, the degree of improvement seen in tropospheric forecast skill could not be entirely explained with these purely forecast experiments. It is hypothesized that the cycling of a better model and assimilation provide improved initial conditions, which result in improved forecasts.
The purpose of this paper is to quantify the influence of turbulence in collision statistics by separately studying the impacts of computational domain sizes, eddy dissipation rates (EDRs), and ...droplet sizes and eventually to develop an accurate parameterization of collision kernels. Direct numerical simulations (DNS) were performed with a relatively wide range of EDRs and Taylor microscale Reynolds numbers . EDR measures the turbulence intensity levels. DNS model studies have simulated homogeneous turbulence in a small domain in the cloud's adiabatic core. Clouds clearly have much larger scales than current DNS can simulate. For this reason, it is emphasized that obtained from current DNS is fundamentally only a measure of the computational domain size for a given EDR and cannot completely describe the physical properties of cloud turbulence. Results show that the collision statistics are independent of the domain sizes and hence of the computational for droplet sizes no bigger than 25 mu m as long as the droplet separation distance, which is on the order of the Kolmogorov scale in real clouds, is resolved. Instead, they are found to be highly correlated with EDRs and droplet sizes, and this correlation is used to formulate an improved parameterization scheme. The new scheme well represents the turbulent geometric collision kernel with a relative uncertainty of 14%. A comparison between different parameterizations is made, and the formulas proposed here are shown to improve the fit to the collision statistics.
A model evaluation approach is proposed in which weather and climate prediction models are analyzed along a Pacific Ocean cross section, from the stratocumulus regions off the coast of California, ...across the shallow convection dominated trade winds, to the deep convection regions of the ITCZ—the Global Energy and Water Cycle Experiment Cloud System Study/Working Group on Numerical Experimentation (GCSS/WGNE) Pacific Cross-Section Intercomparison (GPCI). The main goal of GPCI is to evaluate and help understand and improve the representation of tropical and subtropical cloud processes in weather and climate prediction models. In this paper, a detailed analysis of cloud regime transitions along the cross section from the subtropics to the tropics for the season June–July–August of 1998 is presented. This GPCI study confirms many of the typical weather and climate prediction model problems in the representation of clouds: underestimation of clouds in the stratocumulus regime by most models with the corresponding consequences in terms of shortwave radiation biases; overestimation of clouds by the 40-yr ECMWF Re-Analysis (ERA-40) in the deep tropics (in particular) with the corresponding impact in the outgoing longwave radiation; large spread between the different models in terms of cloud cover, liquid water path and shortwave radiation; significant differences between the models in terms of vertical cross sections of cloud properties (in particular), vertical velocity, and relative humidity. An alternative analysis of cloud cover mean statistics is proposed where sharp gradients in cloud cover along the GPCI transect are taken into account. This analysis shows that the negative cloud bias of some models and ERA-40 in the stratocumulus regions as compared to the first International Satellite Cloud Climatology Project (ISCCP) is associated not only with lower values of cloud cover in these regimes, but also with a stratocumulus-to-cumulus transition that occurs too early along the trade wind Lagrangian trajectory. Histograms of cloud cover along the cross section differ significantly between models. Some models exhibit a quasi-bimodal structure with cloud cover being either very large (close to 100%) or very small, while other models show a more continuous transition. The ISCCP observations suggest that reality is in-between these two extreme examples. These different patterns reflect the diverse nature of the cloud, boundary layer, and convection parameterizations in the participating weather and climate prediction models.
Atmospheric physics is represented in numerical models by parameterizations that use resolved‐scale information to estimate the effects of physical processes on the atmospheric state. Over time, our ...understanding of these processes improves, new techniques are introduced to represent physics in a numerical model, and increased resolution changes the relative importance of different parameterizations within the system. As a result, the physical parameterization packages of numerical weather prediction (NWP) models undergo regular updates as older schemes are replaced with newer ones that offer an improved, and often more complex, depiction of relevant physical processes. Such changes are typically combined with a rebalancing of the physics suite because of strong interactions between parameterization schemes and the presence of compensating errors in the system. In this study, a major update to the package of physical parameterizations used in Canadian operational NWP is introduced. The primary goals of this effort were to improve the global energy budget and to facilitate an increase in the vertical resolution of operational configurations. Both of these objectives were achieved, along with a significant improvement in guidance quality for global and regional prediction systems.
Key Points
An updated suite of parameterizations for atmospheric physics is developed for Canadian NWP
Representation of the global energy budget is improved in the model climatology
Significant improvements in guidance quality are seen in upper‐air, surface, and precipitation scores
Hybridization is increasingly recognized as an integral part of the dynamics of species range expansion and contraction. Thus, it is important to understand the reproductive barriers between ...co-occurring species. Extending previous studies that argued that the rare Eucalyptus risdonii was expanding into the range of the surrounding E. amygdalina by both seed and pollen dispersal, we here investigate the long-term fitness of both species and their hybrids and whether expansion is continuing.
We assessed the survival of phenotypes representing a continuum between the two pure species in a natural hybrid swarm after 29 years, along with seedling recruitment. The performance of pure species as well as of artificial and natural hybrids was also assessed over 28 years in a common garden trial.
In the hybrid zone, E. amygdalina adults showed greater mortality than E. risdonii, and the current seedling cohort is still dominated by E. risdonii phenotypes. Morphologically intermediate individuals appeared to be the least fit. Similar results were observed after growing artificial first-generation and natural hybrids alongside pure species families in a common garden trial. Here, the survival, reproduction, health and growth of the intermediate hybrids were significantly less than those of either pure species, consistent with hybrid inferiority, although this did not manifest until later reproductive ages. Among the variable progeny of natural intermediate hybrids, the most E. risdonii-like phenotypes were the most fit.
This study contributes to the increasing number of reports of hybrid inferiority in Eucalyptus, suggesting that post-zygotic barriers contribute to the maintenance of species integrity even between closely related species. However, with fitness rapidly recovered following backcrossing, it is argued that hybridization can still be an important evolutionary process, in the present case appearing to contribute to the range expansion of the rare E. risdonii in response to climate change.
The goal of this work is to answer the question of whether nonuniformity in the spatial distribution of sizes and positions of cloud droplets and/or variable vertical velocity in a turbulent medium ...can contribute to the broadening of the droplet size distribution. A numerical approach to simulate the growth and trajectory of several tens of thousands of cloud droplets in a turbulent environment whose properties vary from droplet to droplet is used. The finite inertia of particles in a turbulent fluid causes particles to diverge from regions of high vorticity and to converge preferentially in regions of low vorticity, thus creating strong deviations in particle concentration. As a first step, the inertia effect was examined in the context of nongrowing, sedimenting, or nonsedimenting droplets. It was found that statistically significant preferential concentration is possible in conditions typical of cloud droplets in cumulus clouds. In the absence of sedimentation, preferential concentration increases as a function of the Stokes number St. Allowing the droplets to sediment decreases preferential concentration to a degree that increases with the velocity ratio SU. A series of experiments including condensational growth of droplets was then performed. It was found that while the increasing preferential concentration of droplets, as a result of increasing eddy dissipation rate, does result in increases in the instantaneous dispersion of the supersaturation perturbation distribution, the width of the size distribution of droplets, which is a function of the dispersion in the time integral of the supersaturation perturbations, decreases. This result is a consequence of the decrease in decorrelation time of the supersaturation perturbations as the turbulence intensity increases. Comparison of the results herein with the observations made in quasi-adiabatic cloud cores leads one to the conclusion that the microscopic approach, even under the most favorable condition of no turbulence, produces too little broadening to explain the observations.
Abstract
Direct numerical simulations of an evolving turbulent flow field have been performed to explore how turbulence affects the motion and collisions of cloud droplets. Large numbers of droplets ...are tracked through the flow field and their positions, velocities, and collision rates have been found to depend on the eddy dissipation rate of turbulent kinetic energy. The radial distribution function, which is a measure of the preferential concentration of droplets, increases with eddy dissipation rate. When droplets are clustered there is an increased probability of finding two droplets closely separated; thus, there is an increase in the collision kernel. For the flow fields explored in this study, the clustering effect accounts for an increase in the collision kernel of 8%–42%, as compared to the gravitational collision kernel. The spherical collision kernel is also a function of the radial relative velocities among droplets and these velocities increase from 1.008 to 1.488 times the corresponding gravitational value. For an eddy dissipation rate of about 100 cm2 s−3, the turbulent collision kernel is 1.06 times the magnitude of the gravitational value, while for an eddy dissipation rate of 1500 cm2 s−3, this increases to 2.08 times. Therefore, these results demonstrate that turbulence could play an important role in the broadening and evolution of the droplet size distribution and the onset of precipitation.
There is compelling evidence that members of the caspase (interleukin-1beta converting enzyme/CED-3) family of cysteine proteases and the cytotoxic lymphocyte-derived serine protease granzyme B play ...essential roles in mammalian apoptosis. Here we use a novel method employing a positional scanning substrate combinatorial library to rigorously define their individual specificities. The results divide these proteases into three distinct groups and suggest that several have redundant functions. The specificity of caspases 2, 3, and 7 and Caenorhabditis elegans CED-3 (DEXD) suggests that all of these enzymes function to incapacitate essential homeostatic pathways during the effector phase of apoptosis. In contrast, the optimal sequence for caspases 6, 8, and 9 and granzyme B ((I/L/V)EXD) resembles activation sites in effector caspase proenzymes, consistent with a role for these enzymes as upstream components in a proteolytic cascade that amplifies the death signal.
The atmospheric boundary layer mixing height (MH) is an important bulk parameter in air quality (AQ) modelling. Formulating this parameter under statically stable conditions, such as in the Arctic, ...has historically been difficult. In an effort to improve AQ modelling capacity in North America, MH is studied in two geographically distinct areas: the Arctic (Barrow, Alaska) and the southern Great Plains (Lamont, Oklahoma). Observational data from the Atmospheric Radiation Measurement program, Climate Research Facility and numerical weather forecasting data from Environment Canada's Regional Global Environmental Multiscale (GEM15) model have been used in order to examine the suitability of available parameterizations for MH under statically stable conditions and also to compare the level of agreement between observed and modelled MH. The analysis period is 1 October 2011 to 1 October 2012. The observations alone suggest that profile methods are preferred over surface methods in defining MH under statically stable conditions. Surface methods exhibit poorer comparison statistics with observations than profile methods. In addition, the fitted constants for surface methods are site-dependent, precluding their applicability for modelling under general conditions. The comparison of observations and GEM15 MH suggests that although the agreement is acceptable in Lamont, the default model surface method contributes to a consistent overprediction of MH in Barrow in all seasons. An alternative profile method for MH is suggested based on the bulk Richardson number. This method is shown to reduce the model bias in Barrow by a factor of two without affecting model performance in Lamont.