Secondary Microseisms Characterization and Green’s Function Extraction at the Larderello-Travale Geothermal Field (Italy)

Over the past ten years, the cross-correlation of long-time series of ambient seismic noise (ASN) has been widely adopted to extract the surface-wave part of the Green’s Functions (GF). This stochastic procedure relies on the assumption that ASN wave-field is diffuse and stationary. At frequencies <1Hz, the ASN is mainly composed by surface-waves, whose origin is attributed to the sea-wave climate. Consequently, marked directional properties may be observed, which call for accurate investigation about location and temporal evolution of the ASN-sources before attempting any GF retrieval. Within this general context, this thesis is aimed at a thorough investigation about feasibility and robustness of the noise-based methods toward the imaging of complex geological structures at the local (∼10-50km) scale. The study focused on the analysis of an extended (11 months) seismological data set collected at the Larderello-Travale geothermal field (Italy), an area for which the underground geological structures are well-constrained thanks to decades of geothermal exploration. Focusing on the secondary microseism band (SM;f>0.1Hz), I first investigate the spectral features and the kinematic properties of the noise wavefield using beamforming analysis, highlighting a marked variability with time and frequency. For the 0.1-0.3Hz frequency band and during Spring- Summer-time, the SMs waves propagate with high apparent velocities and from well-defined directions, likely associated with ocean-storms in the south- ern hemisphere. Conversely, at frequencies >0.3Hz the distribution of back- azimuths is more scattered, thus indicating that this frequency-band is the most appropriate for the application of stochastic techniques. For this latter frequency interval, I tested two correlation-based methods, acting in the time (NCF) and frequency (modified-SPAC) domains, respectively yielding esti- mates of the group- and phase-velocity dispersions. Velocity data provided by the two methods are markedly discordant; comparison with independent geological and geophysical constraints suggests that NCF results are more robust and reliable.

Starch distribution in pear tree organs in relation to training systems, rootstocks and fruit quality

Starch is the main form in which plants store carbohydrates reserves, both in terms of amounts and distribution among different plant species. Carbohydrates are direct products of photosynthetic activity, and it is well know that yield efficiency and production are directly correlated to the amount of carbohydrates synthesized and how these are distributed among vegetative and reproductive organs. Nowadays, in pear trees, due to the modernization of orchards, through the introduction of new rootstocks and the development of new training systems, the understanding and the development of new approaches regarding the distribution and storage of carbohydrates, are required. The objective of this research work was to study the behavior of carbohydrate reserves, mainly starch, in different pear tree organs and tissues: i.e., fruits, leaves, woody organs, roots and flower buds, at different physiological stages during the season. Starch in fruit is accumulated at early stages, and reached a maximum concentration during the middle phase of fruit development; after that, its degradation begins with a rise in soluble carbohydrates. Moreover, relationships between fruit starch degradation and different fruit traits, soluble sugars and organic acids were established. In woody organs and roots, an interconversion between starch and soluble carbohydrates was observed during the dormancy period that confirms its main function in supporting the growth and development of new tissues during the following spring. Factors as training systems, rootstocks, types of bearing wood, and their position on the canopy, influenced the concentrations of starch and soluble carbohydrates at different sampling dates. Also, environmental conditions and cultural practices must be considered to better explain these results. Thus, a deeper understanding of the dynamics of carbohydrates reserves within the plant could provide relevant information to improve several management practices to increase crop yield efficiency.

Expression and function of 5-hydroxytryptamine 4 receptors in smooth muscle preparations from the duodenum, ileum, and pelvic flexure of horses without gastrointestinal tract disease

OBJECTIVE: To evaluate the expression of the 5-hydroxytryptamine 4 (5-HT4) receptor subtype and investigate the modulating function of those receptors on contractility in intestinal tissues obtained from horses without gastrointestinal tract disease. SAMPLE POPULATION: Smooth muscle preparations from the duodenum, ileum, and pelvic flexure collected immediately after slaughter of 24 horses with no history or signs of gastrointestinal tract disease. PROCEDURES: In isometric organ baths, the contractile activities of smooth muscle preparations in response to 5-hydroxytryptamine and electric field stimulation were assessed; the effect of tegaserod alone or in combination with 5-hydroxytryptamine on contractility of intestinal specimens was also investigated. Presence and distribution of 5-HT4 receptors in intestinal tissues and localization on interstitial cells of Cajal were examined by use of an immunofluorescence technique. RESULTS: Widespread 5-HT4 receptor immunoreactivity was observed in all intestinal smooth muscle layers; 5-HT4 receptors were absent from the myenteric plexus and interstitial cells of Cajal. In electrical field-stimulated tissue preparations of duodenum and pelvic flexure, tegaserod increased the amplitude of smooth muscle contractions in a concentration-dependent manner. Preincubation with tegaserod significantly decreased the basal tone of the 5-HT-evoked contractility in small intestine specimens, compared with the effect of 5-HT alone, thereby confirming that tegaserod was acting as a partial agonist. CONCLUSIONS AND CLINICAL RELEVANCE: In horses, 5-HT4 receptors on smooth muscle cells appear to be involved in the contractile response of the intestinal tract to 5-hydroxytryptamine. Results suggest that tegaserod may be useful for treatment of reduced gastrointestinal tract motility in horses.

