Hydrological effects on gravity and correlations between gravitational variations and level of the Alzette River at the station of Walferdange,Luxembourg

The gravitational effects of water storage variations driven by local precipitation events are modeled for the Walferdange Underground Laboratory for Geodynamics (WULG) in the Grand Duchy of Luxembourg. A modified mass continuity model is implemented, which uses rainfall data from Walferdange as input. In the absence of soil moisture and groundwater level information, the model is empirically parameterized. Model outputs are compared to the gravity time series registered with the Observatory Superconducting Gravimeter CT040 located in the WULG. We find that the model explains 77% of the gravity residuals. In addition, a statistical analysis is carried out to determine the relationship between precipitation, gravity variations and water level changes in the nearby Alzette River. A time delay of 88 ± 34 min between the maximum variation rates of the water level and gravity signal has been calculated. The signals have an admittance of 45 ± 5 cm μGal^?1.     

Applying local Green's functions to study the influence of the crustal structure on hydrological loading displacements

The influence of the elastic Earth properties on seasonal or shorter periodic surface deformations due to atmospheric surface pressure and terrestrial water storage variations is usually modeled by applying a local half-space model or an one dimensional spherical Earth model like PREM from which a unique set of elastic load Love numbers, or alternatively, elastic Green's functions are derived. The first model is valid only if load and observer almost coincide, the second model considers only the response of an average Earth structure. However, for surface loads with horizontal scales less than 2500 km², as for instance, for strong localized hydrological signals associated with heavy precipitation events and river floods, the Earth elastic response becomes very sensitive to inhomogeneities in the Earth crustal structure.We derive a set of local Green's functions defined globally on a 1° × 1° grid for the 3-layer crustal structure TEA12. Local Green's functions show standard deviations of ±12% in the vertical and ±21% in the horizontal directions for distances in the range from 0.1° to 0.5°. By means of Green's function scatter plots, we analyze the dependence of the load response to various crustal rocks and layer thicknesses. The application of local Green's functions instead of a mean global Green's function introduces a variability of 0.5–1.0 mm into the hydrological loading displacements, both in vertical and in horizontal directions. Maximum changes due to the local crustal structures are from −25% to +26% in the vertical and −91% to +55% in the horizontal displacements. In addition, the horizontal displacement can change its direction significantly. The lateral deviations in surface deformation due to local crustal elastic properties are found to be much larger than the differences between various commonly used one-dimensional Earth models.     

Monitoring of Guadalentín valley (southern Spain) through a fast SAR Interferometry method

Ground subsidence is a critical problem in populated areas, which requires attention and routine monitoring. Specifically, subsidence related to the exploitation of aquifer systems is of particular interest in terms of risk for the possible damage to buildings and availability of water resources. We propose to realise a survey of the aquifer systems near the city of Lorca, southern Spain, using Synthetic Aperture Radar Interferometry (InSAR) with a light method necessitating a small number of interferograms rather than analysing hundreds of radar images. We computed nine interferograms with good coherence over the years 2004 and 2005. We analysed them by a two-step filtering process in order to characterise the atmospheric contribution and the ground deformation. We are then able to identify two deforming areas and quantify subsidence rates at 0.20 ± 0.02 mm/day south of Lorca and southeast to the city of Alhama de Murcia with the same efficiency as the studies using sophisticated InSAR methods. Using continuous records of water levels in boreholes, we deduced a mean storage coefficient (S_k) at 4.1 × 10^− 3 which is very close to the previous published values, and coincides in this case to an inelastic deformation mode of all the system. With this value, we determine the groundwater level variations from the measurements of the ground subsidence.     

Short-term changes in the aerosol characteristics at Kharagpur (22°19′N, 87°19′E) during winter

In-situ measurements of number density, size distribution, and mass loading of near-surface aerosols were carried out at Kharagpur, a site on the eastern part of Indo-Gangetic Plains during the winter month of December 2004. The data have been used to investigate wintertime characteristics of aerosols and their effects on the occurrence of haze. The aerosol number density is found to be of the order of 10⁹ m^?3 and mass loading is ~265±70 μg m^?3 (5–8 times that reported from south Indian sites). The diurnal patterns and day-to-day variations in aerosol number density and mass loading are closely associated with atmospheric boundary layer height. During haze events, the number density of submicron particles is found to be 2–5 times higher than that during non-hazy period. This could be attributed to the enhanced concentration of anthropogenic aerosols, low atmospheric boundary layer height/ventilation coefficient and airflow convergence.     

