Seismic tomography from local crustal earthquakes beneath eastern Rif Mountains of Morocco

We applied a tomographic method to image an aseismic strike–slip fault in North Morocco and found that the occurrence of earthquakes is not only controlled by the state of tectonic stress but also by material heterogeneity in the crust. We have constructed an integrated model of seismic, electric, magnetic and heat flow properties across northeastern Morocco primarily based on a tomography inversion of local earthquake arrival times. The seismic images obtained show a pronounced low-velocity zone at 5 km depth parallels to the Nekor fault, coinciding with an anomalously high conductive and low gravity structure, which is interpreted as a fault gouge zone and/or a fluid-filled subsurface rock matrix. Below 10 km depth, a weak positive velocity zone indicates that the fault gouge is stable. The seismicity and the seismic velocity results for the Al-Hoceimas region show that the concentrations of earthquakes are confined in the high velocity area. This anomaly is interpreted to be a brittle and competent layer of the upper crust that sustains seismogenic stress. On the eastern coast line of Morocco, we infer that a high density, high velocity body exists in the shallowest layers of the upper crust, probably formed by Miocene volcanic rocks.     

Development of water management strategies in southern Mesopotamia during the fourth and third millennium B.C.E.

The last two decades witnessed increasing scholarly interest in the history of water management in southern Mesopotamia. Thanks to many geoarchaeological research projects conducted throughout the central and southern Iraqi floodplains, a general understanding of the macrophases of anthropogenic manipulation of this vast hydraulic landscape has been achieved. However, current narratives mostly rely on studies at a regional scale and are based on excessively long chronological phases (often spanning a whole millennium). A finer-tuned analysis at a submillennial scale is needed to better appreciate the dynamics that led to the development of artificial canals and irrigation systems and the creation of harbours in cities and other navigation-related facilities. The Iraqi-Italian QADIS project is addressing this issue through a systematic geoarchaeological investigation in the south-eastern area of the Qadisiyah province. We aim to update the current narrative by analysing case studies involving specific periods of occupation. We performed 17 boreholes to propose a date on the functioning period of the hydraulic works in five selected archaeological sites of this region. This approach allowed us to understand changes in water management strategies in both the short and the medium term (i.e., on a scale of centuries). In this paper, we present the results for the fourth and third millennia B.C.E. This period witnessed a crucial passage from the basic exploitation of natural watercourses for irrigation and occasional navigation to the emergence of the first system of artificial canals and intraurban harbours.     

Surface-wave analysis for static corrections in mineral exploration: A case study from central Sweden

In mineral exploration, increased interest towards deeper mineralizations makes seismic methods attractive. One of the critical steps in seismic processing workflows is the static correction, which is applied to correct the effect of the shallow, highly heterogeneous subsurface layers, and improve the imaging of deeper targets. We showed an effective approach to estimate the statics, based on the analysis of surface waves (groundroll) contained in the seismic reflection data, and we applied it to a legacy seismic line acquired at the iron-oxide mining site of Ludvika in Sweden. We applied surface-wave methods that were originally developed for hydrocarbon exploration, modified as a step-by-step workflow to suit the different geologic context of hard-rock sites. The workflow starts with the detection of sharp lateral variations in the subsurface, the existence of which is common at hard-rock sites. Their location is subsequently used, to ensure that the dispersion curves extracted from the data are not affected by strong lateral variations of the subsurface properties. The dispersion curves are picked automatically, windowing the data and applying a wavefield transform. A pseudo-2D time-average S-wave velocity and time-average P-wave velocity profile are obtained directly from the dispersion curves, after inverting only a reference curve. The time-average P-wave velocity profile is then used for the direct estimation of the one-way traveltime, which provides the static corrections. The resulting P-wave statics from the field data were compared with statics computed through conventional P-wave tomography. Their difference was mostly negligible with more than 91% of the estimations being in agreement with the conventional statics, proving the effectiveness of the proposed workflow. The application of the statics obtained from surface waves provided a stacked section comparable with that obtained by applying tomostatics.     

NanoSr – A New Carbonate Microanalytical Reference Material for In Situ Strontium Isotope Analysis

