SPH calculations of asteroid disruptions: The role of pressure dependent failure models

We present recent improvements of the modeling of the disruption of strength dominated bodies using the Smooth Particle Hydrodynamics (SPH) technique. The improvements include an updated strength model and a friction model, which are successfully tested by a comparison with laboratory experiments. In the modeling of catastrophic disruptions of asteroids, a comparison between old and new strength models shows no significant deviation in the case of targets which are initially non-porous, fully intact and have a homogeneous structure (such as the targets used in the study by Benz and Asphaug, 1999). However, for many cases (e.g. initially partly or fully damaged targets and rubble-pile structures) we find that it is crucial that friction is taken into account and the material has a pressure dependent shear strength. Our investigations of the catastrophic disruption threshold (27, as a function of target properties and target sizes up to a few 100 km show that a fully damaged target modeled without friction has a Q(D)*:, which is significantly (5-10 times) smaller than in the case where friction is included. When the effect of the energy dissipation due to compaction (pore crushing) is taken into account as well, the targets become even stronger (Q(D)*; is increased by a factor of 2-3). On the other hand, cohesion is found to have an negligible effect at large scales and is only important at scales less than or similar to 1 km. Our results show the relative effects of strength, friction and porosity on the outcome of collisions among small (less than or similar to 1000 km) bodies. These results will be used in a future study to improve existing scaling laws for the outcome of collisions (e.g. Leinhardt and Stewart, 2012). (C) 2014 Elsevier Ltd. All rights reserved.

Bayesian models for the analysis of genetic structure when populations are correlated

Motivation: Population allele frequencies are correlated when populations have a shared history or when they exchange genes. Unfortunately, most models for allele frequency and inference about population structure ignore this correlation. Recent analytical results show that among populations, correlations can be very high, which could affect estimates of population genetic structure. In this study, we propose a mixture beta model to characterize the allele frequency distribution among populations. This formulation incorporates the correlation among populations as well as extending the model to data with different clusters of populations. Results: Using simulated data, we show that in general, the mixture model provides a good approximation of the among-population allele frequency distribution and a good estimate of correlation among populations. Results from fitting the mixture model to a dataset of genotypes at 377 autosomal microsatellite loci from human populations indicate high correlation among populations, which may not be appropriate to neglect. Traditional measures of population structure tend to over-estimate the amount of genetic differentiation when correlation is neglected. Inference is performed in a Bayesian framework.

Efficiency Bounds for Estimating Linear Functionals of Nonparametric Regression Models with Endogenous Regressors

Consider a nonparametric regression model Y=mu*(X) + e, where the explanatory variables X are endogenous and e satisfies the conditional moment restriction E[e|W]=0 w.p.1 for instrumental variables W. It is well known that in these models the structural parameter mu* is 'ill-posed' in the sense that the function mapping the data to mu* is not continuous. In this paper, we derive the efficiency bounds for estimating linear functionals E[p(X)mu*(X)] and int_{supp(X)}p(x)mu*(x)dx, where p is a known weight function and supp(X) the support of X, without assuming mu* to be well-posed or even identified.

Cross-Country Evidence on Output Growth Volatility: Nonstationary Variance and GARCH Models

This paper revisits the issue of conditional volatility in real GDP growth rates for Canada, Japan, the United Kingdom, and the United States. Previous studies find high persistence in the volatility. This paper shows that this finding largely reflects a nonstationary variance. Output growth in the four countries became noticeably less volatile over the past few decades. In this paper, we employ the modified ICSS algorithm to detect structural change in the unconditional variance of output growth. One structural break exists in each of the four countries. We then use generalized autoregressive conditional heteroskedasticity (GARCH) specifications modeling output growth and its volatility with and without the break in volatility. The evidence shows that the time-varying variance falls sharply in Canada, Japan, and the U.K. and disappears in the U.S., excess kurtosis vanishes in Canada, Japan, and the U.S. and drops substantially in the U.K., once we incorporate the break in the variance equation of output for the four countries. That is, the integrated GARCH (IGARCH) effect proves spurious and the GARCH model demonstrates misspecification, if researchers neglect a nonstationary unconditional variance.

Interrelationships between the structural and functional aspects of social relationships, comorbidity and health related quality of life in breast cancer survivors

It is widely acknowledged in theoretical and empirical literature that social relationships, comprising of structural measures (social networks) and functional measures (perceived social support) have an undeniable effect on health outcomes. However, the actual mechanism of this effect has yet to be clearly understood or explicated. In addition, comorbidity is found to adversely affect social relationships and health related quality of life (a valued outcome measure in cancer patients and survivors). ^ This cross sectional study uses selected baseline data (N=3088) from the Women's Healthy Eating and Living (WHEL) study. Lisrel 8.72 was used for the latent variable structural equation modeling. Due to the ordinal nature of the data, Weighted Least Squares (WLS) method of estimation using Asymptotic Distribution Free covariance matrices was chosen for this analysis. The primary exogenous predictor variables are Social Networks and Comorbidity; Perceived Social Support is the endogenous predictor variable. Three dimensions of HRQoL, physical, mental and satisfaction with current quality of life were the outcome variables. ^ This study hypothesizes and tests the mechanism and pathways between comorbidity, social relationships and HRQoL using latent variable structural equation modeling. After testing the measurement models of social networks and perceived social support, a structural model hypothesizing associations between the latent exogenous and endogenous variables was tested. The results of the study after listwise deletion (N=2131) mostly confirmed the hypothesized relationships (TLI, CFI >0.95, RMSEA = 0.05, p=0.15). Comorbidity was adversely associated with all three HRQoL outcomes. Strong ties were negatively associated with perceived social support; social network had a strong positive association with perceived social support, which served as a mediator between social networks and HRQoL. Mental health quality of life was the most adversely affected by the predictor variables. ^ This study is a preliminary look at the integration of structural and functional measures of social relationships, comorbidity and three HRQoL indicators using LVSEM. Developing stronger social networks and forming supportive relationships is beneficial for health outcomes such as HRQoL of cancer survivors. Thus, the medical community treating cancer survivors as well as the survivor's social networks need to be informed and cognizant of these possible relationships. ^

The application of latent variable models to the assessment of determinants of HIV risk behavior

Studies on the relationship between psychosocial determinants and HIV risk behaviors have produced little evidence to support hypotheses based on theoretical relationships. One limitation inherent in many articles in the literature is the method of measurement of the determinants and the analytic approach selected. ^ To reduce the misclassification associated with unit scaling of measures specific to internalized homonegativity, I evaluated the psychometric properties of the Reactions to Homosexuality scale in a confirmatory factor analytic framework. In addition, I assessed the measurement invariance of the scale across racial/ethnic classifications in a sample of men who have sex with men. The resulting measure contained eight items loading on three first-order factors. Invariance assessment identified metric and partial strong invariance between racial/ethnic groups in the sample. ^ Application of the updated measure to a structural model allowed for the exploration of direct and indirect effects of internalized homonegativity on unprotected anal intercourse. Pathways identified in the model show that drug and alcohol use at last sexual encounter, the number of sexual partners in the previous three months and sexual compulsivity all contribute directly to risk behavior. Internalized homonegativity reduced the likelihood of exposure to drugs, alcohol or higher numbers of partners. For men who developed compulsive sexual behavior as a coping strategy for internalized homonegativity, there was an increase in the prevalence odds of risk behavior. ^ In the final stage of the analysis, I conducted a latent profile analysis of the items in the updated Reactions to Homosexuality scale. This analysis identified five distinct profiles, which suggested that the construct was not homogeneous in samples of men who have sex with men. Lack of prior consideration of these distinct manifestations of internalized homonegativity may have contributed to the analytic difficulty in identifying a relationship between the trait and high-risk sexual practices. ^

INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

The association between parity and selected structural birth defects

BACKGROUND: Parity is a risk factor in neonatal morbidity and mortality. This dissertation examined the association between first births and selected birth defects. The first aim was to assess the risk of 66 birth defects among first births and third or greater births. The second aim was to determine if maternal race, maternal age, infant sex or infant birth weight modify the association between first births and selected birth defects. METHODS: The Texas Birth Defects Registry provided data for 1999-2009. For the first aim, odds ratios were calculated for each birth defect. For the second aim, analysis was restricted to the ten birth defects significantly associated with first births. Stratified analyses were conducted and interaction terms were added to logistic regression models to assess whether differences in the odds ratios for the effect of first birth were statistically significant across strata. RESULTS: Findings for the first aim showed that first births had significantly increased odds of having an infant with 24 of the 66 birth defects. Third or greater births had significantly increased odds of having four of the 66 birth defects. For the second aim, a number of significant effect modifiers were observed. For patent ductus arteriosis, obstructive urinary defects and gastroschisis, the effect of first births was significantly modified by black or U.S.-born Hispanic mothers. The effect of first birth was also significantly modified among mothers ≥30 years for mitral valve insufficiency, atrial septal defect and congenital hip dislocation. The effect of first births was significantly modified among infants with low birth weight for hypospadias, congenital hip dislocation and gastroschisis. CONCLUSIONS: First births were associated with an elevated risk of 24 categories of birth defects. For some of the birth defects studied, the effect of first birth is modified by maternal age, maternal race and low birth weight. Knowledge of the increased risk for birth defects among women having their first birth allows physicians and midwives to provide better patient care and spur further research into the etiology of associated birth defects. This knowledge may bring about interventions prior to conception in populations most likely to conceive.^

The Gepatschferner from 1850 - 2006, with links to digital elevation models and glacier margins

This theses investigates changes at Gepatschferner in length, area and volume since the last glacier maximum in 1850. Changes are discussed for the following time periods: 1850-1922, 1922-1971, 1971-1997, 1997-2006. Digital elevation models were created for 1850 from geomorphological data and for 1922 and 1971 from historical maps. Existing DEMs for 1997 and 2006 were further analysed. Since 1850 Gepatschferner has retreated by 2 km in length and has lost 32% of its area and 36% of its volume. The rate of loss of volume is increasing faster than the rate of loss of area and losses in the upper regions of the glacier are becoming increasingly more important to overall losses. The largest losses per 50 m elevation increment occur at the tongue. These losses are greatest in the most recent time step studied, 1997-2006, and exceed previous values by 40% and more. The data base includes the glacier margins, elevations models as they have been compiled within the thesis (DEMs of 1997 and 2006 are part of the glacier inventories, length changes are part of the length change data base of the Austrian Alpine Club).

Trioctahedral micas in xenolithic ejecta from recent volcanism of the Somma-Vesuvius

This study reports the first crystal chemical database resulting from a detailed structural investigation of trioctahedral micas found in xenolithic ejecta produced during the AD 1631, 1872 and 1944 eruptions, three explosive episodes of recent volcanic period of Vesuvius volcano (Southern Italy). Three xenolith types were selected: metamorphic/metasomatic skarns, pyrometamorphic/hydrothermally altered nodules and mafic cumulates. They are related to different magma chemistry and effusive styles: from sub-plinian and most evolved (AD 1631 eruption) to violent strombolian with medium evolution degree (AD 1872 eruption) to vulcanian-effusive, least evolved (AD 1944 eruption) event, respectively. Both xenoliths and micas were investigated employing multiple techniques: the xenoliths were characterized by X-ray fluorescence, inductively-coupled plasma-mass spectrometry, optical microscopy, X-ray powder diffraction, and quantitative energy-dispersive microanalysis; the micas were studied by electron probe microanalysis and single crystal X-ray diffraction. The mica-bearing xenoliths show variable texture and mineralogical assemblage, clearly related to their different origin. Based on the major oxide chemistry, only one xenolithic sample falls in the skarn compositional field from the Somma-Vesuvius literature, some fall close to the skarns and cumulate fields, others plot close to the syenite/foidolite/essexite field. A subgroup of the selected ejecta does not fall or approach any of the compositional fields. Trace and rare earth element patterns show some petrological affinity between studied xenoliths and erupted magmas with typical Eu, Ta and Nb negative anomalies. Strongly depleted patterns were detected for the 1631 metamorphic/metasomatic skarns xenoliths. Three distinct mica groups were distinguished: 1) Mg-, Al-rich, low Ti-bearing, low to moderate F-bearing varieties (1631 xenolith), 2) Al-moderate, F- and Mg-rich, Ti-, Fe-poor varieties (1872 xenolith), and 3) Al-, Ti- and Fe-rich, F-poor phases (1944 xenolith). All the analysed mica crystals are 1M polytypes with the expected space group C2/m. Micas from xenoliths of the 1631 Vesuvius eruption are phlogopites characterized by a combination of low extent of oxy-type and variable extent OH-F-substitutions, as testified by the range of F concentration (from ~ 0.20 to 0.80 apfu). Micas from xenoliths of the 1872 Vesuvius eruption exhibit structural peculiarities typical of fluorophlogopites, i.e. OH-F-substitution is predominant. Micas from the xenolith of the 1944 Vesuvius eruption display features typical of oxy-substituted micas. The variability of the crystal chemical features of the studied micas are consistent with the remarkable variation of their host rocks. Micas from 1631 nodules are related to metasomatic, skarn-type environment, deriving from the metamorphosed wall-rocks hosting the magma reservoir. The fluorophlogopites from the 1872 xenoliths testify for strongly dehydrated environmental conditions compared to those of the 1631 and 1944 hosts. Finally, magma storage condition at depth, associated to a decreasing aH2O may have promoted major oxy-type substitutions in 1944 biotites.