A search for atmospheric effects on gravity at different time and space scales

We investigate in this paper various approaches to correct gravity changes for the effect of atmospheric pressure changes. Two specific locations are considered: Strasbourg (France) as mid-latitude station, where regular pressure fronts occur and Djougou (Benin) as equatorial station with large thermally driven S₁ and S₂ waves of planetary extent. We first review the classical approaches based on a constant or frequency-dependent admittance using only local pressure and gravity data. We consider then a model of atmospheric loading and show the barometric admittance in terms of elastic, Newtonian and total load, as a function of the distance from the station. We consider both a 2D pressure model (surface loading) and a 2.5D model, where the density decreases with height (standard atmosphere). Assuming horizontal advection in the atmospheric dynamics, we convert this spatially dependent admittance into a frequency-dependent admittance. Using global pressure data from European Centre for Medium-Range Weather Forecasts (ECMWF) at about 12 km spatial resolution and 3 h sampling, we compute the model-predicted pressure admittance for Djougou and Strasbourg and we simulate the frequency dependence inferred from gravity and pressure observations below 4 cycle per day. A long gravity and pressure data set (1996–2013) from Strasbourg is used to investigate the low frequency part of the pressure admittance while a common 2.5 year data set (August 2010–February 2013) for Strasbourg and Djougou is then analyzed to investigate the high frequency part of the admittance. In both cases, our results are in close agreement with the predictions inferred from an atmospheric 2.5D loading model with a distance–time relationship due to horizontal advection. The frequency dependence of the barometric admittance is explained by the competing contributions of Newtonian attraction and elastic surface deformation according to the distance from the gravimeter. In the far field (low frequencies), the magnitude of the admittance decreases with frequency because of the combined elasticity effect and Newtonian attraction (when the atmosphere is below the horizon) while, on the contrary, in the near field (high frequencies), elasticity becomes negligible and the pressure admittance mainly decreases with increasing frequency because of the decreasing attraction effect of the atmospheric masses inside the cylindrical pressure cell centered on the sensor location of decreasing radius. In the last part, we show that there is variability in time in the pressure admittance for both stations.     

Atmospheric polychlorinated biphenyl (pcb) inputs to a coastal city near the marmara sea

Recent studies have shown that polychlorinated biphenyl (PCB) levels are high in coastal sites. Various types of deposition events, including dry deposition, air–water exchange and wet deposition, were analyzed as part of the study to illustrate the pollution level. The atmospheric levels and deposition mechanisms of 82 PCBs were investigated in Mudanya, a coastal city in Turkey. A total of 175 samples and 112 blanks were collected. Air–water exchange and particle phase dry deposition fluxes were 120 ± 90 ng/m²-d and 125 ± 160 ng/m²-d, respectively. The wet-dry deposition sampler (WDDS) consisted of two reservoirs, wet and dry. Collected with the WDDS, the dry deposition flux was 21 ± 20 ng/m²-d, and the wet deposition fluxes (during the rainy period) were 145 ± 130 ng/m²-d and 195 ± 270 ng/m²-d for the dissolved and particulate phases, respectively. Mass transfer coefficients (MTCs), dry deposition velocities and washout ratios were calculated and compared with the literature values.     

Sea level rise in the north-western part of the Arabian Gulf

Relative sea level variations in the north-western part of the Arabian Gulf have been estimated in the past using no more than 10 to 15 years of observations. In this study, we have almost doubled the period to 28.7 years by examining all available tide gauge data in the area and constructing a mean gauge time-series from seven coastal tide gauges. We found for the period 1979–2007 a relative sea level rise of 2.2 ± 0.5 mm/year. Using the subsidence observed at 6 GPS stations within a radius of 100 km of the tide gauges as an indication of the vertical land motion, the corresponding absolute sea level rise is 1.5 ± 0.8 mm/year that is in agreement with the global estimate of 1.9 ± 0.1 mm/year (Church and White, 2011) for the same studied period. By taking into account the temporal correlations we conclude that previous published results underestimate the true sea level rate uncertainty in this area by a factor of 5–10.     

