The evolution of the high‐elevated depocenters of the northern Sierras Pampeanas (ca. 28˚ SL), Argentine broken foreland,South‐Central Andes: the Pipanaco Basin

The Pipanaco Basin, in the southern margin of the Andean Puna plateau at ca. 28°SL, is one of the largest and highest intermontane basins within the northernmost Argentine broken foreland. With a surface elevation >1000 m above sea level, this basin represents a strategic location to understand the subsidence and subsequent uplift history of high‐elevation depositional surfaces within the distal Andean foreland. However, the stratigraphic record of the Pipanaco Basin is almost entirely within the subsurface, and no geophysical surveys have been conducted in the region. A high‐resolution gravity study has been designed to understand the subsurface basin geometry. This study, together with stratigraphic correlations and flexural and backstripping analysis, suggests that the region was dominated by a regional subsidence episode of ca. 2 km during the Miocene‐Pliocene, followed by basement thrusting and ca. 1–1.5 km of sediment filling within restricted intermontane basin between the Pliocene‐Pleistocene. Based on the present‐day position of the basement top as well as the Neogene‐Present sediment thicknesses across the Sierras Pampeanas, which show slight variations along strike, sediment aggradation is not the most suitable process to account for the increase in the topographic level of the high‐elevation, close‐drainage basins of Argentina. The close correlation between the depth to basement and the mean surface elevations recorded in different swaths indicates that deep‐seated geodynamic process affected the northern Sierras Pampeanas. Seismic tomography, as well as a preliminary comparison between the isostatic and seismic Moho, suggests a buoyant lithosphere beneath the northern Sierras Pampeanas, which might have driven the long‐wavelength rise of this part of the broken foreland after the major phase of deposition in these Andean basins.     

Volcanic ash forecast during the June 2011 Cordón Caulle eruption

We modelled the transport and deposition of ash from the June 2011 eruption from Cordón Caulle volcanic complex, Chile. The modelling strategy, currently under development at the Argentinean Naval Hydrographic Service and National Meteorological Service, couples the weather research and forecasting (WRF/ARW) meteorological model with the FALL3D ash dispersal model. The strategy uses volcanological inputs inferred from satellite imagery, eruption reports and preliminary grain-size data obtained during the first days of the eruption from an Argentinean ash sample collection network. In this sense, the results shown here can be regarded as a quasi-syn-eruptive forecast for the first 16 days of the eruption. Although this article describes the modelling process in the aftermath of the crisis, the strategy was implemented from the beginning of the eruption, and results were made available to the Buenos Aires Volcanic Ash Advisory Centers and other end users. The model predicts ash cloud trajectories, concentration of ash at relevant flight levels, expected deposit thickness and ash accumulation rates at relevant localities. Here, we validate the modelling strategy by comparing results with satellite retrievals and syn-eruptive ground deposit measurements. Results highlight the goodness of the combined WRF/ARW-FALL3D forecasting system and point out the usefulness of coupling both models for short-term forecast of volcanic ash clouds.     

A study of NAO variability and its possible non-linear influences on European surface temperature

 The relationship between European winter temperature spatial and temporal modes of variability and the North Atlantic Oscillation (NAO) has been studied during the period 1852–1997. Temporal modes of variability of the NAO and temperatures are analysed using wavelet transform. Results show that the NAO presents a strong non-stationary behaviour. The most important feature is the existence of a quasi-periodic oscillation, with a period between 6–10 years and maximum amplitude of eight years, during the periods 1842–1868 and 1964–1994. Between 1875 and 1939 the spectra of the NAO is almost white. The possible relationship between the occurrence of extreme events of the NAO and its spectral behaviour has been analysed. The results indicate that quasi-periodic oscillations in the NAO do not lead to more extreme episodes, but rather that an extreme value of the oscillation is more likely to persist for few years. Particularly energetic modes of coherent variability between temperature and NAO are found between 2–6 years for 1857–1879 and 1978–1984, and between 6–10 years from 1961 to 1991. The relationship between the NAO and temperatures as a function of the state of the oscillation has been studied using composites. Empirical evidence has been found suggesting that winter temperatures, in a great part of the study area, do not vary in a linear manner with respect to phase and intensity of the NAO. Regions in the study area differ in sensitivity to changes in the NAO. The spatial patterns of variability of the temperatures are found to be independent of the NAO spectra. Received: 8 April 1999 / Accepted: 19 September 2000     

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

