Chemical and isotopic assessment of surface water–shallow groundwater interaction in the arid Grande river basin,North-Central Chile

ABSTRACT

This paper analyses the composition of surface water and shallow groundwater in the Grande River basin, North-Central Chile, using this information to characterize water interactions. Chemical and isotopic data for surface water and groundwater (7 and 6 sampling locations, respectively) were obtained from three sampling campaigns performed in March–April (autumn), August–September (late winter) and December (early summer) 2012. Precipitation samples were also collected. Data was processed using spatial distribution charts, Piper and Stiff diagrams, and multivariate analysis. In general, the results for each method converge on a high degree of connectivity between surface water and shallow groundwater in the study area. Furthermore, approximately a 10% of groundwater contribution to the surface flow discharge was estimated for a particular reach. This multi-method approach was useful for the characterization of surface water–groundwater interactions in the Grande River basin, and may become a suitable and replicable scheme for studies in arid and semi-arid basins facing similar water management challenges.

Editor D. Koutsoyiannis; Associate editor B. Dewals     

Kinetic energy–rainfall intensity relationship for Central Cebu, Philippines for soil erosion studies

In the study of soil erosion, specifically on detachment of soil particles by raindrop impact, kinetic energy is a commonly suggested indicator of the raindrop's ability to detach soil particles from the soil mass. Since direct measurement of kinetic energy requires sophisticated and costly instruments, the alternative approach is to estimate it from rainfall intensity. The present study aims at establishing a relationship between rainfall intensity and kinetic energy for rainfalls in Central Cebu, Philippines as a preface of a wider regional investigation.

Drop size distributions of rainfalls were measured using the disdrometer RD-80. There are two forms of kinetic energy considered here. One is kinetic energy per unit area per unit time (KE_R, J m⁻² h⁻¹) and the other is kinetic energy per unit area per unit depth (KE, J m⁻² mm⁻¹). Relationships between kinetic energy per unit area per unit time (KE_R) and rainfall intensity (I) were obtained using linear and power relations. The exponential model and the logarithmic model were fitted to the KE–I data to obtain corresponding relationships between kinetic energy per unit area per unit depth of rainfall (KE) and rainfall intensity (I). The equation obtained from the exponential model produced smaller standard error of estimates than the logarithmic model.     

A New Appraisal of Sri Lankan BB Zircon as a Reference Material for LA‐ICP‐MS U‐Pb Geochronology and Lu‐Hf Isotope Tracing

A potential zircon reference material (BB zircon) for laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) U‐Pb geochronology and Hf isotope geochemistry is described. A batch of twenty zircon megacrysts (0.5–1.5 cm³) from Sri Lanka was studied. Within‐grain rare earth element (REE) compositions are largely homogeneous, albeit with some variation seen between fractured and homogeneous domains. Excluding fractured cathodoluminescence bright domains, the variation in U content for all analysed crystals ranged from 227 to 368 μg g^?1 and the average Th/U ratios were between 0.20 and 0.47. The Hf isotope composition (0.56–0.84 g/100 g Hf) is homogeneous within and between the grains – mean ¹⁷⁶Hf/¹⁷⁷Hf of 0.281674 ± 0.000018 (2s). The calculated alpha dose of 0.59 × 10¹⁸ g^?1 for a number of BB grains falls within the trend of previously studied, untreated zircon samples from Sri Lanka. Aliquots of the same crystal (analysed by ID‐TIMS in four different laboratories) gave consistent U‐Pb ages with excellent measurement reproducibility (0.1–0.4% RSD). Interlaboratory assessment (by LA‐ICP‐MS) from individual crystals returned results that are within uncertainty equivalent to the TIMS ages. Finally, we report on within‐ and between‐grain homogeneity of the oxygen isotope systematic of four BB crystals (13.16‰ VSMOW).     

Zero-offset P–S energy in horizontally layered media

