Inter-comparison of three-dimensional models of volcanic plumes

We performed an inter-comparison study of three-dimensional models of volcanic plumes. A set of common volcanological input parameters and meteorological conditions were provided for two kinds of eruptions, representing a weak and a strong eruption column. From the different models, we compared the maximum plume height, neutral buoyancy level (where plume density equals that of the atmosphere), and level of maximum radial spreading of the umbrella cloud. We also compared the vertical profiles of eruption column properties, integrated across cross-sections of the plume (integral variables). Although the models use different numerical procedures and treatments of subgrid turbulence and particle dynamics, the inter-comparison shows qualitatively consistent results. In the weak plume case (mass eruption rate 1.5 × 106 kg s− 1), the vertical profiles of plume properties (e.g., vertical velocity, temperature) are similar among models, especially in the buoyant plume region. Variability among the simulated maximum heights is ~ 20%, whereas neutral buoyancy level and level of maximum radial spreading vary by ~ 10%. Time-averaging of the three-dimensional (3D) flow fields indicates an effective entrainment coefficient around 0.1 in the buoyant plume region, with much lower values in the jet region, which is consistent with findings of small-scale laboratory experiments. On the other hand, the strong plume case (mass eruption rate 1.5 × 109 kg s− 1) shows greater variability in the vertical plume profiles predicted by the different models. Our analysis suggests that the unstable flow dynamics in the strong plume enhances differences in the formulation and numerical solution of the models. This is especially evident in the overshooting top of the plume, which extends a significant portion (~ 1/8) of the maximum plume height. Nonetheless, overall variability in the spreading level and neutral buoyancy level is ~ 20%, whereas that of maximum height is ~ 10%. This inter-comparison study has highlighted the different capabilities of 3D volcanic plume models, and identified key features of weak and strong plumes, including the roles of jet stability, entrainment efficiency, and particle non-equilibrium, which deserve future investigation in field, laboratory, and numerical studies.YJS was partially supported by the ERI Cooperative Research Program and KAKENHI (25750142). The computations of SK-3D were carried out in part on the Earth Simulator at the JAMSTEC and also on the Primergy RX200S6 at the Research Computer System, Kyushu University. AC was partially supported by a grant of the International Research Promotion Office Earthquake Research Institute, the University of Tokyo. AC, TEO and MC were partially supported by the EU-funded project MEDiterranean Supersite Volcanoes (MEDSUV; grant no. 308665). MC acknowledges CINECA award N. HP10BKFD9F (2013) for high performance computing resources and support. AVE acknowledges NSF Postdoctoral Fellowship EAR1250029, a U.S. Geological Survey Mendenhall fellowship, and grant GID 61233 from NASA Ames Supercomputing Center

Who died as a result of the tsunami? – Risk factors of mortality among internally displaced persons in Sri Lanka: a retrospective cohort analysis

BACKGROUND: Describing adverse health effects and identifying vulnerable populations during and after a disaster are important aspects of any disaster relief operation. This study aimed to describe the mortality and related risk factors which affected the displaced population over a period of two and a half months after the 2004 Indian Ocean tsunami in an eastern coastal district of Sri Lanka. METHODS: A cross-sectional household survey was conducted in 13 evacuation camps for internally displaced persons (IDP). Information on all pre-tsunami family members was collected from householders, and all deaths which occurred during the recall period (77 to 80 days starting from the day of the tsunami) were recorded. The distribution of mortality and associated risk factors were analysed. Logistic regression modelling using the generalized estimating equations method was applied in multivariate analysis. RESULTS: Overall mortality rate out of 3,533 individuals from 859 households was 12.9% (446 deaths and 11 missing persons). The majority of the deaths occurred during and immediately after the disaster. A higher mortality was observed among females (17.5% vs. 8.2% for males, p < 0.001), children and the elderly (31.8%, 23.7% and 15.3% for children aged less than 5 years, children aged 5 to 9 years and adults over 50 years, respectively, compared with 7.4% for adults aged 20 to 29 years, p < 0.001). Other risk factors, such as being indoors at the time of the tsunami (13.8% vs. 5.9% outdoors, p < 0.001), the house destruction level (4.6%, 5.5% and 14.2% in increasing order of destruction, p < 0.001) and fishing as an occupation (15.4% vs. 11.2% for other occupations, p < 0.001) were also significantly associated with increased mortality. These correlations remained significant after adjusting for the confounding effects by multivariate analysis. CONCLUSION: A significantly high mortality was observed in women and children among the displaced population in the eastern coastal district of Sri Lanka who were examined by us. Reconstruction activities should take into consideration these changes in population structure

Computer simulation of breast reduction surgery

