Large-scale silicate liquid immiscibility: a possible example from northeastern Brazil

Ultrapotassic, peralkalic silica-saturated plutons (580 Ma) are widespread in the Cachoeirinha-Salgueiro foldbelt, northeastern Brazil. They consist of alkali-feldspar syenites with pyroxenite as co-magmatic inclusions and syn-plutonic or late-stage dikes. Pyroxenite and syenite have the same mineral phases (aegirine-augite, microcline, sphene, apatite, blue amphibole, magnetite), but only in different proportions. Rare inclusions of a “mixed” rock (about 60% syenite+40% pyroxenite in an emulsion-like texture) are also present. Pyroxenes in the three units are all only slightly zoned, silica-saturated and extremely low in Al₂O₃ (0.2–1.4%). Amphiboles are mostly K-rich richterite, characterized by high SiO₂, low Al₂O₃ and TiO₂ contents and low Mg#.

The three rock types have similar REE chondrite-normalized patterns, with negative slopes and lack of Eu anomaly, with the total REE in the pyroxenite greater than that of the syenite. Trace element patterns for the mixed rock are intermediate between those for the pyroxenite and syenite. Major element partitioning between pyroxenite and syenite has the same sense as that one observed between immiscible liquids in volcanic lavas and trace element partitioning is similar to the experimentally determined partition of immiscible liquid pairs.

The rocks have similar high δ¹⁸O values (avg. w.r.+ 8‰_SMOW, corrected from pyroxene), high initial ⁸⁷Sr/⁸⁶Sr ratios (about 0.710), and low ¹⁴³Nd/¹⁴⁴Nd (avg. 0.51104).

Field and geochemical characteristics indicate chemical equilibrium among the three rock types and suggest liquid immiscibility between syenite and pyroxenite, the mixed rock representing the original magma composition.     

Columnar joints in the Patino Formation sandstones, Eastern Paraguay： a dynamic interaction between dyke intrusion, quartz dissolution and cooling-induced fractures

The Patiho Formation sandstones, which crop out in Aregud neighborhood in Eastern Paraguay and show columnar joints near the contact zone with a nephelinite dyke, have as their main characteristics the high pro- portion of syntaxial quartz overgrowth and a porosity originated from different processes, initially by dissolu- tion and later by partial filling and fracturing. Features like the presence of floating grains in the syntaxial cement, the transitional interpenetrative contact between the silica-rich cement and grains as well as the intense fracture porosity are strong indications that the cement has been formed by dissolution and reprecipita- tion of quartz from the framework under the effect of thermal expansion followed by rapid contraction. The increase of the silica-rich cement towards the dyke in association with the orthogonal disposition of the columns relative to dyke walls are indicative that the igneous body may represent the main heat source for the interstitial aqueous solutions previously existing in the sediments. At macroscopic scale, the increasing of inter- nal tensions in the sandstones is responsible for the nucleation of polygons, leading to the individualization of prisms, which are interconnected by a system of joints, formed firstly on isotherm surfaces of low tem- perature and later on successive adjacent planes towards the dyke heat source.     

Rock magnetic analyses as a tool to investigate diversity of drift pumice clasts: An example from Japan's 2021 Fukutoku-Oka-no-Ba (FOB) eruption

The crystallization of nanolites within magma chambers has recently raised a strong interest due to their impact on increasing melt viscosity and triggering magmatic eruptions. In 2021, the Fukutoku-Oka-no-Ba (FOB) underwater eruption produced large quantities of pumices that eventually formed rafts drifting at the surface of the ocean to the East coasts of Japan. Pumices collected along the shore shortly after grounding show various colors, microscopic and Raman analyses made by Yoshida et al. (Island Arc, 31, 1, 2022) revealed the presence of magnetite nanolites in some of them. In this study, we explore the magnetic properties of a batch of pumices of different colors from the FOB eruption, aiming to refine characterization of iron oxide nanolites. We used various analytical techniques such as SEM and FEG-SEM observations, EDS-X analyses, and rock magnetic experiments, including thermomagnetic analyses, hysteresis curves, coercivity analyses and FORC measurements. Our findings reveal that the iron oxides present in the FOB samples are Ti-magnetite, with minor amounts of Mg and Al. The magnetic crystals show a wide range of sizes, from extra small iron oxide nanolites (ESION) in the pumices with the lighter colors, to more bulky grains reaching the micrometer size in some of the dark color samples, significant diffusion is inferred in that case. Consequently, the magnetic characterization of iron oxide crystals within the Fukutoku-Oka-no-Ba pumices reveals varying stages of nucleation, dissolution, growth, and diffusion processes, providing evidences for the heterogeneous state of the magma during the eruption.     

