Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.

A mathematical model to describe fat oxidation kinetics during graded exercise.

PURPOSE: The purpose of this study was to develop a mathematical model (sine model, SIN) to describe fat oxidation kinetics as a function of the relative exercise intensity [% of maximal oxygen uptake (%VO2max)] during graded exercise and to determine the exercise intensity (Fatmax) that elicits maximal fat oxidation (MFO) and the intensity at which the fat oxidation becomes negligible (Fatmin). This model included three independent variables (dilatation, symmetry, and translation) that incorporated primary expected modulations of the curve because of training level or body composition. METHODS: Thirty-two healthy volunteers (17 women and 15 men) performed a graded exercise test on a cycle ergometer, with 3-min stages and 20-W increments. Substrate oxidation rates were determined using indirect calorimetry. SIN was compared with measured values (MV) and with other methods currently used [i.e., the RER method (MRER) and third polynomial curves (P3)]. RESULTS: There was no significant difference in the fitting accuracy between SIN and P3 (P = 0.157), whereas MRER was less precise than SIN (P < 0.001). Fatmax (44 +/- 10% VO2max) and MFO (0.37 +/- 0.16 g x min(-1)) determined using SIN were significantly correlated with MV, P3, and MRER (P < 0.001). The variable of dilatation was correlated with Fatmax, Fatmin, and MFO (r = 0.79, r = 0.67, and r = 0.60, respectively, P < 0.001). CONCLUSIONS: The SIN model presents the same precision as other methods currently used in the determination of Fatmax and MFO but in addition allows calculation of Fatmin. Moreover, the three independent variables are directly related to the main expected modulations of the fat oxidation curve. SIN, therefore, seems to be an appropriate tool in analyzing fat oxidation kinetics obtained during graded exercise.

Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps)

