Tectono-metamorphic evolution of a nappe stack: A case study of the Swiss Alps

Fold-and-thrust belts are prominent structures that occur at the front of compressional orogens. To unravel the tectonic and metamorphic evolution of such complexes, kinematic investigations, quantitative microstructural analysis and geothermometry (calcite–graphite, calcite–dolomite) were performed on carbonate mylonites from thrust faults of the Helvetic nappe stack in Central Switzerland. Paleo-isotherms of peak temperature conditions and cooling stages (fission track) of the nappe pile were reconstructed in a vertical section and linked with the microstructural and kinematic evolution. Mylonitic microstructures suggest that under metamorphic conditions close to peak temperature, strain was highly localized within thrust faults where deformation temperatures spatially continuously increased in both directions, from N to S within each nappe and from top–down in the nappe stack, covering a temperature range of 180–380 °C. Due to the higher metamorphic conditions, thrusting of the lowermost nappe, the Doldenhorn nappe, was accompanied by a much more pronounced nappe internal ductile deformation of carbonaceous rock types than was the case for the overlying Wildhorn- and Gellihorn nappes. Ongoing thrusting brought the Doldenhorn nappe closer to the surface. The associated cooling resulted in a freezing in of the paleo-isotherms of peak metamorphic conditions. Contemporaneous shearing localized in the basal thrust, initially still in the ductile deformation regime and finally as brittle faulting and cataclasis inducing ultimately an inverse metamorphic zonation. With ongoing exhumation and the formation of the Helvetic antiformal nappe stack, a bending of large-scale tectonic structures (thrusts, folds), peak temperature isotherms and cooling isotherms occurred. While this local bending can directly be attributed to active deformation underneath the section investigated up to times of 2–3 ma, a more homogeneous uplift of the entire region is suggested for the very late and still active exhumation stage.

Quantitative Derivation of Effective Evolution Equations for the Dynamics of Bose-Einstein Condensates

This thesis proves certain results concerning an important question in non-equilibrium quantum statistical mechanics which is the derivation of effective evolution equations approximating the dynamics of a system of large number of bosons initially at equilibrium (ground state at very low temperatures). The dynamics of such systems are governed by the time-dependent linear many-body Schroedinger equation from which it is typically difficult to extract useful information due to the number of particles being large. We will study quantitatively (i.e. with explicit bounds on the error) how a suitable one particle non-linear Schroedinger equation arises in the mean field limit as number of particles N → ∞ and how the appropriate corrections to the mean field will provide better approximations of the exact dynamics. In the first part of this thesis we consider the evolution of N bosons, where N is large, with two-body interactions of the form N³ᵝv(Nᵝ⋅), 0≤β≤1. The parameter β measures the strength and the range of interactions. We compare the exact evolution with an approximation which considers the evolution of a mean field coupled with an appropriate description of pair excitations, see [18,19] by Grillakis-Machedon-Margetis. We extend the results for 0 ≤ β < 1/3 in [19, 20] to the case of β < 1/2 and obtain an error bound of the form p(t)/Nᵅ, where α>0 and p(t) is a polynomial, which implies a specific rate of convergence as N → ∞. In the second part, utilizing estimates of the type discussed in the first part, we compare the exact evolution with the mean field approximation in the sense of marginals. We prove that the exact evolution is close to the approximate in trace norm for times of the order o(1)√N compared to log(o(1)N) as obtained in Chen-Lee-Schlein [6] for the Hartree evolution. Estimates of similar type are obtained for stronger interactions as well.

Evolution of mylonitic microfabrics (EMM), a computer application for educational purposes

