A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Psychological pressure in competitive environments: Evidence from a randomized natural experiment

Much like cognitive abilities, emotional skills can have major effects on performance and economic outcomes. This paper studies the behavior of professionalsubjects involved in a dynamic competition in their own natural environment. Thesetting is a penalty shoot-out in soccer where two teams compete in a tournamentframework taking turns in a sequence of five penalty kicks each. As the kicking order is determined by the random outcome of a coin flip, the treatment and control groups are determined via explicit randomization. Therefore, absent any psychological effects, both teams should have the same probability of winning regardless of the kicking order. Yet, we find a systematic first-kicker advantage. Using data on 2,731 penalty kicks from 262 shoot-outs for a three decade period, we find that teams kicking first win the penalty shoot-out 60.5% of the time. A dynamic panel data analysis shows that the psychological mechanism underlying this result arises from the asymmetry in the partial score. As most kicks are scored, kicking first typically means having the opportunity to lead in the partial score, whereas kicking second typically means lagging in the score and having the opportunity to, at most, get even. Having a worse prospect than the opponent hinders subjects' performance.Further, we also find that professionals are self-aware of their own psychological effects. When a recent change in regulations gives winners of the coin toss the chance to choose the kicking order, they rationally react to it by systematically choosing to kick first. A survey of professional players reveals that when asked to explain why they prefer to kick first, they precisely identify the psychological mechanism for which we find empirical support in the data: they want to lead in the score inorder to put pressure on the opponent.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Variance reduction methods for simulation of densities on Wiener space

We develop a general error analysis framework for the Monte Carlo simulationof densities for functionals in Wiener space. We also study variancereduction methods with the help of Malliavin derivatives. For this, wegive some general heuristic principles which are applied to diffusionprocesses. A comparison with kernel density estimates is made.

Tax simulation in the SHP

This paper describes a simulation package designed to estimate the annual income taxes paid by respondents of the Swiss Household Panel (SHP). In Switzerland, the 26 cantons have their own tax system. Additionally, tax levels vary between the over 2000 municipalities and over time. The simulation package takes account of this complexity by building on existing tables on tax levels which are provided by the Swiss Federal Tax Administration Office. Because these are limited to a few types of households and only 812 municipalities, they have to be extended to cover all households and municipalities. A further drawback of these tables is that they neglect several deductions. The tax simulation package fills this gap by taking additionally account of deductions for children, double-earner couples, third pillar and support for dependent persons according to cantonal legislation. The resulting variable on direct taxes not only serves to calculate household income net of taxes, but can also be a variable for analysis by its own account.

Decision making in uncertain and changing environments

We consider an agent who has to repeatedly make choices in an uncertainand changing environment, who has full information of the past, who discountsfuture payoffs, but who has no prior. We provide a learning algorithm thatperforms almost as well as the best of a given finite number of experts orbenchmark strategies and does so at any point in time, provided the agentis sufficiently patient. The key is to find the appropriate degree of forgettingdistant past. Standard learning algorithms that treat recent and distant pastequally do not have the sequential epsilon optimality property.

Carabid beetle assemblages in three environments in the Araucaria humid forest of southern Brazil

Carabid beetle assemblages in three environments in the Araucaria humid forest of southern Brazil. Carabidae is composed mainly by ground-dwelling predator beetles. It is the fourth most diverse group within Coleoptera, but its diversity in the Neotropical region is understudied. Here we describe and analyze the diversity of carabid beetles in a region of subtropical rain forest dominated by Araucaria angustifolia with different landscapes. Three areas were chosen in an environmental integrity gradient: primary forests, secondary forests and old Pinus plantations. Pitfall traps were taken monthly, in a total of 14 samples per area. 1733 adult carabid beetles, belonging to 18 species, were sampled. There were differences in richness and abundance between the sampled areas. The total scores followed the same tendency: primary forests (14 species/747 individuals), secondary forests (13/631) and Pinus forests (10/355). An analysis of similarity shows differences in species composition, for both areas and seasons. Galerita lacordarei was the most abundant species for all samples and seasons. Carabid species show similar responses in accordance with habitat heterogeneity and disturbance. The abundance of Galerita lacordarei was influenced by temperature, for all sampled sites. Environmental changes affect the carabid assemblages and decrease diversity, possibly interfering in local dynamics. Seasonality patterns seem to indicate an increase in individual movement during summer, probably in search of resources. It is suggested that microhabitat patchiness is probably an important factor affecting carabid beetle diversity at small spatial scales.