Timing of the most recent Neoglacial advance and retreat in Potter Cove, Antarctic Peninsula

The timing of the most recent Neoglacial advance in the Antarctic Peninsula is important for establishing global climate teleconnections and providing important post-glacial rebound corrections to gravity-based satellite measurements of ice loss. However, obtaining accurate ages from terrestrial geomorphic and sedimentary indicators of the most recent Neoglacial advance in Antarctica has been hampered by the lack of historical records and the difficulty of dating materials in Antarctica. Here we use a new approach to dating flights of raised beaches in the South Shetland Islands of the northern Antarctic Peninsula to bracket the age of a Neoglacial advance that occurred between 1500 and 1700 AD, broadly synchronous with compilations for the timing of the Little Ice Age in the northern hemisphere. Our approach is based on optically stimulated luminescence of the underside of buried cobbles to obtain the age of beaches previously shown to have been deposited immediately inside and outside the moraines of the most recent Neoglacial advance. In addition, these beaches mark the timing of an apparent change in the rate of isostatic rebound thought to be in response to the same glacial advance within the South Shetland Islands. We use a Maxwell viscoelastic model of glacial-isostatic adjustment (GIA) to determine whether the rates of uplift calculated from the raised beaches are realistic given the limited constraints on the ice advance during this most recent Neoglacial advance. Our rebound model suggests that the subsequent melting of an additional 16-22% increase in the volume of ice within the South Shetland Islands would result in a subsequent uplift rate of 12.5 mm/yr that lasted until 1840 AD resulting in a cumulative uplift of 2.5 m. This uplift rate and magnitude are in close agreement with observed rates and magnitudes calculated from the raised beaches since the most recent Neoglacial advance along the South Shetland Islands and falls within the range of uplift rates from similar settings such as Alaska.

