Origin of peraluminous granites and granodiorites, Iberian massif, Spain: an experimental test of granite petrogenesis

The discrimination between potential source materials involved in the genesis of Iberian granites and granodiorites, as well as the role of mantle-crust interactions, are examined using constraints imposed by melting experiments, melting-assimilation experiments and Sr-Nd isotope systematics. The Sr-Nd isotope relationships indicate the existence of different genetic trends in which juvenile mantle materials are involved by different mechanisms: (1) a source trend, traced by a particular evolution of the pre-Hercynian basement and indicating mantle participation at the time of sedimentation; (2) a set of magmatic trends traced by gabbro-tonalite-enclave-granodiorite associations, implying the incorporation of new mantle material at the time of granite generation. These relationships strongly support a pure crustal origin for the peraluminous leucogranites, derived from partial melting of crustal protoliths, and a hybrid origin for the peraluminous granodiorites. These granodiorites are the most abundant granitic rocks of the Central Iberian zone (CIZ) of the Iberian massif, implying that processes of hybridisation by assimilation and/or magma mixing played an important role in granitoid production during the Hercynian orogeny. These hypotheses have been tested by means of melting and assimilation experiments. Melting experiments in the range 800–900 °C and at pressures of 3, 6, 10 and 15 kbar indicate that: (1) several potential source materials such as Bt-Ms gneisses and metagreywackes are suitable for the production of peraluminous leucogranite melts; (2) the melt compositions are always leucogranitic, regardless of pressure; (3) pressure exerts a strong influence on the fertility of the source: experiments at 3 kbar produce more than 20 vol% of melt, compared with less than 5 vol% of melt produced at 10 and 15 kbar and at the same temperature. The melting-assimilation experiments carried out at 1000 °C and 4, 7 and 10 kbar and using a proportion of 50% gabbro and 50% gneiss give high melt proportions (more than 50 vol.%) and noritic residues. These melts have the composition of leucogranodiorites, and overlap with part of the compositional range of peraluminous granodiorites of the Iberian massif. The generation of more mafic granodiorites may be explained by the incorporation of some residual orthopyroxene to the granodiorite magmas. The low solubility of Fe + Mg prevents the generation of granodiorite melts with more than 3 wt% of MgO + FeO at all crustal pressures. The large volumes of peraluminous, hybrid granodiorites, produced by assimilation of crustal rocks by mantle magmas, imply that an important episode of crustal growth took place during the Late-Palaeozoic Hercynian orogeny in the Iberian massif. Received: 30 June 1998 / Accepted: 27 November 1998     

Frontal dynamics boost primary production in the summer stratified Mediterranean sea

Bio-physical glider measurements from a unique process-oriented experiment in the Eastern Alboran Sea (AlborEx) allowed us to observe the distribution of the deep chlorophyll maximum (DCM) across an intense density front, with a resolution (～ 400 m) suitable for investigating sub-mesoscale dynamics. This front, at the interface between Atlantic and Mediterranean waters, had a sharp density gradient (Δρ ～ 1 kg/m³ in ～ 10 km) and showed imprints of (sub-)mesoscale phenomena on tracer distributions. Specifically, the chlorophyll-a concentration within the DCM showed a disrupted pattern along isopycnal surfaces, with patches bearing a relationship to the stratification (buoyancy frequency) at depths between 30 and 60 m. In order to estimate the primary production (PP) rate within the chlorophyll patches observed at the sub-surface, we applied the Morel and Andrè (J Geophys Res 96:685–698 1991) bio-optical model using the photosynthetic active radiation (PAR) from Argo profiles collected simultaneously with glider data. The highest production was located concurrently with domed isopycnals on the fresh side of the front, suggestive that (sub-)mesoscale upwelling is carrying phytoplankton patches from less to more illuminated levels, with a contemporaneous delivering of nutrients. Integrated estimations of PP (1.3 g C m^?2d^?1) along the glider path are two to four times larger than the estimations obtained from satellite-based algorithms, i.e., derived from the 8-day composite fields extracted over the glider trip path. Despite the differences in spatial and temporal sampling between instruments, the differences in PP estimations are mainly due to the inability of the satellite to measure DCM patches responsible for the high production. The deepest (depth > 60 m) chlorophyll patches are almost unproductive and probably transported passively (subducted) from upper productive layers. Finally, the relationship between primary production and oxygen is also investigated. The logarithm of the primary production in the DCM interior (chlorophyll (Chl) > 0.5 mg/m³) shows a linear negative relationship with the apparent oxygen utilization, confirming that high chlorophyll patches are productive. The slope of this relationship is different for Atlantic, mixed interface waters and Mediterranean waters, suggesting the presence of differences in planktonic communities (whether physiological, population, or community level should be object of further investigation) on the different sides of the front. In addition, the ratio of optical backscatter to Chl is high within the intermediate (mixed) waters, which is suggestive of large phytoplankton cells, and lower within the core of the Atlantic and Mediterranean waters. These observations highlight the relevance of fronts in triggering primary production at DCM level and shaping the characteristic patchiness of the pelagic domain. This gains further relevance considering the inadequacy of optical satellite sensors to observe DCM concentrations at such fine scales.     