The in situ measurement of Sr isotopes in carbonates by MC‐ICP‐MS is limited by the availability of suitable microanalytical reference materials (RMs), which match the samples of interest. Whereas several well‐characterised carbonate reference materials for Sr mass fractions > 1000 µg g^?1 are available, there is a lack of well‐characterised carbonate microanalytical RMs with lower Sr mass fractions. Here, we present a new synthetic carbonate nanopowder RM with a Sr mass fraction of ca. 500 µg g^?1 suitable for microanalytical Sr isotope research (‘NanoSr’). NanoSr was analysed by both solution‐based and in situ techniques. Element mass fractions were determined using EPMA (Ca mass fraction), as well as laser ablation and solution ICP‐MS in different laboratories. The ⁸⁷Sr/⁸⁶Sr ratio was determined by well‐established bulk methods for Sr isotope measurements and is 0.70756 ± 0.00003 (2s). The Sr isotope microhomogeneity of the material was determined by LA‐MC‐ICP‐MS, which resulted in ⁸⁷Sr/⁸⁶Sr ratios of 0.70753 ± 0.00007 (2s) and 0.70757 ± 0.00006 (2s), respectively, in agreement with the solution data within uncertainties. Thus, this new reference material is well suited to monitor and correct microanalytical Sr isotope measurements of low‐Sr, low‐REE carbonate samples. NanoSr is available from the corresponding author.     

An individual-based population dynamic model for estimating biomass yield and nutrient fluxes through an off-shore mussel (Mytilus galloprovincialis) farm

The fluxes of carbon, nitrogen and phosphorus through an off-shore long-line Mytilus galloprovincialis farm during a typical rearing cycle were estimated by combining a simple population dynamic model, based on a new individual model, and a set of field data, concerning the composition of the seston, as well as that of mussel meat and faeces. The individual model, based on an energy budget, was validated against a set of original field data, which were purposely collected from July 2006 to May 2007 in the North-Western Adriatic Sea (Italy) and was further tested using historical data. The model was upscaled to the population level by means of a set of Monte Carlo simulations, which were used for estimating the size structure of the population. The daily fluxes of C, N and P associated with mussel filtration, excretion and faeces and pseudo-faeces production were integrated over the 10-month-long rearing cycle and compared with the total amount of C, N and P removed by harvesting. The results indicate that the individual model compares well with an existing literature model and provides reliable estimations of the growth of mussel specimen over a range of trophic conditions which are typical of the Northern Adriatic Sea coastal area. The results of the budget calculation indicate that, even though the harvest represents a net removal of phosphorus and nitrogen from the ecosystem, the mussel farm increases the retention time of both nutrients in the coastal area, via the deposition of faeces and pseudo-faeces on the sea-bed. In fact, the amount of nitrogen associated with deposition is approximately twice the harvested one and the amount of phosphorus is approximately five times higher. These findings are in qualitative agreement with the results of literature budget and model calculations carried out in a temperate coastal embayment. This agreement suggests that the proper assessment of the overall effect of long-line mussel farming on both the benthic and pelagic ecosystem asks for an integrated modelling approach, which should include the dynamic of early diagenesis processes, as well as of that of nutrients released from the surface sediment.     

Sequence stratigraphy based on high-resolution seismic profiles in the late Pleistocene and Holocene deposits of the Venice area