Exploring the sources of China's challenge to Japan : models of industrial organisation in the motorcycle industry

In view of the recent rise of China, this paper looks into one of the most important yet relatively overlooked ingredients of the Chinese success: industrial organisation. It will examine the case of the motorcycle industry, in which the rise of Chinese manufacturers even disrupted the established dominance of Japanese industry leaders. Adopting the modified version of the global value chain governance framework, this paper shows that the rise of China has been driven by a distinctive arm’s-length model of industrial organisation, which is in sharp contrast to the conventional captive model that has sustained the Japanese leadership.

Refining predictions of iberian plant distribution: lessons from p. nigra and p. sylvestris palaeoecological-based habitat suitability models

In recent years, challenged by the climate scenarios put forward by the IPCC and its potential impact on plant distribution, numerous predictive techniques -including the so called habitat suitability models (HSM)- have been developed. Yet, as the output of the different methods produces different distribution areas, developing validation tools are strong needs to reduce uncertainties. Focused in the Iberian Peninsula, we propose a palaeo-based method to increase the robustness of the HSM, by developing an ecological approach to understand the mismatches between the palaeoecological information and the projections of the HSMs. Here, we present the result of (1) investigating causal relationships between environmental variables and presence of Pinus sylvestris L. and P. nigra Arn. available from the 3rd Spanish Forest Inventory, (2) developing present and past presence-predictions through the MaxEnt model for 6 and 21 kyr BP, and (3) assessing these models through comparisons with biomized palaeoecological data available from the European Pollen Database for the Iberian Peninsula.