GPS constraints on the deformation of Azerbaijan and surrounding regions

In this study, we present new GPS observations in Azerbaijan to provide an improved basis for determining the distribution of crustal deformation throughout the country and surrounding areas. The deformation field in the region has been analyzed with a dense GPS network configuration and a reliable quantification of the ongoing deformation was achieved. Results show that while contraction is dominant over the whole region, it is mostly concentrated on the middle and eastern parts of Caucasus Thrust Fault reaching up to 6.4 ± 0.2 mm/yr and Lesser Caucasus Fault does no accommodate more than 1–2 mm/yr of contraction. New network also clearly substantiates that the West Caspian Fault, which is a continuation of Caucasus Thrust Fault in the south, accommodates right-lateral slip rates of 7.1 ± 0.3 mm/yr in addition to 5.5 ± 0.3 mm/yr contraction rates.     

Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units,Lyon Turin transect,Western Alps)

In the central Western Alps, a combined structural, petrological and ⁴⁰Ar–³⁹Ar geochronological study of the Modane-Aussois and Southern Vanoise units yields important constraints on the timing of deformation and exhumation of the Briançonnais zone. These data help to decipher the respective roles of oceanic subduction, continental subduction and collision in the burial and exhumation of the main units through time. In the Modane-Aussois unit top to the NW thrusting (D1) was followed by top to the east shearing (D2) interpreted by some as normal faulting and by others as backthrusting. Pseudosection calculations imply that D1 deformation occurred at 1.0 ± 0.1 GPa and 350 ± 30 °C. Analysis of chlorite–phengite pairs yield P–T estimates between 0.15 and 0.65 GPa and between 220 and 350 °C for the D2 event. Phengites along the D1 schistosity (sample M80) yields an ⁴⁰Ar–³⁹Ar age of 37.12 ± 0.39 Ma, while D2 phengites yield ages of 35.42 ± 0.38 (sample M173) and 31.60 ± 0.33 Ma (sample M196). It was not possible to test whether these ages are altered by excess argon or not. Our interpretation is that the D1/D2 transition occurred at ∼37 Ma at the beginning of decompression, and that D2 lasted until at least ∼32 Ma. Pseudosection calculation suggests that the Southern Vanoise unit was buried at 1.6 ± 0.2 GPa and 500–540 °C. D1 deformation occurred during exhumation until 0.7–10.5 GPa and 370 ± 30 °C. Published ages suggest that D1 deformation possibly started at ∼50 Ma and lasted until ∼37 Ma. D2 deformations started at P–T conditions close to that recorded in Modane-Aussois unit and lasted until 0.2 ± 0.1 GPa and 280 ± 30 °C at ∼28 Ma. The gap of 0.6 ± 0.3 GPa and 150 ± 130 °C between peak metamorphic conditions in the two units was concealed by thrusting of the South Vanoise unit on top of the Modane-Aussois unit during D1 Deformation. Top to the east deformation (D2) affects both units and is interpreted as backthrusting.Based on these data, we propose a geodynamic reconstruction where the oceanic subduction of the Piedmont unit until ∼50 Ma, is followed by its exhumation at the time of continental subduction of the continental Southern Vanoise unit until ∼45 Ma. The Southern Vanoise is in turn underthrusted by the Modane-Aussois unit until ∼37 Ma (D1). Between 37 and 31 Ma the Modane-Aussois and Southern Vanoise units exhume together during backthrusting to the east (D2). This corresponds to the collision stage and to the activation of the Penninic Thrust. In the ∼50 Ma to ∼31 Ma time period the main thrusts propagated westward as the tectonic context switched from oceanic to continental subduction and finally to collision. During each stage, external units are buried while internal ones are exhumed.     

Seasonal changes in particulate biogenic and lithogenic silica in the upwelling system off Concepción (∼36°S), Chile,and their relationship to fluctuations in marine productivity and continental input