Vertebroplasty: experimental characterization of polymethylmethacrylate bone cement spreading as a function of viscosity, bone porosity, and flow rate

STUDY DESIGN: This is an experimental study on an artificial vertebra model and human cadaveric spine. OBJECTIVE: Characterization of polymethylmethacrylate (PMMA) bone cement distribution in the vertebral body as a function of cement viscosity, bone porosity, and injection speed. Identification of relevant parameters for improved cement flow predictability and leak prevention in vertebroplasty. SUMMARY OF BACKGROUND DATA: Vertebroplasty is an efficient procedure to treat vertebral fractures and stabilize osteoporotic bone in the spine. Severe complications result from bone cement leakage into the spinal canal or the vascular system. Cement viscosity has been identified as an important parameter for leak prevention but the influence of bone structure and injection speed remain obscure. METHODS: An artificial vertebra model based on open porous aluminum foam was used to simulate bone of known porosity. Fifty-six vertebroplasties with 4 different starting viscosity levels and 2 different injection speeds were performed on artificial vertebrae of 3 different porosities. A validation on a human cadaveric spine was executed. The experiments were radiographically monitored and the shape of the cement clouds quantitatively described with the 2 indicators circularity and mean cement spreading distance. RESULTS: An increase in circularity and a decrease in mean cement spreading distance was observed with increasing viscosity, with the most striking change occurring between 50 and 100 Pas. Larger pores resulted in significantly reduced circularity and increased mean cement spreading distance whereas the effect of injection speed on the 2 indicators was not significant. CONCLUSION: Viscosity is the key factor for reducing the risk of PMMA cement leakage and it should be adapted to the degree of osteoporosis encountered in each patient. It may be advisable to opt for a higher starting viscosity but to inject the material at a faster rate.

Subcellular distribution of the insulin-like growth factor (IGF) binding proteins (IGFBPs) 2 and 3 in articular chondrocytes

The insulin-like growth factor (IGF) is a major anabolic regulator in articular cartilage. The IGF-binding proteins (IGFBPs) are increased during osteoarthritis (OA), but the function of the later proteins remains unknown. In general, the IGFBPs are pluripotential effectors capable of IGF regulation and of acting on their own to control key cell functions, including survival and proliferation. The independent functions are often associated with their cell location, and therefore this study explores the distribution of IGFBP-2 and IGFBP-3 in articular chondrocytes. Immunohistochemistry was used to localize IGFBP-2 in normal human articular cartilage. Bovine chondrocytes were used for subcellular fractionation (hypotonic cell lysis) under nonreducing conditions and nuclear purification (centrifugation on sucrose cushions). Cell fraction markers and IGFBPs were assayed in the subcellular fractions by Western immunoblot. The IHC results showed association of IGFBP-2 with chondrocytes, but not with the nuclei. Subcellular fractionation of isolated chondrocytes yielded intact nuclei as assessed at the light microscopic level; the nuclear marker histone H1 was exclusively associated with this fraction. More than 90% of the cytoplasmic marker GAPDH and all the detectable IGFBP-2 were in the cytoplasmic fraction. Immunoreactive IGFBP-3 was found in the cytoplasmic and peri-nuclear/nuclear fractions. Chondrocytes contain intracellular IGFBP-2 and IGFBP-3 but only IGFBP-3 is associated with nuclei. This suggests the hypothesis that the actions of these IGFBPs in articular cartilage extend beyond the classic modulation of IGF receptor action.

IMPROVING FLOW DISTRIBUTION IN SMALL-SCALE WATER-SUPPLY SYSTEMS THROUGH THE USE OF FLOW-REDUCING DISCS AND METHODS FOR ANALYZING THE EFFECTIVENESS OF SOLAR POWERED PUMPING FOR A DRINKING WATER SUPPLY IN RURAL PANAMA

