Impact of irrigated agriculture on groundwater-recharge salinity: a major sustainability concern in semi-arid regions

Intensive irrigated agriculture substantially modifies the hydrological cycle and often has major environmental impacts. The article focuses upon a specific concern—the tendency for progressive long-term increases in the salinity of groundwater recharge derived from irrigated permeable soils and replenishment of unconfined aquifers in more arid regions. This process has received only scant attention in the water-resource literature and has not been considered by agricultural science. This work makes an original contribution by analysing, from scientific principles, how the salinisation of groundwater recharge arises and identifies the factors affecting its severity. If not proactively managed, the process eventually will impact irrigation waterwell salinity, the productivity of agriculture itself, and can even lead to land abandonment. The types of management measure required for mitigation are discussed through three detailed case histories of areas with high-value groundwater-irrigated agriculture (in Spain, Argentina and Pakistan), which provide a long-term perspective on the evolution of the problem over various decades.     

Archaeoseismological evidence of past earthquakes in Rome (fifth to ninth century A.D.) used to quantify dating uncertainties and coseismic damage

The transformation of Rome during the Late Antiquity and the Early Middle Ages has been investigated by archaeologists and historians. Social and political changes are the main aspects which led to a progressive modification of the urban framework; abandonment, spoliation and transformation of buildings are quite diffused as documented by the archaeological literature. The consequence of these practices is a higher vulnerability of the buildings which, from the seismological point of view, played a main role in increasing the effects of seismic shaking. A number of earthquakes have struck Rome during the period of investigation (fifth to ninth century A.D.), known from historical sources: 443, 484–508, 618, 801, 847; in some cases (443, 484–508, 801) damage has been documented. In contrast, the archaeological sources characterise collapse layers and evidence of destruction at different sites with changing and not always conclusive chronological constraints. Consequently, collapse and destruction have been alternatively attributed to the above-mentioned earthquakes. Through a geoarchaeological and stratigraphic analysis of potentially coseismic collapse units, we want (1) to describe the archaeoseismic evidence derived from recent excavations and from the available literature (e.g. Piazza Madonna di Loreto, Piazza Venezia, Palazzo Valentini Crypta Balbi, Colosseo, Basilica Hilariana, Basilica di Santa Petronilla, Santa Maria Antiqua,…); (2) to discuss the chronological problems and the uncertainty of attribution of the collapse units to known historical earthquakes; (3) to discuss the earthquake damage exaggeration due to erroneous attribution of seismic origin to the evidence of destruction derived from archaeological data. Finally, we will infer the role that earthquakes may have had on the development of the urban landscape in the fifth to ninth century A.D.     

Scales and causes of heterogeneity in bars in a large multi‐channel river: Río Paraná, Argentina

To date, published studies of alluvial bar architecture in large rivers have been restricted mostly to case studies of individual bars and single locations. Relatively little is known about how the depositional processes and sedimentary architecture of kilometre‐scale bars vary within a multi‐kilometre reach or over several hundreds of kilometres downstream. This study presents Ground Penetrating Radar and core data from 11, kilometre‐scale bars from the Río Paraná, Argentina. The investigated bars are located between 30 km upstream and 540 km downstream of the Río Paraná – Río Paraguay confluence, where a significant volume of fine‐grained suspended sediment is introduced into the network. Bar‐scale cross‐stratified sets, with lengths and widths up to 600 m and thicknesses up to 12 m, enable the distinction of large river deposits from stacked deposits of smaller rivers, but are only present in half the surface area of the bars. Up to 90% of bar‐scale sets are found on top of finer‐grained ripple‐laminated bar‐trough deposits. Bar‐scale sets make up as much as 58% of the volume of the deposits in small, incipient mid‐channel bars, but this proportion decreases significantly with increasing age and size of the bars. Contrary to what might be expected, a significant proportion of the sedimentary structures found in the Río Paraná is similar in scale to those found in much smaller rivers. In other words, large river deposits are not always characterized by big structures that allow a simple interpretation of river scale. However, the large scale of the depositional units in big rivers causes small‐scale structures, such as ripple sets, to be grouped into thicker cosets, which indicate river scale even when no obvious large‐scale sets are present. The results also show that the composition of bars differs between the studied reaches upstream and downstream of the confluence with the Río Paraguay. Relative to other controls on downstream fining, the tributary input of fine‐grained suspended material from the Río Paraguay causes a marked change in the composition of the bar deposits. Compared to the upstream reaches, the sedimentary architecture of the downstream reaches in the top ca 5 m of mid‐channel bars shows: (i) an increase in the abundance and thickness (up to metre‐scale) of laterally extensive (hundreds of metres) fine‐grained layers; (ii) an increase in the percentage of deposits comprised of ripple sets (to >40% in the upper bar deposits); and (iii) an increase in bar‐trough deposits and a corresponding decrease in bar‐scale cross‐strata (<10%). The thalweg deposits of the Río Paraná are composed of dune sets, even directly downstream from the Río Paraguay where the upper channel deposits are dominantly fine‐grained. Thus, the change in sedimentary facies due to a tributary point‐source of fine‐grained sediment is primarily expressed in the composition of the upper bar deposits.     

