Snow shielding factors for cosmogenic nuclide dating inferred from long-term neutron detector monitoring

The depth-dependent attenuation of the secondary cosmic-ray particle flux due to snow cover and its effects on production rates of cosmogenic nuclides constitutes a potential source of uncertainty for studies conducted in regions characterized by frequent seasonal snow burial. Recent experimental and numerical modelling studies have yielded new constraints on the effect of hydrogen-rich media on the production rates of cosmogenic nuclides by low- and high-energy neutrons (<10(-3) MeV and >10(2) MeV, respectively). Here we present long-term neutron-detector monitoring data from a natural setting that we use to quantify the effect of snow cover on the attenuation of fast neutrons (0.1-10 MeV), which are responsible for the production of Ne-21 from Mg and Cl-36 from K. We use data measured between July 2001 and May 2008 at seven stations located throughout the Ecrins-Pelvoux massif (French Western Alps) and its surroundings, at elevations ranging from 200 to 2500 m a.s.l. From the cosmic-ray fluxes recorded during summer, when snow is absent, we infer an apparent attenuation length of 148 g cm(-2) in the atmosphere at a latitude of similar to 45 degrees N and for altitudes ranging from similar to 200 to 2500 m a.s.l. Using snow water-equivalent (SWE) values obtained through snow-coring campaigns that overlap in time the neutron monitoring for five stations, we show that fast neutrons are much more strongly attenuated in snow than predicted by a conventional mass-shielding formulation and the attenuation length estimated in the atmosphere. We suggest that such strong attenuation results from boundary effects at the atmosphere/snow interface induced by the high efficiency of water as a neutron moderator. Finally, we propose an empirical model that allows calculating snow-shielding correction factors as a function of SWE for studies using Ne-21 and Cl-36 analyses in Mg- and K-rich minerals, respectively. This empirical model is of interest for studies with a focus on cosmic-ray exposure dating, particularly if the target rocks are made up of mafic to ultramafic units where seasonal snow-cover is a common phenomenon.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Age and chemical composition of apatites from five DSDP and ODP Sites in the Pacific, Atlantic and Indian Oceans

Since studies on deep-sea cores were carried out in the early 1990s it has been known that ambient temperature may have a marked affect on apatite fission track annealing. Due to sluggish annealing kinetics, this effect cannot be quantified by laboratory annealing experiments. The unknown amount of low-temperature annealing remains one of the main uncertainties for extracting thermal histories from fission track data, particularly for samples which experienced slow cooling in shallow crustal levels. To further elucidate these uncertainties, we studied volcanogenic sediments from five deep-sea drill cores, that were exposed to maximum temperatures between ~10° and 70°C over geological time scales of ~15-120 Ma. Mean track lengths (MTL) and etch pit diameters (Dpar) of all samples were measured, and the chemical composition of each grain analyzed for age and track length measurements was determined by electron microprobe analysis. Thermal histories of the sampled sites were independently reconstructed, based on vitrinite reflectance measurements and/or 1D numerical modelling. These reconstructions were used to test the most widely used annealing models for their ability to predict low-temperature annealing. Our results show that long-term exposure to temperatures below the temperature range of the nominal apatite fission track partial annealing zone results in track shortening ranging between 4 and 11%. Both chlorine content and Dpar values explain the downhole annealing patterns equally well. Low chlorine apatite from one drill core revealed a systematic relation between Si-content and Dpar value. The question whether Si-substitution in apatite has direct and systematic effects on annealing properties however, cannot be addressed by our data. For samples, which remained at temperatures <30°C, and which are low in chlorine, the Laslett et al. [Laslett G., Green P., Duddy I. and Gleadow A. (1987) Thermal annealing of fission tracks in apatite. Chem. Geol. 65, 1-13] annealing model predicts MTL up to 0.6 µm longer than those actually measured, whereas for apatites with intermediate to high chlorine content, which experienced temperatures >30°C, the predictions of the Laslett et al. (1987) model agree with the measured MTL data within error levels. With few exceptions, predictions by the Ketcham et al. [Ketcham R., Donelick R. and Carlson W. (1999) Variability of apatite fission-track annealing kinetics. III: Extrapolation to geological time scales. Am. Mineral. 84/9, 1235-1255] annealing model are consistent with the measured data for samples which remained at temperatures below ~30°C. For samples which experienced maximum temperatures between ~30 and 70°C, and which are rich in chlorine, the Ketcham et al. (1999) model overestimates track annealing.