`Evolution of mylonitic microfabrics' (EMM) is an interactive Filemaker Pro 3.0 application that documents a series of see-through deformation experiments on polycrystalline norcamphor. The application comprises computer animations, graphics and text explanations designed to give students and researchers insight into the interaction and dynamic nature of small-scale, mylonitic processes like intracrystalline glide, dynamic recrystallization and strain localization (microshearing). EMM shows how mylonitic steady state is achieved at different strain rates and temperatures. First, rotational mechanisms like glide-induced vorticity, subgrain rotation recrystallization and rigid-body rotation bring grains' crystal lattices into orientations that are favorable for intracrystalline glide. In a second stage, selective elimination of grains whose lattices are poorly oriented for glide involves grain boundary migration. This strengthens the texture. Temperature and strain rate affect both the relative activity of different strain accommodation mechanisms and the rate of microfabric change. Steady-state microfabrics are characterized by stable texture, grain size and shape-preferred orientations of grains and domains. This involves the cyclical generation and elimination of dynamically recrystallized grains and microshear zones.

Small strain elastic behaviour of unsaturated soil investigated by bender/extender element testing

The aim of this project was to investigate very small strain elastic behaviour of soils under unsaturated conditions, using bender/extender element (BEE) testing. The behaviour of soils at very small strains has been widely studied under saturated conditions, whereas much less work has been performed on very small strain behaviour under unsaturated conditions. A suction-controlled double wall triaxial apparatus for unsaturated soil testing was modified to incorporate three pairs of BEEs transmitting both shear and compression waves with vertical and horizontal directions of wave transmission and wave polarisation. Various different techniques for measuring wave travel time were investigated in both the time domain and the frequency domain and it was concluded that, at least for the current experimental testing programme, peak-to-first-peak in the time domain was the most reliable technique for determining wave travel time. An experimental test programme was performed on samples of compacted speswhite kaolin clay. Two different forms of compaction were employed (i.e. isotropic and anisotropic). Compacted kaolin soil samples were subjected to constant suction loading and unloading stages at three different values of suction, covering both unsaturated conditions (s= 50kPa and s= 300kPa) and saturated conditions (s=0). Loading and unloading stages were performed at three different values of stress ratio (η=0, η=1 and η=-1 ). In some tests a wetting-drying cycle was performed before or within the loading stage, with the wetting-drying cycles including both wetting-induced swelling and wetting-induced collapse compression. BEE tests were performed at regular intervals throughout all test stages, to measure shear wave velocity Vs and compression wave velocity Vp and hence to determine values of shear modulus G and constrained modulus M. The experimental test programme was designed to investigate how very small strain shear modulus G and constrained modulus M varied with unsaturated state variables, including how anisotropy of these parameters developed either with stress state (stress-induced anisotropy) or with previous straining (strain-induced anisotropy). A new expression has been proposed for the very small strain shear modulus G of an isotropic soil under saturated and unsaturated conditions. This expression relates the variation of G to only mean Bishop’s stress p* and specific volume v, and it converges to a well-established expression for saturated soils as degree of saturation approaches 1. The proposed expression for G is able to predict the variation of G under saturated and unsaturated conditions at least as well as existing expressions from the literature and it is considerably simpler (employing fewer state variables and fewer soil constants). In addition, unlike existing expressions from the literature, the values of soil constants in the proposed new expression can be determined from a saturated test. It appeared that, in the current project at least, any strain-induced anisotropy of very small strain elastic behaviour was relatively modest, with the possible exception of loading in triaxial extension. It was therefore difficult to draw any firm conclusion about evolution of strain-induced anisotropy and whether it depended upon the same aspects of soil fabric as evolution of anisotropy of large strain plastic behaviour. Stress-induced anisotropy of very small strain elastic behaviour was apparent in the experimental test programme. An attempt was made to extend the proposed expression for G to include the effect of stress-induced anisotropy. Interpretation of the experimental results indicated that the value of shear modulus was affected by the values of all three principal Bishop’s stresses (in the direction of wave transmission, the direction of wave polarisation and the third mutually perpendicular direction). However, prediction of stress-induced anisotropy was only partially successful, and it was concluded that the effect of Lode angle was also significant.

Examining the growth and evolution of magnetic fields in clusters of galaxies: a numerical perspective

Magnetic fields are ubiquitous in galaxy cluster atmospheres and have a variety of astrophysical and cosmological consequences. Magnetic fields can contribute to the pressure support of clusters, affect thermal conduction, and modify the evolution of bubbles driven by active galactic nuclei. However, we currently do not fully understand the origin and evolution of these fields throughout cosmic time. Furthermore, we do not have a general understanding of the relationship between magnetic field strength and topology and other cluster properties, such as mass and X-ray luminosity. We can now begin to answer some of these questions using large-scale cosmological magnetohydrodynamic (MHD) simulations of the formation of galaxy clusters including the seeding and growth of magnetic fields. Using large-scale cosmological simulations with the FLASH code combined with a simplified model of the acceleration of cosmic rays responsible for the generation of radio halos, we find that the galaxy cluster frequency distribution and expected number counts of radio halos from upcoming low-frequency sur- veys are strongly dependent on the strength of magnetic fields. Thus, a more complete understanding of the origin and evolution of magnetic fields is necessary to understand and constrain models of diffuse synchrotron emission from clusters. One favored model for generating magnetic fields is through the amplification of weak seed fields in active galactic nuclei (AGN) accretion disks and their subsequent injection into cluster atmospheres via AGN-driven jets and bubbles. However, current large-scale cosmological simulations cannot directly include the physical processes associated with the accretion and feedback processes of AGN or the seeding and merging of the associated SMBHs. Thus, we must include these effects as subgrid models. In order to carefully study the growth of magnetic fields in clusters via AGN-driven outflows, we present a systematic study of SMBH and AGN subgrid models. Using dark-matter only cosmological simulations, we find that many important quantities, such as the relationship between SMBH mass and galactic bulge velocity dispersion and the merger rate of black holes, are highly sensitive to the subgrid model assumptions of SMBHs. In addition, using MHD calculations of an isolated cluster, we find that magnetic field strengths, extent, topology, and relationship to other gas quantities such as temperature and density are also highly dependent on the chosen model of accretion and feedback. We use these systematic studies of SMBHs and AGN inform and constrain our choice of subgrid models, and we use those results to outline a fully cosmological MHD simulation to study the injection and growth of magnetic fields in clusters of galaxies. This simulation will be the first to study the birth and evolution of magnetic fields using a fully closed accretion-feedback cycle, with as few assumptions as possible and a clearer understanding of the effects of the various parameter choices.

Numerical investigations on the strain-adaptive bone remodelling in the periprosthetic femur: Influence of the boundary conditions

Background: There are several numerical investigations on bone remodelling after total hip arthroplasty (THA) on the basis of the finite element analysis (FEA). For such computations certain boundary conditions have to be defined. The authors chose a maximum of three static load situations, usually taken from the gait cycle because this is the most frequent dynamic activity of a patient after THA. Materials and methods: The numerical study presented here investigates whether it is useful to consider only one static load situation of the gait cycle in the FE calculation of the bone remodelling. For this purpose, 5 different loading cases were examined in order to determine their influence on the change in the physiological load distribution within the femur and on the resulting strain-adaptive bone remodelling. First, four different static loading cases at 25%, 45%, 65% and 85% of the gait cycle, respectively, and then the whole gait cycle in a loading regime were examined in order to regard all the different loadings of the cycle in the simulation. Results: The computed evolution of the apparent bone density (ABD) and the calculated mass losses in the periprosthetic femur show that the simulation results are highly dependent on the chosen boundary conditions. Conclusion: These numerical investigations prove that a static load situation is insufficient for representing the whole gait cycle. This causes severe deviations in the FE calculation of the bone remodelling. However, accompanying clinical examinations are necessary to calibrate the bone adaptation law and thus to validate the FE calculations.

Study of the development and evolution of seasonal thermocline and turbulent processes in the northern parts of Persian Gulf using numerical modeling

The Persian Gulf (PG) is a semi-enclosed shallow sea which is connected to open ocean through the Strait of Hormuz. Thermocline as a suddenly decrease of temperature in subsurface layer in water column leading to stratification happens in the PG seasonally. The forcing comprise tide, river inflow, solar radiation, evaporation, northwesterly wind and water exchange with the Oman Sea that influence on this process. In this research, analysis of the field data and a numerical (Princeton Ocean Model, POM) study on the summer thermocline development in the PG are presented. The Mt. Mitchell cruise 1992 salinity and temperature observations show that the thermocline is effectively removed due to strong wind mixing and lower solar radiation in winter but is gradually formed and developed during spring and summer; in fact as a result of an increase in vertical convection through the water in winter, vertical gradient of temperature is decreased and thermocline is effectively removed. Thermocline development that evolves from east to west is studied using numerical simulation and some existing observations. Results show that as the northwesterly wind in winter, at summer transition period, weakens the fresher inflow from Oman Sea, solar radiation increases in this time interval; such these factors have been caused the thermocline to be formed and developed from winter to summer even over the northwestern part of the PG. The model results show that for the more realistic monthly averaged wind experiments the thermocline develops as is indicated by summer observations. The formation of thermocline also seems to decrease the dissolved oxygen in water column due to lack of mixing as a result of induced stratification. Over most of PG the temperature difference between surface and subsurface increases exponentially from March until May. Similar variations for salinity differences are also predicted, although with smaller values than observed. Indeed thermocline development happens more rapidly in the Persian Gulf from spring to summer. Vertical difference of temperature increases to 9 centigrade degrees in some parts of the case study zone from surface to bottom in summer. Correlation coefficients of temperature and salinity between the model results and measurements have been obtained 0.85 and 0.8 respectively. The rate of thermcline development was found to be between 0.1 to 0.2 meter per day in the Persian Gulf during the 6 months from winter to early summer. Also it is resulted from the used model that turbulence kinetic energy increases in the northwestern part of the PG from winter to early summer that could be due to increase in internal waves activities and stability intensified through water column during this time.

Evolution of psychosocial factors at work in a French region

International audience

Evolution of strength and physical properties of ultramafic and carbonate rocks under hydrothermal conditions

Interaction of rocks with fluids can significantly change mineral assemblage and structure. This so-called hydrothermal alteration is ubiquitous in the Earth’s crust. Though the behavior of hydrothermally altered rocks can have planet-scale consequences, such as facilitating oceanic spreading along slow ridge segments and recycling volatiles into the mantle at subduction zones, the mechanisms involved in the hydrothermal alteration are often microscopic. Fluid-rock interactions take place where the fluid and rock meet. Fluid distribution, flux rate and reactive surface area control the efficiency and extent of hydrothermal alteration. Fluid-rock interactions, such as dissolution, precipitation and fluid mediated fracture and frictional sliding lead to changes in porosity and pore structure that feed back into the hydraulic and mechanical behavior of the bulk rock. Examining the nature of this highly coupled system involves coordinating observations of the mineralogy and structure of naturally altered rocks and laboratory investigation of the fine scale mechanisms of transformation under controlled conditions. In this study, I focus on fluid-rock interactions involving two common lithologies, carbonates and ultramafics, in order to elucidate the coupling between mechanical, hydraulic and chemical processes in these rocks. I perform constant strain-rate triaxial deformation and constant-stress creep tests on several suites of samples while monitoring the evolution of sample strain, permeability and physical properties. Subsequent microstructures are analyzed using optical and scanning electron microscopy. This work yields laboratory-based constraints on the extent and mechanisms of water weakening in carbonates and carbonation reactions in ultramafic rocks. I find that inundation with pore fluid thereby reducing permeability. This effect is sensitive to pore fluid saturation with respect to calcium carbonate. Fluid inundation weakens dunites as well. The addition of carbon dioxide to pore fluid enhances compaction and partial recovery of strength compared to pure water samples. Enhanced compaction in CO2-rich fluid samples is not accompanied by enhanced permeability reduction. Analysis of sample microstructures indicates that precipitation of carbonates along fracture surfaces is responsible for the partial restrengthening and channelized dissolution of olivine is responsible for permeability maintenance.

Evolution of factors affecting mechanical olive harvesting

Harvest efficiency is defined as the percentage of fruits harvested by total production. The percentage of fruits harvested is less than 100% when working with trunk shakers to detach olives. It is important to increase the percentage of fruits harvested in order to increase farmer’s income. This objective can be achieved knowing the evolution of the main factors affecting fruit detachment. Fruit removal force (FRF), fruit weight (P) and the ratio between them are important for harvest efficiency. Field trials took place for two years (2013-2014) in Vilariça Valley, northeast Portugal in an olive orchard with ‘Cobrançosa Transmontana’ cultivar. It was adopted a mechanical harvesting system based on a trunk shaker to detach fruits, and an inverted umbrella to collect fruits. Elementary operation times were measured in seconds to evaluate work rates. FRF and P were measured in the ripening period, to evaluate their evolution. In this paper are presented the preliminary results of the ratio FRF (fruit removal force)/fruit weight evolution during the ripening period (P) and the results of the equipment work rate (trees h-1). The ratio FRF/P has predominantly descendant values in the weeks before harvest, from 140 to 80 as a result of a FRF downward variation from 4.9 to 2.94 N and an upward variation of P from 0.0294 to 0.0637 N. The FRF/P ratio stabilizes the decline in the last week of November just before harvesting, registering in some cases a slight increase in consequence of FRF increase higher than P increase (contrary to the tendency of previous weeks). Equipment work rate showed values between 40 and 57 trees h-1, confirming previous results.

El Estribo Volcanic Complex: Evolution from a shield volcano to a cinder cone, Pátzcuaro Lake, Michoacán, México

El Estribo Volcanic Complex (EVC) is located in the northern part of the Michoacán–Guanajuato Volcanic Field within the Trans-Mexican Volcanic Belt (TMVB). El Estribo is located at the southern edge of the E-W Pátzcuaro fault that belongs to the Pátzcuaro-Jarácuaro graben, a western extension of the E-W Morelia–Acambay fault system. Stratigraphy, geochronology, chemistry, and mineral assemblages suggest that the volcanic complex was constructed in two periods separated by a ~ 100 ka volcanic hiatus: a) emission of lava flows that constructed a shield volcano between 126 ka, and b) mixed phreatomagmatic to Strombolian activity that formed a cinder cone ~ 28 ka. The magmas that fed these monogenetic volcanoes were able to use the same feeding system. The cinder cone itself was constructed by Strombolian fallouts and remobilized scoria beds, followed by an erosion period, and by a mixed phreatomagmatic to magmatic phase (Strombolian fallouts ending with lava flows). Soft-sedimentary deformation of beds and impact sags, cross-bedding, as well as pitting and hydrothermal cracks found in particles support the phreatomagmatic phase. The erupted magmas through time ejected basaltic andesitic lava flows (56.21–58.88% SiO2) that built the shield volcano and then basaltic andesitic scoria (57.65–59.05% SiO2) that constructed the cinder cone. Although they used the same feeding system, the geochemical data and the mineral chemistry of the magmas indicate that the shield volcano and the cinder cone were fed by different magma batches erupted thousands of years apart. Therefore, the location of El Estribo Volcanic Complex along an E-W fault that has generated two sector collapses of the shield volcano to the north may be directly linked to this complex redistribution of the magmatic paths to the surface. Our findings show that magmatic feeding systems within monogenetic volcanic fields could be long lived, questioning the classic view of the monogenetic nature of their volcanoes and yielding information about the potential volcanic risk of these settings, usually considered risk-free.

Nonlinear evolution of harmonically forced perturbations on a wingtip vortex

Wingtip vortices are created by flying airplanes due to lift generation. The vortex interaction with the trailing aircraft has sparked researchers’ interest to develop an efficient technique to destroy these vortices. Different models have been used to describe the vortex dynamics and they all show that, under real flight conditions, the most unstable modes produce a very weak amplification. Another linear instability mechanism that can produce high energy gains in short times is due to the non-normality of the system. Recently, it has been shown that these non-normal perturbations also produce this energy growth when they are excited with harmonic forcing functions. In this study, we analyze numerically the nonlinear evolution of a spatially, pointwise and temporally forced perturbation, generated by a synthetic jet at a given radial distance from the vortex core. This type of perturbation is able to produce high energy gains in the perturbed base flow (10^3), and is also a suitable candidate for use in engineering applications. The flow field is solved for using fully nonlinear three-dimensional direct numerical simulation with a spectral multidomain penalty method model. Our novel results show that the nonlinear effects are able to produce locally small bursts of instability that reduce the intensity of the primary vortex.

DIFFUSE-INTERFACE FIELD APPROACH TO MODELING SELF-ASSEMBLY OF HETEROGENEOUS COLLOIDAL SYSTEMS AND RELATED DIPOLE-DIPOLE INTERACTION PHENOMENA

Colloid self-assembly under external control is a new route to fabrication of advanced materials with novel microstructures and appealing functionalities. The kinetic processes of colloidal self-assembly have attracted great interests also because they are similar to many atomic level kinetic processes of materials. In the past decades, rapid technological progresses have been achieved on producing shape-anisotropic, patchy, core-shell structured particles and particles with electric/magnetic charges/dipoles, which greatly enriched the self-assembled structures. Multi-phase carrier liquids offer new route to controlling colloidal self-assembly. Therefore, heterogeneity is the essential characteristics of colloid system, while so far there still lacks a model that is able to efficiently incorporate these possible heterogeneities. This thesis is mainly devoted to development of a model and computational study on the complex colloid system through a diffuse-interface field approach (DIFA), recently developed by Wang et al. This meso-scale model is able to describe arbitrary particle shape and arbitrary charge/dipole distribution on the surface or body of particles. Within the framework of DIFA, a Gibbs-Duhem-type formula is introduced to treat Laplace pressure in multi-liquid-phase colloidal system and it obeys Young-Laplace equation. The model is thus capable to quantitatively study important capillarity related phenomena. Extensive computer simulations are performed to study the fundamental behavior of heterogeneous colloidal system. The role of Laplace pressure is revealed in determining the mechanical equilibrium of shape-anisotropic particles at fluid interfaces. In particular, it is found that the Laplace pressure plays a critical role in maintaining the stability of capillary bridges between close particles, which sheds light on a novel route to in situ firming compact but fragile colloidal microstructures via capillary bridges. Simulation results also show that competition between like-charge repulsion, dipole-dipole interaction and Brownian motion dictates the degree of aggregation of heterogeneously charged particles. Assembly and alignment of particles with magnetic dipoles under external field is studied. Finally, extended studies on the role of dipole-dipole interaction are performed for ferromagnetic and ferroelectric domain phenomena. The results reveal that the internal field generated by dipoles competes with external field to determine the dipole-domain evolution in ferroic materials.

Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.

Assessment of an oral Mycobacterium bovis BCG vaccine and an inactivated M. bovis preparation for wild boar in terms of adverse reactions, vaccine strain survival, and uptake by nontarget species

Wildlife vaccination is increasingly being considered as an option for tuberculosis control. We combined data from laboratory trials and an ongoing field trial to assess the risk of an oral Mycobacterium bovis BCG vaccine and a prototype heat-inactivated Mycobacterium bovis preparation for Eurasian wild boar (Sus scrofa). We studied adverse reactions, BCG survival, BCG excretion, and bait uptake by nontarget species. No adverse reactions were observed after administration of BCG (n = 27) or inactivated M. bovis (n = 21). BCG was not found at necropsy (175 to 300 days postvaccination [n = 27]). No BCG excretion was detected in fecal samples (n = 162) or in urine or nasal, oral, or fecal swab samples at 258 days postvaccination (n = 29). In the field, we found no evidence of loss of BCG viability in baits collected after 36 h (temperature range, 11°C to 41°C). Camera trapping showed that wild boar (39%) and birds (56%) were the most frequent visitors to bait stations (selective feeders). Wild boar activity patterns were nocturnal, while diurnal activities were recorded for all bird species. We found large proportions of chewed capsules (29%) (likely ingestion of the vaccine) and lost baits (39%) (presumably consumed), and the proportion of chewed capsules showed a positive correlation with the presence of wild boar. Both results suggest proper bait consumption (68%). These results indicate that BCG vaccination in wild boar is safe and that, while bait consumption by other species is possible, this can be minimized by using selective cages and strict timing of bait deployment.