As continued global funding and coordination are allocated toward the improvement of access to safe sources of drinking water, alternative solutions may be necessary to expand implementation to remote communities. This report evaluates two technologies used in a small water distribution system in a mountainous region of Panama; solar powered pumping and flow-reducing discs. The two parts of the system function independently, but were both chosen for their ability to mitigate unique issues in the community. The design program NeatWork and flow-reducing discs were evaluated because they are tools taught to Peace Corps Volunteers in Panama. Even when ample water is available, mountainous terrains affect the pressure available throughout a water distribution system. Since the static head in the system only varies with the height of water in the tank, frictional losses from pipes and fittings must be exploited to balance out the inequalities caused by the uneven terrain. Reducing the maximum allowable flow to connections through the installation of flow-reducing discs can help to retain enough residual pressure in the main distribution lines to provide reliable service to all connections. NeatWork was calibrated to measured flow rates by changing the orifice coefficient (θ), resulting in a value of 0.68, which is 10-15% higher than typical values for manufactured flow-reducing discs. NeatWork was used to model various system configurations to determine if a single-sized flow-reducing disc could provide equitable flow rates throughout an entire system. There is a strong correlation between the optimum single-sized flow- reducing disc and the average elevation change throughout a water distribution system; the larger the elevation change across the system, the smaller the recommended uniform orifice size. Renewable energy can jump the infrastructure gap and provide basic services at a fraction of the cost and time required to install transmission lines. Methods for the assessment of solar powered pumping systems as a means for rural water supply are presented and assessed. It was determined that manufacturer provided product specifications can be used to appropriately design a solar pumping system, but care must be taken to ensure that sufficient water can be provided to the system despite variations in solar intensity.

Relaxation in hypertrophic cardiomyopathy and hypertensive heart disease: relations between hypertrophy and diastolic function

AIM To determine the relation between the extent and distribution of left ventricular hypertrophy and the degree of disturbance of regional relaxation and global left ventricular filling. METHODS Regional wall thickness (rWT) was measured in eight myocardial regions in 17 patients with hypertrophic cardiomyopathy, 12 patients with hypertensive heart disease, and 10 age matched normal subjects, and an asymmetry index calculated. Regional relaxation was assessed in these eight regions using regional isovolumetric relaxation time (rIVRT) and early to late peak filling velocity ratio (rE/A) derived from Doppler tissue imaging. Asynchrony of rIVRT was calculated. Doppler left ventricular filling indices were assessed using the isovolumetric relaxation time, the deceleration time of early diastolic filling (E-DT), and the E/A ratio. RESULTS There was a correlation between rWT and both rIVRT and rE/A in the two types of heart disease (hypertrophic cardiomyopathy: r = 0.47, p < 0.0001 for rIVRT; r = -0.20, p < 0.05 for rE/A; hypertensive heart disease: r = 0.21, p < 0.05 for rIVRT; r = -0.30, p = 0.003 for rE/A). The degree of left ventricular asymmetry was related to prolonged E-DT (r = 0. 50, p = 0.001) and increased asynchrony (r = 0.42, p = 0.002) in all patients combined, but not within individual groups. Asynchrony itself was associated with decreased E/A (r = -0.39, p = 0.01) and protracted E-DT (r = 0.69, p < 0.0001) and isovolumetric relaxation time (r = 0.51, p = 0.001) in all patients. These correlations were still significant for E-DT in hypertrophic cardiomyopathy (r = 0.56, p = 0.02) and hypertensive heart disease (r = 0.59, p < 0.05) and for isovolumetric relaxation time in non-obstructive hypertrophic cardiomyopathy (n = 8, r = 0.87, p = 0.005). CONCLUSIONS Non-invasive ultrasonographic examination of the left ventricle shows that in both hypertrophic cardiomyopathy and hypertensive heart disease, the local extent of left ventricular hypertrophy is associated with regional left ventricular relaxation abnormalities. Asymmetrical distribution of left ventricular hypertrophy is indirectly related to global left ventricular early filling abnormalities through regional asynchrony of left ventricular relaxation.

Expression and function of $\beta$-tubulin isotypes in Chinese hamster ovary (CHO) cells

