U–Pb evidence for a polyorogenic evolution of the HP–HT units of the NW Iberian Massif

. A isotope dilution thermal ionisation mass spectrometry U-Pb geochronological study was carried out on the high-pressure and high-temperature units (HP-HT units) overlying the oceanic suture in the Allochthonous Complexes of the NW Iberian Variscan Belt. The rocks investigated are seven granulite- to eclogite-facies paragneisses and one leucosome within mafic high-pressure granulites in the Ordenes and Cabo Ortegal Complexes of NW Spain. U-Pb dating of zircon, monazite, titanite and rutile reveal the presence of a pervasive Early Ordovician metamorphic event at ca. 500-480 Ma and a later Early Devonian event at ca. 400-380 Ma. The U-Pb ages, in conjunction with petrological and structural data, indicate that the high-pressure event recorded by these rocks is Early Ordovician in age. Monazite ages in the paragneisses suggest that peak metamorphic conditions were reached at ca. 500-485 Ma. Subsequently, the rock ensemble underwent exhumation accompanied by partial melting and zircon growth at ca. 485-470 Ma. Melting of mafic granulites was coeval with this latter episode as indicated by zircon crystallisation age in the leucosomes dated at ca. 486 Ma. Based on these data and on the general features of magmatism and metamorphic evolution, it is proposed that this process took place at a convergent plate boundary within a peri-Gondwanan oceanic domain. Monazite, titanite and rutile data in some of the samples studied show evidence of a second metamorphic episode that took place between ca. 400 and 380 Ma (with a peak at ca. 390-385 Ma). This Early Devonian event, at variance with the previous one, was not pervasive, but, rather, was localised in areas of intense Variscan tectonothermal reworking. It is claimed that this later metamorphic event was recorded by the U-Pb system in areas where monazite and titanite growth was enhanced by fluid circulation in highly strained rocks (Variscan shear zones). According to previous structural studies and Ar-Ar dating of fabrics, this Early Devonian episode took place as the HP-HT units were deformed and thrusted upon the ophiolitic units in the early stages of the Variscan collision.     

Evolution of fluids associated with metasedimentary sequences from Chaves (North Portugal)

In order to identify and characterise fluids associated with metamorphic rocks from the Chaves region (North Portugal), fluid inclusions were studied in quartz veinlets, concordant with the main foliation, in graphitic-rich and nongraphitic-rich lithologies from areas with distinct metamorphic grade. The study indicates multiple fluid circulation events with a variety of compositions, broadly within the C–H–O–N–salt system. Primary fluid inclusions in quartz contain low salinity aqueous–carbonic, H₂O–CH₄–N₂–NaCl fluids that were trapped near the peak of regional metamorphism, which occurred during or immediately after D2. The calculated P–T conditions for the western area of Chaves (CW) is P=300–350 MPa and T500 °C, and for the eastern area (CE), P=200–250 MPa and T=400–450 °C. A first generation of secondary fluid inclusions is restricted to discrete cracks at the grain boundaries of quartz and consists of low salinity aqueous–carbonic, H₂O–CO₂–CH₄–N₂–NaCl fluids. P–T conditions from the fluid inclusions indicate that they were trapped during a thermal event, probably related with the emplacement of the two-mica granites.

A second generation of secondary inclusions occurs in intergranular fractures and is characterised by two types of aqueous inclusions. One type is a low salinity, H₂O–NaCl fluid and the second consists of a high salinity, H₂O–NaCl–CaCl₂ fluid. These fluid inclusions are not related to the metamorphic process and have been trapped after D3 at relatively low P (hydrostatic)–T conditions (P<100 MPa and T<300 °C).

Both the early H₂O–CH₄–N₂–NaCl fluids in quartz from the graphitic-rich lithologies and the later H₂O–CO₂–CH₄–N₂–NaCl carbonic fluid in quartz from graphitic-rich and nongraphitic-rich lithologies seem to have a common origin and evolution. They have low salinity, probably resulting from connate waters that were diluted by the water released from mineral dehydration during metamorphism. Their main component is water, but the early H₂O–CH₄–N₂–NaCl fluids are enriched in CH₄ due to interaction with the C-rich host rocks.

