A hybrid optimization scheme for self-potential measurements due to multiple sheet-like bodies in arbitrary 2D resistivity distributions

Self-potential is a passive geophysical method that can be applied in a straightforward manner with minimum requirements in the field. Nonetheless, interpretation of self-potential data is particularly challenging due to the inherited non-uniqueness present in all potential methods. Incorporating information regarding the target of interest can facilitate interpretation and increase the reliability of the final output. In the current paper, a novel method for detecting multiple sheet-like targets is presented. A numerical framework is initially described that simulates sheet-like bodies in an arbitrary 2D resistivity distribution. A scattered field formulation based on finite differences is employed that allows the edges of the sheet to be independent of the grid geometry. A novel analytical solution for two-layered models is derived and subsequently used to validate the accuracy of the proposed numerical scheme. Lastly, a hybrid optimization is proposed that couples linear least-squares with particle-swarm optimization in order to effectively locate the edges of multiple sheet-like bodies. Through numerical and real data, it is proven that the hybrid optimization overcomes local minimal that occurs in complex resistivity distributions and converges substantially faster compared to traditional particle-swarm optimization.     

On the seismic response of a propped r.c. diaphragm wall in a saturated clay

A series of nonlinear dynamic FE simulations have been performed to investigate the seismic performance of a flexible propped retaining wall in a saturated clay. The simulations have been carried out considering different acceleration time histories at the bedrock and two different inelastic soil models: the classical elastoplastic Modified Cam Clay model and the advanced hypoplastic model for clays proposed by Ma?ín (Int J Numer Anal Methods Geomech 29:311–336, 2005) equipped with the intergranular strains extension. The results of the simulations highlight the major role played by the choice of the constitutive model for the soil on the predicted seismic response, in terms of predicted wall displacements and structural loads. In particular, the results show that a key role is played by the model ability to correctly reproduce soil dilatancy as a function of the current stress state and loading history. This has a major impact on the inelastic volumetric deformations accumulated during the undrained seismic shearing and on the development of excess pore water pressures around the excavations.

     

Atmospheric input of nitrogen into the North Sea: ANICE project overview

The aim of the atmospheric nitrogen inputs into the coastal ecosystem (ANICE) project is to improve transport–chemistry models that estimate nitrogen deposition to the sea. To achieve this, experimental and modelling work is being conducted which aims to improve understanding of the processes involved in the chemical transformation, transport and deposition of atmospheric nitrogen compounds. Of particular emphasis within ANICE is the influence of coastal zone processes. Both short episodes with high deposition and chronic nitrogen inputs are considered in the project. The improved transport–chemistry models will be used to assess the atmospheric inputs of nitrogen compounds into the European regional seas (the North Sea is studied as a prototype) and evaluate the impact of various emission reduction strategies on the atmospheric nitrogen loads. Assessment of the impact of atmospheric nitrogen on coastal ecosystems will be based on comparisons of phytoplankton nitrogen requirements, other external nitrogen inputs to the ANICE area of interest and the direct nitrogen fluxes provided by ANICE. Selected results from both the experimental and modelling components are presented here. The experimental results show the large spatial and temporal variability in the concentrations of gaseous nitrogen compounds, and their influences on fluxes. Model calculations show the strong variation of both concentrations and gradients of nitric acid at fetches of up to 25 km. Aerosol concentrations also show high temporal variability and experimental evidence for the reaction between nitric acid and sea salt aerosol is provided by size-segregated aerosol composition measured at both sides of the North Sea. In several occasions throughout the experimental period, air mass back trajectory analysis showed connected flow between the two sampling sites (the Weybourne Atmospheric Observatory on the North Norfolk coast of the UK and Meetpost Noordwijk, a research tower at 9 km off the Dutch coast). Results from the METRAS/SEMA mesoscale chemistry transport model system for one of these cases are presented. Measurements of aerosol and rain chemical composition, using equipment mounted on a commercial ferry, show variations in composition across the North Sea. These measurements have been compared to results obtained with the transport–chemistry model ACDEP which calculates the atmospheric inputs into the whole North Sea area. Finally, the results will be made available for the assessment of the impact of atmospheric nitrogen on coastal ecosystems.     

Origin and emplacement of ultramafic–mafic intrusions in the Erro-Tobbio mantle peridotite (Ligurian Alps, Italy)

