Assessing the vulnerability of blanket peat to climate change using an ensemble of statistical bioclimatic envelope models

We assessed the vulnerability of blanket peat to climate change in Great Britain using an ensemble of 8 bioclimatic envelope models. We used 4 published models that ranged from simple threshold models, based on total annual precipitation, to Generalised Linear Models (GLMs, based on mean annual temperature). In addition, 4 new models were developed which included measures of water deficit as threshold, classification tree, GLM and generalised additive models (GAM). Models that included measures of both hydrological conditions and maximum temperature provided a better fit to the mapped peat area than models based on hydrological variables alone. Under UKCIP02 projections for high (A1F1) and low (B1) greenhouse gas emission scenarios, 7 out of the 8 models showed a decline in the bioclimatic space associated with blanket peat. Eastern regions (Northumbria, North York Moors, Orkney) were shown to be more vulnerable than higher-altitude, western areas (Highlands, Western Isles and Argyle, Bute and The Trossachs). These results suggest a long-term decline in the distribution of actively growing blanket peat, especially under the high emissions scenario, although it is emphasised that existing peatlands may well persist for decades under a changing climate. Observational data from long-term monitoring and manipulation experiments in combination with process-based models are required to explore the nature and magnitude of climate change impacts on these vulnerable areas more fully.

Initiation of deep convection at marginal instability in an ensemble of mesoscale models: a case-study from COPS

The present study investigates the initiation of precipitating deep convection in an ensemble of convection-resolving mesoscale models. Results of eight different model runs from five non-hydrostatic models are compared for a case of the Convective and Orographically-induced Precipitation Study (COPS). An isolated convective cell initiated east of the Black Forest crest in southwest Germany, although convective available potential energy was only moderate and convective inhibition was high. Measurements revealed that, due to the absence of synoptic forcing, convection was initiated by local processes related to the orography. In particular, the lifting by low-level convergence in the planetary boundary layer is assumed to be the dominant process on that day. The models used different configurations as well as different initial and boundary conditions. By comparing the different model performance with each other and with measurements, the processes which need to be well represented to initiate convection at the right place and time are discussed. Besides an accurate specification of the thermodynamic and kinematic fields, the results highlight the role of boundary-layer convergence features for quantitative precipitation forecasts in mountainous terrain.

Comments on “An ensemble cumulus convection parameterization with explicit cloud treatment”

Wagner and Graf (2010) derive a population evolution equation for an ensemble of convective plumes, an analogue with the Lotka–Volterra equation, from the energy equations for convective plumes provided by Arakawa and Schubert (1974). Although their proposal is interesting, as the present note shows, there are some problems with their derivation.

Ensemble prediction for nowcasting with a convection-permitting model—I: description of the system and the impact of radar-derived surface precipitation rates

A key strategy to improve the skill of quantitative predictions of precipitation, as well as hazardous weather such as severe thunderstorms and flash floods is to exploit the use of observations of convective activity (e.g. from radar). In this paper, a convection-permitting ensemble prediction system (EPS) aimed at addressing the problems of forecasting localized weather events with relatively short predictability time scale and based on a 1.5 km grid-length version of the Met Office Unified Model is presented. Particular attention is given to the impact of using predicted observations of radar-derived precipitation intensity in the ensemble transform Kalman filter (ETKF) used within the EPS. Our initial results based on the use of a 24-member ensemble of forecasts for two summer case studies show that the convective-scale EPS produces fairly reliable forecasts of temperature, horizontal winds and relative humidity at 1 h lead time, as evident from the inspection of rank histograms. On the other hand, the rank histograms seem also to show that the EPS generates too much spread for forecasts of (i) surface pressure and (ii) surface precipitation intensity. These may indicate that for (i) the value of surface pressure observation error standard deviation used to generate surface pressure rank histograms is too large and for (ii) may be the result of non-Gaussian precipitation observation errors. However, further investigations are needed to better understand these findings. Finally, the inclusion of predicted observations of precipitation from radar in the 24-member EPS considered in this paper does not seem to improve the 1-h lead time forecast skill.