Environmental issues in urban groundwater systems: a multidisciplinary study of the Paranhos and Salgueiros spring waters,Porto (NW Portugal)

This study presents the results of a multidisciplinary approach, using hydrogeochemical, isotopic and ecotoxicological analyses, performed to assess the nature and suitability for use of Paranhos and Salgueiros spring waters (Porto city, NW Portugal). Based on the surface activities located along the course of the springs, 23 water samples were collected. All the samples were analysed for major element concentrations. The isotopic techniques employed included δ²H, δ¹⁸O and ³H. Standard acute bioassays with Daphnia magna were also performed. The hydrogeochemical analyses showed a nitrate and sulphate-enriched composition for these groundwaters, resulting mainly from urban drainage and sewer leakage. In the ecotoxicological analyses, no significant mortality was observed in any of the tests performed. The results obtained in this study suggest that Porto urban groundwater could be suitable for irrigation uses.     

ERORUN-STAFOR: A collaborative observatory for the multidisciplinary study of the critical zone processes in a tropical volcanic watershed including a Tropical Montane Cloud Forest

Tropical volcanic islands are biodiversity hotspots where the Critical Zone (CZ) still remains poorly studied. In such steep topographic environments associated with extreme climatic events (cyclones), deployment and maintenance of monitoring equipment is highly challenging. While a few Critical Zone Observatories (CZOS) are located in tropical volcanic regions, none of them includes a Tropical Montane Cloud Forest (TMCF) at the watershed scale. We present here the dataset of the first observatory from the French network of critical zone observatories (OZCAR) located in an insular tropical and volcanic context, integrating a ‘Tropical Montane Cloud Forest’: The ERORUN-STAFOR observatory. This collaborative observatory is located in the northern part of La Réunion island (Indian Ocean) within the 45.0 km² watershed of Rivière des Pluies (i.e., Rainfall river) which hosts the TMCF of Plaines des Fougères, one of the best preserved natural habitats in La Réunion Island. Since 2014, the ERORUN-STAFOR monitoring in collaboration with local partners collected a multidisciplinary dataset with a constant improvement of the instrumentation over time. At the watershed scale and in its vicinity, the ERORUN-STAFOR observatory includes 10 measurement stations covering the upstream, midstream and downstream part of the watershed. The stations record a total of 48 different variables through continuous (sensors) or periodic (sampling) monitoring. The dataset consists of continuous time series variables related to (i) meteorology, including precipitation, air temperature, relative humidity, wind speed and direction, net radiation, atmospheric pressure, cloud water flux, irradiance, leaf wetness and soil temperature, (ii) hydrology, including water level and temperature, discharge and electrical conductivity (EC) of stream, (iii) hydrogeology, including (ground)water level, water temperature and EC in two piezometers and one horizontally drilled groundwater gallery completed by soil moisture measurements under the canopy. The dataset is completed by periodic time series variables related to (iv) hydrogeochemistry, including field parameters and water analysis results. The periodic sampling survey provides chemical and isotopic compositions of rainfall, groundwater, and stream water at different locations of this watershed. The ERORUN-STAFOR monitoring dataset extends from 2014 to 2022 with an acquisition frequency from 10 min to hourly for the sensor variables and from weekly to monthly frequency for the sampling. Despite the frequent maintenance of the monitoring sites, several data gaps exist due to the remote location of some sites and instrument destruction by cyclones. Preliminary results show that the Rivière des Pluies watershed is characterized by high annual precipitation (>3000 mm y⁻¹) and a fast hydrologic response to precipitation (≈2 h basin lag time). The long-term evolution of the deep groundwater recharge is mainly driven by the occurrence of cyclone events with a seasonal groundwater response. The water chemical results support existing hydrogeological conceptual models suggesting a deep infiltration of the upstream infiltrated rainfall. The TMCF of Plaine des Fougères shows a high water storage capacity (>2000% for the Bryophytes) that makes this one a significant input of water to groundwater recharge which still needs to be quantified. This observatory is a unique research site in an insular volcanic tropical environment offering three windows of observation for the study of critical zone processes through upstream-midstream-downstream measurements sites. This high-resolution dataset is valuable to assess the response of volcanic tropical watersheds and aquifers at both event and long-term scales (i.e., global change). It will also provide insights in the hydrogeological conceptual model of volcanic islands, including the significant role of the TMCFs in the recharge processes as well as the watershed hydrosedimentary responses to extreme climatic events and their respective evolution under changing climatic conditions. All data sets are available at https://doi.org/10.5281/zenodo.7983138 .     

