Surface-water quality for the Salí River watershed in NW Argentina

Spatial and temporal variations of Na, Mg, Si, K, Ca, SO₄^2-, Cl^-, HCO₃^-, Ag, Al, As, Au, B, Ba, Be, Bi, Br, Ce, Cd, Co, Cr, Cs, Cu, Dy, Er, Eu, F, Fe, Ga, Ge, Gd, Hf, Hg, Ho, I, La, Li, Mn, Mo, Nb, Nd, Ni, P, Pb, Pr, Rb, Sb, Sm, Sn, Sr, Ta, Tb, Te, Th, Ti, Tl, U, V, W, Y, Yb, Zn, and Zr were monitored through 37 sampling stations to determine the main aspects influencing the surface-water quality in the Salí River watershed (Tucumán Province, NW Argentina). The influence of the regional geological setting on water chemistry allows to distinguish three sub-basins. The interaction with sedimentary rocks was found to be dominant in the northern and central-eastern sub-basin as well as in the southern sub-basin, whereas the metamorphic-granitic basement of the Pampean Ranges was noted in the central sub-basin. In addition, anthropogenic activities affect the spatial variation of K, P, Mn, Rb, and Pb as well as dissolved oxygen concentrations and Eh. Temporal water-quality variation is related to the spatial distribution of precipitation and to the seasonal character of the main local industries (sugar cane, alcohol, citrus), increasing P and K concentrations and decreasing dissolved oxygen concentration and Eh in winter. Cl^-, Na, SO₄^2-, Al, As, B, Fe, Mn, Se and U concentrations exceed the regulated drinking-water thresholds at several sampling stations.     

Separate lead isotope analyses of leachate and residue rock fractions: implications for metal source tracing in ore deposit studies

Separate lead isotope analyses of leachate and residue fractions are applied to a broad spectrum of rocks commonly investigated in metallogenic studies. Resulting data highlight a systematic behavior of leachate and residue fractions with respect to lead isotope compositions, which essentially depends on the mineralogical composition of the rock. Granitoid and high-grade metamorphic rocks have residue compositions virtually identical to common lead. In contrast, low-grade metasedimentary rocks may have residue compositions swamped by radiogenic lead of leach-resistant zircons. Mafic magmatic rocks have residues that are often more radiogenic than leachates, depending on the ratio of leach-refractory zircons to common lead in the residual fraction of these rocks. Separate leachate and residue analyses of source rocks provide two lead isotope end members whose mixture may represent lead with the appropriate ore fluid composition. Our leaching experiments indicate that hot acid solutions (and by inference hydrothermal fluids) preferentially leach radiogenic lead from medium- to high-grade metamorphic and granitoid rocks, whereas they preferentially leach common lead from low-grade metasedimentary and mafic magmatic rocks. The method presented in this study provides a reliable alternative to other methods (i.e., age-correction of bulk-rock compositions) for the determination of the common lead signature of felsic to intermediate magmatic rocks. This may be preferable to age-corrected bulk-rock analyses, where ages to apply for corrections of bulk-rock data are not known or where moderately to highly altered rocks must be used. Case studies of orogenic gold and MVT districts of Peru (Pataz and San Vicente, respectively) show that separate leachate and residue lead isotope analyses carried out systematically on whole rocks allow a more thorough evaluation of metal source reservoirs than does the standard method of age-corrected or uncorrected bulk-rock analyses.     

Predictive GIS-Based Model of Rockfall Activity in Mountain Cliffs

Rockfall susceptibility has been analysed in mountain cliffs of the Cantabrian Range, North Spain. The main aim of this analysis has been to build a predictive model of rockfall activity from a low number of environmental and geological variables. The rockfall activity has been quantified in a GIS. The cartographic information used shows the spatial distribution of all the recent talus screes as well as their associated source areas in the rock-slopes. The area relation At/Ar (recent talus scree polygon/source basins) in the rock slopes has been used as the rockfall activity indicator. This relation has been validated in 50 pilot rock-slopes and compared with the relation number of recent rock fragments/source basin, obtained from field work. The environmental factors causing rockfall depend on the rock slope situation, and these are: altitude and sun radiation on the rock cliff. The geological factors considered are: lithology, relative position of the main discontinuities with respect to the topographic surface and two morphologic parameters: the roughness and slope gradient. A logistic regression analysis has been applied to a population of 442 limestone and quartzite rock cliffs. The dependent variable is the rockfall activity indicator, which allows the definition of two classes of rock cliff units: low and high activity. The independent variables are altitude, sun radiation (equinox radiation, summer solstice radiation, winter solstice radiation), slope roughness, slope gradient,anisotropy and lithology. Results suggest that it is possible tobuild a valid cartographic predictive model for rockfall activity in mountain rock cliffs from a limited number of easily obtainable variables. The method is especially applicable in massive rock slopes or in regions with uniform rock mass characteristics.     

