Bismuth and cobalt minerals as indicators of stringer zones to massive sulphide deposits, Iberian Pyrite Belt

The stringer zones and commonly the interaction zone at the base of the massive sulphide mounds in the Iberian Pyrite Belt contain bismuth and cobalt minerals that are not found in the overlying massive sulphides. These are fairly rare cobalt sulphoarsenides (cobaltite, alloclasite, galucodot) that were formed at the beginning of the massive sulphide genesis, and fairly common bismuth sulphides (bismuthinite, hammarite, wittichenite, cosalite, kobellite, joseite, etc.), including species rare at world scale (nuffieldite, giessenite, jaskolskiite) that were deposited from last stage high-temperature (> 300 °C) copper-bearing fluids containing Bi (Te, Se). The last stage fluids precipitated chalcopyrite containing Cu, Bi, Te, (Se) sulphosalts at the base of the sulphide mound to form a high cupriferous zone. Their interaction with the massive sulphides is reflected by the formation of an exchange zone, a few metres thick, showing chalcopyrite disease textures, at the base of the mound; this zone forms the upper limit of potentially economic copper enrichment and of bismuth minerals. Gold is undoubtedly in part, if not totally, related to this last phase. The bismuth concentrations being equivalent in the massive sulphides and the stringers, the presence of bismuth minerals in the stringer zones results from high-temperature conditions combined with a rarity of galena, which impedes absorption of available Bi. The distribution of these bismuth minerals provides a basic mineralogical zoning in the stringer zone, with a deep, low-aS₂ zone containing native bismuth and tellurides and a shallow, higher-aS₂ zone in contact with the massive ore sensu stricto and containing complex bismuth sulphides. These results make it possible to distinguish between sulphide veinlets belonging to stockwork zones of massive orebodies and veinlets of an ambiguous nature, and provide mineralogical criteria for the proximity of copper-rich zones. They enrich the very complex mineralogy of the Iberian Pyrite Belt.     

Test of the vibrational modelling for the λ-type transitions: Application to the α-β quartz transition

Vibrational modelling is at the present time the only known way to predict the heat capacities of the Earth's mantle minerals at high-pressure and high-temperature. To test the validity of this method for λ-type transitions, we have applied it to the α-β quartz transition (T ₀=846±1 K). Raman spectra of quartz were recorded up to 900 K. Measured frequency shifts of the α-quartz Raman modes were then used in conjunction with available high-pressure Raman data to calculate intrinsic mode anharmonicity, through the parameter a _i=(?Lnv_i/?T)_v. Vibrational modelling of the heat capacity at constant volume, Cv, and at constant pressure, Cp, including anharmonic corrections deduced from the a _i parameters, are compared to experimental data. Taking into account the soft-mode associated to the α-β quartz transition, the model reproduces the excess of Cp related to the transition. Then, this study confirms that detecting a soft-mode from vibrational data allows one to predict λ-type transitions. However, when modelling the thermodynamic properties, the contribution of a soft-mode cannot be established from spectroscopic data. Therefore, one needs first to determine this contribution in order to predict the heat capacities of Earth's mantle minerals displaying λ-type transitions. In α-quartz, this contribution has been determined as 0.007% of the total number of the optic modes in the model of the density of states.     

Present-day trends of vertical ground motion along the coast lines

Vertical ground motion (VGM) rates stand as crucial information, either for predicting the impact of the actual sea level rise along low-lying coasts or refining geodynamic problems. Because present day VGM rates have a magnitude smaller than 10 mm/yr, they remain challenging to quantify and often elusive. We focus on the quantification of global-scale VGM rates in order to identify global or regional trends. We computed VGM rates by combining tide gauges records and local satellite altimetry, which yield a new dataset of 634 VGM rates. We further compare this database to previous studies that use geodetic techniques and tide gauges records in order to evaluate the consistency of both our results and previous ones. The magnitudes differ by less than 5 mm/yr, and similar subsidence and uplift general tendencies appear. Even if the asset of our database stands in the greater number of sites, the combination of all studies, each with different pros and cons, yields a hybrid dataset that makes our attempt to extract VGM trends more robust than any other, independent study. Fennoscandia, the West coast of North America, and the eastern coast of Australia are uplifting, while the eastern coast of North America, the British Isles and Western Europe, the eastern Mediterranean Sea, Japan, and the western coast of Australia are subsiding. Glacial Isostatic Adjustment (GIA) is expected to provide a major contribution to the present-day signal. Aside from Fennoscandia, observed VGM often depart from the GIA model predictions of Peltier (2004). This either results from an underestimate of the model predictions or from the influence of other processes: indeed, the influence of the geodynamic setting appears in particular along the coasts of western North America or Japan, where the alternation of transform faults and subduction zones makes it possible to assign contrasted behaviours to the local geodynamic context. Local mechanisms like anthropogenic processes or sediment compaction, also contribute to VGM. This remains true for the critical cases of Venice, the Gulf of Mexico, the Ganges delta, and the Maldives, which are particularly exposed to the current sea level rise.     