From the early H₂O–CH₄–N₂–NaCl to the later aqueous–carbonic H₂O–CO₂–CH₄–N₂–NaCl fluids, there is an enrichment in CO₂ that is more significant for the fluids associated with nongraphitic-rich lithologies.

The aqueous–carbonic fluids, enriched in H₂O and CH₄, are primarily associated with graphitic-rich lithologies. However, the aqueous–carbonic CO₂-rich fluids were found in both graphitic and nongraphitic-rich units from both the CW and CE studied areas, which are of medium and low metamorphic grade, respectively.     

Andesite and dacite genesis via contrasting processes: the geology and geochemistry of El Valle Volcano,Panama

The easternmost stratovolcano along the Central American arc is El Valle volcano, Panama. Several andesitic and dacitic lava flows, which range in age 5–10 Ma, are termed the old group. After a long period of quiescence (approximately 3.4 Ma), volcanic activity resumed approximately 1.55 Ma with the emplacement of dacitic domes and the deposition of dacitic pyroclastic flows 0.9–0.2 Ma. These are referred to as the young group. All of the samples analyzed are calc-alkaline andesites and dacites. The mineralogy of the two groups is distinct; two pyroxenes occur in the old-group rocks but are commonly absent in the young group. In contrast, amphibole has been found only in the young-group samples. Several disequilibrium features have been observed in the minerals (e.g., oscillatory zoning within clinopyroxenes). These disequilibrium textures appear to be more prevalent among the old- as compared with the young-group samples and are most likely the result of magma-mixing, assimilation, and/or polybaric crystallization. Mass-balance fractionation models for major and trace elements were successful in relating samples from the old group but failed to show a relationship among the young-group rocks or between the old- and young-group volcanics. We believe that the old-group volcanics were derived through differentiation processes from basaltic magmas generated within the mantlewedge. The young group, however, does not appear to be related to more primitive magmas by differentiation. The young-group samples cannot be related by fractionation including realistic amounts of amphibole. Distinctive geochemical features of the young group, including La/Yb ratios〉15, Yb〈1, Sr/Y〉150, and Y〈6, suggest that these rocks were derived from the partial melting of the subducted lithosphere. These characteristics can be explained by the partial melting of a source with residual garnet and amphibole. Dacitic material with the geochemical characteristics of subducted-lithosphere melting is generated apparently only where relatively hot crust is subducted, based on recent work. The young dacite-genesis at El Valle volcano is related to the subduction of relatively hot lithosphere.     

Rutile R632 – A New Natural Reference Material for U‐Pb and Zr Determination

A new natural rutile reference material is presented, suitable for U‐Pb dating and Zr‐in‐rutile thermometry by microbeam methods. U‐Pb dating of rutile R632 using laser ablation ICP‐MS with both magnetic sector field and quadrupole instruments as well as isotope dilution‐thermal ionisation mass spectrometry yielded a concordia age of 496 ± 2 Ma. The high U content (> 300 μg g^?1) enabled measurement of high‐precision U‐Pb ages despite its young age. The sample was found to have a Zr content of 4294 ± 196 μg g^?1, which makes it an excellent complementary reference material for Zr‐in‐rutile thermometry. Individual rutile grains have homogeneous compositions of a number of other trace elements including V, Cr, Fe, Nb, Mo, Sn, Sb, Hf, Ta and W. This newly characterised material significantly expands the range of available rutile reference materials relevant for age and temperature determinations.     

Genetic and chemical relationships between berthierine,chlorite and cordierite in nodules associated to granitic pegmatites of Sierra Albarrana (Iberian Massif,Spain)

