A hydrogeological conceptual model of the Suio hydrothermal area (central Italy)

A hydrogeological conceptual model has been developed that describes the hydrothermal system of Suio Terme (central Italy). The studied area is located along the peri-Tyrrhenian zone of the central Apennines, between the Mesozoic and Cenozoic carbonate platform sequences of the Aurunci Mountains and the volcanic sequences of the Roccamonfina. A multi-disciplinary approach was followed, using new hydrogeological surveys, the interpretation of stratigraphic logs of boreholes and water wells, and geophysical data—seismic sections, shear-wave velocity (Vs) crustal model and gravimetric model. The collected information allowed for construction of a conceptual hydrogeological model and characterization of the hydrothermal system. The Suio hydrothermal system is strongly influenced by the Eastern Aurunci hydrostructure. Along the southeastern side, the top of the hydrostructure sinks to ?1,000 m relative to sea level via a series of normal faults which give origin to the Garigliano graben. Geological and hydrogeological data strongly suggest the propagation and mixing of hot fluids, with cold waters coming from the shallow karst circuit. The aquitard distribution, the normal tectonic displacements and the fracturing of the karst hydrostructure strongly influence the hydrothermal basin. Carbon dioxide and other gasses play a key role in the whole circuit, facilitating the development of the hydrothermal system. The current level of knowledge suggests that the origin of the Suio hydrothermalism is the result of interaction between the carbonate reservoir of the Eastern Aurunci Mountains and the hot and deep crust of this peri-Tyrrhenian sector, where the Roccamonfina volcano represents the shallowest expression.     

Sensitivity of simulated flow fields and bathymetries in meandering channels to the choice of a morphodynamic model

Morphodynamic models are used by river practitioners and scientists to simulate geomorphic change in natural and artificial river channels. It has long been recognized that these models are sensitive to the choice of parameter values, and proper calibration is now common practice. This paper investigates the less recognized impact of the choice of the model itself. All morphodynamic models purport to simulate the same flow and sediment dynamics, often relying on the same governing equations. Yet in solving these equations, the models have different underlying assumptions, for example regarding spatial discretization, turbulence, sediment inflow, lateral friction, and bed load transport. These differences are not always considered by the average model user, who might expect similar predictions from calibrated models. Here, a series of numerical simulations in meandering channels was undertaken to test whether six morphodynamic codes (BASEMENT, CCHE‐2D, NAYS, SSIIM‐1, TELEMAC‐2D and TELEMAC‐3D) would yield significantly different equilibrium bathymetries if subjected to identical, initial flow conditions. We found that, despite producing moderately similar velocity patterns on a fixed‐flat bed (regression coefficient r of 0.77 ± 0.20), the codes disagree substantially with respect to simulated bathymetries (r = 0.49 ± 0.31). We relate these discrepancies to differences in the codes' assumptions. Results were configuration specific, i.e. codes that perform well for a given channel configuration do not necessarily perform well with higher or lower sinuosity configurations. Finally, limited correlation is found between accuracy and code complexity; the inclusion of algorithms that explicitly account for the effects of local bed slope and channel curvature effects on transport magnitude and direction does not guarantee accuracy. The range of solutions obtained from the evaluated codes emphasizes the need for carefully considering the choice of code. We recommend the creation of a central repository providing universal validation cases and documentation of recognized fluvial codes in commonly studied fluvial settings. Copyright © 2016 John Wiley & Sons, Ltd.     

Feeding habits of Solea senegalensis in earthen ponds in Sado estuary

The senegale sole, Solea senegalensis Kaup, 1858, is a commercially valuable fish and it is beginning to be cultured in a semi-intensive way in fish farms in southern European countries. The present study was initiated to investigate the diet of S. senegalensis under semi-extensive conditions. The feeding habits of sole were studied in two earthen ponds of a fish farm in the Sado estuary. In one earthen pond, artificial fish ration was given. S. senegalensis feeds on few prey items, its diet is mainly composed of insect larvae (Chironomus salinarus) and polychaeta (Hediste diversicolor). The diet composition of this species suggests feeding specialization, by consuming mainly annelids and insect larvae and by avoiding other items, extremely abundant in the environment, such as gastropods. In the water reservoir where ration was given, some fish consumed simultaneously benthic organism together with ration. However, benthic organisms seem to be the most important component of S. senegalensis diet.     

Remarkable variability in dyke features at the Vicuña Pampa Volcanic Complex,Southern Central Andes

