Magmatic-dominated fluid evolution in the Jurassic Nambija gold skarn deposits (southeastern Ecuador)

The Jurassic (approximately 145 Ma) Nambija oxidized gold skarns are hosted by the Triassic volcanosedimentary Piuntza unit in the sub-Andean zone of southeastern Ecuador. The skarns consist dominantly of granditic garnet (Ad_20–98) with subordinate pyroxene (Di_46–92Hd_17–42Jo_0–19) and epidote and are spatially associated with porphyritic quartz-diorite to granodiorite intrusions. Endoskarn is developed at the intrusion margins and grades inwards into a potassic alteration zone. Exoskarn has an outer K- and Na-enriched zone in the volcanosedimentary unit. Gold mineralization is associated with the weakly developed retrograde alteration of the exoskarn and occurs mainly in sulfide-poor vugs and milky quartz veins and veinlets in association with hematite. Fluid inclusion data for the main part of the prograde stage indicate the coexistence of high-temperature (500°C to >600°C), high-salinity (up to 65 wt.% eq. NaCl), and moderate- to low-salinity aqueous-carbonic fluids interpreted to have been trapped at pressures around 100–120 MPa, corresponding to about 4-km depth. Lower-temperature (510–300°C) and moderate- to low-salinity (23–2 wt.% eq. NaCl) aqueous fluids are recorded in garnet and epidote of the end of the prograde stage. The microthermometric data (Th from 513°C to 318°C and salinity from 1.0 to 23 wt.% eq. NaCl) and δ¹⁸O values between 6.2‰ and 11.5‰ for gold-bearing milky quartz from the retrograde stage suggest that the ore-forming fluid was dominantly magmatic. Pressures during the early retrograde stage were in the range of 50–100 MPa, in line with the evidence for CO₂ effervescence and probable local boiling. The dominance of magmatic low-saline to moderately saline oxidizing fluids during the retrograde stage is consistent with the depth of the skarn system, which could have delayed the ingression of external fluids until relatively low temperatures were reached. The resulting low water-to-rock ratios explain the weak retrograde alteration and the compositional variability of chlorite, essentially controlled by host rock compositions. Gold was precipitated at this stage as a result of cooling and pH increase related to CO₂ effervescence, which both result in destabilization of gold-bearing chloride complexes. Significant ingression of external fluids took place after gold deposition only, as recorded by δ¹⁸O values of 0.4‰ to 6.2‰ for fluids depositing quartz (below 350°C) in sulfide-rich barren veins. Low-temperature (<300°C) meteoric fluids (δ¹⁸O_water between −10.0‰ and −2.0‰) are responsible for the precipitation of late comb quartz and calcite in cavities and veins and indicate mixing with cooler fluids of higher salinities (about 100°C and 25 wt.% eq. NaCl). The latter are similar to low-temperature fluids (202–74.5°C) with δ¹⁸O values of −0.5‰ to 3.1‰ and salinities in the range of 21.1 to 17.3 wt.% eq. CaCl₂, trapped in calcite of late veins and interpreted as basinal brines. Nambija represents a deep equivalent of the oxidized gold skarn class, the presence of CO₂ in the fluids being partly a consequence of the relatively deep setting at about 4-km depth. As in other Au-bearing skarn deposits, not only the prograde stage but also the gold-precipitating retrograde stage is dominated by fluids of magmatic origin.     

Camouflaged Galactic cosmic microwave background polarization foregrounds: total and polarized contributions of the kinetic Sunyaev–Zeldovich effect

We consider the role of the Galactic kinetic Sunyaev–Zeldovich (SZ) effect as a cosmic microwave background (CMB) polarization foreground. While the Galactic thermal SZ effect has previously been studied and discarded as a potential CMB foreground, we find that the kinetic SZ effect is dominant in the Galactic case. We analyse the detectability of the kinetic SZ effect by means of an optimally matched filter technique applied to a simulation of an ideal observation. We obtain no detection, getting a signal-to-noise ratio of 0.1, thereby demonstrating that the kinetic SZ effect can also safely be ignored as a CMB foreground. However, we provide maps of the expected signal for inclusion in future high-precision data processing. Furthermore, we rule out the significant contamination of the polarized CMB signal by second scattering of Galactic kinetic SZ photons, since we show that the scattering of the CMB quadrupole photons by Galactic electrons is a stronger effect than the SZ second scattering, and has already been shown to produce no significant polarized contamination.     