We analyzed the temporal and vertical distribution of biogenic (BSi) and lithogenic (LSi) silica, and diatom abundance in the upwelling center off Concepción, Chile, from April 2004 to May 2005. Measurements were performed at the FONDAP COPAS Time Series Station 18 (36°30.8′S, 73°07.7′W; 88 m water depth), and were combined with primary production estimates and river runoff data to assess the relationships between water column BSi and primary production, and between LSi and river runoff. Throughout the sampling period, water-column-integrated (0–80 m) BSi averaged 252±287 mmol m⁻², and was about six times higher than average LSi (44±30 mmol m⁻²). The highest water column BSi observed during the upwelling season (786±281 mmol m⁻²) coincided with increments in total diatom abundance, and high integrated chlorophyll a concentration and primary production. In contrast, LSi was nearly two times higher in winter (85±43 mmol m⁻²) than the annual average, in agreement with the period of substantial discharges from the Itata and Bio-Bio rivers. The observed temporal patterns in BSi and LSi are coincident with primary production-related factors and riverine outflow, respectively, suggesting that the BSi and LSi pools are separate. With respect to the vertical distribution in the water column, most of the BSi and diatoms were found in surface waters (0–30 m depth), whereas LSi was most abundant at depth. Our study attempts to make an inventory of both BSi and LSi in the water column off Concepción, and gives the present-day background information necessary to assess potential future changes in the hydrological cycle that, in turn, may induce modifications in the Si path from the watersheds to the ocean.     

Seismo-meteo-electromagnetic phenomena observed during a 5-year interval around the 2011 Tohoku earthquake

The purpose of the paper is to try to find ULF electromagnetic precursors to the 2011 Tohoku earthquake (EQ), on the basis of extensive investigation of radiations in vertical component of the magnetic field or with a large ratio of the vertical to full horizontal component. Nighttime records have been analyzed of three Japanese fluxgate magnetometers located in a distance from 300 to 1300 km from the epicenter of the main shock, and the frequency range from 10 to 150 mHz was used for the analysis. We have applied wavelet analyses to improve the detection of pulsed signals. All obtained scalograms have been averaged over the nighttime interval from 01 h to 05 h JST and flattened by means of multiplication by square of frequency. The sequence of spectra thus obtained has been compared with the evolution of seismicity, which has resulted in that the radiation in the vertical component has been detected. It exhibits seasonal variations with winter maxima, but it increases further by approaching the moment of the EQ and decreases after that. This radiation seems to be correlated with atmospheric parameters – air humidity, temperature, vapor pressure and rainfall. So, we consider that this radiation cannot be caused by subsurface sources, but its possible sources can be atmospheric discharges. The evolution of this phenomenon can be explained by a seasonal variation of atmospheric parameters and also its variations under the influence of injection of gas from the focal zone of a forthcoming EQ.     

Crustal loading in vertical GPS time series in Fennoscandia

We compare time series of vertical position from GPS with modelled vertical deformation caused by variation in continental water storage, variation in the level of the Baltic Sea, and variation in atmospheric pressure. Monthly time series are used. The effect of continental water storage was calculated from three different global models. The effect of non-tidal variation in Baltic Sea level was calculated using tide gauge observations along the coasts. Atmospheric loading was computed from a numerical weather model. The loading time series are then compared with three different GPS time series at seven stations in Fennoscandia. A more detailed analysis is computed at three coastal stations. When the monthly GPS time series are corrected using the load models, their root-mean-square scatter shows an improvement between 40 and 0%, depending on the site and on the GPS solution. The modelled load effect shows a markedly seasonal pattern of 15 mm peak-to-peak, of which the uncorrected GPS time series reproduce between 60 and 0%.     

Cenozoic changes in atmospheric lead recorded in central Pacific ferromanganese crusts

The possible sources of pre-anthropogenic Pb contributed to the world's oceans have been the focus of considerable study. The role of eolian dust versus riverine inputs has been of particular interest. With better calibration of isotopic records from central Pacific ferromanganese crusts using Os isotope stratigraphy it now appears that deep water Pb isotopic compositions were effectively homogeneous over a distance of 5000 km for the past 80 Myr. The composition shifted slightly from high ²⁰⁶Pb/²⁰⁴Pb ratios in the range of 18.87 ± 0.02 before 65 Ma to lower values of 18.62 ± 0.02 by 45 Ma and then gradually increased again very slightly to the present day ratio of 18.67 ± 0.02. The regional homogeneity provides evidence of a dominant well-mixed atmospheric source the most likely candidate for which is volcanic aerosols contributed either directly or as soluble condensates on eolian dust. The slight shift in Pb isotope composition of deep waters in the central Pacific between 65 and 45 Ma may be the result of a regional- or perhaps global-scale change in the sources of volcanic exhalations and volcanic activity caused by an increase in the importance of melting and assimilation of older continental crustal components over the Cenozoic.     