The Erro-Tobbio peridotites (Voltri Massif, Ligurian Alps) represent subcontinental lithospheric mantle tectonically exhumed during Permo–Mesozoic extension of the Europe–Adria lithosphere. Previous studies have shown that exhumation started during Permian times, and occurred along kilometer-scale lithospheric shear zones which enhanced progressive deformation and recrystallization from spinel- to plagioclase-facies conditions. Ongoing field and petrologic investigations have revealed that the peridotites experienced, during uplift, a composite history of diffuse melt migration and multiple episodes of ultramafic–mafic intrusions. In this paper we present the results of field, structural and petrologic–geochemical investigations into a sector of the Erro-Tobbio peridotite unit that preserves well this multiple intrusion history. Melt impregnation in the peridotites is evidenced by significant plagioclase enrichment and crystallization of unstrained orthopyroxene replacing kinked mantle olivine and clinopyroxene; impregnating melts were thus opx-saturated. Melt–rock interaction caused chemical changes in mantle minerals (e.g. Al decrease and REE increase in cpx; Ti and Cr# enrichment in spinel). Nevertheless, clinopyroxenes still exhibit LREE depletion (Ce_N/Sm_N = 0.006–0.011), indicating a depleted signature for the percolating melts. Melt impregnation was thus related to diffuse porous flow migration of depleted MORB-type melt fractions that modified their compositions towards opx saturation by mantle–melt interaction during ascent. The impregnated peridotites are intruded by a hectometer-scale stratified cumulate body, mostly consisting of troctolites and plagioclase wehrlites, showing gradational, interfingered contacts with the host mantle rocks. Subsequent intrusion events are revealed by the occurrence of olivine gabbros as decameter-wide lenses, variably thick (centimeter- to meter-scale) dykes and thin dykelets, which crosscut both the peridotite foliation and the magmatic layering in the cumulates. Overall, major and trace element compositions of minerals in the intrusives indicate that they represent variably differentiated cumulus products crystallized from rather primitive N-MORB-type aggregated melts. Slightly more evolved compositions are shown by olivine gabbros, relative to the troctolites and plagioclase wehrlites of the cumulate body. Mineral chemistry features (e.g. the Fo–An correlation and high Na, Ti, Mg# in cpx) indicate that the studied intrusive rocks crystallized at moderate pressure conditions (3–5 kbar, i.e. 9–15 km depth). Our study thus points to a progressive transition from porous flow melt migration to emplacement of magmas in fractures, presumably related to progressive change of lithospheric mantle rheology during extension-related uplift and cooling.     

Water-table and discharge changes associated with the 2016–2017 seismic sequence in central Italy: hydrogeological data and a conceptual model for fractured carbonate aquifers

A seismic sequence in central Italy from August 2016 to January 2017 affected groundwater dynamics in fractured carbonate aquifers. Changes in spring discharge, water-table position, and streamflow were recorded for several months following nine Mw 5.0–6.5 seismic events. Data from 22 measurement sites, located within 100 km of the epicentral zones, were analyzed. The intensity of the induced changes were correlated with seismic magnitude and distance to epicenters. The additional post-seismic discharge from rivers and springs was found to be higher than 9 m³/s, totaling more than 0.1 km³ of groundwater release over 6 months. This huge and unexpected contribution increased streamflow in narrow mountainous valleys to previously unmeasured peak values. Analogously to the L’Aquila 2009 post-earthquake phenomenon, these hydrogeological changes might reflect an increase of bulk hydraulic conductivity at the aquifer scale, which would increase hydraulic heads in the discharge zones and lower them in some recharge areas. The observed changes may also be partly due to other mechanisms, such as shaking and/or squeezing effects related to intense subsidence in the core of the affected area, where effects had maximum extent, or breaching of hydraulic barriers.

     

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.     

Impact cratering in H2O‐bearing targets on Mars: Thermal field under craters as starting conditions for hydrothermal activity

Abstract– We present a case modeling study of impact crater formation in H₂O‐bearing targets. The main goal of this work was to investigate the postimpact thermal state of the rock layers modified in the formation of hypervelocity impact craters. We present model results for a target consisting of a mixture of H₂O‐ice and rock, assuming an ice/water content variable with depth. Our model results, combined with results from previous work using dry targets, indicate that for craters larger than about 30 km in diameter, the onset of postimpact hydrothermal circulation is characterized by two stages: first, the formation of a mostly dry, hot central uplift followed by water beginning to flow in and circulate through the initially dry and hot uplifted crustal rocks. The postimpact thermal field in the periphery of the crater is dependent on crater size: in midsize craters, 30–50 km in diameter, crater walls are not strongly heated in the impact event, and even though ice present in the rock may initially be heated enough to melt, overall temperatures in the rock remain below melting, undermining the development of a crater‐wide hydrothermal circulation. In large craters (with diameters more than 100 km or so), the region underneath the crater floor and walls is heated well above the melting point of ice, thus facilitating the onset of an extended hydrothermal circulation. These results provide preliminary constraints in characterizing the many water‐related features, both morphologic and spectroscopic, that high‐resolution images of Mars are now detecting within many Martian craters.