Crustal xenoliths from Calbuco Volcano,Andean Southern Volcanic Zone: implications for crustal composition and magma-crust interaction

Crustal xenoliths in the 1961 andesite flow of Calbuco Volcano, in the southern Southern Volcanic Zone (SSVZ) of the Andes, consist predominantly of pyroxene granulites and hornblende gabbronorites. The granulites contain plagioclase+pyroxene+magnetite±amphibole, and have pronounced granoblastic textures. Small amounts of relict amphibole surrounded by pyroxene-plagioclase-magnetite-glass symplectites are found in some specimens. These and similar textures in the gabbronorites are interpreted as evidence of dehydration melting. Mineral and bulk rock geochemical data indicate that the granulites are derived from an incompatible trace element depleted basaltic protolith that underwent two stages of metamorphism: a moderate pressure, high temperature stage accompanied by melting and melt extraction from some samples, followed by thermal metamorphism after entrainment in the Calbuco andesite lavas. High _Nd ^T values (+4.0 to +8.6), Nd-isotope model ages of 1.7–2.0 Ga, and trace element characteristics like chondrite normalized La/Yb< and La/Nb1 indicate that the protoliths were oceanic basalts. Similar oceanic metabasalts of greenschist to amphibolite facies are found in the Paleozoic metamorphic belt that underlies the Chilean coastal ranges. Mineral and bulk rock compositions of the gabbronorite xenoliths indicate that they are cognate, crystallizing from the basaltic andesite magma at Calbuco. Crystallization pressures for the gabbros based on total Al contents in amphibole are 6–8 kbar. These pressures point to middle to lower crustal storage of the Calbuco magma. Neither granulite nor gabbro xenoliths have the appropriate geochemical characteristics to be contaminants of Calbuco andesites, although an ancient sedimentary contaminant is indicated by the lava compositions. The presence of oceanic metabasaltic xenoliths, together with the sedimentary isotopic imprint, suggests that the lower crust beneath the volcano is analogous to the coastal metamorphic belt, which is an accretionary complex of intercalated basalts and sediments that formed along the Paleozoic Gondwanan margin. If this is the case, the geochemical composition of the lower and middle crust beneath the SSVZ is significantly different from that of most recent SSVZ volcanic rocks.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Screening model for volatile pollutants in dual porosity soils

This paper develops mass fraction models for transport and fate of agricultural pollutants in structured two-region soils. Mass fraction index models, based on a semi-infinite domain solution, are derived that describe leaching at depth, vapor losses through soil surface, absorption, and degradation in the dynamic- and stagnant-water soil regions. The models predict that leaching is the result of the combined effect of the upward vapor-phase transport relative to downward advection, residence time relative to half-life, dispersion, and lateral diffusive mass transfer. Simulations show that leached fraction of volatile compounds does not always decrease monotonically with increased residence time relative to the pollutant half-life, as a result of complex interactions among the different physical and biochemical processes. The results show that leaching, volatilization, and degradation losses can be affected significantly by lateral diffusive mass transfer into immobile-water regions and advection relative to dispersion (i.e. Peclet number) in the mobile-water regions. It is shown that solute diffusion into the immobile phase and subsequent biochemical decay reduces leaching and vapor losses through soil surface. Potential use of the modified leaching index for the screening of selected pesticides is illustrated for different soil textures and infiltration rates. The analysis may be useful to the management of pesticides and the design of landfills.     

Groundwater contamination with arsenic and other trace elements in an area of the pampa,province of Córdoba,Argentina

A geochemical study of groundwater of the pampa in the province of Córdoba, Argentina, was performed; the area covered approximately 10,000 km².Physical-chemical parameters, dissolved solids, and seven trace elements were determined in 60 selected water samples. Systematic and accurate measurements of arsenic, flourine, and vanadium were performed for the first time. Three trace element contaminants not reported earlier were found: an important one, selenium, and two others of less known effects, uranium and molybdenum.Eighty-four percent of the water analyzed showed arsenic contents over 0.05 mg/L, maximum contaminant level established by the U.S. Environmental Protection Agency (1982). The frequency distribution of trace elements was analyzed, and its fit to the lognormal distribution was proved by means of the Pearson and Kolmogorov-Smirnov test; the geographic distribution of the seven trace elements was mapped and its correlation with the anion-cation composition of the water was studied.The maximum arsenic, fluorine, vanadium, and uranium contents were found in the western part of the area under study, in waters containing dominant alkali metals in the cation composition. Maximum selenium and antimony contents were found in the eastern part of the area, while molybdenum distribution does not show any relationship with the other two groups. In addition, the geographic distribution of the trace elements seems to be related to the subsurface structure, which has been inferred using interactive digital analysis of Landsat imagery. The movements of the subsoil have disturbed surface and subsurface drainage influencing the water salinity and trace element contents.In order to investigate the origin of the contamination, 54 loess samples were collected in wells at depths ranging from the surface down to the water table. This loess, which has a high proportion of volcanic components, mainly rhyolitic glass, exhibits a chemical composition corresponding to that of a dacite.The loess and the volcanic glass show anomalous contents of all contaminant trace elements, mainly arsenic and selenium. For this reason loess is considered the most important contamination source in the groundwater under study.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