The existence of non‐zero reflectivity for zero‐offset PS waves in horizontally layered media is discussed. Field measurements have suggested the occurrence of this phenomenon. We show that, in some cases, if anisotropy is considered in the problem, this reflectivity can be predicted. By using an approximated formulation to determine the coefficient of reflection in arbitrarily anisotropic media, it is clear that, for some elastic arrangements of the subsurface, reflected energy is associated with the converted wave resulting from normal P‐wave incidence.     

Ecotoxicological evaluation of dredged sediments from Rodrigo de Freitas Lagoon (Rio de Janeiro State,Brazil) using bioassays with earthworms and collembolans

The toxicity and metal bioavailability were studied in dredged sediments from Rodrigo de Freitas Lagoon (Rio de Janeiro, RJ, Brazil) using acute and avoidance tests with Eisenia andrei, and reproduction tests with Folsomia candida. The sediment was mixed with an artificial soil, and two natural soils (ferralsol and chernosol—representative Brazilian tropical soils) to obtain the following doses: 1, 3, 6, 12, 18, 24 and 30%. Total metal concentrations were determined in the sediment to support the interpretation of ecotoxicological data. Metal concentrations in the mixtures were in agreement with the threshold limits established by Brazilian law. However, significant avoidance responses were found on doses ≥?3% and were the most sensitive endpoint. Earthworm mortality found in artificial soil mixtures (LC50?=?3.9) suggests higher toxicity levels than those obtained in ferralsol (LC50?=?7.6%) and chernosol (11.0%) treatments. Earthworm mortality, avoidance responses and collembolan reproduction levels found in ferralsol mixtures (LC50?=?9.2; avoidance EC50?=?2.3%; reproduction EC50?=?2.8%) were higher compared to chernosol treatments (LC50?=?11.0%; avoidance EC50?=?4.3%; reproduction EC50?=?4.9%). The reduction of toxicity levels in chernosol mixtures is probably due to the abundance of expansive clay minerals in chernosols with capacity of adsorbing metals and other xenobiotic substances from soil pore water, decreasing metal bioavailability. Finally, threshold limits defined by Brazilian legislation for soil quality and land disposal of dredged sediments are not sufficient to prevent noxious effects on soil fauna and should be complemented with a preliminary ecotoxicological evaluation.     

Wake Response to an Ocean-Feedback Mechanism: Madeira Island Case Study

We focus on an island wake episode that occurred in the Madeira Archipelago region of the north-east Atlantic at $32.5^{\circ }\mathrm{N}, 17^{\circ }\mathrm{W}$ . The Weather Research and Forecasting numerical model was used in a (one-way) downscaling mode, considering initial and boundary conditions from the European Centre for Medium-range Weather Forecasts system. The current literature emphasizes adiabatic effects on the dynamical aspects of atmospheric wakes. Changes in mountain height and consequently its relation to the atmospheric inversion layer should explain the shift in wake regimes, from a ‘strong-wake’ to ‘weak-wake’ scenario. Nevertheless, changes in sea-surface temperature variability in the lee of an island can induce similar regime shifts because of exposure to stronger solar radiation. Increase in evaporation contributes to the enhancement of convection and thus to the uplift of the stratified atmospheric layer above the critical height, with subsequent internal gravity wave activity.     

Long-term abundance patterns of macroalgae in relation to environmental variables in the Tagus Estuary (Portugal)