Effect of the addition of granitic powder to an acidic soil from Galicia (NW Spain) in comparison with lime

The impact of diffuse pollution, agricultural land use and climate change on the long-term response of subsurface–surface water quality is not well understood, but is a prerequisite for evaluation of water management options. The goal of this study is to model geochemical evolution of water chemistry from the infiltration through soil into the unsaturated zone, transport through bedrocks and granular aquifers to a river in order to identify zones of steep concentration gradients and high dynamics under transient flow conditions. A numerical model was constructed comprising a 2-D 1,500 m × 150 m vertical cross-section of typical sedimentary rock formations, a glacio-fluvial quaternary gravel aquifer in the valley and soil layers. The model coupled saturated/un-saturated flow and reactive transport under steady state and transient conditions. Geochemical interactions, include intra-aqueous kinetic reactions of oxygen with dissolved organic matter, as well as kinetics of carbonate dissolution/precipitation. This model section was chosen to provide insight in to the principal processes and time scales affecting water chemistry along different flow paths. The numerical simulator MIN3P was used, a finite volume program for variably saturated subsurface flow and multi-component reactive transport. The results show that subsurface water residence times range from approximately 2 to 2,000 years. Different zones are to be expected with respect to the development of mineral equilibria; namely, purely atmospherically influenced, as well as open and closed system carbonate dissolution. Short-term responses to daily averaged changes in precipitation, however, are only visible to some extent in the shallower and near-river parts of flow system and solute loads. This can most likely be explained by directional changes in flow paths, indicating that equilibrium geochemical condition predominate at the hillslope scale, i.e. water quality depends on transport pathways rather than on kinetic effects. The extent of reducing conditions is controlled by the presence of organic-rich layers (i.e. peat deposits), the dissolution kinetics of aquifer organic matter and the subsequent mixing with oxygenated water by hydrodynamic dispersion.     

Extracting water-related features using reflectance data and principal component analysis of Landsat images

This study aimed to map water features using a Landsat image rather than traditional land cover. We involved the original bands, spectral indices and principal components (PCs) of a principal component analysis (PCA) as input data, and performed random forest (RF) and support vector machine (SVM) classification with water, saturated soil and non-water categories. The aim was to compare the efficiency of the results based on various input data. Original bands provided 93% overall accuracy (OA) and bands 4–5–7 were the most informative in this analysis. Except for MNDWI (modified normalized differenced water index, with 98% OA), the performance of all water indices was between 60 and 70% (OA). The PCA-based approach conducted on the original bands resulted in the most accurate identification of all classes (with only 1% error in the case of water bodies). We therefore show that both water bodies and saturated soils can be identified successfully using this approach.     

Coastal Occupations in Tierra del Fuego,Southernmost South America: A Late Holocene Hunter‐Gatherer Context at Marazzi 2

This paper focuses on the late Holocene occupation of hunter‐gatherers at the Marazzi 2 site located on the northwestern steppe of Tierra del Fuego, Chile. Our aim is to understand stratigraphy, formation processes, and pedogenesis with respect to human occupation over the last 3000 years. Based on archaeological excavations on a fluvial terrace of the Torcido River, we integrate soil micromorphology, mineralogy, geochemistry, magnetic susceptibility, and geomorphology, as well as the micro‐ and macrofrequency distributions of archaeological material. A micro‐taphonomical perspective is also applied to anthropogenic components recorded in sedimentary thin sections. We discuss various events in the interplay between soil development and human occupation through time. Marazzi 2 was witness to aggradation and pedogenesis with an episode of surface stability coincident with a phase of more intense human occupation at about 860 BC. Taphonomic analysis suggests that there are multiple occupation events mixed together, probably by biomechanical processes, resulting in a palimpsest with poor temporal resolution.     