The Mont Collon mafic complex is one of the best preserved examples of the Early Permian magmatism in the Central Alps, related to the intra-continental collapse of the Variscan belt. It mostly consists (> 95 vol.%) of ol+hy-nonnative plagioclase-wehrlites, olivine- and cpx-gabbros with cumulitic structures, crosscut by acid dikes. Pegmatitic gabbros, troctolites and anorthosites outcrop locally. A well-preserved cumulative, sequence is exposed in the Dents de Bertol area (center of intrusion). PT-calculations indicate that this layered magma chamber emplaced at mid-crustal levels at about 0.5 GPa and 1100 degrees C. The Mont Collon cumulitic rocks record little magmatic differentiation, as illustrated by the restricted range of clinopyroxene mg-number (Mg#(cpx)=83-89). Whole-rock incompatible trace-element contents (e.g. Nb, Zr, Ba) vary largely and without correlation with major-element composition. These features are characteristic of an in-situ crystallization process with variable amounts of interstitial liquid L trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distribution in the latter is homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in-situ crystallization equation successfully duplicated the REE concentrations in cumulitic minerals of all rock facies of the intrusion. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 20%, relative to the least evolved facies of the intrusion. L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. However, the in-situ crystallization model reaches its limitations with rock containing high modal content of REE-bearing minerals (i.e. zircon), such as pegmatitic gabbros. Dikes of anorthositic composition, locally crosscutting the layered lithologies, evidence that the Mont Collon rocks evolved in open system with mixing of intercumulus liquids of different origins and possibly contrasting compositions. The proposed model is not able to resolve these complex open systems, but migrating liquids could be partly responsible for the observed dispersion of points in some correlation diagrams. Absence of significant differentiation with recurrent lithologies in the cumulitic pile of Dents de Bertol points to an efficiently convective magma chamber, with possible periodic replenishment, (c) 2005 Elsevier B.V. All rights reserved.

Genetic- or transforming growth factor-beta 1-induced changes in epidermal peroxisome proliferator-activated receptor beta/delta expression dictate wound repair kinetics.

Advances in wound care are of great importance in clinical injury management. In this respect, the nuclear receptor peroxisome proliferator-activated receptor (PPAR)beta/delta occupies a unique position at the intersection of diverse inflammatory or anti-inflammatory signals that influence wound repair. This study shows how changes in PPARbeta/delta expression have a profound effect on wound healing. Using two different in vivo models based on topical application of recombinant transforming growth factor (TGF)-beta1 and ablation of the Smad3 gene, we show that prolonged expression and activity of PPARbeta/delta accelerate wound closure. The results reveal a dual role of TGF-beta1 as a chemoattractant of inflammatory cells and repressor of inflammation-induced PPARbeta/delta expression. Also, they provide insight into the so far reported paradoxical effects of the application of exogenous TGF-beta1 at wound sites.

Normal kinetics of intestinal glucose absorption in the absence of GLUT2: evidence for a transport pathway requiring glucose phosphorylation and transfer into the endoplasmic reticulum.

Glucose is absorbed through the intestine by a transepithelial transport system initiated at the apical membrane by the cotransporter SGLT-1; intracellular glucose is then assumed to diffuse across the basolateral membrane through GLUT2. Here, we evaluated the impact of GLUT2 gene inactivation on this transepithelial transport process. We report that the kinetics of transepithelial glucose transport, as assessed in oral glucose tolerance tests, was identical in the presence or absence of GLUT2; that the transport was transcellular because it could be inhibited by the SGLT-1 inhibitor phlorizin, and that it could not be explained by overexpression of another known glucose transporter. By using an isolated intestine perfusion system, we demonstrated that the rate of transepithelial transport was similar in control and GLUT2(-/-) intestine and that it was increased to the same extent by cAMP in both situations. However, in the absence, but not in the presence, of GLUT2, the transport was inhibited dose-dependently by the glucose-6-phosphate translocase inhibitor S4048. Furthermore, whereas transport of [(14)C]glucose proceeded with the same kinetics in control and GLUT2(-/-) intestine, [(14)C]3-O-methylglucose was transported in intestine of control but not of mutant mice. Together our data demonstrate the existence of a transepithelial glucose transport system in GLUT2(-/-) intestine that requires glucose phosphorylation and transfer of glucose-6-phosphate into the endoplasmic reticulum. Glucose may then be released out of the cells by a membrane traffic-based pathway similar to the one we previously described in GLUT2-null hepatocytes.

Effects of gastric bypass and gastric banding on glucose kinetics and gut hormone release.

BACKGROUND: Bariatric surgery markedly improves glucose homeostasis in patients with type 2 diabetes even before any significant weight loss is achieved. Procedures that involve bypassing the proximal small bowel, such as Roux-en-Y gastric bypass (RYGBP), are more efficient than gastric restriction procedures such as gastric banding (GB). OBJECTIVE: To evaluate the effects of RYGBP and GB on postprandial glucose kinetics and gastro-intestinal hormone secretion after an oral glucose load. METHODS AND PROCEDURES: This study was a cross-sectional comparison among non-diabetic, weight-stable women who had undergone RYGBP (n = 8) between 9 and 48 months earlier or GB (n = 6) from 25 to 85 months earlier, and weight- and age-matched control subjects (n = 8). The women were studied over 4 h following ingestion of an oral glucose load. Total glucose and meal glucose kinetics were assessed using glucose tracers and plasma insulin, and gut hormone concentrations were simultaneously monitored. RESULTS: Patients who had undergone RYGBP showed a a more rapid appearance of exogenous glucose in the systemic circulation and a shorter duration of postprandial hyperglycemia than patients who had undergone GB and C. The response in RYGBP patients was characterized by early and accentuated insulin response, enhanced postprandial levels of glucagon-like peptide-1 (GLP-1) and polypeptide YY (PYY), and greater postprandial suppression of ghrelin. DISCUSSION: These findings indicate that RYGBP is associated with alterations in glucose kinetics and glucoregulatory hormone secretion. These alterations are probably secondary to the anatomic rearrangement of the foregut, given the fact that they are not observed after GB. Increased PYY and GLP-1 concentrations and enhanced ghrelin suppression are compatible with reduced food intake after RYGBP.

Repeated injections of 131I-rituximab show patient-specific stable biodistribution and tissue kinetics.

PURPOSE: It is generally assumed that the biodistribution and pharmacokinetics of radiolabelled antibodies remain similar between dosimetric and therapeutic injections in radioimmunotherapy. However, circulation half-lives of unlabelled rituximab have been reported to increase progressively after the weekly injections of standard therapy doses. The aim of this study was to evaluate the evolution of the pharmacokinetics of repeated 131I-rituximab injections during treatment with unlabelled rituximab in patients with non-Hodgkin's lymphoma (NHL). METHODS: Patients received standard weekly therapy with rituximab (375 mg/m2) for 4 weeks and a fifth injection at 7 or 8 weeks. Each patient had three additional injections of 185 MBq 131I-rituximab in either treatment weeks 1, 3 and 7 (two patients) or weeks 2, 4 and 8 (two patients). The 12 radiolabelled antibody injections were followed by three whole-body (WB) scintigraphic studies during 1 week and blood sampling on the same occasions. Additional WB scans were performed after 2 and 4 weeks post 131I-rituximab injection prior to the second and third injections, respectively. RESULTS: A single exponential radioactivity decrease for WB, liver, spleen, kidneys and heart was observed. Biodistribution and half-lives were patient specific, and without significant change after the second or third injection compared with the first one. Blood T(1/2)beta, calculated from the sequential blood samples and fitted to a bi-exponential curve, was similar to the T(1/2) of heart and liver but shorter than that of WB and kidneys. Effective radiation dose calculated from attenuation-corrected WB scans and blood using Mirdose3.1 was 0.53+0.05 mSv/MBq (range 0.48-0.59 mSv/MBq). Radiation dose was highest for spleen and kidneys, followed by heart and liver. CONCLUSION: These results show that the biodistribution and tissue kinetics of 131I-rituximab, while specific to each patient, remained constant during unlabelled antibody therapy. RIT radiation doses can therefore be reliably extrapolated from a preceding dosimetry study.

Ambulatory measurement of ankle kinetics for clinical applications.

This study aimed to design and validate the measurement of ankle kinetics (force, moment, and power) during consecutive gait cycles and in the field using an ambulatory system. An ambulatory system consisting of plantar pressure insole and inertial sensors (3D gyroscopes and 3D accelerometers) on foot and shank was used. To test this system, 12 patients and 10 healthy elderly subjects wore shoes embedding this system and walked many times across a gait lab including a force-plate surrounded by seven cameras considered as the reference system. Then, the participants walked two 50-meter trials where only the ambulatory system was used. Ankle force components and sagittal moment of ankle measured by ambulatory system showed correlation coefficient (R) and normalized RMS error (NRMSE) of more than 0.94 and less than 13% in comparison with the references system for both patients and healthy subjects. Transverse moment of ankle and ankle power showed R>0.85 and NRMSE<23%. These parameters also showed high repeatability (CMC>0.7). In contrast, the ankle coronal moment of ankle demonstrated high error and lower repeatability. Except for ankle coronal moment, the kinetic features obtained by the ambulatory system could distinguish the patients with ankle osteoarthritis from healthy subjects when measured in 50-meter trials. The proposed ambulatory system can be easily accessible in most clinics and could assess main ankle kinetics quantities with acceptable error and repeatability for clinical evaluations. This system is therefore suggested for field measurement in clinical applications.

Are oxygen uptake kinetics modified when using a respiratory snorkel?

PURPOSE: The aim of this study was to compare VO2 kinetics during constant power cycle exercise measured using a conventional facemask (CM) or a respiratory snorkel (RS) designed for breath-by-breath analysis in swimming. METHODS: VO2 kinetics parameters-obtained using CM or RS, in randomized counterbalanced order-were compared in 10 trained triathletes performing two submaximal heavy-intensity cycling square-wave transitions. These VO2 kinetics parameters (ie, time delay: td1, td2; time constant: τ1, τ2; amplitude: A1, A2, for the primary phase and slow component, respectively) were modeled using a double exponential function. In the case of the RS data, this model incorporated an individually determined snorkel delay (ISD). RESULTS: Only td1 (8.9 ± 3.0 vs 13.8 ± 1.8 s, P < .01) differed between CM and RS, whereas all other parameters were not different (τ1 = 24.7 ± 7.6 vs 21.1 ± 6.3 s; A1 = 39.4 ± 5.3 vs 36.8 ± 5.1 mL x min(-1) x kg(-1); td2 = 107.5 ± 87.4 vs 183.5 ± 75.9 s; A2' (relevant slow component amplitude) = 2.6 ± 2.4 vs 3.1 ± 2.6 mL x min(-1) x kg(-1) for CM and RS, respectively). CONCLUSIONS: Although there can be a small mixture of breaths allowed by the volume of the snorkel in the transition to exercise, this does not appear to significantly influence the results. Therefore, given the use of an ISD, the RS is a valid instrument for the determination of VO2 kinetics within submaximal exercise.

Differences in whole-body fat oxidation kinetics between cycling and running.

This study aimed to quantitatively describe and compare whole-body fat oxidation kinetics in cycling and running using a sinusoidal mathematical model (SIN). Thirteen moderately trained individuals (7 men and 6 women) performed two graded exercise tests, with 3-min stages and 1 km h(-1) (or 20 W) increment, on a treadmill and on a cycle ergometer. Fat oxidation rates were determined using indirect calorimetry and plotted as a function of exercise intensity. The SIN model, which includes three independent variables (dilatation, symmetry and translation) that account for main quantitative characteristics of kinetics, provided a mathematical description of fat oxidation kinetics and allowed for determination of the intensity (Fat(max)) that elicits maximal fat oxidation (MFO). While the mean fat oxidation kinetics in cycling formed a symmetric parabolic curve, the mean kinetics during running was characterized by a greater dilatation (i.e., widening of the curve, P < 0.001) and a rightward asymmetry (i.e., shift of the peak of the curve to higher intensities, P = 0.01). Fat(max) was significantly higher in running compared with cycling (P < 0.001), whereas MFO was not significantly different between modes of exercise (P = 0.36). This study showed that the whole-body fat oxidation kinetics during running was characterized by a greater dilatation and a rightward asymmetry compared with cycling. The greater dilatation may be mainly related to the larger muscle mass involved in running while the rightward asymmetry may be induced by the specific type of muscle contraction.

Enriched basaltic andesites from mid-crustal fractional crystallization, recharge, and assimilation (Pilavo Volcano, Western Cordillera of Ecuador)

The origin of andesite is an important issue in petrology because andesite is the main eruptive product at convergent margins, corresponds to the average crustal composition and is often associated with major Cu-Au mineralization. In this study we present petrographic, mineralogical, geochemical and isotopic data for basaltic andesites of the latest Pleistocene Pilavo volcano, one of the most frontal volcanoes of the Ecuadorian Quaternary arc, situated upon thick (30-50 km) mafic crust composed of accreted Cretaceous oceanic plateau rocks and overlying mafic to intermediate Late Cretaceous-Late Tertiary magmatic arcs. The Pilavo rocks are basaltic andesites (54-57 center dot 5 wt % SiO(2)) with a tholeiitic affinity as opposed to the typical calc-alkaline high-silica andesites and dacites (SiO(2) 59-66 wt %) of other frontal arc volcanoes of Ecuador (e.g. Pichincha, Pululahua). They have much higher incompatible element contents (e.g. Sr 650-1350 ppm, Ba 650-1800 ppm, Zr 100-225 ppm, Th 5-25 ppm, La 15-65 ppm) and Th/La ratios (0 center dot 28-0 center dot 36) than Pichincha and Pululahua, and more primitive Sr ((87)Sr/(86)Sr similar to 0 center dot 7038-0 center dot 7039) and Nd (epsilon(Nd) similar to +5 center dot 5 to +6 center dot 1) isotopic signatures. Pilavo andesites have geochemical affinities with modern and recent high-MgO andesites (e.g. low-silica adakites, Setouchi sanukites) and, especially, with Archean sanukitoids, for both of which incompatible element enrichments are believed to result from interactions of slab melts with peridotitic mantle. Petrographic, mineral chemistry, bulk-rock geochemical and isotopic data indicate that the Pilavo magmatic rocks have evolved through three main stages: (1) generation of a basaltic magma in the mantle wedge region by flux melting induced by slab-derived fluids (aqueous, supercritical or melts); (2) high-pressure differentiation of the basaltic melt (at the mantle-crust boundary or at lower crustal levels) through sustained fractionation of olivine and clinopyroxene, leading to hydrous, high-alumina basaltic andesite melts with a tholeiitic affinity, enriched in incompatible elements and strongly impoverished in Ni and Cr; (3) establishment of one or more mid-crustal magma storage reservoirs in which the magmas evolved through dominant amphibole and clinopyroxene (but no plagioclase) fractionation accompanied by assimilation of the modified plutonic roots of the arc and recharge by incoming batches of more primitive magma from depth. The latter process has resulted in strongly increasing incompatible element concentrations in the Pilavo basaltic andesites, coupled with slightly increasing crustal isotopic signatures and a shift towards a more calc-alkaline affinity. Our data show that, although ultimately originating from the slab, incompatible element abundances in arc andesites with primitive isotopic signatures can be significantly enhanced by intra-crustal processes within a thick juvenile mafic crust, thus providing an additional process for the generation of enriched andesites.