Ensemble prediction for nowcasting with a convection-permitting model - II: forecast error statistics

A 24-member ensemble of 1-h high-resolution forecasts over the Southern United Kingdom is used to study short-range forecast error statistics. The initial conditions are found from perturbations from an ensemble transform Kalman filter. Forecasts from this system are assumed to lie within the bounds of forecast error of an operational forecast system. Although noisy, this system is capable of producing physically reasonable statistics which are analysed and compared to statistics implied from a variational assimilation system. The variances for temperature errors for instance show structures that reflect convective activity. Some variables, notably potential temperature and specific humidity perturbations, have autocorrelation functions that deviate from 3-D isotropy at the convective-scale (horizontal scales less than 10 km). Other variables, notably the velocity potential for horizontal divergence perturbations, maintain 3-D isotropy at all scales. Geostrophic and hydrostatic balances are studied by examining correlations between terms in the divergence and vertical momentum equations respectively. Both balances are found to decay as the horizontal scale decreases. It is estimated that geostrophic balance becomes less important at scales smaller than 75 km, and hydrostatic balance becomes less important at scales smaller than 35 km, although more work is required to validate these findings. The implications of these results for high-resolution data assimilation are discussed.

A method for merging flow-dependent forecast error statistics from an ensemble with static statistics for use in high resolution variational data assimilation

The background error covariance matrix, B, is often used in variational data assimilation for numerical weather prediction as a static and hence poor approximation to the fully dynamic forecast error covariance matrix, Pf. In this paper the concept of an Ensemble Reduced Rank Kalman Filter (EnRRKF) is outlined. In the EnRRKF the forecast error statistics in a subspace defined by an ensemble of states forecast by the dynamic model are found. These statistics are merged in a formal way with the static statistics, which apply in the remainder of the space. The combined statistics may then be used in a variational data assimilation setting. It is hoped that the nonlinear error growth of small-scale weather systems will be accurately captured by the EnRRKF, to produce accurate analyses and ultimately improved forecasts of extreme events.

Metal cluster expansion by addition of low valent group 14 reagents: III. Reaction of bis[bis(trimethylsilyl)methyl]tin with M3(CO)12 or its derivatives (M = Fe, Ru, or Os): the crystal and molecular structures of M3(μ-SnR2)2(CO)10(R = CH(SiMe3)2; M = Ru or Os)

The stannylene [SnR2] (R = CH(SiMe3)2) reacts in different ways with the three dodecacarbonyls of the iron triad: [Fe3(CO)12] gives [Fe2(CO)8(μ-SnR2)], [Ru3(CO)12] gives the planar pentametallic cluster [Ru3(CO)10(μ-SnR2)2], for which a full structural analysis is reported, while [Os3(CO)12] fails to react. Different products are also obtained from three nitrile derivatives: [Fe3-(CO)11(MeCN)] gives [Fe2(CO)6(μ-SnR2)2], which has a structure significantly different from that of known Fe2Sn2 clusters, [Ru3(CO)10(MeCN)2] gives the pentametallic cluster described above, while [Os3(CO)10(MeCN)2] gives the isostructural osmium analogue, which shows the unusual feature of a CO group bridging two osmium atoms.

Metal cluster expansion by addition of low-valent group 4B reagents: crystal and molecular structure of [Os3SnH2(CO)10{CH(SiMe3)2}2]; a compound with hydrogen bridging osmium and tin

Cluster expansion of [Os3H2(CO)10] with [SnR2][R = CH(SiMe3)2] take place in high yield to give [Os3SnH2(CO)10R2], the first closed triosmium–main-group metal cluster to be structurally characterized; a novel feature is the presence of a hydrogen atom bridging the tin atom and one of the osmium atoms.

Targeted observations of a polar low in the Norwegian sea

