Groundwater evolution and mean water age inferred from hydrochemical and isotopic tracers in a tropical confined aquifer

The coastal confined aquifer in the Gulf of Urabá (Colombia) is an important water source for the banana agro‐industry as well as for urban and rural communities. However, the main processes controlling recharge and mixing in the aquifer are still poorly understood. Hydrochemical analyses and stable isotope monitoring were conducted to (a) determine groundwater recharge origin, mean groundwater age, and the main processes governing groundwater chemistry and the potential mixing of marine water and the influence of diffusive processes from the two surrounding aquitard layers. Hydrochemical data indicate that the main processes affecting the dissolved chemical composition include cation exchange, dissolution of carbonated and CO₂, and silicate weathering. δ¹⁸O and δ²H compositions combined with ¹⁴C data highlight the differences in climatic conditions between the recharge zone and the confined section of the aquifer, which is close to the Atlantic Ocean. Groundwater samples with ¹⁴C ages from recent to 28,300 years BP show a depleted isotopic trend ranging from ?6.43‰ to ?9.14‰ in δ¹⁸O and from ?43.2‰ to ?65.7‰ in δ²H. The most depleted δ¹⁸O and δ²H compositions suggest a cooler recharge climate than the current conditions (corresponding to the last glacial period of the late Pleistocene). Depleted δ¹³C values in the total dissolved inorganic carbon indicate the existence of organic material oxidation processes within the geologic formation. These results can be used or transferred to enhance groundwater modelling efforts in other confined coastal aquifers of South America where scarcity of long‐term monitoring data limits water resources planification under a changing climate.     

Runoff initiation in a preserved semiarid Caatinga small watershed,Northeastern Brazil

This study analyses some hydrological driving forces and their interrelation with surface‐flow initiation in a semiarid Caatinga basin (12 km²), Northeastern Brazil. During the analysis period (2005 – 2014), 118 events with precipitation higher than 10 mm were monitored, providing 45 events with runoff, 25 with negligible runoff and 49 without runoff. To verify the dominant processes, 179 on‐site measurements of saturated hydraulic conductivity (Ksat) were conducted. The results showed that annual runoff coefficient lay below 0.5% and discharge at the outlet has only occurred four days per annum on average, providing an insight to the surface‐water scarcity of the Caatinga biome. The most relevant variables to explain runoff initiation were total precipitation and maximum 60‐min rainfall intensity (I₆₀). Runoff always occurred when rainfall surpassed 31 mm, but it never occurred for rainfall below 14 mm or for I₆₀ below 12 mm h^?1. The fact that the duration of the critical intensity is similar to the basin concentration time (65 min) and that the infiltration threshold value approaches the river‐bank saturated hydraulic conductivity support the assumption that Hortonian runoff prevails. However, none of the analysed variables (total or precedent precipitation, soil moisture content, rainfall intensities or rainfall duration) has been able to explain the runoff initiation in all monitored events: the best criteria, e.g. failed to explain 27% of the events. It is possible that surface‐flow initiation in the Caatinga biome is strongly influenced by the root‐system dynamics, which changes macro‐porosity status and, therefore, initial abstraction. Copyright © 2016 John Wiley & Sons, Ltd.     

A simple water balance model adapted for soil water repellency: application on Portuguese burned and unburned eucalypt stands

Soil water repellency can impact soil hydrology, overland flow generation and associated soil losses. However, current hydrological models do not take it into account, which creates a challenge in repellency‐prone regions. This work focused on the adaptation for soil water repellency of a daily water balance model. Repellency is estimated from soil moisture content using site‐specific empirical relations and used to limit maximum soil moisture. This model was developed and tested using approximately 2 years of data from one long‐unburned and two recently burned eucalypt plantations in northern Portugal, all of which showed strong seasonal soil water repellency cycles. Results indicated important improvements for the burned plantations, with the Nash–Sutcliffe efficiency increasing from ?0.55 and ?0.49 to 0.55 and 0.65. For the unburned site, model performance was already good without the modification and efficiency only improved slightly from 0.71 to 0.74, mostly due to the better simulation of delayed soil wetting after dry periods. Results suggested that even a simple approach to simulate soil water repellency can markedly improve the performance of hydrological models in eucalypt forests, especially after fire. Copyright © 2015 John Wiley & Sons, Ltd.     

Beach erosion and recovery during consecutive storms at a steep‐sloping,meso‐tidal beach

This study analyses beach morphological change during six consecutive storms acting on the meso‐tidal Faro Beach (south Portugal) between 15 December 2009 and 7 January 2010. Morphological change of the sub‐aerial beach profile was monitored through frequent topographic surveys across 11 transects. Measurements of the surf/swash zone dimensions, nearshore bar dynamics, and wave run‐up were extracted from time averaged and timestack coastal images, and wave and tidal data were obtained from offshore stations. All the information combined suggests that during consecutive storm events, the antecedent morphological state can initially be the dominant controlling factor of beach response; while the hydrodynamic forcing, and especially the tide and surge levels, become more important during the later stages of a storm period. The dataset also reveals the dynamic nature of steep‐sloping beaches, since sub‐aerial beach volume reductions up to 30 m³/m were followed by intertidal area recovery (–2 < z < 3 m) with rates reaching ~10 m³/m. However, the observed cumulative dune erosion and profile pivoting imply that storms, even of regular intensity, can have a dramatic impact when they occur in groups. Nearshore bars seemed to respond to temporal scales more related to storm sequences than to individual events. The formation of a prominent crescentic offshore bar at ~200 m from the shoreline appeared to reverse the previous offshore migration trend of the inner bar, which was gradually shifted close to the seaward swash zone boundary. The partially understood nearshore bar processes appeared to be critical for storm wave attenuation in the surf zone; and were considered mainly responsible for the poor interpretation of the observed beach behaviour on the grounds of standard, non‐dimensional, morphological parameters. Copyright © 2011 John Wiley & Sons, Ltd.     