Impact And Explosive Loads On Concrete Buildings Using Shell And Beam Type Elements

The threat of impact or explosive loads is regrettably a scenario to be taken into account in the design of lifeline or critical civilian buildings. These are often made of concrete and not specifically designed for military threats. Numerical simulation of such cases may be undertaken with the aid of state of the art explicit dynamic codes, however several difficult challenges are inherent to such models: the material modeling for the concrete anisotropic failure, consideration of reinforcement bars and important structural details, adequate modeling of pressure waves from explosions in complex geometries, and efficient solution to models of complete buildings which can realistically assess failure modes. In this work we employ LS-Dyna for calculation, with Lagrangian finite elements and explicit time integration. Reinforced concrete may be represented in a fairly accurate fashion with recent models such as CSCM model [1] and segregated rebars constrained within the continuum mesh. However, such models cannot be realistically employed for complete models of large buildings, due to limitations of time and computer resources. The use of structural beam and shell elements for this purpose would be the obvious solution, with much lower computational cost. However, this modeling requires careful calibration in order to reproduce adequately the highly nonlinear response of structural concrete members, including bending with and without compression, cracking or plastic crushing, plastic deformation of reinforcement, erosion of vanished elements etc. The main objective of this work is to provide a strategy for modeling such scenarios based on structural elements, using available material models for structural elements [2] and techniques to include the reinforcement in a realistic way. These models are calibrated against fully three-dimensional models and shown to be accurate enough. At the same time they provide the basis for realistic simulation of impact and explosion on full-scale buildings