The sequence-stratigraphic investigation by Very High-Resolution (VHR) seismic profiles allowed recognition of the detailed architecture of the late Pleistocene and Holocene succession of the Venice area. In this way deposits previously known by the analyses of scattered cores, mainly taken along the lagoon margin and the littoral strips, have been correlated at regional scale including the near offshore sector and the result has pointed out the lateral variability of the stratal architecture. Late Pleistocene deposits consist of an aggrading floodplain and fluvial channel fills accumulated during decreasing eustatic sea level, and they are coeval with offlapping forced regressive marine wedges in the Central Adriatic basin. The Holocene sequence is composed of three main seismic units separated by major stratal surfaces. Unit 1 (up to 9 m thick) is formed by channelized deposits separated by areas showing sub-horizontal and hummocky reflectors, and is bounded at the base by a surface that records prolonged conditions of subaerial exposure and at the top by a flatter surface resulting from erosion by marine processes. Deposits of Unit 1 are interpreted as estuarine and distributary channel fills, and back-barrier strata. Unit 2 is well distinguishable from Unit 1 only in the offshore area and at the barrier island bounding the Venice Lagoon, and is composed of a prograding marine wedge (up to 10 m thick) that interacts laterally with ebb tidal deltas. Unit 3 consists of a tidal channel complex and inlet deposits, which testify the evolution of the lagoon area. Tidal channels are entrenched in the lagoon mud flat (coeval with Units 1–2) and cut the Pleistocene–Holocene boundary in several places.Following current sequence-stratigraphic concepts, the Holocene sequence is composed of a paralic transgressive systems tract (TST) (Unit 1) overlying a sequence boundary (the Pleistocene–Holocene boundary) and overlain by a marine highstand systems tract (HST) (Unit 2) in seaward locations and by highstand lagoonal deposits landwards. TST and HST are separated by a downlap surface that is amalgamated with a wave ravinement surface in several places. Unit 3 is coeval with the upper part of Unit 2, and its development has been favoured by human interventions, which led to a transgression limited to the lagoon area.Local factors during the deposition, i.e. subsidence, sediment supply, physiography, and current/wave regimes, led to a significant lateral variability in the architecture of the Holocene sequence, as evidenced by the extreme thickness variation of the TST along both depositional strike and dip. The HST, instead, shows less pronounced strike variations in the stratal architecture. Also, present data clearly evidence that the human impact has a great relevance in influencing the late Holocene sedimentation.     

Positive fault inversion triggering 'intrusive diapirism': an analogue modelling perspective

Understanding the formation and the development of salt structures is very important especially because they are of significant economical interest for hydrocarbon trapping and for long-term storage of radioactive waste and energy reserves. Generally, the activity of normal faults developed in extensional regimes is considered the most efficient mechanism for salt diapirs. The results of analogue models reported in this paper suggest a new triggering mechanism for the rise of salt structures during basin inversion. This mechanism relates the localization of ductile diapirs to early normal faults only after their inversion during later shortening. In this case, diapiric growth is related to the strong dip-slip reactivation component along the fault extruding the silicone-simulating salt upward. Some natural cases, in which the timing and the mechanism of diapiric growth is not clear, can be re-interpreted in the light of these analogue model results.     

Star formation in dwarf irregular galaxies: preliminary results

In order to study the history of star formation in dwarf irregular galaxies we have started a project of deep CCD observations of irregulars in the Local Group. For the theoretical interpretation of the colour-magnitude diagrams resulting from these observations we present a new method of numerical simulations based on stellar evolutionary tracks. The preliminary results concerning the galaxies WLM and Sextans B are summarized.     

Semi-analytical solution for viscous Stokes flow in two eccentrically nested spheres

Applying the infinite Prandtl number approximation, a semi-analytical solution for computing 2-D axisymmetric viscous Stokes flow in a model consisting of two eccentrically nested spheres of different viscosities is derived. Since numerical codes based on spectral or finite techniques for modelling mantle flow in a spherical geometry in the presence of lateral viscosity variation are becoming more and more popular, reliable examples for testing and validating such codes are extremely useful. The eccentrically nested sphere solution was used to test a numerical algorithm based on a mixed spherical-harmonic finite-element formulation of the Stokes problem, and good agreement was obtained.     

A new project to monitor land subsidence in the northern Venice coastland (Italy)

The paper deals with the implementation of a levelling and Global Positioning System (GPS) network to control land subsidence in the coastal area north of the Venice Lagoon. About 480 km of levelling lines with 527 benchmarks, 45 of which suited for Differential GPS measurements, were established in 2004. A complete survey of the net was carried out in 2004 soon after its establishment. The 2004 records have been compared with previous scattered data obtained by the use of levelling surveys, DGPS and SAR interferometry. The results show a trend in land settlement that increases from the lagoon margin to the north and jeopardize the Venice coastland. Groundwater withdrawals for domestic, agricultural, and health spas uses, peat oxidation of reclaimed marshlands for farming, natural consolidation of the Holocene deposits, and tectonics of the pre-Quaternary basement are the causes of land subsidence in the study area. Since most of the area lies below the mean sea level and on account of the expected sea level rise due to global change, a detailed monitoring of land displacements in the near future will be of paramount importance to plan necessary works for coastland protection.