Modeling asteroid collisions and impact processes

As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.

Natural variations at position 93 of the invariant Vα24-Jα18 α chain of human iNKT-cell TCRs strongly impact on CD1d binding

Human invariant natural killer T (NKT) cell TCRs bind to CD1d via an "invariant" Vα24-Jα18 chain (iNKTα) paired to semi-invariant Vβ11 chains (iNKTβ). Single-amino acid variations at position 93 (p93) of iNKTα, immediately upstream of the "invariant" CDR3α region, have been reported in a substantial proportion of human iNKT-cell clones (4-30%). Although p93, a serine in most human iNKT-cell TCRs, makes no contact with CD1d, it could affect CD1d binding by altering the conformation of the crucial CDR3α loop. By generating recombinant refolded iNKT-cell TCRs, we show that natural single-nucleotide variations in iNKTα, translating to serine, threonine, asparagine or isoleucine at p93, exert a powerful effect on CD1d binding, with up to 28-fold differences in affinity between these variants. This effect was observed with CD1d loaded with either the artificial α-galactosylceramide antigens KRN7000 or OCH, or the endogenous glycolipid β-galactosylceramide, and its importance for autoreactive recognition of endogenous lipids was demonstrated by the binding of variant iNKT-cell TCR tetramers to cell surface expressed CD1d. The serine-containing variant showed the strongest CD1d binding, offering an explanation for its predominance in vivo. Complementary molecular dynamics modeling studies were consistent with an impact of p93 on the conformation of the CDR3α loop.

Modeling of red blood cell life-spans in hematologically normal populations

Despite the impact of red blood cell (RBC) Life-spans in some disease areas such as diabetes or anemia of chronic kidney disease, there is no consensus on how to quantitatively best describe the process. Several models have been proposed to explain the elimination process of RBCs: random destruction process, homogeneous life-span model, or a series of 4-transit compartment model. The aim of this work was to explore the different models that have been proposed in literature, and modifications to those. The impact of choosing the right model on future outcomes prediction--in the above mentioned areas--was also investigated. Both data from indirect (clinical data) and direct life-span measurement (biotin-labeled data) methods were analyzed using non-linear mixed effects models. Analysis showed that: (1) predictions from non-steady state data will depend on the RBC model chosen; (2) the transit compartment model, which considers variation in life-span in the RBC population, better describes RBC survival data than the random destruction or homogenous life-span models; and (3) the additional incorporation of random destruction patterns, although improving the description of the RBC survival data, does not appear to provide a marked improvement when describing clinical data.

A model of osteoporosis impact in Switzerland 2000-2020

The aim of our study was to develop a modeling framework suitable to quantify the incidence, absolute number and economic impact of osteoporosis-attributable hip, vertebral and distal forearm fractures, with a particular focus on change over time, and with application to the situation in Switzerland from 2000 to 2020. A Markov process model was developed and analyzed by Monte Carlo simulation. A demographic scenario provided by the Swiss Federal Statistical Office and various Swiss and international data sources were used as model inputs. Demographic and epidemiologic input parameters were reproduced correctly, confirming the internal validity of the model. The proportion of the Swiss population aged 50 years or over will rise from 33.3% in 2000 to 41.3% in 2020. At the total population level, osteoporosis-attributable incidence will rise from 1.16 to 1.54 per 1,000 person-years in the case of hip fracture, from 3.28 to 4.18 per 1,000 person-years in the case of radiographic vertebral fracture, and from 0.59 to 0.70 per 1,000 person-years in the case of distal forearm fracture. Osteoporosis-attributable hip fracture numbers will rise from 8,375 to 11,353, vertebral fracture numbers will rise from 23,584 to 30,883, and distal forearm fracture numbers will rise from 4,209 to 5,186. Population-level osteoporosis-related direct medical inpatient costs per year will rise from 713.4 million Swiss francs (CHF) to CHF946.2 million. These figures correspond to 1.6% and 2.2% of Swiss health care expenditures in 2000. The modeling framework described can be applied to a wide variety of settings. It can be used to assess the impact of new prevention, diagnostic and treatment strategies. In Switzerland incidences of osteoporotic hip, vertebral and distal forearm fracture will rise by 33%, 27%, and 19%, respectively, between 2000 and 2020, if current prevention and treatment patterns are maintained. Corresponding absolute fracture numbers will rise by 36%, 31%, and 23%. Related direct medical inpatient costs are predicted to increase by 33%; however, this estimate is subject to uncertainty due to limited availability of input data.