High resolution mapping of Earth tide response based on GPS data in Japan

We observe the Earth tidal fields at diurnal and semi-diurnal periods using Kinematic Precise Point Positioning (KPPP) GPS analysis. Our KPPP GPS solutions compare well with super-conducting gravimeter (SG) observations and a theoretical Earth tidal model, that includes both ocean tide loading model and body tides. We make a high resolution map of the observed Earth tidal response fields using the Japanese GEONET GPS network which consists of 1200 sites. We find that: (1) the average phase of GPS data lags 0.11 ± 0.04° from our theoretical Earth tidal model, (2) the average amplitude ratio between GPS and the theoretical Earth tidal model is 1.007 ± 0.003, (3) the amplitude in the Kyushu district is about 1.0–1.5 ± 0.3% larger than in the Hokkaido district, and (4) the amplitude at the Japan Sea side is about 0.5 ± 0.2% larger than that at the Pacific Ocean side. These results suggest that we may be able to place constraints on Earth structure using GPS-derived tidal information.     

Selected trace metals (As,Cd and Hg) distribution and contamination in the coastal wetland sediment of the northern Beibu Gulf,South China Sea

Contamination with As, Cd and Hg, their spatial and temporal distribution are reported from the coastal wetland sediments of the northern Beibu Gulf, South China Sea. The content of As, Cd, Hg and TOC in surface sediments is 8.1 ± 5.8 μg g^?1, 0.08 ± 0.14 μg g^?1, 0.034 ± 0.028 μg g^?1 and 0.45 ± 0.39%, respectively. The mean sedimentation rates are 0.93–1.37 cm year^?1 during 1920s to 2008 determined by ²¹⁰Pb and ¹³⁷Cs dating in three cores. The vertical profiles of As, Cd and Hg content in the cores retrieved from Qin and Nanliu River estuaries show increasing trends during 1985–2008 due to anthropogenic impact caused by local economic development. Locally the surface sediments have potential ecological risk of As to benthos according to the NOAA sediment quality guidelines.     

Empirical relationship between leaf wax n-alkane δD and altitude in the Wuyi,Shennongjia and Tianshan Mountains,China: Implications for paleoaltimetry

Estimating past elevation not only provides evidence for vertical movements of the Earth's lithosphere, but also increases our understanding of interactions between tectonics, relief and climate in geological history. Development of biomarker hydrogen isotope-based paleoaltimetry techniques that can be applied to a wide range of sample types is therefore of continuing importance. Here we present leaf wax-derived n-alkane δD (δD_wax) values along three soil altitudinal transects, at different latitudes, in the Wuyi, Shennongjia and Tianshan Mountains in China, to investigate δD_wax gradients and the apparent fractionation between leaf wax and precipitation (ε_wax-p).We find that soil δD_wax track altitudinal variations of precipitation δD along the three transects that span variable environment conditions and vertical vegetation spectra. An empirical δD_wax-altitude relation is therefore established in which the average δD_wax lapse rate of ? 2.27 ± 0.38‰/100 m is suitable for predicting relative paleoelevation change (relative uplift). The application of this empirical gradient is restricted to phases in the mountain uplift stage when the atmospheric circulation had not distinctly changed and to when the climate was not arid. An empirical δD_wax–latitude–altitude formula is also calculated: δD_wax = 3.483LAT ? 0.0227ALT ? 261.5, which gives the preliminary spatial distribution pattern of δD_wax in modern China.Mean value of ε_wax-p in the extreme humid Wuyi Mountains is quite negative (? 154‰), compared to the humid Shennongjia (? 129‰) and the arid (but with abundant summer precipitation) Tianshan Mountains (? 130‰), which suggests aridity or water availability in the growing season is the primary factor controlling soil/sediment ε_wax-p. Along the Tianshan transects, values of ε_wax-p are speculated to be constant with altitude; while along the Wuyi and Shennongjia transects, ε_wax-p are also constant at the low-mid altitudes, but become slightly more negative at high altitudes which could be attributed to overestimates of precipitation δD or the vegetation shift to grass/conifer.Additionally, a reversal of altitude effect in the vertical variation of δD_wax was found in the alpine zone of the Tianshan Mountains, which might be caused by atmospheric circulation change with altitude. This implies that the paleo-circulation pattern and its changes should also be evaluated when stable isotope-based paleoaltimetry is applied.     