On the inclusion of ground-based gravity measurements to the calibration process of a global rainfall-discharge reservoir model: case of the Durzon karst system (Larzac, southern France)

This work examines the relevance of the inclusion of ground-based gravity data in the calibration process of a global rainfall-discharge reservoir model. The analysis is performed for the Durzon karst system (Larzac, France). The first part of the study focuses on the hydrological interpretation of the ground-based gravity measurements. The second part of the study investigates further the information content of the gravity data with respect to water storage dynamics modelling. The gravity-derived information is found unable to either reduce equifinality of the single-objective, discharge-based model calibration process or enhance model performance through assimilation.     

Historical evidence of climate change impact on drought outlook in river basins: analysis of annual maximum drought severities through daily SPI definitions

In the context of major outcomes of a steadily changing climate, extreme climatic conditions and the associated events in various forms of weather-related natural disasters, e.g. droughts, floods, and heat waves, are more frequently experienced on the global scale in recent years. In support of this argument, there are adequate numbers of explicit signals over such a persistent outlook, which is greatly illustrated by historical data and observations. This study, which is mainly oriented to investigating the drought behaviour in Thracian, Aegean and Mediterranean transects of Turkey's major river basins, is actually inspired by the foreseen potential of using annual maximum drought severity series (based on drought definition through the standardized precipitation index (SPI)) within a framework that resembles the use of flood discharge directly from flow measurements in a river basin. To this end, a series of spatial analyses were employed to identify different aspects of flood appearance in the study extent, including trend views on annual average drought severity series, shifts in the starting time of the annually most severe flood periods, and changes in spatial coverage views of average drought conditions under different drought severity categories. The framework of the analytical approaches depends greatly on validated international datasets and open-source computational algorithms. The results from the analyses that were conducted in two consecutive periods of 1958–1980 and 1981–2004 revealed that Turkey's western and southern river basin systems seemed to have experienced quite different behaviours between the two periods in terms of drought severity magnitudes, drought durations and annual occurrence times.

     

Determination of the 87Rb decay constant

The decay constant ⁸⁷Rb has been redetermined by measuring the amount of radiogenic ⁸⁷Sr produced over a period of 19 years, in 20 g samples of purified RbClO₄, using isotope dilution techniques. The rubidium sample was spiked with ⁸⁴Sr and the nanogram quantities of strontium separated by coprecipitation with Ba(NO₃)₂. Analyses were carried out on a 25cm, 90° sector mass spectrometer equipped with a Spiraltron electron multiplier. Measurement of three independent ratios permitted continuous monitoring of the ion beam fractionation. The average of nine determinations gives a value for the decay constant of $\text{1.419(±0.012) × 10}{}^{\mathrm{?11}}\text{yr}{}^{\mathrm{?1}}\text{(2σ)}$. [$\text{τ}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{= 4.89(±0.04) × 10}{}^{10}\text{yr}$.]     

A stochastic methodology for risk assessment of a large earthquake when a long time has elapsed

We propose a stochastic methodology for risk assessment of a large earthquake when a long time has elapsed from the last large seismic event. We state an approximate probability distribution for the occurrence time of the next large earthquake, by knowing that the last large seismic event occurred a long time ago. We prove that, under reasonable conditions, such a distribution is exponential with a rate depending on the asymptotic slope of the cumulative intensity function corresponding to a nonhomogeneous Poisson process. As it is not possible to obtain an empirical cumulative distribution function of the waiting time for the next large earthquake, an estimator of its cumulative distribution function based on existing data is derived. We conduct a simulation study for detecting scenario in which the proposed methodology would perform well. Finally, a real-world data analysis is carried out to illustrate its potential applications, including a homogeneity test for the times between earthquakes.     

Salinisation de l’aquifère libre de la Chaouia côtière (Azemmour-Tnine Chtouka), Maroc

L’aquifère libre de la Chaouia côtière constitue un exemple des aquifères les plus exploités au Maroc. Ce travail est consacré à l’étude des processus de la salinisation des eaux souterraines par l’analyse physico-chimique de 39 puits répartis dans la zone. Deux types de faciès ont été révélés, l’un est chloruré-sodique dans la frange côtière ; l’autre faciès est de type chloruré bicarbonaté-calcique caractérisant les eaux exploitées dans les calcaires marneux du Crétacé. L’influence marine (aérosols et intrusion marine), la dissolution/précipitation de la roche aquifère et l’infiltration des eaux d’irrigation, sont parmi les causes de l’augmentation de la salinité des eaux souterraines, en plus de l’exploitation excessive de l’eau souterraine.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.