Consistency of the activity–composition models of Holland,Green, and Powell (2018) with experiments on natural and synthetic compositions: A comparative study

The recent publication of new activity–composition models by Holland, Green and Powell (2018; Journal of Petrology 59 : 881–900), with a melt model calibrated for source compositions ranging from peridotitic to granitic, opens the door to the modelling of multiple petrogenetic processes at supersolidus conditions in which the composition of the melt phase changes considerably, without having to change the melt model. This melt model is also the first one using the internally consistent thermodynamic databases published by T. Holland and R. Powell that contains TiO₂ and Fe₂O₃, further expanding the application of this model to more realistic geological scenarios. The accompanying mineral models are also the first in containing some minor elements, like TiO₂ in garnet and K₂O in clinopyroxene. Consequently, it is relevant to test the applicability of these new models to a large P–T–X range of conditions before they can be used in full. Thermodynamic calculations made with the software Perple_X using these models were compared to experimental results, namely the modal proportions and the composition of the melt and several mineral solution phases. The experiments chosen for the comparative study covered a wide range of source compositions (from mafic to felsic), pressure (from 0.3 to 2.1 GPa), temperature (from 700 to 1,150°C) and total and added water content (structural water: 0.15–1.48 wt%; added water: 0–8 wt%; total water: 0.15%–8.15%). The results indicate that the extended melt model reproduces well the composition of the experimental melts, with an inverse correlation between component amount and fit: the best match is found for SiO₂ (−0.8% on average) and the worst match is found for those elements with the lowest amounts, TiO₂ and MgO (+241% and +235% on average, respectively; values indicate calculated minus experimental, times 100 and divided by experimental). The TiO₂ content in the melt model increases dramatically with increasing pressure, from +90% for P < 1.5 GPa to +593% for P > 1.5 GPa. No comparison was made on the Fe₂O₃ content, as the published iron contents of the experimental melts were always reported as FeO_t. In some cases, there is a substantial mismatch in the modal proportions between experiments and calculations, with the reactant phases less abundant and product phases more abundant in the calculations, an effect that is attributed to kinetic effects in the experiments and to the selected clinoamphibole model. Finally, the extended melt model was compared to the tonalitic melt model of Green et al. (2016; Journal of Metamorphic Geology, 34 : 845–869). Both melt models produce very similar results for SiO₂, Al₂O₃, Na₂O and K₂O, with slightly better results for the tonalitic melt model in FeO_t and MgO and for the extended melt model in CaO. No comparison is made on TiO₂ because the tonalitic melt model does not include this component. In summary, the new activity–composition models represent a significant contribution to thermodynamic calculations on the evolution of siliceous magmas where their composition, temperature and pressure changes substantially.     

Pliocene lahar deposits in the Coastal Cordillera of central Chile: Implications for uplift, avalanche deposits, and porphyry copper systems in the Main Andean Cordillera

Lahar deposits occur within a shallow marine sedimentary succession of the Pliocene La Cueva Formation in the Coastal Cordillera of central Chile (33°40′–34°15′S). Provenance studies of the abundant volcanic material in the lahar deposits suggest that they derive from denudation by mass wasting of Oligocene–Miocene volcanic rocks on the western slopes of the Main Andean Cordillera at the same latitude. Pliocene rock debris deposits preserved in the region of El Teniente (34°S) and scattered along the westernmost part of the Andes of central Chile indicate catastrophic erosive events related to the rapid uplift of the cordilleran block. This rock debris was deposited by avalanches and transformed further downslope into lahars by dilution with stream water. Lahars were channeled along the ancient drainage system that reached a shallow Pliocene sea at the site of the present Coastal Cordillera. The exceedingly rapid exhumation of active porphyry systems during the Early Pliocene in this part of the Andes may have played a role in affecting hydrothermal processes, brecciation, and diatreme formation at the porphyry systems of El Teniente and Río Blanco–Los Bronces.     

The high hydraulic conductivity of three wooded tropical peat swamps in northeast Peru: measurements and implications for hydrological function

The form and functioning of peatlands depend strongly on their hydrological status, but there are few data available on the hydraulic properties of tropical peatlands. In particular, the saturated hydraulic conductivity (K) has not previously been measured in neotropical peatlands. Piezometer slug tests were used to measure K at two depths (50 and 90 cm) in three contrasting forested peatlands in the Peruvian Amazon: Quistococha, San Jorge and Buena Vista. Measured K at 50 cm depth varies between 0.00032 and 0.11 cm s^?1, and at 90 cm, it varies between 0.00027 and 0.057 cm s^?1. Measurements of K taken from different areas of Quistococha showed that spatial heterogeneity accounts for ~20% of the within‐site variance and that depth is a good predictor of K. However, K did not vary significantly with depth at Buena Vista and San Jorge. Statistical analysis showed that ~18% of the variance in the K data can be explained by between‐site differences. Simulations using a simple hydrological model suggest that the relatively high K values could lead to lowering of the water table by >10 cm within ~48 m of the peatland edge for domed peatlands, if subjected to a drought lasting 30 days. However, under current climatic conditions, even with high K, peatlands would be unable to shed the large amount of water entering the system via rainfall through subsurface flow alone. We conclude that most of the water leaves these peatlands via overland flow and/or evapotranspiration. Copyright © 2013 John Wiley & Sons, Ltd.     