Evidences for multiple remagnetization of Proterozoic dykes from Iguerda inlier (Anti-Atlas Belt,Southern Morocco)

No paleomagnetic data exist for Paleo-Mesoproterozoic times of the West African Craton (WAC). Therefore, paleogeographic reconstructions for such old geological times are difficult to constrain. Gaps on the sedimentary record and intense remagnetizations are the major problems that paleomagnetic studies come across. Recent geochronological results for dyke swarms that intrude several Proterozoic inliers of WAC in the Anti-Atlas Belt (southern Morocco) revealed ages between Paleoproterozoic and early Neoproterozoic, opening for the first time a window of opportunity to conduct paleomagnetic studies and tentatively infer about the paleoposition of WAC during Proterozoic. On this scope we conducted a paleomagnetic study on seven Proterozoic dykes of the Iguerda inlier. The meaning of the obtained paleomagnetic directions was evaluated by rock magnetic and mineral analyses, complemented by petrographic observations. Our samples record the presence of a complex history of remagnetization, mostly assigned to several Phanerozoic thermal/chemical events, in particular to the late stages of Pan African orogeny (s.l.), to the Late Carboniferous Variscan orogeny, and even to more recent events. The recognized remagnetization processes are related to widespread metamorphic events under greenschist facies followed by low-temperature oxidation, both responsible for the formation of new magnetic phases, like magnetite and hematite. These events obliterated the primary (magmatic) thermo-remanent magnetization and promoted multiple remagnetizations of the dykes, thermally and chemically. For only one dyke the presence of primary magnetization is possible to infer, though not to confirm, and would place WAC at an equatorial position around 1750 Ma.     

Anomalous Diffuse CO2 Emissions at the Masaya Volcano (Nicaragua) Related to Seismic-Volcanic Unrest

Anomalous changes in the diffuse emission of carbon dioxide within the Masaya caldera have been observed before two seismic events that occurred at 10 and 30 km from the observation site. Their epicenters are located, respectively, south of Managua in Las Colinas (4.3 magnitude) and the Xiloa caldera (3.6 magnitude), in 2002 and 2003, recorded by the geochemical station located at El Comalito, Masaya volcano (Nicaragua). Anomalous increases were observed, which occurred around 50 and 8 days before the main seismic event that took place in Las Colinas, and 4 days before the seismic swarm at the Xiloa caldera, with a maximum CO₂ efflux of 9.3 and 10.7 kg m^?2 day^?1, respectively. The anomalous CO₂ efflux increases remained after filtering with multiple regression analysis was applied to the CO₂ efflux time series, which indicated that atmospheric variables, during the first 4 months, explained 23 % CO₂ variability, whereas, during the rest of the time series, CO₂ efflux values are poorly controlled with only 6 %. The observed anomalies of the diffuse CO₂ emission rate might be related to pressure changes within the volcanic–hydrothermal system and/or to geostructural changes in the crust due to stress/strain changes caused before and during the earthquakes’ formation, and seem not to be related to the activity of the main crater of Masaya volcano.     

Physical properties of morphological units on Comet 9P/Tempel 1 derived from near-IR Deep Impact spectra

In this paper we analyze near-infrared thermal emission spectra of the spatially resolved nucleus of Comet 9P/Tempel 1 obtained by the NASA spacecraft Deep Impact. Maps of spectral reddening, the product X^′ between the beaming function and directional emissivity, as well as surface temperature are constructed. Thermophysical modeling is used to estimate the degree of small scale surface roughness and thermal inertia by detailed reproduction of the empirical temperature map. Mie and Hapke theories are used in combination with numerically calculated beaming functions to analyze the X^′ map and place constraints on composition and grain size of the surface material. We show that it is absolutely mandatory to include small scale surface roughness in thermophysical modeling of this object, since the resulting self heating is vital for reproducing the measured temperatures. A small scale self heating parameter in the range 0.6?ξ?0.75 is common, but smoother areas where 0.2?ξ?0.3 are also found. Contrary to models neglecting small scale surface roughness, we find that the thermal inertia of Comet 9P/Tempel 1 generally is high (1000-3000 J m⁻² K⁻¹ s^−1/2), although it may be substantially lower (40-380 J m⁻² K⁻¹ s^−1/2) in specific areas. We obtain a disk-averaged reddening of 3.5% kÅ⁻¹, with statistically significant local variations around that value on a ±1.0% kÅ−1 level. Vast regions appear covered by small (∼0.1 μm) highly absorbing grains such as carbon or iron-rich silicates. Other regions appear dominated by somewhat larger (∼0.5 μm) and/or less absorbing grains such as troilite or magnesium-rich silicates. Surface variations in reddening, roughness, thermal inertia, composition and/or grain size are moderately to strongly correlated to the locations of morphological units on the surface. The existence of morphological units with differing physical properties may be primordial, hence reflecting a diversity in the building block cometesimals, or resulting from evolutionary processes.