Phosphorus dynamics in water and sediments in urbanized and non-urbanized rivers in Southern Brazil

Phosphorus flux models show that the removal tends to be a common feature in polluted estuarine systems whereas the release of P to the adjacent coastal area occurs in pristine environments. This study analyzes the distribution of P in water and sediments along six rivers in the south coast of Brazil. Three rivers located inside protected areas were considered non-polluted and used as a control of urbanization. The other three, situated within urbanized areas, were considered as having highly elevated concentrations of P. Results showed a different behavior of P in water and sediments located in urbanized and non-urbanized areas. The concentrations of dissolved organic (P-org) and inorganic (P-inorg) phosphorus in water, and the total phosphorus (P-tot) and polyphosphate (P-poly) in sediments where higher in the urbanized rivers compared to the non-urbanized ones. Both P of punctual origin and of diffuse origin contributed to the maintenance of elevated concentrations and disturbed the natural fluxes along the polluted rivers. The minimum and maximum concentrations in urbanized areas varied from 0.39 to 12.45 (microM) for P-org and 0.00 to 5.92 (microM) for P-inorg in water, and in sediments from 89.90 to 808.16 (microM g(-1)) for P-tot and 0.00 to 76.51 (microM g(-1)) for P-poly. In non-urbanized areas concentrations in water varied from 0.22 to 1.20 (microM) for P-inorg and 4.43 to 5.56 (microM) for P-org, and in sediments from 45.91 to 652.26 (microM g(-1)) for P-tot and 0.00 to 8.61 (microM g(-1)) for P-poly. Using a hierarchical sampling design and a simple model of variation of P (K(d) model) it was possible to demonstrate that urbanized and non-urbanized estuaries may act as sinks or sources of P. The variation of P analysis in different points along each of the six rivers showed that release and removal areas may occur within a same river independently of its urbanization.     

Assessment of open boundary conditions on the elliptic formulation of the mild-slope equation

The performance of open boundaries in a finite differences scheme of the elliptic mild-slope equation is assessed. The wave propagation results show that lowest order parabolic radiation boundary conditions, unlike sponge layers combined with first order radiation boundary conditions, are an efficient alternative to first order radiation boundary conditions in order to improve the accuracy of the numerical solution of the problem.     

Paleoenvironmental reconstruction of a Late Quaternary lagoon system in southern Brazil (Jaguaruna region,Santa Catarina state) based on multi‐proxy analysis

This paper provides a paleoenvironmental reconstruction of a Late Quaternary lagoon system in the Jaguaruna region of Santa Catarina state, southern Brazil. Integrated results of bulk sedimentary organic matter characterization (δ¹³C, δ¹⁵N and C/N), microfossil (pollen and diatom) and grain‐size analysis from three shallow cores (～2.5 m depth) allowed us to propose an evolving paleogeographic scenario in this coastal region for the last ca. 5500 cal a BP. The lagoonal system in this area was more extensive during the mid‐Holocene than today, with a gradual and continuous lagoon–sea disconnection until the present. We add to the debate regarding relative sea‐level (RSL) variations for the Brazilian coast during the Holocene and discuss the importance of sedimentary dynamics for interpreting changes in coastal ecosystems. The multi‐proxy analysis suggests that changes in coastal ecosystems could be directly related to local sedimentary processes, which are not necessarily linked to RSL fluctuations and/or to climatic variations. Copyright © 2011 John Wiley & Sons, Ltd.     

The U-WEDGE Transformation Method for Multivariate Geostatistical Simulation

To speed up multivariate geostatistical simulation it is common to transform the set of attributes into spatially uncorrelated factors that can be simulated independently. Spatial decorrelation methods are usually based on the diagonalisation of the variance/covariance and semivariogram matrices of the set of attributes for a chosen family of lag spacings. These matrices are symmetric and there are several efficient methods for the approximate joint diagonalisation of a family of symmetric matrices. One of these is the uniformly weighted exhaustive diagonalisation with Gauss iterations (U-WEDGE) method. In contrast to the method of minimum/maximum autocorrelation factors (MAF), where a two structure linear model of coregionalisation is assumed, U-WEDGE can be applied directly to the set of experimental semivariogram matrices without having to place restrictions on the number of structures in the linear model of coregionalisation, thus removing one of the restrictions placed on the subsequent modelling of the spatial structure of the factors. We use an iron-ore data set to illustrate the method and present a comparison between the simulated attributes obtained from U-WEDGE and MAF with the full co-simulation of the attributes.