Numerical modeling of storm surges in the coast of Mozambique: the cases of tropical cyclones Bonita (1996) and Lisette (1997)

The coast of Mozambique is often affected by storms, particularly tropical cyclones during summer or sometimes midlatitude systems in the southern part. Storm surges combined with high freshwater discharge can drive huge coastal floods, affecting both urban and rural areas. To improve the knowledge about the impact of storm surges in the coast of Mozambique, this study presents the first attempt to model this phenomenon through the implementation of the Princeton Ocean Model (POM) in the Southwestern Indian Ocean domain (SWIO; 2–32°S, 28–85°E) using a regular grid with 1/6° of spatial resolution and 36 sigma levels. The simulation was performed for the period 1979–2010, and the most interesting events of surges were related to tropical cyclones Bonita (1996) and Lisette (1997) that occurred in the Mozambique Channel. The results showed that the model represented well the amplitude and phase of principal lunar and solar tidal constituents, as well as it captured the spatial pattern and magnitudes of SST with slight positive bias in summer and negative bias in winter months. In terms of SSH, the model underestimated the presence of mesoscale eddies, mainly in the Mozambique Channel. Our results also showed that the atmospheric sea level pressure had a significant contribution to storm heights during the landfall of the tropical cyclones Bonita (1996) and Lisette (1997) in the coast of Mozambique contributing with about 20 and 16% of the total surge height for each case, respectively, surpassing the contribution of the tide-surge nonlinear interactions by a factor of 2.     

A Method to Estimate the Hydraulic Conductivity of the Ground by TRT Analysis

The knowledge of hydraulic properties of aquifers is important in many engineering applications. Careful design of ground‐coupled heat exchangers requires that the hydraulic characteristics and thermal properties of the aquifer must be well understood. Knowledge of groundwater flow rate and aquifer thermal properties is the basis for proper design of such plants. Different methods have been developed in order to estimate hydraulic conductivity by evaluating the transport of various tracers (chemical, heat etc.); thermal response testing (TRT) is a specific type of heat tracer that allows including the hydraulic properties in an effective thermal conductivity value. Starting from these considerations, an expeditious, graphical method was proposed to estimate the hydraulic conductivity of the aquifer, using TRT data and plausible assumption. Suggested method, which is not yet verified or proven to be reliable, should be encouraging further studies and development in this direction.     

SlugIn 1.0: A Free Tool for Automated Slug Test Analysis

The correct characterization of aquifer parameters is essential for water‐supply and water‐quality investigations. Slug tests are widely used for these purposes. While free software is available to interpret slug tests, some codes are not user‐friendly, or do not include a wide range of methods to interpret the results, or do not include automatic, inverse solutions to the test data. The private sector has also generated several good programs to interpret slug test data, but they are not free of charge. The computer program SlugIn 1.0 is available online for free download, and is demonstrated to aid in the analysis of slug tests to estimate hydraulic parameters. The program provides an easy‐to‐use Graphical User Interface. SlugIn 1.0 incorporates automated parameter estimation and facilitates the visualization of several interpretations of the same test. It incorporates solutions for confined and unconfined aquifers, partially penetrating wells, skin effects, shape factor, anisotropy, high hydraulic conductivity formations and the Mace test for large‐diameter wells. It is available in English and Spanish and can be downloaded from the web site of the Geological Survey of Spain. Two field examples are presented to illustrate how the software operates.     

Storm surge frequency reduction in Venice under climate change