A developing polar low is targeted with dropsonde observations to improve the forecast of its landfall. Accurately forecasting a polar low's strength and location remains a challenge; polar lows form over the ocean in poorly observed regions, therefore initial condition errors may contribute significantly to forecast error. The targeted polar low formed in the Norwegian Sea on 3 March 2008, during the Norwegian IPY-THORPEX field campaign. Two flights, six hours apart, released dense networks of dropsondes into a sensitive region covering the polar low and Arctic front to its west. The impact of the targeted observations is assessed using the limited-area Met Office Unified Model and three-dimensional variational (3D-Var) data assimilation scheme. Forecasts were verified using ECMWF analysis data, which show good agreement with both dropsonde data from a flight through the mature polar low, and 10 m QuikSCAT winds. The impact of the targeted data moved southwards with the polar low as it developed and then hit the Norwegian coast after 24 hours. The results show that the forecast of the polar low is sensitive to the initial conditions; targeted observations from the first flight did not improve the forecast, but those from the second flight clearly improved the forecast polar low position and intensity. However, caution should be applied to attributing the forecast improvement to the assimilation of the targeted observations from a single case-study, especially in this case as the forecast improvement is moderate relative to the spread from an operational ensemble forecast

Sensory and instrumental analyses of volatiles generated during the extrusion cooking of oat flours

Three batches of oats were extruded under four combinations of process temperature (150 or 180 °C) and process moisture (14.5 and 18%). Two of the extrudates were evaluated by a sensory panel, and three were analyzed by GC-MS. Maillard reaction products, such as pyrazines, pyrroles, furans, and sulfur-containing compounds, were found in the most severely processed extrudates (high-temperature, low-moisture). These extrudates were also described by the assessors as having toasted cereal attributes. Lipid degradation products, such as alkanals, 2-alkenals, and 2,4-alkadienals, were found at much higher levels in the extrudates of the oat flour that had been debranned. It contained lower protein and fiber levels than the others and showed increased lipase activity. Extrudates from these samples also had significantly lower levels of Maillard reaction products that correlated, in the sensory analysis, with terms such as stale oil and oatmeal. Linoleic acid was added to a fourth oat flour to simulate the result of increased lipase activity, and GC-MS analysis showed both an increase in lipid degradation products and a decrease in Maillard reaction products.

Performance of palm- based fat blends with a low saturated fat content in puff pastry

Four fat blends based on palm fractions in combination with high oleic sunflower oil (HOSF) with a relatively low saturated fatty acid content (29.2±0.85%, i.e. less than 50% of that of butter) were prepared. The saturated fat was located in different triacylglycerols (TAG) structures in each blend. Principal saturated TAG were derived from palm stearin (POs, containing tripalmitoyl glycerol - PPP), palm mid fraction (PMF, containing 1,3-dipalmitoyl-2-oleoyl glycerol - POP) and interesterified PMF (inPMF, containing PPP, POP and rac-1,2-dipalmitoyl-3-oleoyl glycerol - PPO). Thus, in blend 1, composed of POs and HOSF, the saturates resided principally in PPP. In blend 2, composed of POs, PMF and HOSF, the principal saturate-containing TAG were PPP and POP. Blend 3, composed of inPMF and HOSF, was similar to blend 2 except that the disaturated TAG comprised a 2:1 mixture of PPO:POP. Finally, blend 4, a mixture of PMF and HOSF, had saturates present mainly as POP. The physical properties and the functionality of blends, as shortenings for puff pastry laminated in a warm bakery environment (20-30°C), were compared with each other, and with butter. Puff pastry prepared with blend 1 (POs:HOSF 29:71) and blend 4 (PMF:HOSF 41:59), was very hard; blend 2 (POs:PMF:HOSF 13:19:68) was most similar to butter in the compressibility of the baked product and it performed well in an independent baking trial; blend 3 (inPMF:HOSF 40:60) gave a product that required a higher force for compression than butter.

Trends in Austral jet position in ensembles of high- and low-top CMIP5 models