Chlorite and berthierine occur through alteration of cordierite within enclaves of metamorphic rocks transformed by the Sierra Albarrana pegmatites. The coexistence of both phyllosilicates allows us to study their stability relationships and to compare their chemical compositions. Samples showing incipient replacement of cordierite by small cryptocrystalline aggregates can be identified by X-ray diffraction (XRD) as berthierine with small quantities of chlorite. Electron Microprobe (EMP) analyses give mixed compositions of berthierine and cordierite. Samples with extensive replacement of cordierite by aggregates show similar characteristics to those with incipient replacement, but some small crystals are present. The last type of sample shows complete replacement of cordierite by crystals showing optical properties of chlorite and EMPA compositions coherent with chlorite or berthierine. Their XRD pattern corresponds to chlorite and their high resolution transmission electron microscopy (HRTEM) images only show perfect sequences of 14 Å lattice fringes. The cryptocrystalline aggregates of the samples with incipient and extensive replacement present coexistent areas of 14 Å and 7 Å lattice fringes that are intergrown at different levels: (1) large areas (> 1 m) of 7 Å layers; (2) packets of 7 Å layers between 14 Å layer areas, with visible 7 Å to 14 Å lateral changes; (3) random mixed-layers 7 Å/14 Å. Chlorite is the final stable product of alteration of cordierite, with berthierine as an intermediate metastable phase. Energy dispersive X-ray spectrometry microanalyses of 14 Å, 7 Å and (14+7) Å areas show lack of systematic differences in chemical compositions between both phyllosilicates which may be considered as true polymorphs.     

Regional retrograde alteration of sub-greenschist facies chlorite to smectite

Dioctahedral smectite is present as a retrograde alteration product of chlorite in Permian-Triassic red slates of the Malaguide Complex in Sierra de Espuña (Betic Cordillera). Mineral assemblages and textures, illite crystallinity indices, and fluid inclusion data indicate sub-greenschist facies conditions that reached at least 180°C in the higher-grade tectonic unit of the Malaguide Complex, preceding formation of smectite. Smectite, having K as the dominant interlayer cation, occurs ubiquitously intercalated with trioctahedral chlorite as thin packets of layers and as individual layers that commonly change to chlorite along layers. Although some chlorite is typically homogeneous and trioctahedral, much chlorite shows signs of alteration and has compositions corresponding to different degrees of smectite contaimination. The incompatibility of metamorphic grade with the occurrence of smectite, the general association of chlorite and smectite, and the textural relations collectively show that dioctahedral smectite is derived through replacement of trioctahedral chlorite. Such replacement occurs on a regional basis and demonstrates that caution must be used in interpreting the occurrence of smectite in pelites as being due to prograde processes. Alteration of trioctahedral chlorite under oxidizing conditions due to introduction of phreatic water after uplift of the Betic Cordillera is proposed as the cause of formation of smectite.     

Hydrogeology and sustainable future groundwater abstraction from the Agua Verde aquifer in the Atacama Desert,northern Chile

The hyper-arid conditions prevailing in Agua Verde aquifer in northern Chile make this system the most important water source for nearby towns and mining industries. Due to the growing demand for water in this region, recharge is investigated along with the impact of intense pumping activity in this aquifer. A conceptual model of the hydrogeological system is developed and implemented into a two-dimensional groundwater-flow numerical model. To assess the impact of climate change and groundwater extraction, several scenarios are simulated considering variations in both aquifer recharge and withdrawals. The estimated average groundwater lateral recharge from Precordillera (pre-mountain range) is about 4,482 m³/day. The scenarios that consider an increase of water withdrawal show a non-sustainable groundwater consumption leading to an over-exploitation of the resource, because the outflows surpasses inflows, causing storage depletion. The greater the depletion, the larger the impact of recharge reduction caused by the considered future climate change. This result indicates that the combined effects of such factors may have a severe impact on groundwater availability as found in other groundwater-dependent regions located in arid environments. Furthermore, the scenarios that consider a reduction of the extraction flow rate show that it may be possible to partially alleviate the damage already caused to the aquifer by the continuous extractions since 1974, and it can partially counteract climate change impacts on future groundwater availability caused by a decrease in precipitation (and so in recharge), if the desalination plant in Taltal increases its capacity.     

Climatology of extratropical cyclones over the South American–southern oceans sector