Late Pleistocene dust deposition in the Patagonian steppe - extending and refining the paleoenvironmental and tephrochronological record from Laguna Potrok Aike back to 55 ka

Paleoenvironmental records extending well into the last glacial period are scarce in the steppe regions of southern South America. Here, we present a continuous record for the past 55 ka from the maar lake Laguna Potrok Aike (51°58′ S, 70°23′ W, southern Patagonia, Argentina). Previous studies on a sedimentary core from a lake level terrace near the northern margin of the lake covered parts of Oxygen Isotope Stage (OIS) 3 (59–29 ka) whereas a second core from the centre of the basin comprised the last 16 ka. Tephrostratigraphical constraints and OSL ages from a third core located below the lake level terrace provide the crucial piece to close the gap between the previous coring sites. High-resolution XRF and magnetic susceptibility as well as grain size data indicate a positive hydrological balance alongside with relatively high aeolian activity during the glacial which is contemporaneous with increased dust fluxes in Antarctica. This is therefore the first evidence for contemporaneity of aeolian deposition in both the target area (Antarctica) and in the major source area of Patagonia. During the Holocene climatic conditions driving sediment deposition seem to have been more variable and less dominated by wind compared to glacial times. The identification of a minor lake level lowering at approximately 4 cal ka BP allows to refine earlier paleoenvironmental reconstructions for the Holocene. Within error margins the OSL ages are consistent with published radiocarbon-dated records offering hence a valuable tool for further studies of the sediments from Laguna Potrok Aike. The new chronology confirms the age of three tephra layers up to now only found in Laguna Potrok Aike sediments and ascribed to OIS 3.     

Imaging a shallow salt diapir using ambient seismic vibrations beneath the densely built-up city area of Hamburg,Northern Germany

Salt diapirs are common features of sedimentary basins. If close to the surface, they can bear a significant hazard due to possible dissolution sinkholes, karst formation and collapse dolines or their influence on ground water chemistry. We investigate the potential of ambient vibration techniques to map the 3-D roof morphology of shallow salt diapirs. Horizontal-to-vertical (H/V) spectral peaks are derived at more than 900 positions above a shallow diapir beneath the city area of Hamburg, Germany, and are used to infer the depth of the first strong impedance contrast. In addition, 15 small-scale array measurements are conducted at different positions in order to compute frequency-dependent phase velocities of Rayleigh waves between 0.5 and 25 Hz. The dispersion curves are inverted together with the H/V peak frequency to obtain shear-wave velocity profiles. Additionally, we compare the morphology derived from H/V and array measurements to borehole lithology and a gravity-based 3-D model of the salt diapir. Both methods give consistent results in agreement with major features indicated by the independent data. An important result is that H/V and array measurements are better suited to identify weathered gypsum caprocks or gypsum floaters, while gravity-derived models better sample the interface between sediments and homogeneous salt. We further investigate qualitatively the influence of the 3-D subsurface topography of the salt diapir on the validity of local 1-D inversion results from ambient vibration dispersion curve inversion.     

Regional gravity modeling in terms of spherical base functions

This article provides a survey on modern methods of regional gravity field modeling on the sphere. Starting with the classical theory of spherical harmonics, we outline the transition towards space-localizing methods such as spherical splines and wavelets. Special emphasis is given to the relations among these methods, which all involve radial base functions. Moreover, we provide extensive applications of these methods and numerical results from real space-borne data of recent satellite gravity missions, namely the Challenging Minisatellite Payload (CHAMP) and the Gravity Recovery and Climate Experiment (GRACE). We also derive high-resolution gravity field models by effectively combining space-borne and surface measurements using a new weighted level-combination concept. In addition, we outline and apply a strategy for constructing spatio-temporal fields from regional data sets spanning different observation periods.     