Trends in the position of the DJF Austral jet have been analysed for multi-model ensemble simulations of a subset of high- and low-top models for the periods 1960-2000, 2000-2050, and 2050-2098 under the CMIP5 historical, RCP4.5, and RCP8.5 scenarios. Comparison with ERA-Interim, CFSR and the NCEP/NCAR reanalysis shows that the DJF and annual mean jet positions in CMIP5 models are equatorward of reanalyses for the 1979-2006 mean. Under the RCP8.5 scenario, the mean jet position in the high-top models moves 3 degrees poleward of its 1860-1900 position by 2098, compared to just over 2 degrees for the low-top models. Changes in jet position are linked to changes in the meridional temperature gradient. Compared to low-top models, the high-top models predict greater warming in the tropical upper troposphere due to increased greenhouse gases for all periods considered: up to 0.28 K/decade more in the period 2050-2098 under the RCP8.5 scenario. Larger polar lower-stratospheric cooling is seen in high-top models: -1.64 K/decade compared to -1.40 K/decade in the period 1960-2000, mainly in response to ozone depletion, and -0.41 K/decade compared to -0.12 K/decade in the period 2050-2098, mainly in response to increases in greenhouse gases. Analysis suggests that there may be a linear relationship between the trend in jet position and meridional temperature gradient, even under strong forcing. There were no clear indications of an approach to a geometric limit on the absolute magnitude of the poleward shift by 2100.

Investigating global tropical cyclone activity with a hierarchy of AGCMs: the role of model resolution

The ability to run General Circulation Models (GCMs) at ever-higher horizontal resolutions has meant that tropical cyclone simulations are increasingly credible. A hierarchy of atmosphere-only GCMs, based on the Hadley Centre Global Environmental Model (HadGEM1), with horizontal resolution increasing from approximately 270km to 60km (at 50N), is used to systematically investigate the impact of spatial resolution on the simulation of global tropical cyclone activity, independent of model formulation. Tropical cyclones are extracted from ensemble simulations and reanalyses of comparable resolutions using a feature-tracking algorithm. Resolution is critical for simulating storm intensity and convergence to observed storm intensities is not achieved with the model hierarchy. Resolution is less critical for simulating the annual number of tropical cyclones and their geographical distribution, which are well captured at resolutions of 135km or higher, particularly for Northern Hemisphere basins. Simulating the interannual variability of storm occurrence requires resolutions of 100km or higher; however, the level of skill is basin dependent. Higher resolution GCMs are increasingly able to capture the interannual variability of the large-scale environmental conditions that contribute to tropical cyclogenesis. Different environmental factors contribute to the interannual variability of tropical cyclones in the different basins: in the North Atlantic basin the vertical wind shear, potential intensity and low-level absolute vorticity are dominant, while in the North Pacific basins mid-level relative humidity and low-level absolute vorticity are dominant. Model resolution is crucial for a realistic simulation of tropical cyclone behaviour, and high-resolution GCMs are found to be valuable tools for investigating the global location and frequency of tropical cyclones.

Piezotolerant small-colony variants with increased thermotolerance, antibiotic susceptibility, and low invasiveness in a clonal Staphylococcus aureus population

Following a pressure treatment of a clonal Staphylococcus aureus culture with 400 MPa for 30 min, piezotolerant variants were isolated. Among 21 randomly selected survivors, 9 were piezotolerant and all formed small colonies on several agar media. The majority of the isolates showed increased thermotolerance, impaired growth, and reduced antibiotic resistance compared to the wild type. However, several nonpiezotolerant isolates also demonstrated impaired growth and the small-colony phenotype. In agglutination tests for the detection of protein A and fibrinogen, the piezotolerant variants showed weaker agglutination reactions than the wild type and the other isolates. All variants also showed defective production of the typical S. aureus golden color, a characteristic which has previously been linked with virulence. They were also less able to invade intestinal epithelial cells than the wild type. These S. aureus variants showed phenotypic similarities to previously isolated Listeria monocytogenes piezotolerant mutants that contained mutations in ctsR. Because of these similarities, possible alterations in the ctsR hypermutable regions of the S. aureus variants were investigated through amplified fragment length polymorphism analysis. No mutations were identified, and subsequently we sequenced the ctsR and hrcA genes of three representative variants, finding no mutations. This work demonstrates that S. aureus probably possesses a strategy resulting in an abundance of multiple-stressresistant variants within clonal populations. This strategy, however, seems to involve genes and regulatory mechanisms different from those previously reported for L. monocytogenes. We are in the process of identifying these mechanisms.