Salt and nitrate exports from the sprinkler-irrigated Malfarás creek watershed (Ebro river valley,Spain) during 2010

Irrigated agriculture is a clear source of non-point pollution by salts and nitrogen species. The impact of such pollution should be quantified according to specific cases. The case of the Malfarás creek basin, a sprinkler irrigation district located in the semiarid Ebro valley in northeast Spain, has been evaluated. The main crops in the district were corn, barley and alfalfa, occupying 93 % of the irrigated area. The fate of water, salts and nutrients was evaluated by a daily water balance developed at a field scale for the natural year 2010. The yearly data of the whole set of 101 irrigated fields plus the non-irrigated area compared to the measured drainage produced a basin water balance with a low degree of error. The basin consumed 90 % of the total water input of which 68 % was used for crop evapotranspiration and the rest was lost due to non-productive uses. 16 % of the incoming water left the irrigation area as drainage water. The irrigated area was responsible for 87 % of the drainage. The average volume of drained water was 152 mm year^?1 for the whole basin area. The irrigated area drained 183 mm year^?1. The basin exported 473 kg of salt per hectare during 2010. This value was the lowest of the sprinkler irrigation areas in the Ebro valley, mainly due to the lower soil salinity. All the crops except barley received a nitrogen surplus of 10–50 % above their needs. The extra nitrogen entered the water cycle increasing the nitrate concentration in the aquifer water (150 mg L^?1) and drainage water (98 mg L^?1). In 2010 the mass of nitrogen exported by drainage was 49 kg per irrigated hectare. This value is too high for this type of irrigation system and implies that 17 % of nitrogen applied as a fertilizer was lost to drainage water. The key to decreasing the nitrogen leaching and pollution that it causes could be appropriate time-controlled fertigation along with better irrigation scheduling.     

Applications géodésiques du système DORIS à l'Institut géographique national

IGN is in charge of the installation and maintenance of the DORIS orbit determination network. More recently, in collaboration with JPL, precise geodetic computations were performed. The goal of this paper is to recall the various historic contributions of IGN to the DORIS system in their international context and then to describe a new estimation technique developed for a multi-satellite mode, making full profit of a better modeling for satellites and ground clocks as well as tropospheric correction parameters. Derived geodetic results demonstrate a precision in the order of 1 cm for station positions. To cite this article: P. Willis et al., C. R. Geoscience 337 (2005).     

Preliminary data on chemical changes in the Aral Sea during low-level periods from the last 9000 years

The withdrawal of the Aral Sea tributaries (Amu and Syr Daria) for cultures has led to significant falls of its level and an important increase in its salinity. During the Holocene, a succession of low and high water inputs occurred. Silty deposits correspond to the high levels and carbonates to the low levels. This study makes a distinction between the Syr Daria and the Amu Daria water inputs during low-level periods by using mineralogical and chemical compositions of the carbonates deposits. Waters from the Syr Daria are more sulphatic and have a low iron content in comparison with that of the Amu Daria. The Syr Daria was the major tributary around 7500, 4956 and 970 yr?BP, whereas around 6200 and 3610 yr?BP, inflow also from the Amu Daria is observed. To cite this article: L. Le Callonnec et al., C. R. Geoscience 337 (2005).     

Holocene relative sea level changes in the Västervik-Gamlebyviken region on the southeast coast of Sweden,southern Baltic Sea

We reconstruct the Holocene shore displacement of the Västervik-Gamlebyviken area on the southeast coast of Sweden, characterised by a maritime cultural landscape and archaeological significance since the Mesolithic. Sediment cores were retrieved from four lake basins that have been raised above sea level due to the postglacial land uplift and eustatic sea level changes after the melting of the Fennoscandian Ice Sheet. The cores were radiocarbon dated and analysed for loss on ignition and diatoms. The isolation thresholds of the basins were determined using LiDAR data. The results provide evidence for the initiation of the first Littorina Sea transgression in this area at 8.5 thousand calibrated years before present (cal. ka BP). A relative sea level rise by ∼7 m a.s.l. is recorded between 8.0 and 7.5 cal. ka BP with a highstand at ∼22 m a.s.l. between 7.5 and 6.2 cal. ka BP. These phases coincide with the second and third Littorina Sea transgressions, respectively, in the Blekinge area, southern Sweden and are consistent with the final deglaciation of North America. After 6.2 cal. ka BP, the relative sea level dropped below 22 m a.s.l., and remained at ∼20 m a.s.l. until 4.6 cal. ka BP coinciding with the fourth Littorina Sea transgression in Blekinge. From 4.6 to 4.2 cal. ka BP, the shore displacement shows a regression rate of 10 mm a⁻¹ followed by a slowdown with a mean value of 4.6 mm a⁻¹ until 1.6 cal. ka BP, when the relative sea level dropped below 3.3 m a.s.l. The Middle to Late Holocene highstand and other periods of minor sea level transgressions and/or higher salinity between 6.2 and 1.7 cal. ka BP are attributed to a combination of warmer climate and higher inflow of saline waters in the southern Baltic Sea due to stronger westerlies, caused by variations in the North Atlantic atmospheric patterns.     

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.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