Dykes at the Vicuña Pampa Volcanic Complex, which are mostly basaltic (trachy)‐andesite and (trachy)‐andesite, are exposed at the base and along the walls of a large depression resulting from intense degradation. Dykes intruding stiff layers (lavas, plugs and necks) are thin, mostly dip >60° and have coherent textures, whereas dykes intruding more compliant materials (breccias and conglomerates) tend to be thicker, have lower dips and have coherent, brecciated or mixed textures (coherent and brecciated textural domains in a single or compound dyke). Single dykes with brecciated and mixed textures are only found intruding near‐surface units. Dykes with mixed textures always have sharp contacts between domains. Dykes with sinuous domain contacts and enclaves of one domain inside the other are interpreted as resulting from dyke arrest, partial cooling and reinjection of new magma. Dykes with straight domain contacts are considered to be compound dykes, with a new dyke intruding along the margins of an older, solidified one.     

Thermal expansion and dehydroxylation of phengite micas

Phengite samples (2M ₁ and 3T politypes) and a synthetic end-member muscovite specimen were studied by in situ high-temperature synchrotron radiation X-ray diffraction. The measured volume thermal expansion of 2M ₁ phengite (<α _V> ≈ 36.6 × 10⁻⁶ K⁻¹) was systematically greater than <α _V> of the 3T polytype (≈33.3 × 10⁻⁶ K⁻¹). A positive linear correlation between the average thermal expansion on (001) plane and the mean tetrahedral rotation angle at ambient condition is proposed on the ground of new measurements and literature data. Dehydroxylation processes were observed in 2M ₁, starting at 1,000 K in 3T at 800 and 945 K in synthetic muscovite. Rietveld refinements allowed a determination of structural variations upon heating of phengite samples and their dehydroxylate phases. The phengite structure expands by regularizing the tetrahedral sheet and by reducing the bond length differences between the outer and inner coordination shell of the interlayer site. The dehydroxylate phase derived from 2M ₁ is characterized by fivefold polyhedra in the low temperature form as a consequence of two OH groups reacting to form H₂O + O (residual). The dehydroxylate exhibits an increase of the cation–cation distances along the M–Or–M bonds with respect to low-temperature phengite structures. For the 3T phase, we were unable to achieve completion of dehydroxylation. The refined structural model of the dehydroxylate phase shows two hydroxyl sites, but at a short distance from one another. This result suggests that the dehydroxylation reaction did not proceed to completion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Combined geological and gravimetric mapping and modelling for an improved understanding of observed high-resolution gravity variations: a case study for the Global Geodynamics Project (GGP) station Moxa, Germany

Geodynamic observatories around the globe continuously monitor signals like gravity, tilt and strain as a function of time. However, global signals are often masked by local effects, caused by the direct surroundings of the station, including the local geological setting. This link is well established for superconducting gravimeters (SG) that observe the gravity field variations at very high resolution. An enhancement of the SG time series by the application of local correction is exemplarily shown here in a very practicable procedure for the Geodynamic Observatory Moxa, Germany. We show how the combination of geological and gravimetrical mapping and modelling around Moxa results in a significant correction of the original gravity data, as can be proven by comparison with the satellite-derived gravity field. Detailed geological mapping of the observatory surroundings, including measurements of fold axes, foliations and joints, was the basis of the present study. The fold structure in the area of interest was interpreted by geometrical 3D modelling. The complete representation of different rock types in space also provides a better understanding of the local hydro-geological situation. Using a combination of the 3D geological model with the high-resolution Bouguer map of the observatory surroundings, we developed a 3D density model. This model enables the correction of small gravity effects caused by mass variations close to the SG, in particular local hydro-geological mass changes. Thus, the procedure presented here, based on a close connection of geology and gravimetry, is a powerful tool for the reduction of local gravity effects on SG recordings. It should be applicable to SG stations worldwide, where similar hydrologically driven mass changes can be assumed.     

Marbles and carbonate rocks from central Morocco: a petrographic,mineralogical and geochemical study