Fly-through movie for the reconstructed folds in the onset region of the Petermann Glacier, North Greenland

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier.

(Table 1) Moisture content of buried snowpack deposit samples from Pearse Valley, McMurdo Dry Valleys

Buried snowpack deposits are found within the McMurdo Dry Valleys of Antarctica, which offers the opportunity to study these layered structures of sand and ice within a polar desert environment. Four discrete buried snowpacks are studied within Pearse Valley, Antarctica, through in situ observations, sample analyses, O-H isotope measurements and numerical modelling of snowpack stability and evolution. The buried snowpack deposits evolve throughout the year and undergo deposition, melt, refreeze, and sublimation. We demonstrate how the deposition and subsequent burial of snow can preserve the snowpacks in the Dry Valleys. The modelled lifetimes of the buried snowpacks are dependent upon subsurface stratigraphy but are typically less than one year if the lag thickness is less than c. 7 cm and snow thickness is less than c. 10 cm, indicating that some of the Antarctic buried snowpacks form annually. Buried snowpacks in the Antarctic polar desert may serve as analogues for similar deposits on Mars and may be applicable to observations of the north polar erg, buried ice at the Mars Phoenix landing site, and observations of buried ice throughout the martian Arctic. Numerical modelling suggests that seasonal snows and subsequent burial are not required to preserve the snow and ice on Mars.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.