Likely effects of climate change on groundwater availability in a Mediterranean region of Southeastern Spain

Groundwater resources are typically the main fresh water source in arid and semi‐arid regions. Natural recharge of aquifers is mainly based on precipitation; however, only heavy precipitation events (HPEs) are expected to produce appreciable aquifer recharge in these environments. In this work, we used daily precipitation and monthly water level time series from different locations over a Mediterranean region of Southeastern Spain to identify the critical threshold value to define HPEs that lead to appreciable aquifer recharge in this region. Wavelet and trend analyses were used to study the changes in the temporal distribution of the chosen HPEs (≥20 mm day^?1) over the observed period 1953–2012 and its projected evolution by using 18 downscaled climate projections over the projected period 2040–2099. The used precipitation time series were grouped in 10 clusters according to similarities between them assessed by using Pearson correlations. Results showed that the critical HPE threshold for the study area is 20 mm day^?1. Wavelet analysis showed that observed significant seasonal and annual peaks in global wavelet spectrum in the first sub‐period (1953–1982) are no longer significant in the second sub‐period (1983–2012) in the major part of the ten clusters. This change is because of the reduction of the mean HPEs number, which showed a negative trend over the observed period in nine clusters and was significant in five of them. However, the mean size of HPEs showed a positive trend in six clusters. A similar tendency of change is expected over the projected period. The expected reduction of the mean HPEs number is two times higher under the high climate scenario (RCP8.5) than under the moderate scenario (RCP4.5). The mean size of these events is expected to increase under the two scenarios. The groundwater availability will be affected by the reduction of HPE number which will increase the length of no aquifer recharge periods (NARP) accentuating the groundwater drought in the region. Copyright © 2016 John Wiley & Sons, Ltd.     

Detection and mapping of burnt areas from time series of MODIS-derived NDVI data in a Mediterranean region

Moderate resolution remote sensing data, as provided by MODIS, can be used to detect and map active or past wildfires from daily records of suitable combinations of reflectance bands. The objective of the present work was to develop and test simple algorithms and variations for automatic or semiautomatic detection of burnt areas from time series data of MODIS biweekly vegetation indices for a Mediterranean region. MODIS-derived NDVI 250m time series data for the Valencia region, East Spain, were subjected to a two-step process for the detection of candidate burnt areas, and the results compared with available fire event records from the Valencia Regional Government. For each pixel and date in the data series, a model was fitted to both the previous and posterior time series data. Combining drops between two consecutive points and 1-year average drops, we used discrepancies or jumps between the pre and post models to identify seed pixels, and then delimitated fire scars for each potential wildfire using an extension algorithm from the seed pixels. The resulting maps of the detected burnt areas showed a very good agreement with the perimeters registered in the database of fire records used as reference. Overall accuracies and indices of agreement were very high, and omission and commission errors were similar or lower than in previous studies that used automatic or semiautomatic fire scar detection based on remote sensing. This supports the effectiveness of the method for detecting and mapping burnt areas in the Mediterranean region.     

Mixing, sediment accumulation and focusing using 210Pb and 137Cs

The Portil lagoon is a natural freshwater reservoir located at the southwest of Spain, near the coast. In its surroundings an important tourist complex has been developed since the earlier 1970s. This has resulted in increased loads of nutrients, specially during summer months. In order to evaluate the impact from humans on the lagoon caused by anthropogenic activities, we have determined vertical profiles of excess ²¹⁰Pb and ¹³⁷Cs in a sediment core taken from its bottom. Vertical profile of excess ²¹⁰Pb seems to indicate that the sediment core was affected by large-scale mixing processes. Nevertheless this possibility was ruled out based on the shape of the ¹³⁷Cs vertical distribution after modelling the ¹³⁷Cs profile considering large-scale mixing. The chronology developed from ²¹⁰Pb and ¹³⁷Cs has allowed us to calculate the temporal evolution of the mass sedimentation rate during the last century. This mass sedimentation rate increased in 1973: from 0.08(2) g cm^-2 a^-1 in the period 1900-1973 to 0.17(5) g cm^-2 a^-1 in the period 1973-1995. This increase could be related to land-movements in the drainage area and to permitted sewage inputs into the lagoon, both with its origin in the growing/running of the surrounding tourist complex.     

13C12C ratios of pleistocene mummified remains from beringia

During the Quaternary glacial episodes, when sea level was considerably lower, Asia and North America were linked by large extensions of circumarctic land (Beringia), which remained unglaciated. This land mass served not only as a biogeographical bridge for plants, animals, and humans, but also supported a biome very different from present tundra or boreal coniferous forests, which was dominated by steppes and a rich mammalian megafauna. Carbon stable isotope ratios of Beringian late Pleistocene mummified remains of bison, equids, mammoth, caribou, musk-ox, moose, woolly rhino, and other undetermined species, found preserved in permafrost, indicate that these megaherbivores fed exclusively on C₃ plants, and that C₄ grasses were not differentially ingested by bison, as previously suggested. Paleoclimatic constraints probably prevented the formation of a warm-season (C₄) guild during the later part of the growing season in the steppes of Beringia during the last glaciation.