Groundwater salinization in the Azores archipelago (Portugal)

Groundwater salinization in coastal regions causes severe constraints to water supply and economic losses to society worldwide. In the Azores archipelago, groundwater abstraction in wells drilled in coastal aquifers is very important for water supply, and quality problems have been reported. Therefore, a groundwater chemistry dataset from wells was compiled to study groundwater salinization in these aquifers. Waters are mainly of the Na–Cl type, presenting a slightly acidic to slightly alkaline character, with a pH between 5.63 and 8.50 (median 7.40). Electrical conductivity measurements range from 127 to 9,670 μS/cm (median 862), suggesting highly variable mineralization, with higher values observed on Santa Maria, São Miguel, Pico, Graciosa, and São Jorge islands. The major-ion composition reflects the contribution of seawater to the groundwater compositional evolution, which is essentially explained by seawater intrusion into wells. In many samples, exchange reactions of Na⁺ + K⁺ for Ca²⁺ + Mg²⁺ are associated with salinization. The seawater fraction in groundwater composition reaches a maximum of 22.5%. These results are extremely challenging to water managers in the Azores because failure to comply with national water quality regulations and with European Union groundwater directive requirements often occurs.     

Changes in benthic fish assemblages as a consequence of coastal works in a coastal lagoon: The Mar Menor (Spain, Western Mediterranean)

The benthic fish assemblage of the Mar Menor consisted of 37 species. Dominant species are: Gobius cobitis, Lipophrys pavo and Tripterygion tripteronotus on infralittoral rocks; Pomatoschistus marmoratus, Callionymus pussillus, Callionymus risso and Solea vulgaris on sandy bottoms and Gobius niger, Syngnathus abaster, Hippocampus ramulosus and Symphodus cinereus on Cymodocea nodosa-Caulerpa prolifera mixed beds. From 1985 to 1989 tourist development has led to the creation of new beaches and the installation of artificial rocky structures for retaining sediments. Dredging for the extraction of sand and subsequent pumping altered sediment characteristics causing a real stress leading to the substitution of typical sandy bottoms communities with Cymodocea nodosa by Caulerpa prolifera communities on mud. Soft bottom fish assemblages responded to changes in vegetation cover and substratum characteristics mainly changing the species composition, while artificial hard substrata contain a similar fish community than natural ones, harbouring even richer and more diverse assemblages. This positive effect of breakwaters should not obscure their likely negative effects on hydrodynamics and the subsequent changes of sediment quality and vegetation cover on the breakwaters' area of influence.     

Assessing rheological properties of fluid mud samples through tuning fork data

Ocean Dynamics - At the nautical bottom approach, part of the fluid mud layers can be included in the available depth if they present favorable rheology. As it is difficult to perform in situ...     

A hybrid fast algorithm for first arrivals tomography

A hybrid algorithm, combining Monte-Carlo optimization with simultaneous iterative reconstructive technique (SIRT) tomography, is used to invert first arrival traveltimes from seismic data for building a velocity model. Stochastic algorithms may localize a point around the global minimum of the misfit function but are not suitable for identifying the precise solution. On the other hand, a tomographic model reconstruction, based on a local linearization, will only be successful if an initial model already close to the best solution is available. To overcome these problems, in the method proposed here, a first model obtained using a classical Monte Carlo-based optimization is used as a good initial guess for starting the local search with the SIRT tomographic reconstruction. In the forward problem, the first-break times are calculated by solving the eikonal equation through a velocity model with a fast finite-difference method instead of the traditional slow ray-tracing technique. In addition, for the SIRT tomography the seismic energy from sources to receivers is propagated by applying a fast Fresnel volume approach which when combined with turning rays can handle models with both positive and negative velocity gradients. The performance of this two-step optimization scheme has been tested on synthetic and field data for building a geologically plausible velocity model.This is an efficient and fast search mechanism, which permits insertion of geophysical, geological and geodynamic a priori constraints into the grid model and ray path is completed avoided. Extension of the technique to 3D data and also to the solution of 'static correction' problems is easily feasible.     

Summertime moisture transport over Southeastern South America and extratropical cyclones behavior during inter-El Niño events

The impact of the inter-El Niño (EN) variability on the moisture availability over Southeastern South America (SESA) is investigated. Also, an automatic tracking scheme was used to analyze the extratropical cyclones properties (system density - SD and central pressure - CP) in this region. During the austral summer period from 1977–2000, the differences for the upper-level wave train anomaly composites seem to determine the rainfall composite differences. In fact, the positive rainfall anomalies over most of the SESA domain during the strong EN events are explained by an upper-level cyclonic center over the tropics and an anticyclonic center over the eastern subtropical area. This pattern seems to contribute to upward vertical motion at 500 hPa and reinforcement of the meridional moisture transport from the equatorial Atlantic Ocean and western Amazon basin to the SESA region. These features may contribute to the positive SD and negative CP anomalies explaining part of the positive rainfall anomalies found there. On the other hand, negative rainfall anomalies are located in the northern part of SESA for the weak EN years when compared to those for the strong events. Also, positive anomalies are found in the southern part, albeit less intense. It was associated with the weakening of the meridional moisture transport from the tropics to the SESA that seems have to contributed with smaller SD and CP anomalies over the most part of subtropics, when compared to the strong EN years.