On the automatic parameter calibration of a hypoplastic soil model

Acta Geotechnica - This paper presents an approach for the automatic parameter calibration (AC) of a hypoplastic constitutive soil model. The calibration software developed in this work simplifies...     

Determination of flow exits in the Soultz borehole GPK2 by using the brine displacement method

The article presents a method to derive synthetic flow logs from low rate injection tests in wells, which are otherwise not accessible for logging. The tests were carried out in the 5 km deep geothermal well GPK2 in the Enhanced Geothermal System (EGS) project of Soultz-sous-Forêts (France) and were aimed at characterizing the status of the well before chemical treatments. The method consists of injecting water at constant flow rates into the well initially filled with brine. The wellhead pressure record of this test is transformed into log of the flow velocity versus depth by a technique described in the paper. Applying this method to borehole GPK2, three outlets could be detected and quantified within the cased and uncased section of this borehole. The outlet in the cased section at 3860 m absorbed 15% of the total flow. This depth corresponds to a known casing restriction, which makes the borehole inaccessible for logging operations. The two other outlets are at the casing shoe at 4420 m (15%) and at 4670 m (70%).     

A state-dependent hyperelastic-plastic constitutive model considering shear-induced particle breakage in granular soils

Acta Geotechnica - An elastic–plastic constitutive model considering particle breakage for simulation of crushable granular soils behavior is proposed. In the model, elastic strain rates are...     

Improved estimates of mean sea level changes in the German Bight over the last 166 years

In this paper, mean sea level changes in the German Bight, the south-eastern part of the North Sea, are analysed. Records from 13 tide gauges covering the entire German North Sea coastline and the period from 1843 to 2008 have been used to derive high quality relative mean sea level time series. Changes in mean sea level are assessed using non-linear smoothing techniques and linear trend estimations for different time spans. Time series from individual tide gauges are analysed and then ‘virtual station’ time series are constructed (by combining the individual records) which are representative of the German Bight and the southern and eastern regions of the Bight. An accelerated sea level rise is detected for a period at the end of the nineteenth century and for another one covering the last decades. The results show that there are regional differences in sea level changes along the coastline. Higher rates of relative sea level rise are detected for the eastern part of the German Bight in comparison to the southern part. This is most likely due to different rates of vertical land movement. In addition, different temporal behaviour of sea level change is found in the German Bight compared to wider regional and global changes, highlighting the urgent need to derive reliable regional sea level projections for coastal planning strategies.     

A probabilistic approach for the classification of earthquakes as ��triggered�� or ��not triggered��

The occurrence time of earthquakes can be anticipated or delayed by external phenomena that induce strain energy changes on the faults. ??Anticipated?? earthquakes are generally called ??triggered??; however, it can be controversial to label a specific earthquake as such, mostly because of the stochastic nature of earthquake occurrence and of the large uncertainties usually associated to stress modelling. Here we introduce a combined statistical and physical approach to quantify the probability that a given earthquake was triggered by a given stress-inducing phenomenon. As an example, we consider an earthquake that was likely triggered by a natural event: the M?=?6.2 13 Jan 1976 Kópasker earthquake on the Grímsey lineament (Tj?rnes Fracture Zone, Iceland), which occurred about 3?weeks after a large dike injection in the nearby Krafla fissure swarm. By using Coulomb stress calculations and the rate-and-state earthquake nucleation theory, we calculate the likelihood of the earthquake in a scenario that contains only the tectonic background and excludes the dike and in a scenario that includes the dike but excludes the background. Applying the Bayes?? theorem, we obtain that the probability that the earthquake was indeed triggered by the dike, rather than purely due to the accumulation of tectonic strain, is about 60 to 90?%. This methodology allows us to assign quantitative probabilities to different scenarios and can help in classifying earthquakes as triggered or not triggered by natural or human-induced changes of stress in the crust.