A climatology of extratropical cyclones is presented. Extratropical cyclones, their main characteristics and their predominant tracks, as well as their interannual variability, affect weather in South America. For that purpose, a storm track database has been compiled by applying a cyclone tracking scheme to six-hourly sea level pressure fields, available from the National Center for Environmental Prediction–National Center for Atmospheric Research reanalyses II for the 1979–2003 period. The spatial distribution of the cyclogenesis frequency shows two main centers: one around Northern Argentina, Uruguay, and Southern Brazil in all seasons and the other near to the North Antarctic Peninsula. The lifetime of extratropical cyclones in the South American sector exhibits small seasonality, being typically of the order of 3.0 days during most of the year and slightly higher (3.5 days) in austral summer. The distance travelled by the cyclones formed in the South American sector tends to be smaller than the total paths found in other areas of the Southern Hemisphere. A k-mean clustering technique is used to summarize the analysis of the 25-year climatology of cyclone tracks. Three clusters were found: one storm-track cluster in Northeast Argentina; a second one west of the Andes Cordillera; and a third cluster located to the north of the Antarctic Peninsula (around the Weddell Sea). The influence of the Antarctic Oscillation (AAO) in the variability of extratropical cyclones is explored, and some signals of the impacts of the variability of the AAO can be observed in the position of the extratropical cyclones around 40°S, while the impacts on the intensity is detected around 55°S.     

A copula-based multivariate analysis of Canadian RCM projected changes to flood characteristics for northeastern Canada

In the present work, climate change impacts on three spring (March–June) flood characteristics, i.e. peak, volume and duration, for 21 northeast Canadian basins are evaluated, based on Canadian regional climate model (CRCM) simulations. Conventional univariate frequency analysis for each flood characteristic and copula based bivariate frequency analysis for mutually correlated pairs of flood characteristics (i.e. peak–volume, peak–duration and volume–duration) are carried out. While univariate analysis is focused on return levels of selected return periods (5-, 20- and 50-year), the bivariate analysis is focused on the joint occurrence probabilities P1 and P2 of the three pairs of flood characteristics, where P1 is the probability of any one characteristic in a pair exceeding its threshold and P2 is the probability of both characteristics in a pair exceeding their respective thresholds at the same time. The performance of CRCM is assessed by comparing ERA40 (the European Centre for Medium-Range Weather Forecasts 40-year reanalysis) driven CRCM simulated flood statistics and univariate and bivariate frequency analysis results for the current 1970–1999 period with those observed at selected 16 gauging stations for the same time period. The Generalized Extreme Value distribution is selected as the marginal distribution for flood characteristics and the Clayton copula for developing bivariate distribution functions. The CRCM performs well in simulating mean, standard deviation, and 5-, 20- and 50-year return levels of flood characteristics. The joint occurrence probabilities are also simulated well by the CRCM. A five-member ensemble of the CRCM simulated streamflow for the current (1970–1999) and future (2041–2070) periods, driven by five different members of a Canadian Global Climate Model ensemble, are used in the assessment of projected changes, where future simulations correspond to A2 scenario. The results of projected changes, in general, indicate increases in the marginal values, i.e. return levels of flood characteristics, and the joint occurrence probabilities P1 and P2. It is found that the future marginal values of flood characteristics and P1 and P2 values corresponding to longer return periods will be affected more by anthropogenic climate change than those corresponding to shorter return periods but the former ones are subjected to higher uncertainties.     

Modelling of Ocean Shallow Convection under Fast Ice in the Arctic

Ocean convection in the Antarctic has been studied many times and has been revealed to be responsible for ice-cover reduction. In the Arctic, proof of that phenomenon has not been documented. It is believed that this phenomenon happens on a smaller scale in the Arctic when local circulation of deep warmer water melts and slows ice production. An example of this is the North Water (NOW) polynya in northern Baffin Bay. A polynya is an area of open water in an otherwise ice-covered area. As ice forms under the fast ice near the boundary of the polynya, ocean salts (brine) are ejected from the newly formed ice. This water, which has an increased concentration of salt, sinks and is replaced by warmer water from below, and this slows ice formation. In our study a coupled one-dimensional thermodynamic snow–fast ice model incorporating ocean heat flux input via a shallow convection model was used. Ice thickness was calculated using a thermodynamic model that included a current-induced entrainment model and a convection model to account for brine rejection during ice growth. Atmospheric observations from Grise Fiord and Thule and ocean profiles around the NOW polynya near these sites were used as input to the model. This purely thermodynamic study enables us to obtain ice thickness values that can be compared with qualitative observations. This modelling study compares two sites related to the NOW polynya. The results indicate that the shallow convection model simulates the reduction of fast ice near Thule but not near Grise Fiord.