Background: Plastic surgery of the breast, particularly breast reduction, is considered difficult. It can become a challenge for a less experienced surgeon to understand exactly what to do when facing a particular type of breast and how to avoid unsatisfactory results. Methods: The goal of this study was to create a computer model of the breast that provides a basis for the simulation of breast surgery, particularly breast reduction. The reconstruction of elastic parameters is based on observations of the breast with the patient in different positions. Results: It is shown that several measurements with the patient in different positions allow one to choose the parameters of the model and determine the elastic coefficients of the breast and the skin. The geometry of the breast before and after surgery is simulated. A qualitative study of the incision parameters’ influence on the final geometry of the breast is presented. Conclusion: The developed methodology and software allow one to estimate the form of the breast after the surgery by knowing its form before surgery and taking into consideration the parameters of incision applied by the surgeon at the time of surgery. The described approach can be used for the qualitative and quantitative study of breast reduction surgery with a satisfactory result. Level of Evidence: V (This journal requires that authors assign a level of evidence to each article. For a full description of these Evidence-Based Medicine ratings, please refer to the Table of Contents or the online Instructions to Authors http://www.springer.com/00266.

Comparison of Gene Expression Profiles in Chromate Transformed BEAS-2B Cells

Hexavalent chromium [Cr(VI)] is a potent human carcinogen. Occupational exposure has been associated with increased risk of respiratory cancer. Multiple mechanisms have been shown to contribute to Cr(VI) induced carcinogenesis, including DNA damage, genomic instability, and epigenetic modulation, however, the molecular mechanism and downstream genes mediating chromium's carcinogenicity remain to be elucidated.We established chromate transformed cell lines by chronic exposure of normal human bronchial epithelial BEAS-2B cells to low doses of Cr(VI) followed by anchorage-independent growth. These transformed cell lines not only exhibited consistent morphological changes but also acquired altered and distinct gene expression patterns compared with normal BEAS-2B cells and control cell lines (untreated) that arose spontaneously in soft agar. Interestingly, the gene expression profiles of six Cr(VI) transformed cell lines were remarkably similar to each other yet differed significantly from that of either control cell lines or normal BEAS-2B cells. A total of 409 differentially expressed genes were identified in Cr(VI) transformed cells compared to control cells. Genes related to cell-to-cell junction were upregulated in all Cr(VI) transformed cells, while genes associated with the interaction between cells and their extracellular matrices were down-regulated. Additionally, expression of genes involved in cell proliferation and apoptosis were also changed.This study is the first to report gene expression profiling of Cr(VI) transformed cells. The gene expression changes across individual chromate exposed clones were remarkably similar to each other but differed significantly from the gene expression found in anchorage-independent clones that arose spontaneously. Our analysis identified many novel gene expression changes that may contribute to chromate induced cell transformation, and collectively this type of information will provide a better understanding of the mechanism underlying chromate carcinogenicity

Precision measurement of the top quark mass from dilepton events at CDF II

We report a measurement of the top quark mass, M_t, in the dilepton decay channel of $t\bar{t}\to b\ell'^{+}\nu_{\ell'}\bar{b}\ell^{-}\bar{\nu}_{\ell}$ using an integrated luminosity of 1.0 fb^{-1} of p\bar{p} collisions collected with the CDF II detector. We apply a method that convolutes a leading-order matrix element with detector resolution functions to form event-by-event likelihoods; we have enhanced the leading-order description to describe the effects of initial-state radiation. The joint likelihood is the product of the likelihoods from 78 candidate events in this sample, which yields a measurement of M_{t} = 164.5 \pm 3.9(\textrm{stat.}) \pm 3.9(\textrm{syst.}) \mathrm{GeV}/c^2, the most precise measurement of M_t in the dilepton channel.Comment: 7 pages, 2 figures, version includes changes made prior to publication by journa

Measurement of the Ratios of Branching Fractions B(Bs -> Ds pi pi pi) / B(Bd -> Dd pi pi pi) and B(Bs -> Ds pi) / B(Bd -> Dd pi)

Using 355 pb^-1 of data collected by the CDF II detector in \ppbar collisions at sqrt{s} = 1.96 TeV at the Fermilab Tevatron, we study the fully reconstructed hadronic decays B -> D pi and B -> D pi pi pi. We present the first measurement of the ratio of branching fractions B(Bs -> Ds pi pi pi) / B(Bd -> Dd pi pi pi) = 1.05 pm 0.10 (stat) pm 0.22 (syst). We also update our measurement of B(Bs -> Ds pi) / B(Bd -> Dd pi) to 1.13 pm 0.08 (stat) pm 0.23 (syst) improving the statistical uncertainty by more than a factor of two. We find B(Bs -> Ds pi) = [3.8 pm 0.3 (stat) pm 1.3 (syst)] \times 10^{-3} and B(Bs -> Ds pi pi pi) = [8.4 pm 0.8 (stat) pm 3.2 (syst)] \times 10^{-3}.Comment: 7 pages, 2 figure