(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ~34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved under¬standing of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Dynamic Effects On Culverts for High Speed Trains

Culverts are very common in recent railway lines. Wild life corridors and drainage conducts often fall in this category of partially buried structures. Their dynamic behavior has received far less attention than other structures such as bridges but its large number makes that study an interesting challenge from the point of view of safety and savings. In this paper a complete study of a culvert, including on-site measurements as well as numerical modelling, will be presented. The structure belongs to the high speed railway line linking Segovia and Valladolid, in Spain. The line was opened to traffic in 2004. Its dimensions (3x3m) are the most frequent along the line. Other factors such as reduced overburden (0.6m) and an almost right angle with the track axis make it an interesting example to extract generalized conclusions. On site measurements have been performed in the structure recording the dynamic response at selected points of the structure during the passage of high speed trains at speeds ranging between 200 and 300km/h. The measurements by themselves provide a good insight into the main features of the dynamic behaviour of the structure. A 3D finite element model of the structure, representing its key features was also studied as it allows further understanding of the dynamic response to the train loads . In the paper the discrepancies between predicted and measured vibration levels will be analyzed and some advices on numerical modelling will be proposed

Study of bond stress-slip relationship and radial dilation in prestressed concrete

This paper presents two test procedures for evaluating the bond stress–slip and the slip–radial dilation relationships when the prestressing force is transmitted by releasing the steel (wire or strand) in precast prestressed elements. The bond stress–slip relationship is obtained with short length specimens, to guarantee uniform bond stress, for three depths of the wire indentation (shallow, medium and deep). An analytical model for bond stress–slip relationship is proposed and compared with the experimental results. The model is also compared with the experimental results of other researchers. Since numerical models for studying bond-splitting problems in prestressed concrete require experimental data about dilatancy angle (radial dilation), a test procedure is proposed to evaluate these parameters. The obtained values of the radial dilation are compared with the prior estimated by numerical modelling and good agreement is reached

Some minor problems with B.E.M.

As it is well known B.E.M. works efficiently in the treatment of a bread class of potential and elasticity problems. In this paper we present the results of several runs established with linear elements in plane potential theory and treating the importance of singularities and the pattern and size of elements used in the boundary discretization.

Deformabilidad de escolleras empleadas en banquetas de cimentación de cajones portuarios

En esta tesis doctoral se presenta una investigación sobre el comportamiento deformacional de las escolleras empleadas en banquetas de cimentación de obras portuarias de cajones. El trabajo aborda el estudio de la deformabilidad de escolleras portuarias combinando (i) investigación mediante ensayos de laboratorio; (ii) análisis del comportamiento in situ de las banquetas de escolleras y (iii) cálculos realizados con modelos numéricos. Se expone en primer lugar la investigación experimental realizada en el Laboratorio de Geotecnia del Centro de Estudios y Experimentación de Obras Públicas (CEDEX) para estudiar la deformabilidad de las escolleras mediante ensayos a gran escala, habida cuenta las grandes dimensiones de las partículas de escollera. Se ha tratado de establecer una metodología de ensayo que reproduzca las solicitaciones de las escolleras colocadas en banquetas de cimentación de obras portuarias. Asimismo, se ha hecho una interpretación exhaustiva de los resultados de los ensayos con el fin de establecer unos valores que caractericen la deformabilidad de las escolleras analizadas. Es posible dar un intervalo de valores de la compresibilidad de las escolleras portuarias que, dada la escasez de literatura existente, constituyen unos valores de referencia. Asimismo, se ha propuesto una metodología para para estimar la deformabilidad de escolleras portuarias in situ. La información disponible ha permitido realizar estudios de la deformabilidad in situ en dos muelles españoles con semejanzas estructurales y constructivas. La interpretación conjunta de los resultados ha sugerido unos valores de deformabilidad in situ. Conviene destacar que la práctica ausencia de rangos de valores de compresibilidad in situ para estos rellenos empleados en obras portuarias pone de manifiesto la importancia de los resultados obtenidos. Evidencias de diferencias de comportamiento de las escolleras empleadas en banquetas de cimentaci ón de cajones portuarios en laboratorio e in situ han sido documentadas. La evaluación conjunta del comportamiento tenso-deformacional de las escolleras en laboratorio e in situ ha estimulado la búsqueda de una correlación entre la compresibilidad de las escolleras en ambos escenarios. Finalmente, se ha elaborado un modelo numérico con la formulación matemática del método sincrético (Perucho (2004, 2008)) que supone una opción interesante para evaluar la deformabilidad de los rellenos granulares. En la práctica, el empleo del modelo sincrético requiere la determinación de unos microparámetros. La disponibilidad de numerosos resultados de laboratorio realizados en las escolleras portuarias ha permitido calibrar el modelo realizado. De esta manera, se dispone de una herramienta de cálculo para evaluar la deformabilidad de los relleno granulares con un método numérico. The focus of this Thesis is to explore the deformational behavior of large rock fill materials used as rock mattress foundations for gravity caissons structures. The determination of the compressibility of large granular media focuses on (i) laboratory testing, (ii) in situ performance analysis of rock mattress foundations for caissons, and (iii) numerical modelling. First, the results of the large-scale laboratory research program, conducted at the Geotechnical Laboratory for the Center for Studies and Experimentation for Public Works (CEDEX), to determine the deformability of large rock fill materials is presented. The testing procedure was specifically designed to reproduce the loading sequence of in situ rubble mound foundations. A thoughtful analysis of the laboratory testing results suggests a range of compressibility for large granular media. The lack of currently available information regarding large rock fill deformability places a certain emphasis on the results of the testing program. Second, the results of this research includes a procedure for evaluating in situ rock fill deformational behavior. Data, collected from monitoring two caisson-type quays in Spain, provides information to study in situ rock mattress foundations. Careful interpretation of in situ data reveals a range of deformability of rock mattress foundations in caisson-type quays. Based upon a review of available literature, assessments on the behavior of rock mattress foundations for caissons using in situ analysis are quite limited. The data from this research are likely to contribute to the knowledge of the in situ behavior of rock mattress foundations for caissons. Additionally, findings indicate an appreciable variation between the laboratory and the in situ behaviour of materials from rock mattress foundations for caissons. Dissimilarities between laboratory and in situ moduli of deformation are examined in detail. Correlations between laboratory and in situ values are made. Finally, numerical modeling, based upon the research of Perucho (2004, 2008), is presented to predict the deformation behavior of large granular media. The determination of microparameters that control macropropierties requires extensive calibration effort. The calibration process was carried out using the results of large-scale laboratory testing available from previous analysis. The presented numerical method is both versatile and attractive as it reasonably predicts the compressibility of large rock fill materials.

Vibro-acoustic analysis of spacecraft structures with thin air layers

Esta Tesis presenta un estudio sobre el comportamiento vibroacústico de estructuras espaciales que incluyen capas de aire delgadas, así como sobre su modelización numérica. Las capas de aire pueden constituir un elemento fundamental en estos sistemas, como paneles solares plegados, que se consideran el caso de estudio en este trabajo. Para evaluar la influencia de las capas de aire en la respuesta dinámica del sistema se presenta el uso de modelos unidimensionales. La modelización de estos sistemas se estudia para los rangos de baja y alta frecuencia. En el rango de baja frecuencia se propone un conjunto de estrategias de simulación basadas en técnicas numéricas que se utilizan habitualmente en la industria aeroespacial para facilitar la aplicación de los resultados de la Tesis en los modelos numéricos actuales. Los resultados muestran el importante papel de las capas de aire en la respuesta del sistema. El uso de modelos basados en elementos finitos o de contorno para estos elementos proporciona resultados equivalentes aunque la aplicabilidad de estos últimos puede estar condicionada por la geometría del problema. Se estudia asimismo el uso del Análisis Estadístico de la Energía (SEA) para estos elementos. Una de las estrategias de simulación propuestas, que incluye una formulación energética para el aire que rodea a la estructura, se propone como estimador preliminar de la respuesta del sistema y sus frecuencias propias. Para el rango de alta frecuencia, se estudia la influencia de la definición del propio modelo SEA. Se presenta el uso de técnicas de reducción para determinar una matriz de pérdidas SEA reducida para definiciones incompletas del sistema (si algún elemento que interactúa con el resto no se incluye en el modelo). Esta nueva matriz tiene en cuenta la contribución de las subestructuras que no se consideran parte del modelo y que suelen ignorarse en el procedimiento habitual para reducir el tamaño del mismo. Esta matriz permite también analizar sistemas que incluyen algún componente con problemas de accesibilidad para medir su respuesta. Respecto a la determinación de los factores de pérdidas del sistema, se presenta una metodología que permite abordar casos en los que el método usual, el Método de Inyección de Potencia (PIM), no puede usarse. Se presenta un conjunto de métodos basados en la técnicas de optimización y de actualización de modelos para casos en los que no se puede medir la respuesta de todos los elementos del sistema y también para casos en los que no todos los elementos pueden ser excitados, abarcando un conjunto de casos más amplio que el abordable con el PIM. Para ambos rangos de frecuencia se presentan diferentes casos de análisis: modelos numéricos para validar los métodos propuestos y un panel solar plegado como caso experimental que pone de manifiesto la aplicación práctica de los métodos presentados en la Tesis. ABSTRACT This Thesis presents an study on the vibro-acoustic behaviour of spacecraft structures with thin air layers and their numerical modelling. The air layers can play a key role in these systems as solar wings in folded configuration that constitute the study case for this Thesis. A method based on one-dimensional models is presented to assess the influence of the air layers in the dynamic response of the system. The modelling of such systems is studied for low and high frequency ranges. In the low frequency range a set of modelling strategies are proposed based on numerical techniques used in the industry to facilitate the application of the results in the current numerical models. Results show the active role of the air layers in the system response and their great level of influence. The modelling of these elements by means of Finite Elements (FE) and Boundary Elements (BE) provide equivalent results although the applicability of BE models can be conditioned by the geometry of the problem. The use of Statistical Energy Analysis (SEA) for these systems is also presented. Good results on the system response are found for models involving SEA beyond the usual applicability limit. A simulation strategy, involving energetic formulation for the surrounding fluid is proposed as fast preliminary approach for the system response and the coupled eigenfrequencies. For the high frequency range, the influence of the definition of the SEA model is presented. Reduction techniques are used to determine a Reduced SEA Loss Matrix if the system definition is not complete and some elements, which interact with the rest, are not included. This new matrix takes into account the contribution of the subsystems not considered that are neglected in the usual approach for decreasing the size of the model. It also allows the analysis of systems with accessibility restrictions on some element in order to measure its response. Regarding the determination of the loss factors of a system, a methodology is presented for cases in which the usual Power Injection Method (PIM) can not be applied. A set of methods are presented for cases in which not all the subsystem responses can be measured or not all the subsystems can be excited, as solar wings in folded configuration. These methods, based on error minimising and model updating techniques can be used to calculate the system loss factors in a set of cases wider than the PIM’s. For both frequency ranges, different test problems are analysed: Numerical models are studied to validate the methods proposed; an experimental case consisting in an actual solar wing is studied on both frequency ranges to highlight the industrial application of the new methods presented in the Thesis.

Determination of the JWL constants for ANFO and emulsion explosives from cylinder test data

The Jones-Wilkins-Lee (JWL) equation of state parameters for ANFO and emulsion-type explosives have been obtained from cylinder test expansion measurements. The calculation method comprises a new radial expansion function, with a non-zero initial velocity at the onset of the expansion in order to comply with a positive Gurney energy at unit relative volume, as the isentropic expansion from the CJ state predicts. The equations reflecting the CJ state conditions and the measured expansion energy were solved for the JWL parameters by a non-linear least squares scheme. The JWL parameters of thirteen ANFO and emulsion type explosives have been determined in this way from their cylinder test expansion data. The results were evaluated through numerical modelling of the tests with the LS-DYNA hydrocode; the expansion histories from the modelling were compared with the measured ones, and excellent agreement was found.

ROCK OUTCROPS IN THE CANADIAN SHIELD: AN INVESTIGATION OF CONTAMINANT TRANSPORT FROM SURFACE SOURCES IN FRACTURED ROCK AQUIFERS

Flow, recharge and transport dynamics in fractured rock aquifers with low lying rock outcrops is a largely unexplored area of study in hydrogeology. The purpose of this thesis is to examine these topics in an agricultural area in Eastern Ontario. The study consists of a regional scale groundwater quality study, an infiltration experiment that considers bacteria transport from the ground surface to a well, and a numerical modelling study that tests the parameters that affect surface infiltration of a tracer from a rock outcrop to a deeper horizontal fracture. In the water quality study, approximately 65% of the samples contained total coliform, 16% contained E. coli, and 1% contained nitrate-N at greater than 5 mg/L. Occurrence of E. coli increased when considering seasonality, where wells were drilled on rock outcrops, and for shallow well intervals. Nitrate-N did not occur above the Guidelines for Canadian Drinking Water Quality (Health Canada, 2012) of 10 mg/L. Rapid arrival times were observed in the infiltration study for both the microspheres (30 minutes) and a dye tracer (45 minutes) in a well approximately 6.0 m in horizontal and 2.8 m in vertical distance from the tracer source. Transport velocities were approximately 38.9 m/day for the dye tracer and 115.2 m/day for the colloidal tracer. Results of the model runs indicate that overburden can provide an effective protective layer to transport in fractures, that high groundwater velocities occur in larger fracture apertures and higher gradients dilute tracer concentrations, and that lower groundwater velocities occur with smaller fracture apertures and lower gradients result in elevated tracer concentrations. Lower rainfall rates, larger fracture apertures, early tracer time, larger gradients, and lower water levels maintained unsaturated conditions for longer time periods such that tracer transport was delayed until saturated conditions were attained. The overall heterogeneity of this aquifer environment creates a source water protection conundrum where the water quality is generally good, while transport can occur very quickly in proximity to rock outcrops and in areas with limited overburden.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 2: Characteristics and petrogenesis of the high-Ti latite and high-Ti and low-Ti voluminous quartz latite eruptives

As a result of their relative concentration towards the respective Atlantic margins, the silicic eruptives of the Parana (Brazil)-Etendeka large igneous province are disproportionately abundant in the Etendeka of Namibia. The NW Etendeka silicic units, dated at similar to132 Ma, occupy the upper stratigraphic levels of the volcanic sequences, restricted to the coastal zone, and comprise three latites and five quartz latites (QL). The large-volume Fria QL is the only low-Ti type. Its trace element and isotopic signatures indicate massive crustal input. The remaining NW Etendeka silicic units are enigmatic high-Ti types, geochemically different from low-Ti types. They exhibit chemical affinities with the temporally overlapping Khumib high-Ti basalt (see Ewart et al. Part 1) and high crystallization temperatures (greater than or equal to980 to 1120degreesC) inferred from augite and pigeonite phenocrysts, both consistent with their evolution from a mafic source. Geochemically, the high-Ti units define three groups, thought genetically related. We test whether these represent independent liquid lines of descent from a common high-Ti mafic parent. Although the recognition of latites reduces the apparent silica gap, difficulty is encountered in fractional crystallization models by the large volumes of two QL units. Numerical modelling does, however, support large-scale open-system fractional crystallization, assimilation of silicic to basaltic materials, and magma mixing, but cannot entirely exclude partial melting processes within the temporally active extensional environment. The fractional crystallization and mixing signatures add to the complexity of these enigmatic and controversial silicic magmas. The existence, however, of temporally and spatially overlapping high-Ti basalts is, in our view, not coincidental and the high-Ti character of the silicic magmas ultimately reflects a mantle signature.