The 2011 submarine volcanic eruption in El Hierro (Canary Islands)

Forty years after the Teneguía Volcano (La Palma, 1971), a submarine eruption took place off the town of La Restinga, south of El Hierro, the smallest and youngest island of the Canarian Archipelago. Precursors allowed an early detection of the event and its approximate location, suggesting it was submarine. Uncertainties derived from insufficient scientific information available to the authorities during the eruption, leading to disproportionate civil protection measures, which had an impact on the island's economy—based primarily on tourism—while residents experienced extra fear and distress.     

Influence of Mussel Culture on the Vertical Export of Phytoplankton Carbon in a Coastal Upwelling Embayment (Ría de Vigo,NW Iberia)

The goal of this paper is to find out whether suspended mussel culture affects the vertical fluxes of biogenic particles in the Ría de Vigo on a seasonal scale. With this aim, vertical fluxes of particulate organic carbon (POC) and the magnitude and composition of vertical export of phytoplankton carbon (Cphyto) collected in sediment traps were examined by comparing data obtained inside a mussel farming area (RaS) with those found at a reference station (ReS) not affected by mussels. Our results indicate that mussel farming has a strong impact on sedimentation fluxes under the rafts, not only increasing POC flux but also altering the magnitude and composition of Cphyto fluxes. Average POC flux at RaS (2564?±?1936 mg m^?2 day^?1) was four times higher than at ReS (731?±?276 mg m^?2 day^?1), and much of this increase was due to biodeposit fluxes (Cbiodep) which accounted for large proportion of POC flux (35–60 %). Indeed, because of this high Cbiodep flux, only a small proportion of the POC flux was due to Cphyto flux (3–12 %). At the same time, we observed an increased sedimentation of phytoplankton cells at RaS that could be explained by a combination of mechanisms: less energetic hydrodynamic conditions under mussel rafts, ballast effect by sinking mussel feces, and diatom aggregates. Moreover, mussel farming also altered the quality of the Cphyto flux by removing part of the predatory pressure of zooplankton and thus matching diatom composition in water column and sediment traps.     

Trace elements contamination of agricultural soils affected by sulphide exploitation (Iberian Pyrite Belt,Sw Spain)

Agricultural soils of the Riotinto mining area (Iberian Pyrite Belt) have been studied to assess the degree of pollution by trace elements as a consequence of the extraction and treatment of sulphides. Fifteen soil samples were collected and analysed by ICP-OES and INAA for 51 elements. Chemical analyses showed an As–Cu–Pb–Zn association related with the mineralisation of the Iberian Pyrite Belt. Concentrations were 19–994 mg kg⁻¹ for As, 41–4,890 mg kg⁻¹ for Pb, 95–897 mg kg⁻¹ for Zn and of 27–1,160 mg kg⁻¹ for Cu. Most of the samples displayed concentrations of these elements higher than the 90th percentile of the corresponding geological dominium, which suggests an anthropogenic input besides the bedrock influence. Samples collected from sediments were more contaminated than leptosols because they were polluted by leachates or by mining spills coming from the waste rock piles. The weathering of the bedrock is responsible for high concentrations in Co, Cr and Ni, but an anthropogenic input, such as wind-blown dust, seems to be indicative of the high content of As, Cu, Pb and Zn in leptosols. The metal partitioning patterns show that most trace elements are associated with Fe amorphous oxy-hydroxides, or take part of the residual fraction. According to the results obtained, the following mobility sequence is proposed for major and minor elements: Mn, Pb, Cd, > Zn, Cu > Ni > As > Fe > Cr. The high mobility of Pb, Cu and Zn involve an environmental risk in this area, even in soils where the concentrations are not so high.