Dynamic analysis of road vehicle-bridge systems under turbulent wind by means of Finite Element Models

When an automobile passes over a bridge dynamic effects are produced in vehicle and structure. In addition, the bridge itself moves when exposed to the wind inducing dynamic effects on the vehicle that have to be considered. The main objective of this work is to understand the influence of the different parameters concerning the vehicle, the bridge, the road roughness or the wind in the comfort and safety of the vehicles when crossing bridges. Non linear finite element models are used for structures and multibody dynamic models are employed for vehicles. The interaction between the vehicle and the bridge is considered by contact methods. Road roughness is described by the power spectral density (PSD) proposed by the ISO 8608. To consider that the profiles under right and left wheels are different but not independent, the hypotheses of homogeneity and isotropy are assumed. To generate the wind velocity history along the road the Sandia method is employed. The global problem is solved by means of the finite element method. First the methodology for modelling the interaction is verified in a benchmark. Following, the case of a vehicle running along a rigid road and subjected to the action of the turbulent wind is analyzed and the road roughness is incorporated in a following step. Finally the flexibility of the bridge is added to the model by making the vehicle run over the structure. The application of this methodology will allow to understand the influence of the different parameters in the comfort and safety of road vehicles crossing wind exposed bridges. Those results will help to recommend measures to make the traffic over bridges more reliable without affecting the structural integrity of the viaduct