Modeling the stratospheric warming following the Mt. Pinatubo eruption: uncertainties in aerosol extinctions

In terms of atmospheric impact, the volcanic eruption of Mt. Pinatubo (1991) is the best characterized large eruption on record. We investigate here the model-derived stratospheric warming following the Pinatubo eruption as derived from SAGE II extinction data including recent improvements in the processing algorithm. This method, termed SAGE_4λ, makes use of the four wavelengths (385, 452, 525 and 1024 nm) of the SAGE II data when available, and uses a data-filling procedure in the opacity-induced "gap" regions. Using SAGE_4λ, we derived aerosol size distributions that properly reproduce extinction coefficients also at much longer wavelengths. This provides a good basis for calculating the absorption of terrestrial infrared radiation and the resulting stratospheric heating. However, we also show that the use of this data set in a global chemistry–climate model (CCM) still leads to stronger aerosol-induced stratospheric heating than observed, with temperatures partly even higher than the already too high values found by many models in recent general circulation model (GCM) and CCM intercomparisons. This suggests that the overestimation of the stratospheric warming after the Pinatubo eruption may not be ascribed to an insufficient observational database but instead to using outdated data sets, to deficiencies in the implementation of the forcing data, or to radiative or dynamical model artifacts. Conversely, the SAGE_4λ approach reduces the infrared absorption in the tropical tropopause region, resulting in a significantly better agreement with the post-volcanic temperature record at these altitudes.

Modeling the heterogeneous ice and gas coma of Comet 103P/Hartley 2

The spectacular images of Comet 103P/Hartley 2 recorded by the Medium Resolution Instrument (MRI) and High Resolution Instrument (HRI) on board of the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) spacecraft, as the Deep Impact extended mission, revealed that its bi-lobed very active nucleus outgasses volatiles heterogeneously. Indeed, CO2 is the primary driver of activity by dragging out chunks of pure ice out of the nucleus from the sub-solar lobe that appear to be the main source of water in Hartley 2's coma by sublimating slowly as they go away from the nucleus. However, water vapor is released by direct sublimation of the nucleus at the waist without any significant amount of either CO2 or icy grains. The coma structure for a comet with such areas of diverse chemistry differs from the usual models where gases are produced in a homogeneous way from the surface. We use the fully kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J. 685, 659-677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J. 732, 104-120) applied to Comet 103P/Hartley 2 including sublimating icy grains to reproduce the observations made by EPOXI and ground-based measurements. A realistic bi-lobed nucleus with a succession of active areas with different chemistry was included in the model enabling us to study in details the coma of Hartley 2. The different gas production rates from each area were found by fitting the spectra computed using a line-by-line non-LTE radiative transfer model to the HRI observations. The presence of icy grains with long lifetimes, which are pushed anti-sunward by radiation pressure, explains the observed OH asymmetry with enhancement on the night side of the coma.

Weber fractions and time-order errors for long and short durations: Implications for modeling

Using a weighted up-down procedure, in each of eight conditions 28 participants compared durations of auditory (noise bursts) or visual (LED flashes) intervals; filled or unfilled with 3-ms markers; with or without feedback. Standards (Sts) were 100 and 1000 ms, and the ISI 900 ms. Intermixedly, presentation orders were St-Comparison (Co) and Co-St. TOEs were positive for St=100-ms and negative for St=1000 ms. Weber fractions (WFs, JND/St) were lowered by feedback. For visual-filled and visual-empty, WFs were highest for St=100 ms. For auditory-filled and visual-empty, St interacted with Order: lowest WFs occurred for St-Co with St=1000 ms, but for Co-St with St=100 ms. Lowest average WFs occurred with St-Co for visual-filled, but with Co-St for visual-empty. The results refute the generalization of better discrimination with St-Co than with Co-St (”type-B effect”), and support the notion of sensation weighting: flexibly differential impact weights of the compared durations in generating the response.