Architecture of fault zones determined from outcrop, cores, 3-D seismic tomography and geostatistical modeling: example from the Albalá Granitic Pluton, SW Iberian Variscan Massif

The 3-D seismic tomographic data are used together with field, core and well log structural information to determine the detailed 3-D architecture of fault zones in a granitic massif of volume 500×575×168 m at Mina Ratones area in the Albalá Granitic Pluton. To facilitate the integration of the different data, geostatistical simulation algorithms are applied to interpolate the relatively sparse structural (hard) control data conditioned to abundant but indirect 3-D (soft) seismic tomographic data. To effectively integrate geologic and tomographic data, 3-D migration of the velocity model from the time domain into the depth domain was essential. The resulting 3-D model constitutes an image of the fault zone architecture within the granitic massif that honours hard and soft data and provides an evaluation of the spatial variability of structural heterogeneities based on the computation of 3-D experimental variograms of Fracture Index (fault intensity) data. This probabilistic quantitative 3-D model of spatially heterogeneous fault zones is suitable for subsequent fluid flow simulations. The modeled image of the 3-D fault distribution is consistent with the fault architecture in the Mina Ratones area, which basically consists of two families of subvertical structures with NNE–SSW and ENE–WSW trends that displaces the surfaces of low-angle faults (North Fault) and follows their seismically detected staircase geometry. These brittle structures cut two subvertical dykes (27 and 27′ Dykes) with a NNE–SSW to N–S trend. The faults present high FI (FI>12) adjacent bands of irregular geometry in detail that intersect in space delimiting rhombohedral blocks of relatively less fractured granite (FI<6). Both structural domains likely correspond with the protolith and the damaged zone/fault core in the widely accepted model for fault zone architecture. Therefore, the construction of 3-D grids of the FI in granitic areas affected by brittle tectonics permits the quantitative structural characterization of the rock massif.     

Validation of Landslide Susceptibility Maps; Examples and Applications from a Case Study in Northern Spain

A procedure for validating landslide susceptibility maps wasapplied in a study area in northern Spain and the results obtained compared. Validationwas used to carry out sensitivity analysis for individual variables and combinationsof variables. The validity of different map-making methods was tested, as well as theutility of different types of Favourability Functions. The results obtained show thatvalidation is essential to determine the predictive value of susceptibility maps. Italso helps to better select the most suitable function and significant variables, thus improving the efficiency of the mapping process. Validation based on a temporal strategy makes it possible to derive hazard maps from susceptibility maps.     

Strontium isotope evolution of Late Permian and Triassic seawater