A calderas field at the Marifil Formation, new volcanogenic interpretation, Norpatagonian Massif, Argentina

Major volcanogenic structures within the Marifil Formation allowed to determine the location of a claderas field that spreads along 400 kilometers. This field has at least three large calderas with diameters near to 100 kilometers, recognizable by regional distribution of associated dikes and rhyolite lava flows present at the caldera border that separates a monotonous ignimbritic plateau from a depressed (150 meters) inner volcanic collapsed caldera. The Marifil ignimbritic plateau cover more than 50,000 square kilometers with thicknesses that reach 800 meters. Associated with these plateau ignimbrites there is an important fluorite mineralization.     

An Albian–Cenomanian unconformity in the northern Andes: Evidence and tectonic significance

The Villeta Group of Colombia and equivalent stratigraphic units of Venezuela and Ecuador comprise marine sequences ranging from Albian to Santonian in age. Deposition of the Villeta Group was presumed to take place entirely in quiet tectonic conditions in a passive margin setting that occupied NW South America. From a large database of 2D/3D seismic, well, surface geology, and biostratigraphic data, we present evidence for intra-Villeta (mostly late Albian–Cenomanian) deformation in parts of the Upper Magdalena Valley and Eastern Cordillera of Colombia, controlled by transpressional fault reactivation, produced by transpressional fault reactivation and thrusting that resulted in an angular unconformity. This event has been largely unnoticed in the literature, but previously scattered evidence supports our observations, suggesting regionally extensive tectonism. Published fission-track age determinations and other geologic evidence from Colombia and Venezuela suggest significant uplifts around 80–100 m.y., which may reflect changes in the subduction regime, with compressional deformation in certain regions and extensional deformation in others. A late Albian onset of compressional deformation along the Colombian and Peruvian segments of the Andes may be related to the opening of the South Atlantic Ocean at equatorial latitudes. Identification of tectonic activity with development of an unconformity in intra-Villeta times provides new insights into understanding the evolution of the Upper Magdalena Valley and adjacent areas of Colombia and western Venezuela and creates new possibilities for hydrocarbon exploration, with additional trapping phases, better reservoir preservation by early migration and secondary porosity, and ultimately facies changes with stratigraphic potential.     

High-pressure behaviour along the jadeite NaAlSi2O6–aegirine NaFeSi2O6 solid solution up to 10 GPa

The volume variation as a function of pressure along the jadeite–aegirine solid solution was determined at room temperature up to pressures between 6.5 and 9.7 GPa by single-crystal X-ray diffraction. The unit-cell volumes collected at room pressure for the different compositions indicate a slight deviation from linearity along the join. The pressure–volume data have been fitted using a third-order Birch-Murnaghan equation of state (BM3-EoS). The bulk modulus, K _T0, varies from 134.0(7) GPa for pure jadeite to 116.1(5) GPa for pure aegirine. Its evolution with composition along the join is not linear and can be described by the following second order polynomial:

Qanat is not a hazard

Comments are presented on the article by Abbasnejad et al. (Environ Earth Sci 75:1306, 2016) dealing with qanat and hazards in Iran. My rebuttal starts from the direct attribution of the described hazards to qanat and addresses the importance in the correct use of terminology for geological hazards. All of the problems pointed out in Iran (subsidence, sinkholes, pollution) have, to me, an anthropogenic origin and cannot be directly ascribed to qanat. Eventually, I present some additional elements to highlight the remarkable importance of qanat systems and their influence on the development of similar underground structures in many countries of the Mediterranean Basin. This remarks the relevance of qanat as cultural heritage sites and the need for their preservation and valorization.     

Uncertainty quantification of overpressure buildup through inverse modeling of compaction processes in sedimentary basins

This study illustrates a procedure conducive to a preliminary risk analysis of overpressure development in sedimentary basins characterized by alternating depositional events of sandstone and shale layers. The approach rests on two key elements: (1) forward modeling of fluid flow and compaction, and (2) application of a model-complexity reduction technique based on a generalized polynomial chaos expansion (gPCE). The forward model considers a one-dimensional vertical compaction processes. The gPCE model is then used in an inverse modeling context to obtain efficient model parameter estimation and uncertainty quantification. The methodology is applied to two field settings considered in previous literature works, i.e. the Venture Field (Scotian Shelf, Canada) and the Navarin Basin (Bering Sea, Alaska, USA), relying on available porosity and pressure information for model calibration. It is found that the best result is obtained when porosity and pressure data are considered jointly in the model calibration procedure. Uncertainty propagation from unknown input parameters to model outputs, such as pore pressure vertical distribution, is investigated and quantified. This modeling strategy enables one to quantify the relative importance of key phenomena governing the feedback between sediment compaction and fluid flow processes and driving the buildup of fluid overpressure in stratified sedimentary basins characterized by the presence of low-permeability layers. The results here illustrated (1) allow for diagnosis of the critical role played by the parameters of quantitative formulations linking porosity and permeability in compacted shales and (2) provide an explicit and detailed quantification of the effects of their uncertainty in field settings.     

Activity and kinematic behaviour of deep-seated landslides from PS-InSAR displacement rate measurements

Large landslides and deep-seated gravitational slope deformations (DSGSD) represent an important geo-hazard in relation to the deformation of large structures and infrastructures and to the associated secondary landslides. DSGSD movements, although slow (from a few millimetres to several centimetres per year), can continue for very long periods, producing large cumulative displacements and undergoing partial or complete reactivation. Therefore, it is important to map the activity of such phenomena at a regional scale. Ground surface displacements at DSGSD typically range close to the detection limit of monitoring equipment but are suitable for synthetic aperture radar (SAR) interferometry. In this paper, permanent scatterers (PSInSAR?) and SqueeSAR? techniques are used to analyse the activity of 133 DSGSD, in the Central Italian Alps. Statistical indicators for assigning a degree of activity to slope movements from displacement rates are discussed together with methods for analysing the movement and activity distribution within each landslide. In order to assess if a landslide is active or not, with a certain degree of reliability, three indicators are considered as optimal: the mean displacement rate, the activity index (ratio of active PS, displacement rate larger than standard deviation, overall PS) and the nearest neighbor ratio, which allows to describe the degree of clustering of the PS data. According to these criteria, 66% of the phenomena are classified as active in the monitored period 1992–2009. Finally, a new methodology for the use of SAR interferometry data to attain a classification of landslide kinematic behaviour is presented. This methodology is based on the interpretation of longitudinal ground surface displacement rate profiles in the light of numerical simulations of simplified failure geometries. The most common kinematic behaviour is rotational, amounting to 41 DSGSDs, corresponding to the 62.1% of the active phenomena.