The availability of isotype specific antisera for $\beta$-tubulin, coupled with genetic and biochemical analysis, has allowed the determination of $\beta$-tubulin isotype expression and distribution in Chinese hamster ovary (CHO) cells. Using genetic manipulations involving selection for colcemid resistance followed by reversion and reselection for drug resistance, we have succeeded in isolating cell lines that exhibit three major and one minor $\beta$-tubulin spots by two-dimensional gel electrophoresis. In concert with isotype specific antibodies, analysis of these mutants demonstrates that CHO cells express two copies of isotype I, at least one copy of isotype IV, and very small amounts of isotype V. Their stoichiometry is approximately 1:1:0.7:0.2. All three isotypes assemble into both cytoplasmic and spindle microtubules, and are similar in their responses to cold, colcemid, and calcium induced depolymerization. They have comparable turnover rates and are equally sensitive to depression of synthesis upon colchicine treatment. These results suggest that $\beta$-tubulin isotypes are used interchangeably to assemble microtubule structures in CHO cells. However, of 18 colcemid resistant mutants with a demonstrable alteration in $\beta$-tubulin, all were found to have the alteration in isotype I, thus leaving open the possibility that subtle differences in isotype properties may exist. Under various conditions of the cell growth, the relative proportion of each expressed isotype does not significantly seem to change except in the early G1 phase of the cell cycle. At this time the synthesis of isotype V increases more than two fold relative to isotype I and IV, while at the same time, total $\beta$-tubulin synthesis is decreased about 60-70%. ^

Cell surface galactosyltransferase structure and function during mesenchymal cell migration

In order to more fully understand the function of surface GalTase on mesenchymal cells, anti-GalTase IgG was used to (a) examine the role of surface GalTase during mouse mesenchymal cell migration on laminin and fibronectin; (b) define the plasma membrane distribution of GalTase by indirect immunofluorescence on migrating cells; (c) quantitate the level of surface GalTase on migrating cells; and (d) determine whether GalTase is associated with the cytoskeleton.^ Results show that anti-GalTase IgG was able to inhibit migration (48-80% as compared to basal rate) when cells were migrating on laminin-containing matrices. Monovalent Fab fragments inhibited migration on laminin by 90% after 4 hours. On the other hand, anti-GalTase IgG had no effect on cells migrating on fibronectin. This illustrates the substrate specificity of GalTase mediated-migration. When anti-GalTase IgG was used to localize surface GalTase on cells migratory on laminin, the enzyme was restricted to the leading and trailing edges of the cell. Assays indicate that GalTase is elevated approximately 3-fold when cells are migrating on laminin-containing matrices as compared to migratory cells on plastic or fibronectin, or as compared to stationary cells on any substrate. Laminin appears to recruit GalTase from preexisting intracellular pools to the growing lamellipodia.^ Double-label indirect immunofluorescence studies indicate that there is an apparent co-localization between some of the surface GalTase and some actin filaments. This relationship was explored by extracting cells prelabeled with anti-GalTase IgG and quantitated by radiolabeled second antibodies. Results show that 79% of the surface GalTase is associated with the cytoskeleton (as judged by detergent insolubility) when monovalent antibodies (Fab) are used. However virtually all (80-100%) of the surface GalTase can be induced to associate with the cytoskeleton when cross-linked with bivalent antibodies. Furthermore, when cells in suspension are incubated with divalent antibodies, an additional 66% of the surface GalTase can be induced to associate with the cytoskeleton. The elevated levels of surface GalTase detectable on cells migrating on laminin also appear to be associated with the cytoskeleton.^ Several lines of evidence suggest that GalTase is associated with F-actin. Data suggest that laminin induces the expression of surface GalTase to the growing lamellipodia where it becomes associated with the cytoskeleton leading to cell spreading and migration. (Abstract shortened with permission of author.) ^

Cytochrome P450 in mammalian brain: Purification, anatomical distribution and regulation of the cytochrome P450 monooxygenase system in rat brain