Radiogenic heat production in the lithosphere of Sulu ultrahigh-pressure metamorphic belt

The Chinese Continental Scientific Drilling (CCSD) project is located at the Sulu ultrahigh-pressure metamorphic (UHPM) belt. It offers a unique opportunity for studying the radiogenic heat production of both shallower and deeper rocks. Based on the concentrations of radiogenic elements U, Th and K on 349 samples from main hole of CCSD (CCSD MH), pilot holes and exposures, we determined radiogenic heat productions of all major rock types in the Sulu UHPM belt. Results show the mean values of orthogneiss and paragneiss are respectively 1.65 ± 0.81 and 1.24 ± 0.61 µW m^? 3. Due to different composition and grade of retrogressive metamorphism, the eclogites display significant scatter in radiogenic heat production, ranging from 0.01 to 2.85 µW m^? 3, with a mean of 0.44 ± 0.55 µW m^? 3. The radiogenic heat production in ultramafic rocks also varies within a large range of 0.02 to 1.76 µW m^? 3, and the average turns out to be 0.18 ± 0.31 µW m^? 3. Based on the measurements and crustal petrologic model, the vertical distribution model of heat production in Sulu crust is established. The resulting mean heat production (0.76 µW m^? 3) contributes 24 mW m^? 2 to the surface heat flow. 1-D thermal model indicates that the temperature at the Moho reaches above 750 °C, and the thermal thickness of the lithosphere is ~ 75 km, in good agreement with the geophysical results. The high teat flow (~ 75 mW m^? 2) together with thin lithosphere presents strong support for the extension events during the late Cretaceous and Cenozoic.     

Melting and thermal expansion in the Fe–FeO system at high pressure

Melting in the Fe–FeO system was investigated at pressures up to 93 GPa using synchrotron X-ray diffraction (XRD) and a laser heated diamond anvil cell (DAC). The criteria for melting were the disappearance of reflections associated with one of the end-member phases upon raising the temperature above the eutectic and the reappearance of those reflections on dropping the temperature below the eutectic. The Fe–FeO system is a simple eutectic at 50 GPa and remains eutectic to at least 93 GPa. The eutectic temperature was bound at several pressure points between 19 and 93 GPa, and in some cases the liquidus temperature was also determined. The eutectic temperature rises rapidly with pressure closely following the melting curve of pure Fe. A detailed phase diagram at 50 GPa is presented; the eutectic temperature is 2500 ± 150 K and the eutectic composition is bound between 7.6 ± 1.0 and 9.5 ± 1.0 wt.% O. The coefficient of thermal expansion of FeO is a strong function of volume and decreases with pressure according to a simple power law.     

Miocene magmatism and tectonics within the Peri-Alboran orogen (western Mediterranean)