Biosphere/lithosphere interface : microbially-mediated CaCO3 precipitation in hypersaline and freshwater environments

Abstract: Microbial mats very efficiently cycle elements, such as C, 0, N, S and H, which makes them key players of redox processes at the biosphere-lithosphere interface. They are characterized by high metabolic activities and high turnover rates (production and consumption) of biomass, which mainly consists of cell material and of extracellular organic matter (EOM). The EOM forms a matrix, embedding the microbial cells and fulfilling various functions within the microbial mat, including: mat attachment to surfaces; creation of micro-domains within the mat; physical stabilization under hy- drodynamic stress and the protection of the cells in multiple other stress conditions. EOM mainly consists of polysaccharides, amino acids, and a variety of chemical func-tional groups {e.g., -C00H, - SH -OH). These groups strongly bind cations such as Ca2+ and Mg2+ and thus exert a strong control on carbonate mineral formation within the microbial mat. A feedback mechanism between community metabolisms, their prod¬ucts, and the surrounding physicochemical microenvironment thus influences the de¬gree of carbonate saturation favoring either carbonate precipitation or dissolution. We investigated the driving forces and mechanisms of microbialite formation in the Sari ne River, FR, Switzerland, the hypersaline lake, Big Pond, Bahamas and in labo¬ratory experiments. The two fundamentally different natural systems allowed us to compare the geochemical conditions and microbial metabolisms, necessary for car¬bonate formation in microbial mats. Although carbonates are oversaturated in both environments, precipitation does not occur on physicochemical substrates (i.e. out¬side the microbial mats). In the Sarine a high crystal nucleation threshold exceeds the carbonate saturation, despite the high carbonate alkalinity in the water column. Cyanobacterial photosynthesis strongly locally enhances the carbonate alkalinity, whereas the EOM attract and immobilize calcium, which increases the saturation state and finally leads to carbonate precipitation within the EOM (in this case the cyanobacterial sheath) as nucleation template. In Big Pond, the presence of calcium- chelating anions (i.e. sulfate) and EOM, as well as the presence of magnesium, lowers the calcium activity in the water column and mat, and thus inhibits carbonate pre¬cipitation. Coupled with other heterotrophic metabolisms, sulfate reduction uses the EOM as carbon source, degrading it. The resulting EOM consumption creates alkalin¬ity, releases calcium and consumes sulfate in mat-micro domains, which leads to the formation of carbonate layers at the top of the microbial mat. Résumé: Interface biosphère/lithosphère: médiation microbienne de la précipitation de CaC03 dans des environnements en eaux douces et hypersalines Les tapis microbiens engendrent une circulation très efficace des éléments, tels que C, 0, N, S et H, ce qui en fait des acteurs clé pour les processus d'oxydoréduction à l'inter¬face biosphère-lithosphère. Ils sont caractérisés par des taux élevés d'activité méta¬bolique, ainsi que par la production et la consommation de biomasse, principalement constituée de cellules microbiennes et de matière organique extracellulaire (MOE). Dans un tapis microbien, les cellules microbiennes sont enveloppées par une matrice de MOE qui a différentes fonctions dont l'attachement du tapis aux surfaces, la créa¬tion de micro-domaines dans le tapis, la stabilisation physique en situation de stress hydrodynamique, et la protection des cellules dans de multiples autres conditions de stress. La MOE se compose principalement de polysaccharides, d'acides aminés, et d'une variété de groupes fonctionnels chimiques (par exemple, COOH, -SH et -OH). Ces groupes se lient fortement aux cations, tels que Ca2+ et Mg2+, et exercent ainsi un contrôle fort sur la formation de CaC03 dans le tapis microbien. Un mécanisme de rétroaction, entre les métabolismes de la communauté microbienne, leurs produits, et le microenvironnement physico-chimique, influence le degré de saturation de car¬bonate, favorisant soit leur précipitation, soit leur dissolution. Nous avons étudié le moteur et les mécanismes de minéralisation dans des tapis de la Sarine, FR, Suisse et du lac hypersalin, Big Pond, aux Bahamas, ainsi que durant des expériences en laboratoire. Les deux systèmes naturels, fondamentalement dif¬férents, nous ont permis de comparer les conditions géochimiques et les métabolis¬mes nécessaires à la formation des carbonates dans des tapis microbiens. Bien que les carbonates soient sursaturés dans les deux environnements, la précipitation ne se produit pas sur des substrats physico-chimiques (en dehors du tapis microbien). Dans la Sarine, malgré un taux d'alcalinité élevé, les valeurs de seuil pour la nucléa- tion de carbonates sont plus hautes que la saturation du carbonate. La photosynthèse cyanobactérienne augmente localement l'alcalinité, alors que la MOE attire et immo¬bilise le calcium, ce qui augmente l'état de saturation et conduit finalement à la pré¬cipitation des carbonates, en utilisant la MOE comme substrat de nucléation. À Big Pond, la présence de chélateurs de calcium, notamment les anions (p.ex. le sulfate) et la MOE, ainsi que la présence de magnésium, réduit l'activité du calcium et inhibe en conséquence la précipitation des carbonates. Couplée avec d'autres métabolismes hétérotrophes, la réduction des sulfates utilise la MOE comme source de carbone, en la dégradant. Cette consommation de MOE crée l'alcalinité, consomme des sulfates et libère du calcium dans des micro-domaines, conduisant à la formation de couches de carbonates dans le haut du tapis microbien.