The cytochrome P450 (P450) monooxygenase system plays a major role in metabolizing a wide variety of xenobiotic as well as endogenous compounds. In performing this function, it serves to protect the body from foreign substances. However, in a number of cases, P450 activates procarcinogens to cause harm. In most animals, the highest level of activity is found in the liver. Virtually all tissues demonstrate P450 activity, though, and the role of the P450 monooxygenase system in these other organs is not well understood. In this project I have studied the P450 system in rat brain; purifying NADPH-cytochrome P450 reductase (reductase) from that tissue. In addition, I have examined the distribution and regulation of expression of reductase and P450 in various anatomical regions of the rat brain.^ NADPH-cytochrome P450 reductase was purified to apparent homogeneity and cytochrome P450 partially purified from whole rat brain. Purified reductase from brain was identical to liver P450 reductase by SDS-PAGE and Western blot techniques. Kinetic studies utilizing cerebral P450 reductase reveal Km values in close agreement with those determined with enzyme purified from rat liver. Moreover, the brain P450 reductase was able to function successfully in a reconstituted microsomal system with partially purified brain cytochrome P450 and with purified hepatic P4501A1 as measured by 7-ethoxycoumarin and 7-ethoxyresorufin O-deethylation. These results indicate that the reductase and P450 components may interact to form a competent drug metabolism system in brain tissue.^ Since the brain is not a homogeneous organ, dependent upon the well orchestrated interaction of numerous parts, pathology in one nucleus may have a large impact upon its overall function. Hence, the anatomical distribution of the P450 monooxygenase system in brain is important in elucidating its function in that organ. Related to this is the regulation of P450 expression in brain. In order to study these issues female rats--both ovariectomized and not--were treated with a number of xenobiotic compounds and sex steroids. The brains from these animals were dissected into 8 discrete regions and the presence and relative level of message for P4502D and reductase determined using polymerase chain reaction. Results of this study indicate the presence of mRNA for reductase and P4502D isoforms throughout the rat brain. In addition, quantitative PCR has allowed the determination of factors affecting the expression of message for these enzymes. ^

Structure and function study of prostaglandin H synthase

Prostaglandin H synthase (PGHS) is a key enzyme in biosynthesis of prostaglandins, thromboxane, and prostacyclin. It has two activities, cyclooxygenase and peroxidase. "PGHS" means PGHS-1. A current hypothesis considers the cyclooxygenase reaction to be a free radical chain reaction, initiated by interaction of the synthase peroxidase with hydroperoxides leading to the production of a tyrosyl free radical. According to this hypothesis, tyrosyl residue(s) may play a key role in the cyclooxygenase reaction. Tetranitromethane (TNM) can relatively selectively nitrate tyrosines at pH 8.0. The effect of TNM on both cyclooxygenase activity and peroxidase activity has been examined: reaction of the synthase holoenzyme with TNM at pH 8.0 led to inactivation of both activities, with the cyclooxygenase activity being lost rapidly and completely, while the peroxidase activity was lost more slowly. Indomethacin, a non-steroidal anti-inflammatory agent, can protect the synthase from the inactivation of TNM. Amino acid analyses indicated that a loss of tyrosine and formation of nitrotyrosine residues occurred during reaction with TNM, and that TNM-reacted holoenzyme with $<$10% residual cyclooxygenase activity had about 2.0 nitrotyrosine/subunit.^ PGH synthase is known to be an endoplasmic reticulum membrane-associated protein. Antibodies directed at particular PGHS peptide segments and indirect immunofluorescence have been used to characterize the membrane topology of crucial portions of PGHS. PGHS was expressed in COS-1 cells transfected with the appropriate cDNA. Stably-transfected human endothelial cells were also used for the topology study. The cells were treated with streptolysin-O, which selectively permeabilizes the plasma membrane, or with saponin to achieve general membrane disruption, before incubation with the antipeptide antibodies. Bound antipeptide antibody was stained by FITC-labelled secondary antibody and visualized by fluorescence microscopy. With the antipeptide antibodies against residues 51-66, 156-170 or 377-390, there was a significant reticular and perinuclear pattern of staining in cells permeabilized with saponin but not in cells permeabilized with SLO alone. Antibodies directed against the endogenous C-terminal peptide or against residues 271-284 produced staining in cells permeabilized with saponin, and also in a lower, but significant fraction of cells permeabilized with SLO. Similar results were obtained when COS-1 cells expressing recombinant PGHS with a viral reporter peptide inserted at the C-terminus were stained with antibody against the reporter epitope.^ The PGHS C-terminal sequence is similar to that of the consensus KDEL ER retention signal. The potential function of the PGHS C-terminus segment in ER retention was examined by mutating this segment and analyzing the subcellular distribution of the mutants expressed in COS-1 cells. None of the mutants had an altered subcellular distribution, although some had greatly diminished the enzyme activities. (Abstract shortened by UMI.) ^

Are shrubs really a sign of declining ecosystem function? Disentangling the myths and truths of woody encroachment in Australia