The aim of this paper concerns Miocene igneous activity in the Alboran Sea and Peri-Alboran area (northern Morocco, western Algeria and Betic Cordilleras in Spain), considering its age and its location with regard to major tectonics structures.We have compiled previous K-Ar isotopic ages of lavas and plutonic boulders and intrusives with an error of ±1σ and completed this set by a new K-Ar isotopic age for andesitic tuffites from Alboran Island. Geochemistry of most of these samples has been considered after previous analyses completed with new data for Spain magmatism. These two sets of data allow us to place the magmatic activity within the regional stratigraphy and tectonics and their chronological framework of the three major tectonic phases of the Maghrebian orogen, at 17 Ma (Burdigalian), 15 Ma (Langhian) and 9 Ma (Tortonian). Petro-geochemical characteristics are compared through time and geographical locations. A major goal of this coupled approach is to help the elaboration of possible geodynamical processes.As an application, we present the case study of the Dellys, Djinet and Thenia region (east of Algiers) where the successive magmatic events between 19.4 ± 1 and 11.6 ± 0.5 Ma are closely related to the local tectonics and sedimentation.The Peri-Alboran igneous activity is placed in a multidisciplinary framework. Timing of activity is defined according to the ages of the neighbouring sedimentary units and the K-Ar ages of igneous rocks.In Spain, the Cabo de Gata-Carboneras magmatic province displays late Oligocene and early Miocene leucogranitic dikes, dated from 24.8 ± 1.3 to 18.1 ± 1.2 Ma; three following andesitic to rhyolitic events took place around 15.1 ± 0.8 to 14.0 ± 0.7 Ma, 11.8 ± 0.6 to 9.4 ± 0.4 Ma, 8.8 ± 0.4 to 7.9 ± 0.4 Ma; this last event displays also granitic rocks. Lamproitic magmas dated between 8.4 ± 0.4 and 6.76 ± 0.04 Ma were emplaced after the Tortonian phase.In Morocco, after the complex building of the Ras Tarf volcanic edifice, major calc-alkaline to shoshonitic volcanoes were built between 9.0 ± 0.5 and 4.8 ± 0.5 Ma, in particular the large Gourougou volcanic complex. Near Oujda, volcanic activity of alkaline affinity leads to multiple emissions of basalts throughout Pliocene times until the beginning of Pleistocene, between 6.2 ± 0.3 and 1.5 ± 0.1 Ma.In the Alboran domain, an age of 19.7 ± 0.8 Ma is reported (this study) for the andesitic tuffites that form the emergent part of the Alboran Island. This age is comparable to that of the Algerian tuffites and cherts “silexites” and the Malaga ones in Spain. Younger activity, completely separated from the previous one, forms the low-K basaltic andesitic dikes from Alboran Island, dated between 9.1 ± 0.5 and 7.5 ± 0.3 Ma. Along the Alboran Ridge both low-K and high-K andesites to dacites were emitted in the estimated range of 10.7–8.7 Ma. Low-K and high-K andesites to dacites sampled at ODP sites 977 and 978 into the East Alboran Basin, are dated between 12.1 ± 0.2 and 9.3 ± 0.1 Ma.We propose to relate with the Trans-Alboran lineament only the post-Tortonian igneous activity younger than 9 Ma.     

The oxidation state of sulfur in magmatic fluids

Sulfur compounds in volcanic gases are responsible for the global cooling after explosive eruptions and they probably controlled the early evolution of the Earth's atmosphere. We have therefore studied the oxidation state of sulfur in aqueous fluids under the pressure and temperature conditions and oxygen fugacities typical for magma chambers (0.5–3 kbar, 650–950 °C, Ni–NiO to Re–ReO₂ buffer conditions). Sulfur speciation was determined by Raman spectroscopy of quenched fluids trapped as inclusions in quartz. Our results show that sulfur in hydrothermal fluids and volcanic gases is much more oxidized than previously thought and in particular, some explosive eruptions may release a significant fraction of sulfur as SO₃ or its hydrated forms. In the pressure range from 500 to 2000 bar, the equilibrium constant K₁ of the reaction 2H₂S + 3O₂ = 2SO₂ + 2H₂O in aqueous fluids can be described by lnK₁ = ?(57.1 ± 7.1) + (173,480 ± 7592)T^? 1, where T is temperature in Kelvin. The equilibrium constant K₂ for the reaction SO₂ + ½O₂ = SO₃ in aqueous fluids, where SO₃ may include hydrated forms, such as H₂SO₄, was found to be strongly pressure dependent, with lnK₂ = ?(5.2 ± 5.7) + (19,243 ± 5993)T^? 1 at 1500 bar; lnK₂ = ?(11.1 ± 1.3) + (25,383 ± 1371)T^? 1 at 2000 bar and lnK₂ = ?(22.1 ± 2.2) + (37,082 ± 2248)T^? 1 at 2500 bar. Our data imply that volcanoes may directly inject hexavalent sulfur in the form of H₂SO₄ into the atmosphere, not only on Earth, but possibly also on Venus and on Mars, when it was still tectonically active. Remote measurements from satellites may have underestimated the sulfur yield of some recent eruptions. Moreover, the mechanisms of the interaction of volcanic gases with the stratosphere need to be reconsidered.