The ⁸⁷Sr/⁸⁶Sr values based on brachiopods and conodonts define a nearly continuous record for the Late Permian and Triassic intervals. Minor gaps in measurements exist only for the uppermost Brahmanian, lower part of the Upper Olenekian, and Middle Norian, and only sparse data are available for the Late Permian. These 219 measurements include 67 brachiopods and 114 conodont samples from the Tethyan realm as well as 37 brachiopods and one conodont sample from the mid-European Middle Triassic Muschelkalk Sea. The Late Permian/Lower Triassic interval is characterized by a steep 1.3 × 10⁻³ rise, from 0.7070 at the base of the Dzhulfian to 0.7082 in the late Olenekian, a rate of change comparable to that in the Cenozoic. In the mid-Triassic (Anisian and Ladinian), the isotope values fall to 0.7075, followed again by a rise to 0.7081 in the Middle/Late Norian. The ⁸⁷Sr/⁸⁶Sr values decline again in the Late Norian (Sevatian) and Rhaetian to 0.7076.The sharp rise in the ⁸⁷Sr/⁸⁶Sr values during the Late Permian/Early Triassic was coincident with widespread clastic sedimentation. Because of the paucity of tectonic uplifts, the enhanced erosion may have been due to intermittent humid phases, during mainly an arid interval, coupled with the absence of a dense protective land plant cover following the mass extinction during the latest Permian. The apex of the ⁸⁷Sr/⁸⁶Sr curve at the Olenekian/Anisian boundary coincides with cessation of the large-scale clastic sedimentation and also marks the final recovery of land vegetation, as indicated by the renewed onset of coal formation in the Middle Triassic. The rising ⁸⁷Sr/⁸⁶Sr values from the Middle Carnian to the Late Norian coincide with the uplift and erosion of the Cimmeride-Indosinian orogens marking the closure of the Palaeotethys. The subsequent Rhaetian decline that continues into Jurassic (Pliensbachian/Toarcian boundary), on the other hand, coincides with the opening of the Vardar Ocean and its eastern continuation in the Izmir-Ankara Ophiolitic Belt.Samples from the Upper Muschelkalk are more radiogenic than the global trend. This may reflect separation of the basin from the open ocean. Due to strong meteoric influx from a large land mass in the north, the Germanic Basin became increasing brackish up section in the north and east, but because of the high evaporation rates, the salt content was not much reduced in the southern and central basin where a rich, but increasingly endemic, marine fauna survived.     

Impact of waves on the sea drag: measurements in the Baltic sea and a model interpretation

Measurements from the Baltic Sea and a wind-over-wave coupled model are used to study the wave impact on the sea drag. The study has been carried out for different wave conditions, namely a pure wind-sea, following-swell/ mixed sea and cross-swell/ mixed sea. Measurements reveal the fact that the sea drag is dependent on the sea-state. In stationary conditions and in the absence of severe cross-swell, swell reduces drag compared to wind-sea at the same wind speed. The cross-swell enhances the drag as compared to the following-swell case and the magnitude of the drag coefficient is increased with increasing the angle of swell propagation to the wind. It is shown that the agreement between the model results and measurements is good for pure wind-sea and stationary mixed-sea cases. Discrepancies occur at light winds, where most of the data represent pure swell conditions. During these pure swell conditions the data are characterized by a large variation of the drag coefficient. The variation is caused by mesoscale variability in the stress co-spectra, wind-cross-swell effects and nonstationarity in the wave and wind fields not represented in the model.     

In-situ Formation of Light-Absorbing Organic Matter in Cloud Water

Current climate models seem to underestimate the flux of solar energy absorbed by the global troposphere. All of these models are constrained with the assumption that cloud droplets consist of pure water. Here we demonstrate in a simple laboratory experiment that aromatic hydroxy-acids which are found in continental fine aerosol can react with hydroxyl radicals under typical conditions prevalent in cloud water influenced by biomass burning. The reactions yield colored organic species which do absorb solar radiation. We also suggest that the products of such reactions may be humic-like substances whose presence in continental aerosol has been confirmed but their source mechanisms are still much sought after. We also attempt to give a first order estimate of the enhancement of water absorption at a visible wavelength under atmospheric conditions.     

Effects of early season defoliation ofAnthyllis cytisoides(a Mediterranean browse species) on further herbivore attack

To determine if defoliation of a woody plant affects foraging by folivorous insects, we examined the infection rate (number of leaves damaged per total number of leaves sampled on marked stems) ofAnthyllis cytisoidesunder three experimental treatments: 10, 50 and 90% plant defoliation. Observations were made for three age classes, established by trunk base perimeter (equal to or lower than 11 cm, between 11 and 20 cm, larger than 20 cm). Plants respond positively to artificial defoliation by increasing total vegetative length of the stem and total inflorescence length. This response is most evident in young individuals.Response to herbivory was measured as overall infection rate and also as infection rate by different feeding guilds—chewing, mining, or sucking insects. We found that increased defoliation elicited increased resistance of leaves to insect attack. This was particularly evident in young plants. Different insect guilds respond in different ways. Attack by chewing insects declines with defoliation for all plant age classes; only sucking insects which feed on the oldest plants reduce feeding rate with plant defoliation. Finally, mining insects present the opposite trend in young and senescent plants.     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.