Petrographic, mineralogical, and stable isotopes (δ¹³C, δ¹⁸O values) compositions were used to characterise marbles and sedimentary carbonate rocks from central Morocco, which are considered to be a likely source of ornamental and building material from Roman time to the present day. This new data set was used in the frame of an archaeometric provenance study on Roman artefacts from the town of Thamusida (Kenitra, north Morocco), to assess the potential employment of these rocks for the manufacture of the archaeological materials. A representative set of samples from marbles and other carbonate rocks (limestone, dolostone) were collected in several quarries and outcrops in the Moroccan Meseta, in a region extending from the Meknes–Khenifra alignment to the Atlantic Ocean. All the samples were studied using a petrographic, mineralogical and geochemical methods. The petrographic and minerological investigations (optical microscopy, electron microscopy, X-ray diffraction) allowed to group the carbonate rocks in limestones, foliated limestone, diagenetic breccias and dolostone. The limestones could be further grouped as mudstones, wackestones–packstones, crinoid grainstones, oolitic grainstone and floatstones. Textural differences allowed to define marbles varieties. The stable carbon and oxygen isotope composition proved to be quite useful in the discrimination of marble sources, with apparently less discriminatory potential for carbonate rocks.     

A physical analogue model to analyse interactions between tensile stresses and dissolution in carbonate slopes

A physical analogue model was developed to analyse the relationship between tensional states in rock masses, seepage and karst processes. Use was made of an experimental apparatus consisting of two hydraulic circuits realized by drilling two holes into each of two blocks of sampled limestone from the central Apennines (Italy). A static load (208.85 kg) was applied to one of the blocks in order to elicit tensile stresses within it. Physical and chemical monitoring data showed that the main process involved was temperature-dependent CaCO₃ dissolution. This process was more marked in the loaded block circuit, as Ca⁺⁺ concentration in circulated water reached 54 mg/L, whereas only 28 mg/L was reached in the unloaded one. The interaction between load and dissolution caused the observed opening of microcracks, as confirmed by further increase of water loss and by dilution in the loaded block circuit, resulting in a decrease of Ca⁺⁺ concentration. These findings were confirmed by recording additional water losses after increasing the load up to 445.05 kg. A finite difference numerical model showed that tensile stresses (max 20 kPa) within the loaded block were clustered at the intersection of the main joints with the flowpaths, thus representing points of preferential and accelerated dissolution.This study is part of a research project on analogical and numerical modelling for geological risk mitigation being conducted at CERI, Research Centre for Geological Risks of the University of Rome “La Sapienza” (Valmontone, Roma, Italy; http://www.ceri@uniroma1.it).     

Modern analogues for Miocene to Pleistocene alkali basaltic phreatomagmatic fields in the Pannonian Basin: “soft-substrate” to “combined” aquifer controlled phreatomagmatism in intraplate volcanic fields Research Article

The Pannonian Basin (Central Europe) hosts numerous alkali basaltic volcanic fields in an area similar to 200 000 km². These volcanic fields were formed in an approximate time span of 8 million years producing smallvolume volcanoes typically considered to be monogenetic. Polycyclic monogenetic volcanic complexes are also common in each field however. The original morphology of volcanic landforms, especially phreatomagmatic volcanoes, is commonly modified. by erosion, commonly aided by tectonic uplift. The phreatomagmatic volcanoes eroded to the level of their sub-surface architecture expose crater to conduit filling as well as diatreme facies of pyroclastic rock assemblages. Uncertainties due to the strong erosion influenced by tectonic uplifts, fast and broad climatic changes, vegetation cover variations, and rapidly changing fluvio-lacustrine events in the past 8 million years in the Pannonian Basin have created a need to reconstruct and visualise the paleoenvironment into which the monogenetic volcanoes erupted. Here phreatomagmatic volcanic fields of the Miocene to Pleistocene western Hungarian alkali basaltic province have been selected and compared with modern phreatomagmatic fields. It has been concluded that the Auckland Volcanic Field (AVF) in New Zealand could be viewed as a prime modern analogue for the western Hungarian phreatomagmatic fields by sharing similarities in their pyroclastic successions textures such as pyroclast morphology, type, juvenile particle ratio to accidental lithics. Beside the AVF two other, morphologically more modified volcanic fields (Pali Aike, Argentina and Jeju, Korea) show similar features to the western Hungarian examples, highlighting issues such as preservation potential of pyroclastic successions of phreatomagmatic volcanoes.     

Dispersive Fluxes Within and Over a Real Urban Canopy: A Large-Eddy Simulation Study

Boundary-Layer Meteorology - Large-eddy simulations (LES) are conducted to study the transport of momentum and passive scalar within and over a real urban canopy in the City of Boston, USA. This...