Characterization and Application of Educational Virtual Environments for Science Teaching in Secondary School

This work describes the characteristics of a representative set of seven different virtual laboratories (VLs) aimed for science teaching in secondary school. For this purpose, a 27-item evaluation model that facilitates the characterization of the VLs was prepared. The model takes into account the gaming features, the overall usability, and also the potential to induce scientific literacy. Five of the seven VLs were then tested with two larger and highly heterogenic groups of students, and in two different contexts – biotechnology and physics, respectively. It is described how the VLs were received by the students, taking into account both their motivation and their self-reported learning outcome. In some cases, students’ approach to work with the VLs was recorded digitally, and analyzed qualitatively. In general, the students enjoyed the VL activities, and claimed that they learned from them. Yet, more investigation is required to address the effectiveness of these tools for significant learning.

PENELOPE-2008: A Code System for Monte Carlo Simulation of Electron and Photon Transport

The computer code system PENELOPE (version 2008) performs Monte Carlo simulation of coupledelectron-photon transport in arbitrary materials for a wide energy range, from a few hundred eV toabout 1 GeV. Photon transport is simulated by means of the standard, detailed simulation scheme.Electron and positron histories are generated on the basis of a mixed procedure, which combinesdetailed simulation of hard events with condensed simulation of soft interactions. A geometry packagecalled PENGEOM permits the generation of random electron-photon showers in material systemsconsisting of homogeneous bodies limited by quadric surfaces, i.e., planes, spheres, cylinders, etc. Thisreport is intended not only to serve as a manual of the PENELOPE code system, but also to provide theuser with the necessary information to understand the details of the Monte Carlo algorithm.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

Stratigraphy, sedimentology and depositional environments of the Permian to uppermost Cretaceous Batain Group, eastern-Oman

Permian to Late Cretaceous allochthonous sedimentary and volcanic rocks exposed in the Batain area (eastern Oman Margin) have received comparably little attention in the past. They largely were considered as part of the Hamrat Duru Group (Hawasina Complex) of the northern Oman Mountains. Structural, kinematic and biostratigraphic results from our mapping campaign in the Batain area have now revealed, that emplacement of these units occurred in a WNW direction during latest Cretaceous/Early Paleogene time. This clearly contrasts with previous models that postulated a S-ward directed obduction in Campanian times such as recorded from the Hawasina Complex and Semail Ophiolite in the Oman Mountains. We herewith establish the `'Batain Group'' comprising all Permian to Late Cretaceous allochthonous units in the Batain Area. These are: 1.) the Permian Qarari Formation deposited in the toe of a slope setting; 2.) the Late Permian to late Liassic Al Jil Formation comprising periplatform detritus and very coarse breccias; 3.) the Scythian to Norian Matbat Formation formed by slope deposits; 4.) the Early Jurassic to early Oxfordian Guwayza Formation with high energy platform detritus; 5.) the Mid-Jurassic to earliest Cretaceous Ruwaydah Formation seamount; and 6.) the Oxfordian to Santonian Wahrah Formation, mainly radiolarites; and 7.) the Santonian to latest Maastrichtian Fayah Formation built by flysch-type sediments. These sedimentary and volcanic rocks represent deposits of the former ``Batain basin'' off eastern-Oman, destroyed by compressional tectonics at the Cretaceous/Paleogene transition. For tectono-stratigraphic reasons the Batain Group does not form part of the Hawasina Complex.