The Removal of Dissolved Metals by Hydroxysulphate Precipitates during Oxidation and Neutralization of Acid Mine Waters, Iberian Pyrite Belt

This study examines the removal of dissolved metals during the oxidation and neutralization of five acid mine drainage (AMD) waters from La Zarza, Lomero, Esperanza, Corta Atalaya and Poderosa mines (Iberian Pyrite Belt, Huelva, Spain). These waters were selected to cover the spectrum of pH (2.2–3.5) and chemical composition (e.g., 319–2,103 mg/L Fe; 2.85–33.3 g/L SO₄⁼) of the IPB mine waters. The experiments were conducted in the laboratory to simulate the geochemical evolution previously recognized in the field. This evolution includes two stages: (1) oxidation of dissolved Fe(II) followed by hydrolysis and precipitation of Fe(III), and (2) progressive pH increase during mixing with fresh waters. Fe(III) precipitates at pH < 3.5 (stages 1 and 2) in the form of schwertmannite, whereas Al precipitates during stage 2 at pH 5.0 in the form of several hydroxysulphates of variable composition (hydrobasaluminite, basaluminite, aluminite). During these stages, trace elements are totally or partially sorbed and/or coprecipitated at different rates depending basically on pH, as well as on the activity of the SO₄⁼ anion (which determines the speciation of metals). The general trend for the metals which are chiefly present as aqueous free cations (Pb²⁺, Zn²⁺, Cu²⁺, Cd²⁺, Mn²⁺, Co²⁺, Ni²⁺) is a progressive sorption at increasing pH. On the other hand, As and V (mainly present as anionic species) are completely scavenged during the oxidation stage at pH < 3.5. In waters with high activities (> 10⁻¹) of the SO ₄⁼ ion, some elements like Al, Zn, Cd, Pb and U can also form anionic bisulphate complexes and be significantly sorbed at pH < 5. The removal rates at pH 7.0 range from around 100% for As, V, Cu and U, and 60–80% for Pb, to less than 20% for Zn, Co, Ni and Mn. These processes of metal removal represent a significant mechanism of natural attenuation in the IPB.     

Stability analysis of tunnels driven in jointed rocks by means of a homogenized limit analysis approach

Conceived as a potential alternative to the classical design methods employed for analyzing the stability of underground structures driven in jointed rocks, the homogenization approach stems from the heuristic idea that, from a macroscopic point of view, a rock mass cut by a network of joints may be perceived as a homogenized continuum. The strength properties of the latter can be theoretically obtained from the failure conditions of its individual constituents: rock matrix and joint interfaces. At the material level, the limit analysis reasoning is used in the context of homogenization to formulate the homogenized strength criterion of a jointed rock mass in the particular situation of a single set of parallel joints. As it could be expected, the obtained closed‐form expressions show the strength anisotropy induced by joint preferential orientation. The support functions (π functions) associated with the homogenized strength criterion are also determined in both plane strain and three‐dimensional cases. This criterion is then applied to the investigation of stability analysis of a tunnel excavated in a jointed rock mass. Upper bounds estimated of the stability factor are derived from the implementation of the kinematic approach directly on the homogenized underground structure. Finally, the approach is applied to analyze and discuss the collapse of the Pinheiros subway station (São Paulo, Brazil). Copyright © 2014 John Wiley & Sons, Ltd.     

Climatic vs. tectonic signals in a continental extensional basin (Teruel,NE Spain) from stable isotope (δ18O) and sequence stratigraphical evolution

Alluvial fans and shallow carbonate lakes interfered in the Teruel half‐graben during the Late Miocene–Pliocene. Tectonic influence is recorded in alluvial and lacustrine–palustrine successions, with long‐term climate changes being recorded in detail in the isotopic signatures of carbonates. Episodes of tectonic activity induced alluvial fan progradation and lake retraction in the whole basin. Three lacustrine stages have been identified, which support the idea that climate also exerted an important control on sedimentation. The transition between stages 1 and 2 occurred during a tectonically calm episode due to an increase in aridity in the Early Turolian; small fans with source areas next to the lake margin prograded, inducing lake‐shore retraction. The transition from stage 2 to 3 was caused by the superimposition of increasing tectonic activity and aridity effects. Our study demonstrates that discrimination of allogenic factors controlling sedimentation in continental closed basins is possible using sequence stratigraphy in combination with other techniques such as geochemistry of carbonates.