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.

Use of high-resolution geophysical data to characterize heterogeneous aquifers: influence of data integration method on hydrological predictions

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. The investigation of exactly how much benefit can be brought by geophysical data in terms of its effect on hydrological predictions, however, has received considerably less attention in the literature. Here, we examine the potential hydrological benefits brought by a recently introduced simulated annealing (SA) conditional stochastic simulation method designed for the assimilation of diverse hydrogeophysical data sets. We consider the specific case of integrating crosshole ground-penetrating radar (GPR) and borehole porosity log data to characterize the porosity distribution in saturated heterogeneous aquifers. In many cases, porosity is linked to hydraulic conductivity and thus to flow and transport behavior. To perform our evaluation, we first generate a number of synthetic porosity fields exhibiting varying degrees of spatial continuity and structural complexity. Next, we simulate the collection of crosshole GPR data between several boreholes in these fields, and the collection of porosity log data at the borehole locations. The inverted GPR data, together with the porosity logs, are then used to reconstruct the porosity field using the SA-based method, along with a number of other more elementary approaches. Assuming that the grid-cell-scale relationship between porosity and hydraulic conductivity is unique and known, the porosity realizations are then used in groundwater flow and contaminant transport simulations to assess the benefits and limitations of the different approaches.

Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

What interventions facilitate client progress through the assimilation model? A task analysis of interventions in the psychodynamic treatment of depression.

A variation of task analysis was used to build an empirical model of how therapists may facilitate client assimilation process, described in the Assimilation of Problematic Experiences Scale. A rational model was specified and considered in light of an analysis of therapist in-session performances (N = 117) drawn from six inpatient therapies for depression. The therapist interventions were measured by the Comprehensive Psychotherapeutic Interventions Rating Scale. Consistent with the rational model, confronting interventions were particularly useful in helping clients elaborate insight. However, rather than there being a small number of progress-related interventions at lower levels of assimilation, therapists' use of interventions was broader than hypothesized and drew from a wide range of therapeutic approaches. Concerning the higher levels of assimilation, there was insufficient data to allow an analysis of the therapist's progress-related interventions.