Increased tidal levels and storm surges related to climate change are projected to result in extremely adverse effects on coastal regions. Predictions of such extreme and small-scale events, however, are exceedingly challenging, even for relatively short time horizons. Here we use data from observations, ERA-40 re-analysis, climate scenario simulations, and a simple feature model to find that the frequency of extreme storm surge events affecting Venice is projected to decrease by about 30% by the end of the twenty-first century. In addition, through a trend assessment based on tidal observations we found a reduction in extreme tidal levels. Extrapolating the current +17 cm/century sea level trend, our results suggest that the frequency of extreme tides in Venice might largely remain unaltered under the projected twenty-first century climate simulations.     

Assessment of calanchi and rill–interrill erosion susceptibilities using terrain analysis and geostochastics: A case study in the Oltrepo Pavese,Northern Apennines,Italy

Soil erosion is one of the most important environmental problems distributed worldwide. In the last decades, numerous studies have been published on the assessment of soil erosion and the related processes and forms using empirical, conceptual and physically based models. For the prediction of the spatial distribution, more and more sophisticated stochastic modelling approaches have been proposed – especially on smaller spatial scales such as river basins. In this work, we apply a maximum entropy model (MaxEnt) to evaluate badlands (calanchi) and rill–interrill (sheet erosion) areas in the Oltrepo Pavese (Northern Apennines, Italy). The aim of the work is to assess the important environmental predictors that influence calanchi and rill–interrill erosion at the regional scale. We used 13 topographic parameters derived from a 12 m digital elevation model (TanDEM-X) and data on the lithology and land use. Additional information about the vegetation is introduced through the normalized difference vegetation index based on remotely sensed data (ASTER images). The results are presented in the form of susceptibility maps showing the spatial distribution of the occurrence probability for calanchi and rill–interrill erosion. For the validation of the MaxEnt model results, a support vector machine approach was applied. The models show reliable results and highlight several locations of the study area that are potentially prone to future soil erosion. Thus, coping and mitigation strategies may be developed to prevent or fight the soil erosion phenomenon under consideration. © 2020 John Wiley & Sons, Ltd.     

Misinterpretation of the Kenessey method for the determination of the runoff coefficient: a review

The calculation of runoff and infiltration parameters is a fundamental task for water budget estimation and water resources management. The “Kenessey method” allows these values to be calculated and its use is widespread among water scientists and practitioners. Due to language barriers and inefficient translation in past decades, the method has undergone some misinterpretations and is now far from its original usefulness and purpose. The use of the modified method can give unreal results for specific climatic zones. The purpose of this paper is to point out these methodological variations and to return to the original approach, improving it with modern technology.     

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.     

Methods for estimating the velocities of the Brazil Current in the pre-salt reservoir area off southeast Brazil (23∘ S–26∘ S)

The Brazil Current (BC) is likely the least observed and investigated subtropical western boundary current in the world. This study proposes a simple and systematic methodology to estimate quasi-synoptic cross-sectional speeds of the BC within the Santos Basin (23^° S–26^° S) based on the dynamic method using several combinations of data: Conductivity, temperature, and depth (CTD), temperature profiles, CTD and vessel-mounted Acoustic Doppler Current Profiler (VMADCP), and temperature profiles and VMADCP. All of the geostrophic estimates agree well with lowered Acoustic Doppler Current Profiler (LADCP) velocity observations and yield volume transports of -5.56 ±1.31 and 2.50 ±1.01 Sv for the BC and the Intermediate Western Boundary Current (IWBC), respectively. The LADCP data revealed that the BC flows southwestward and is ～100 km wide, 500 m deep, and has a volume transport of approximately -5.75 ±1.53 Sv and a maximum speed of 0.59 m s^?1. Underneath the BC, the IWBC flows northeastward and has a vertical extent of approximately 1,300 m, a width of ～60 km, a maximum velocity of ～0.22 m s?¹, and a volume transport of 4.11 ± 2.01 Sv. Our analysis indicates that in the absence of the observed velocities, the isopycnal (σ ₀) of 26.82 kg m^?3 (～500 dbar) is an adequate level of no motion for use in geostrophic calculations. Additionally, a simple linear relationship between the temperature and the specific volume anomaly can be used for a reliable first estimate of the BC-IWBC system in temperature-only transects.