Impact of dams, dam removal and dam-related river engineering structures on sediment connectivity and channel morphology of the Fugnitz and the Kaja Rivers

In terms of changing flow and sediment regimes of rivers, dams are often regarded as the most dominant form of human impact on fluvial systems. Dams can decrease the flux of water and sediments leading to channel changes such as upstream aggradation and downstream degradation. The opposite effects occur when dams are removed. Channel degradation often requires further intervention in terms of river bed and bank protection works. The situation evolves more complex in river systems that are impacted by a series of dams due to feedback processes between the different system compartments. A number of studies have recently investigated geomorphic systems using connectivity approaches to improve the understanding of geomorphic system response to change. This paper presents a case study investigating the impact of dam construction, dam removal and dam-related river bed and bank protection measures on the sediment connectivity and channel morphology of the Fugnitz and the Kaja Rivers using a combination of DEM analyses, field surveys and landscape evolution modelling. For both river systems the results revealed low sediment connectivity accompanied by a fine river bed sediment facies in river sections upstream of active dams and of removed dams with protection measures. Contrarily, high sediment connectivity which was accompanied by a coarse river bed sediment facies was observed in river sections either located downstream of active dams or of removed dams with upstream protection. In terms of channel changes, significant channel degradation was examined at locations downstream of active dams and of removed dams. Channel bed and bank protection measures prevent erosion and channel slope recovery after dam removal. Landscape evolution modeling revealed a complex geomorphic response to dam construction and dam removal as sediment output rates and therefore geomorphic processes have been shown to act in a non-linear manner. These insights are deemed to have major implications for river management and conservation, as quality and state of riverine habitats are determined by channel morphology and river bed sediment composition.

Impact of development of job control and task-related stressors on resigned job attitude.

With more experience in the labor market, some job characteristics increase, some decrease. For example, among young employees who just entered the labor market, job control may initially be low but increase with more routine and experience. Job control is a job resource that is valued in itself and is positively associated with job satisfaction; but job control also helps dealing with stressors at work. There is little research on correlated changes, but the existing evidence suggests a joint development over time. However, even less is known about the relevance of such changes for employees. Usually, research tends to use mean levels to predict mean levels in outcomes but development in job control and stressors may be as relevant for job satisfaction as having a certain level in those job characteristics. Job satisfaction typically is regarded as a positive attitude towards one’s work. What has received less attention is that some employees may lower their expectations if their job situation does not reflect their needs, resulting in a resigned attitude towards one’s job. The present study investigates the development of job control and task-related stressors over ten years and tests the predictive value of changes in job control and task-related stressors for resigned attitude towards one’s job. We used data from a Swiss panel study (N=356) ranging over ten years. Job control, task-related stressors (an index consisting of time pressure, concentration demands, performance constraints, interruptions, and uncertainty about tasks), and resigned attitude towards one’s job were assessed in 1998, 1999, 2001, and 2008. Latent growth modeling revealed that growth rates of job control and task-related stressors were not correlated with one another. We predicted resigned attitude towards one’s job in 2008 a) by initial levels, and b) by changes in job control and stressors, controlling for resigned attitude in 1998. There was some prediction by initial levels (job control: β = -.15, p < .05; task-related stressors: β = .12, p = .06). However, as expected, changes in control and stressors predicted resigned attitude much better, with β = -.37, p < .001, for changes in job control, and β = .31, p < .001, for changes in task-related stressors. Our data confirm the importance of having low levels of task-related stressors and higher levels of job control for job attitudes. However, development in these job characteristics seems even more important than initial levels.