Since European settlement, there has been a dramatic increase in the density, cover and distribution of woody plants in former grassland and open woodland. There is a widespread belief that shrub encroachment is synonymous with declines in ecosystem functions, and often it is associated with landscape degradation or desertification. Indeed, this decline in ecosystem functioning is considered to be driven largely by the presence of the shrubs themselves. This prevailing paradigm has been the basis for an extensive program of shrub removal, based on the view that it is necessary to reinstate the original open woodland or grassland structure from which shrublands are thought to have been derived. We review existing scientific evidence, particularly focussed on eastern Australia, to question the notion that shrub encroachment leads to declines in ecosystem functions. We then summarise this scientific evidence into two conceptual models aimed at optimising landscape management to maximise the services provided by shrub-encroached areas. The first model seeks to reconcile the apparent conflicts between the patch- and landscape-level effects of shrubs. The second model identifies the ecosystem services derived from different stages of shrub encroachment. We also examined six ecosystem services provided by shrublands (biodiversity, soil C, hydrology, nutrient provision, grass growth and soil fertility) by using published and unpublished data. We demonstrated the following: (1) shrub effects on ecosystems are strongly scale-, species- and environment-dependent and, therefore, no standardised management should be applied to every case; (2) overgrazing dampens the generally positive effect of shrubs, leading to the misleading relationship between encroachment and degradation; (3) woody encroachment per se does not hinder any of the functions or services described above, rather it enhances many of them; (4) no single shrub-encroachment state (including grasslands without shrubs) will maximise all services; rather, the provision of ecosystem goods and services by shrublands requires a mixture of different states; and (5) there has been little rigorous assessment of the long-term effectiveness of removal and no evidence that this improves land condition in most cases. Our review provides the basis for an improved, scientifically based understanding and management of shrublands, so as to balance the competing goals of providing functional habitats, maintaining soil processes and sustaining pastoral livelihoods.

Measurement of the Z/ѵ* boson transverse momentum distribution in pp collisions at √s = 7 TeV with the ATLAS detector

This paper describes a measurement of the Z/ѵ* boson transverse momentum spectrum using ATLAS proton-proton collision data at a centre-of-mass energy of √s = 7 TeV at the LHC. The measurement is performed in the Z/ѵ* → e+e− and Z/ѵ* → μ+μ− channels, using data corresponding to an integrated luminosity of 4.7 fb−1. Normalized differential cross sections as a function of the Z/ѵ* boson transverse momentum are measured for transverse momenta up to 800 GeV. The measurement is performed inclusively for Z/ѵ* rapidities up to 2.4, as well as in three rapidity bins. The channel results are combined, compared to perturbative and resummed QCD calculations and used to constrain the parton shower parameters of Monte Carlo generators.

Temperature Prediction Model for Bone Drilling Based on Density Distribution and In Vivo Experiments for Minimally Invasive Robotic Cochlear Implantation

Surgical robots have been proposed ex vivo to drill precise holes in the temporal bone for minimally invasive cochlear implantation. The main risk of the procedure is damage of the facial nerve due to mechanical interaction or due to temperature elevation during the drilling process. To evaluate the thermal risk of the drilling process, a simplified model is proposed which aims to enable an assessment of risk posed to the facial nerve for a given set of constant process parameters for different mastoid bone densities. The model uses the bone density distribution along the drilling trajectory in the mastoid bone to calculate a time dependent heat production function at the tip of the drill bit. Using a time dependent moving point source Green's function, the heat equation can be solved at a certain point in space so that the resulting temperatures can be calculated over time. The model was calibrated and initially verified with in vivo temperature data. The data was collected in minimally invasive robotic drilling of 12 holes in four different sheep. The sheep were anesthetized and the temperature elevations were measured with a thermocouple which was inserted in a previously drilled hole next to the planned drilling trajectory. Bone density distributions were extracted from pre-operative CT data by averaging Hounsfield values over the drill bit diameter. Post-operative [Formula: see text]CT data was used to verify the drilling accuracy of the trajectories. The comparison of measured and calculated temperatures shows a very good match for both heating and cooling phases. The average prediction error of the maximum temperature was less than 0.7 °C and the average root mean square error was approximately 0.5 °C. To analyze potential thermal damage, the model was used to calculate temperature profiles and cumulative equivalent minutes at 43 °C at a minimal distance to the facial nerve. For the selected drilling parameters, temperature elevation profiles and cumulative equivalent minutes suggest that thermal elevation of this minimally invasive cochlear implantation surgery may pose a risk to the facial nerve, especially in sclerotic or high density mastoid bones. Optimized drilling parameters need to be evaluated and the model could be used for future risk evaluation.