Seasonal and inter-annual patterns of macroalgal abundance in a Tagus Estuary oyster reef are described. Macroalgal abundance was estimated as canopy percent cover by three permanent point intercept transects over a 7-year period. Four categories were defined, corresponding to bare substrate and three different macroalgal functional-form groups: (1) ULVA, foliose group, included Ulva spp.; (2) GRACIL, terete corticated macrophyte group, included only Gracilaria gracilis; and (3) FILAM, small (<10 cm) filamentous group, including eight species. A canonical correspondence analysis (CCA) showed that: (1) ULVA were associated with long and hot days, being usually dominant during spring and especially summer; (2) FILAM were associated with mild temperatures and relatively long days, abundant in spring but showed frequent peaks in summer; and (3) GRACIL were also favoured by spring season, although associated to lower temperature and less daylight hours than FILAM. GRACIL and FILAM were present throughout the year. On the contrary, ULVA were absent or with low cover during colder periods. A negative correlation between GRACIL and FILAM seems to indicate competition between the two categories. The applied models explained 23.3% of the temporal variance in category abundance. Rainfall negatively affected macroalgal cover, as indicated by the positive correlation between rainfall and bare substrate. Our conclusions are in agreement with previous studies that consider algae as excellent environmental integrators, even on a small scale, due to a strong link between the macroalgal communities and relevant environmental variables. It is also relevant that this study used open-access databases of environmental variables, which open up new possibilities for mining existing data resources in new ways. Due to large inter-annual variability, long-term studies are essential to understand population dynamics in estuarine phytobenthic communities.     

Influence of upwelling and river runoff interaction on phytoplankton assemblages in a Middle Galician Ria and Comparison with northern and southern rias (NW Iberian Peninsula)

The first oceanographic research (hydrography, nutrient salts, chlorophyll, primary production and phytoplankton assemblages) in a Middle Galician Ria was carried out in Corme-Laxe during 2001, just a year before the Prestige oil spill, being the only reference to evaluate eventual changes in the phytoplankton community. Due to the small size of this ria (6.5 km²), oceanographic processes were driven by the continental water supplied by Anllons River during the wet season (20–30 m³ s⁻¹ in winter), and the strong oceanic influence from the nearby shelf during the dry season. The annual cycle showed a spring bloom with high levels of chlorophyll (up to 14 μg Chl-a L⁻¹) and primary production (3 g C m⁻² d⁻¹) and a summer upwelling bloom (up to 8 μg Chl-a L⁻¹ and 10 g C m⁻² d⁻¹) where the proximity of the Galician upwelling core (<13.5 °C at sea surface) favors the input of upwelled seawater (up to 9 μM of nitrate and silicate) to the bottom ria layer, even during summer stratification events (primary production around 2 g C m⁻² d⁻¹). Thus, phytoplankton assemblages form a “continuum” from spring to autumn with a predominance of diatoms and overlapping species between consecutive periods; only in autumn dinoflagellates and flagellates characterized the phytoplankton community. In the Middle Rias as Corme-Laxe, the nutrient values, Chl-a, primary production and phytoplankton abundance for productive periods were higher than those reported for the Northern (Ria of A Coruña) and Southern Rias (Ria of Arousa) for year 2001; this suggests the importance of the hydrographic events occurring in the zone of maximum upwelling intensity of the Western Iberian Shelf, where a lack of annual cycles studies exists.     

Distinguishing the Regional Atmospheric Controls on Precipitation Isotopic Variability in the Central-Southeast Portion of Brazil

Precipitation isotope ratios (O and H) record the history of water phase transitions and fractionation processes during moisture transport and rainfall formation. Here, we evaluated the isotopic composition of precipitation over the central-southeastern region of Brazil at different timescales. Monthly isotopic compositions were associated with classical effects (rainfall amount, seasonality, and continentality), demonstrating the importance of vapor recirculation processes and different regional atmospheric systems (South American Convergence Zone-SACZ and Cold Fronts-CF). While moisture recycling and regional atmospheric processes may also be observed on a daily timescale, classical effects such as the amount effect were not strongly correlated (δ¹⁸O-precipitation rate r ≤ –0.37). Daily variability revealed specific climatic features, such as δ¹⁸O depleted values (~ –6‰ to –8‰) during the wet season were associated with strong convective activity and large moisture availability. Daily isotopic analysis revealed the role of different moisture sources and transport effects. Isotope ratios combined with d-excess explain how atmospheric recirculation processes interact with convective activity during rainfall formation processes. Our findings provide a new understanding of rainfall sampling timescales and highlight the importance of water isotopes to decipher key hydrometeorological processes in a complex spatial and temporal context in central-southeastern Brazil.