CO2 emissions and heat flow through soil,fumaroles, and steam heated mud pools at the Reykjanes geothermal area,SW Iceland

Carbon dioxide emissions and heat flow through soil, steam vents and fractures, and steam heated mud pools were determined in the Reykjanes geothermal area, SW Iceland. Soil diffuse degassing of CO₂ was quantified by soil flux measurements on a 600 m by 375 m rectangular grid using a portable closed chamber soil flux meter and the resulting data were analyzed by both a graphical statistical method and sequential Gaussian simulations. The soil temperature was measured in each node of the grid and used to evaluate the heat flow. The heat flow data were also analyzed by sequential Gaussian simulations. Heat flow from steam vents and fractures was determined by quantifying the amount of steam emitted from the vents by direct measurements of steam flow rate. The heat loss from the steam heated mud pools was determined by quantifying the rate of heat loss from the pools by evaporation, convection, and radiation. The steam flow rate into the pools was calculated from the observed heat loss from the pools, assuming that steam flow was the only mechanism of heat transport into the pool. The CO₂ emissions from the steam vents and mud pools were determined by multiplying the steam flow rate from the respective sources by the representative CO₂ concentration of steam in the Reykjanes area. The observed rates of CO₂ emissions through soil, steam vents, and steam heated mud pools amounted to 13.5 ± 1.7, 0.23 ± 0.05, and 0.13 ± 0.03 tons per day, respectively. The heat flow through soil, steam vents, and mud pools was 16.9 ± 1.4, 2.2 ± 0.4, and 1.2  ± 0.1 MW, respectively. Heat loss from the geothermal reservoir, inferred from the CO₂ emissions through the soil amounts to 130 ± 16 MW of thermal energy. The discrepancy between the observed heat loss and the heat loss inferred from the CO₂ emissions is attributed to steam condensation in the subsurface due to interactions with cold ground water. These results demonstrate that soil diffuse degassing can be a more reliable proxy for heat loss from geothermal systems than soil temperatures. The soil diffuse degassing at Reykjanes appears to be strongly controlled by the local tectonics. The observed diffuse degassing defines 3–5 elongated N–S trending zones (000–020°). The orientation of the diffuse degassing structures at Reykjanes is consistent with reported trends of right lateral strike slip faults in the area. The natural CO₂ emissions from Reykjanes under the current low-production conditions are about 16% of the expected emissions from a 100 MWe power plant, which has recently been commissioned at Reykjanes.     

Physicochemical parameters of crystallization of melts in intermediate suprasubduction chambers (by the example of Tolbachik and Ichinskii Volcanoes,Kamchatka Peninsula)

In study of plagioclases, amphiboles, and melt inclusions, we have determined the physicochemical parameters of crystallization of melts in the intermediate suprasubduction chambers of volcanoes representing different types of subduction magmatism on the Kamchatka Peninsula: the young basaltic systems of Tolbachik Volcano (Klyuchevskaya group) and ancient Ichinskii Volcano (Sredinnyi Ridge) with alternating basaltic and felsic eruptions. For Tolbachik Volcano, we have found that plagioclase lapilli formed from basaltic melts at 1075-1115 °C and low (< 1 kbar) pressures at depths of 2-3 km. Andesite minerals crystallized within a wider range of temperatures and pressures (1220-1020 °C and 3.3-1.6 kbar) in an intermediate chamber at depths of < 10 km. The melts were generated in basaltic magma chambers (detected well by geophysical methods at depths of 18-20 km) with minimum temperatures of ~ 1290 °C. For Ichinskii Volcano, three levels of intermediate chambers are distinguished. Andesites formed at depths of < 23 km at < 1225 °C. Dacitic melts were generated from an intermediate chamber (14 km) at 1135-1045 °C as a result of differentiation of andesitic magmas. Dacites formed in the uppermost horizons (9-3 km) at 1130-1030 °C. Despite the similarity between differentiation processes in the intermediate chambers of the Kamchatka volcanoes, each volcano is characterized by specific magmatism. The lavas of basaltic volcanoes (Tolbachik) and those of andesitic volcanoes (Ichinskii) differ in genesis and differentiation.     

Hummocky moraines in the Namaras and Susam Valleys,Central Taurids,SW Turkey

Late Quaternary glacial features have been found in the Central Taurid Mountains, in U-shaped valleys located at an altitude of more than 2000 m and surrounded by mountain ranges reaching 2850 m. No present day active glaciers exist in this area where the snowline elevation lies at about 3500 m. The Namaras Valley and its tributary Susam Valley are characterized by coarse loose material forming chaotic knob-and-kettle topography. Mounds, 1–10 m high and 10–30 m wide, are separated by 5–30 m wide, several meters deep, irregular depressions. The upper surfaces of the mounds are covered by angular to subangular limestone pebbles and blocks and internal sediments show a typical diamicton appearance with pebbles suspended in a muddy to sandy matrix. These chaotic structures are interpreted as hummocky disintegration moraines from former active glaciers. Successive cross-valley morainic ridges, up to 200 m high and several hundreds of meters long, limit the down-valley extension of these hummocks, and are interpreted as ice-marginal moraines. In the tributary Susam Valley, part of the coarse loose material forms a 200–250 m long and 90–120 m wide tongue-shaped structure with successive arcuate ridges and furrows at its down-valley reach. This structure, which is connected upward to a talus slope and perched cirque, ressembles the morphology of a periglacial rockglacier but is interpreted as the disintegration moraine controlled by small periodic retreat and readvance of the last active ice-front in this region.     

Magmatic-vapor expansion and the formation of high-sulfidation gold deposits: Chemical controls on alteration and mineralization

Large bulk-tonnage high-sulfidation gold deposits, such as Yanacocha, Peru, are the surface expression of structurally-controlled lode gold deposits, such as El Indio, Chile. Both formed in active andesite–dacite volcanic terranes. Fluid inclusion, stable isotope and geologic data show that lode deposits formed within 1500 m of the paleo-surface as a consequence of the expansion of low-salinity, low-density magmatic vapor with very limited, if any, groundwater mixing. They are characterized by an initial ‘Sulfate’ Stage of advanced argillic wallrock alteration ± alunite commonly with intense silicification followed by a ‘Sulfide’ Stage — a succession of discrete sulfide–sulfosalt veins that may be ore grade in gold and silver. Fluid inclusions in quartz formed during wallrock alteration have homogenization temperatures between 100 and over 500 °C and preserve a record of a vapor-rich environment.Recent data for El Indio and similar deposits show that at the commencement of the Sulfide Stage, ‘condensation’ of Cu–As–S sulfosalt melts with trace concentrations of Sb, Te, Bi, Ag and Au occurred at > 600 °C following pyrite deposition. Euhedral quartz crystals were simultaneously deposited from the vapor phase during crystallization of the vapor-saturated melt occurs to Fe-tennantite with progressive non-equilibrium fractionation of heavy metals between melt-vapor and solid. Vugs containing a range of sulfides, sulfosalts and gold record the changing composition of the vapor.Published fluid inclusion and mineralogical data are reviewed in the context of geological relationships to establish boundary conditions through which to trace the expansion of magmatic vapor from source to surface and consequent alteration and mineralization. Initially heat loss from the vapor is high resulting in the formation of acid condensate permeating through the wallrock. This Sulfate Stage alteration effectively isolates the expansion of magmatic vapor in subsurface fracture arrays from any external contemporary hydrothermal activity. Subsequent fracturing is localized by the embrittled wallrock to provide high-permeability fracture arrays that constrain vapor expansion with minimization of heat loss. The Sulfide Stage vein sequence is then a consequence of destabilization of metal-vapor species in response to depressurization and decrease in vapor density.The geology, mineralogy, fluid inclusion and stable isotope data and geothermometry for high-sulfidation, bulk-tonnage and lode deposits are quite different from those for epithermal gold–silver deposits such as McLaughlin, California that formed near-surface in groundwater-dominated hydrothermal systems where magmatic fluid has been diluted to less than about 30%. High sulfidation gold deposits are better termed ‘Solfataric Gold Deposits’ to emphasize this distinction. The magmatic-vapor expansion hypothesis also applies to the phenomenology of acidic geothermal systems in active volcanic systems and equivalent magmatic-vapor discharges on the flanks of submarine volcanoes.     

Composition and exhalation flux of gases from mud volcanoes in Taiwan

Many mud volcanoes are distributed along the tectonic sutures in southern Taiwan and can be divided into five zones based on their relative positions in different tectonic domains. Most active mud volcanoes are exhaling methane-dominated gases. Nevertheless, some gases show unusual carbon dioxide-dominated and/or nitrogen-excess compositions. This implies that there are multiple sources for the gas compositions of mud volcanoes in Taiwan.For better understanding the total amount of exhalation gases and its flux, the gas flow and compositions were continuously measured in the interval of two minutes at Chung-lun (CL) bubbling mud pool for a few months. The major compositions of gases exhaling from this site were 75~90% of CO₂ and 5~12% of CH₄. The amount of gases exhaling from the mud pool can be estimated to be about 1.4 ton/year for CH₄ and 28 ton/year for CO₂, respectively. The preliminary results of exhaling gas flux from the major vents of representative active mud volcanoes, yielded an estimated total CH₄ output of the mud volcanoes in Taiwan of ca. 29 ton/year during quiescent period.     

Geochemistry,provenance, and metamorphic evolution of Gabal Samra Neoproterozoic metapelites,Sinai, Egypt

Metapelites are exposed at Wadi Ba’ba, east of Abu Zenima city; represent the northwestern extension of the Fieran-Solaf Metamorphic Complex, Sinai Peninsula, Egypt. The metapelites are characterized by qtz + pl (An_24–28) + bt + grt ± crd ± sil mineral assemblage, indicating upper amphibolite facies with peak metamorphic conditions of 700 °C and pressures of 7 kbar, as determined by conventional geothermobarometeric methods. This resulted in incipient migmatization, forms patches of leucosomes and melanosomes. Geochemical investigation indicates that the precursor sediments of the metapelites had been deposited as immature Fe-rich shales from source materials of dominantly intermediate composition. Source area exhibited weak to moderate chemical weathering in a tectonically active continental marginal basin within a continental-arc system. A strong shallow-dipping foliation, characterizing the metapelites, was folded around an open antiform with sub-horizontal south plunging hinge.Phase equilibria calculations in the KFMASH system indicate that the peak metamorphic conditions formed at 730–750 °C and 6.8–7.9 kbar. This was followed by a retrogression formed at 770–785 °C and 3.9–4.5 kbar. Hence, this implies an isothermal decompression and rapid exhumation of the metapelites from depth (25–29 km) in the lower crustal level at peak conditions, continuous to include shallow to middle crustal level (14–17 km), at overprint retrograde conditions. Subsequent isobaric cooling took place at 720–750 °C and 3.6–4.5 kbar. The resulting isothermal decompression followed by isobaric cooling clockwise P–T path of the metapelites is more likely, in which the high-temperatures attained maximum conditions during isothermal decompression were enhanced by heat flux, due to the presence of an active magmatic arc that formed on top of subducting young lithosphere. This is supported by a moderate geothermal gradient of 27–43 °C/km and dating compatibility of the Sinai granitoids and the metamorphic complexes. The P–T path segment records the tectonothermal histories of crustal thickening as a result of the East and West Gondwana collision at the metamorphic peak. This was subsequent by extensional and crustal thinning with syn-metamorphic magmatic intrusions, during P–T path retrogression, which resulted in the final assembly of the Arabian–Nubian Shield during Neoproterozoic.     

Glaciations and paleoclimate of Mount Erciyes,central Turkey,since the Last Glacial Maximum,inferred from 36Cl cosmogenic dating and glacier modeling

Forty-four boulders from moraines in two glacial valleys of Mount Erciyes (38.53°N, 35.45°E, 3917 m), central Turkey, dated with cosmogenic chlorine-36 (³⁶Cl), indicate four periods of glacial activity in the past 22 ka (1 ka = 1000 calendar years). Last Glacial Maximum (LGM) glaciers were the most extensive, reaching 6 km in length and descending to an altitude of 2150 m above sea level. These glaciers started retreating 21.3 ± 0.9 ka (1σ) ago. They readvanced and retreated by 14.6 ± 1.2 ka ago (Lateglacial), and again by 9.3 ± 0.5 ka ago (Early Holocene). The latest advance took place 3.8 ± 0.4 ka ago (Late Holocene). Using glacier modeling together with paleoclimate proxy data from the region, we reconstructed the paleoclimate at these four discrete times. The results show that LGM climate was 8–11 °C colder than today and moisture levels were somewhat similar to modern values, with a range between 20% more and 25% less than today. The analysis of Lateglacial advance suggests that the climate was colder by 4.5–6.4 °C based on up to 1.5 times wetter conditions. The Early Holocene was 2.1–4.9 °C colder and up to twice as wet as today, while the Late Holocene was 2.4–3 °C colder and its precipitation amounts approached to similar conditions as today. Our paleoclimate reconstructions show a general trend of warming for the last 22 ka, and an increase of moisture until Early Holocene, and a decrease after that time. The recent glacier terminates at 3450 m on the northwest side of the mountain. It is a remnant from the last advance (possibly during the Little Ice Age). Repeated measurements of glacier length between 1902 and 2008 reveal a retreat rate of 4.2 m per year, which corresponds to a warming rate of 0.9–1.2 °C per century.     

Monitoring groundwater flow and chemical and isotopic composition at a demonstration site for carbon dioxide storage in a depleted natural gas reservoir

Between March 2008 and August 2009, 65,445 tonnes of ∼75 mol% CO₂ gas were injected in a depleted natural gas reservoir approximately 2000 m below surface at the Otway project site in Victoria, Australia. Groundwater flow and composition were monitored biannually in two overlying aquifers between June 2006 and March 2011, spanning the pre-, syn- and post-injection periods. The shallower (∼0–100 m), unconfined, porous and karstic aquifer of the Port Campbell Limestone and the deeper (∼600–900 m), confined and porous aquifer of the Dilwyn Formation contain valuable fresh to brackish water resources. Groundwater levels in either aquifer have not been affected by the drilling, pumping and injection activities that were taking place, or by the rainfall increase during the project. In terms of groundwater composition, the Port Campbell Limestone groundwater is brackish (electrical conductivity = 801–3900 μS cm⁻¹), cool (temperature = 12.9–22.5 °C), and near-neutral (pH = 6.62–7.45), whilst the Dilwyn Aquifer groundwater is fresher (electrical conductivity = 505–1473 μS cm⁻¹), warmer (temperature = 42.5–48.5 °C), and more alkaline (pH = 7.43–9.35). Carbonate dissolution, evapotranspiration and cation exchange control the composition of the groundwaters. Comparing the chemical and isotopic composition of the groundwaters collected before, during and after injection shows no statistically significant changes; even if they were statistically significant, they are mostly not consistent with those expected if CO₂ addition had taken place. The monitoring program reveals no impact on the groundwater resources attributable to the C storage demonstration project.     

Past glaciation in the tropics

Tropical glaciers are considered along three meridional profiles, for the Australasian sector, Africa and the Americas. Evaluated are the annual mean freezing level 0 °C, modern equilibrium line altitude MEL, and past equilibrium line altitude PEL. The calculation of 0 °C is based on a 1958–1997 global data set; MEL refers to estimates concerning the first half of the 20th century; and the timing of the PEL is not generally known. 0 °C stands around 4000–5000 m, with lower levels in the outer tropics. MEL is reached in the Australasian sector on four mountains, and in Africa on three mountains, near the Equator, near or above 0 °C. In the American cordilleras many peaks are still glaciated, above 0 °C, but in the arid southern tropical Andes even summits above 6000 m do not reach MEL. The PEL stands between 3000 m and 5000 m, high in the equatorial zone, but highest in the arid southern tropical Andes. The height difference between MEL and PEL is of order 1000 m, with regional differences. Deglaciation dates range between 15,000 and 8000 years BP, with later timing towards the higher elevations. This synopsis suggests priority targets for further research: morphological mapping and age determination in the High Atlas of Morocco; timing of deglaciation in the High Semyen of Ethiopia and the Altos de Cuchumatanes of Guatemala; exploration in the highlands of Lesotho and of the Dominican Republic; and glaciation in the arid southern tropical Andes.     

Biogenic components of the Baltic Sea sediments

The contents of biogenic components in 1511 samples of the Baltic Sea sediments (depth range 0–5 cm) are studied, and maps of their distribution are compiled. The sediments contain < 13.03% C_org, < 1.33% N, < 9.0% SiO_2am, < 5.0% CaCO₃, and < 1.45% P. The maximum and elevated contents of components are found in the mud of the sea deeps. The more fraction < 0.01 mm the sediments contain, the higher are the contents of components. Four facies types of carbonaceous mud, precursors of shales, have been recognized: (1) shallow-water (lagoon) lime sapropel, (2) carbonaceous mud of the shallow-water Gulf of Finland, (3) carbonaceous mud of the middle-depth Baltic Sea, and (4) laminated carbonaceous metal-bearing mud. The latter type of mud is strongly enriched in manganese and ore-forming trace elements, which points to its formation in the stagnant environment. In composition the Baltic Sea mud is similar to petroliferous mudstones of the Bazhenov Formation in West Siberia and to ancient black shales.     

Source characteristics of the Yutian earthquake in 2008 from inversion of the co-seismic deformation field mapped by InSAR

On 21 March 2008, an Ms7.3 earthquake occurred at Yutian County, Xinjiang Uygur Autonomous Region, which is in the same year as 2008 Mw 7.9 Wenchuan earthquake. These two earthquakes both took place in the Bayar Har block, while Yutian earthquake is located in the west edge and Wenchuan earthquake is in the east. The research on source characteristics of Yutian earthquake can serve to better understand Wenchuan earthquake mechanism. We attempt to reveal the features of the causative fault of Yutian shock and its co-seismic deformation field by a sensitivity-based iterative fitting (SBIF) method. Our work is based on analysis and interpretation to high-resolution satellite (Quickbird) images as well as D-InSAR data from the satellite Envisat ASAR, in conjunction with the analysis of seismicity, focal mechanism solutions and active tectonics in this region. The result shows that the 22 km long, nearly NS trending surface rupture zone by this event lies on a range-front alluvial platform in the Qira County. It is characterized by distinct linear traces and a simple structure with 1–3 m-wide individual seams and maximum 6.5 m width of a collapse fracture. Along the rupture zone are seen many secondary fractures and fault-bounded blocks by collapse, exhibiting remarkable extension. The co-seismic deformation affected a big range 100 km × 40 km. D-InSAR analysis indicates that the interferometric deformation field is dominated by extensional faulting with a small strike-slip component. Along the causative fault, the western wall fell down and the eastern wall, that is the active unit, rose up, both with westerly vergence. The maximum subsidence displacement is ～2.6 m in the LOS, and the maximum uplift is 1.2 m. The maximum relative vertical dislocation reaches 4.1 m, which is 10 km distant from the starting rupture point to south. The 42 km-long seismogenic fault in the subsurface extends in NS direction as an arc, and it dipping angle changes from 70° near the surface to 52° at depth ～10 km. The slip on the fault plane is concentrated in the depth range 0–8 km, forming a belt of length 30 km along strike on the fault plane. There are three areas of concentrating slip, in which the largest slip is 10.5 m located at the area 10 km distant from the initial point of the rupture.     

Paleomagnetic analysis of the Marwar Supergroup,Rajasthan, India and proposed interbasinal correlations

The Marwar Supergroup refers to a 1000–2000 m thick marine and coastal sequence that covers a vast area of Rajasthan in NW–India. The Marwar Basin uncomformably overlies the ∼750–770 Ma rocks of the Malani Igneous Suite and is therefore considered Late Neoproterozoic to Early Cambrian in age. Upper Vindhyan basinal sediments (Bhander and Rewa Groups), exposed in the east and separated by the Aravalli–Delhi Fold Belt, have long been assumed to coeval with the Marwar Supergroup. Recent studies based on detrital zircon populations of the Marwar and Upper Vindhyan sequences show some similarity in the older populations, but the Vindhyan sequence shows no zircons younger than 1000 Ma whereas samples taken from the Marwar Basin show distinctly younger zircons. This observation led to speculation that the Upper Vindhyan and Marwar sequences did not develop coevally.While there are alternative explanations for why the two basins may differ in their detrital zircon populations, paleomagnetic studies may provide independent evidence for differences/similarities between the assumed coeval basins. We have collected samples in the Marwar Basin and present the paleomagnetic results. Previous paleomagnetic studies of Marwar basinal sediments were misinterpreted as being indistinguishable from the Upper Vindhyan sequence. The vast majority of our samples show directional characteristics similar to the previously published studies. We interpret these results to be a recent overprint. A small subset of hematite-bearing rocks from the Jodhpur Formation (basal Marwar) exhibit directional data (Dec = 89° Inc = −1° α95 = 9°) that are distinct from the Upper Vindhyan pole and may offer additional support for temporally distinct episodes of sedimentation in these proximal regions. A VGP based upon our directional data is reported at 1°S 344°E (dp = 5°, dm = 9°). We conclude that the Marwar Supergroup developed near the close of the Ediacaran Period and is part of a larger group of sedimentary basins that include the Huqf Supergroup (Oman), the Salt-Range (Pakistan), the Krol–Tal belt (Himalayas) and perhaps the Molo Supergroup (Madagascar).     

Quaternary hinterland evolution of the active Banda Arc: Surface uplift and neotectonic deformation recorded by coral terraces at Kisar,Indonesia

Coral terrace surveys and U-series ages of coral yield a surface uplift rate of ∼0.5 m/ka for Kisar Island, which is an emergent island in the hinterland of the active Banda arc–continent collision. Based on this rate, Kisar first emerged from the ocean as recently as ∼450 ka. These uplifted terraces are gently warped in a pattern of east–west striking folds. These folds are strike parallel to more developed thrust-related folds of similar wavelength imaged by a seismic reflection profile just offshore. This deformation shows that the emergence of Kisar is influenced by forearc closure along the south-dipping Kisar Thrust. However, the pinnacle shape of Kisar and the protrusion of its metamorphic rocks through the forearc basin sediments also suggest a component of extrusion along shear zones or active doming.Coral encrusts the island coast in many locations over 100 m above sea level. Terrace morphology and coral ages are best explained by recognizing major surfaces as mostly growth terraces and minor terraces as mostly erosional into older terraces. All reliable and referable coral U-series ages determined by MC-ICP-MS correlate with marine isotope stage (MIS) 5e (118–128 ka). The only unaltered coral samples are found below 6 m elevation; however an unaltered Tridacna (giant clam) shell in growth position at 95 m elevation yields a U-series age of 195 ± 31 ka, which corresponds to MIS 7. This age agrees with the best-fit uplift model for the island. Loose deposits of unaltered coral fragments found at elevations between 8 and 20 m yield U-series ages of <100 years and may represent paleotsunami deposits from previously undocumented tectonic activity in the region.     

Duration of high-pressure metamorphism and cooling during the intraplate Petermann Orogeny

The Ediacaran–Cambrian Petermann Orogeny, central Australia, is an exceptional example of intraplate orogenesis. It involved sub-eclogite facies metamorphism and extreme basin inversion during the exhumation of Musgrave Province basement from beneath the formerly contiguous Centralian Superbasin. Sensitive High Resolution Ion Microprobe (SHRIMP) U–Pb geochronology of zircon, titanite and rutile, along with Ti-in-zircon thermometry from meta-igneous samples, have been used to determine the timing and duration of high-pressure metamorphism and subsequent cooling associated with this orogenic event. Peak metamorphic temperatures of 720–760 °C were attained at 544 ± 7 Ma (U–Pb zircon). Subsequent cooling to 600–660 °C by ~ 521 Ma occurred at a rate of ~ 2.6–7.0 °C Myr^− 1, as recorded by the closure of Pb diffusion in titanite. Further cooling to 585–560 °C by 498–472 Ma occurred at a rate of 0.9–4.8 °C Myr^− 1, as recorded by Pb closure in rutile. The duration of tectonism was long-lived (> 40 Myr) across the central and western parts of the orogenic system, and deformation occurred in a comparatively warm and weak portion of crust, characterised by regional thermal gradients of 17–26 °C km^− 1. This proposed duration of tectonism is much longer than that permitted by a shear heating mechanism, which requires an exceptionally short duration of tectonism, and additionally, an overall cold lithosphere characterised by geothermal gradients of ~ 9 °C km^− 1.     

Tectonics and thermal structure of western Satpura,India

Increased seismicity and occurrences of hot springs having surface temperature of 36–58 °C are observed in the central part of India (74–81° E, 20–25° N), where the NE trending Middle Proterozoic Aravalli Mobile Belt meets the ENE trending Satpura Mobile Belt. Earlier Deep Seismic Sounding (DSS) studies along Thuadara-Sendhwa-Sindad profile in the area has showed Mesozoic Sediments up to around 4 km depth covered by Deccan Trap and the Moho depth with a boundary velocity (P_n) of 8.2 km/s. In the present study, surface heat flow of 48 ± 4 mW m^?2 has been estimated based on P_n velocity, which agrees with the value of heat flow of 52 ± 4 mW m^?2 based on Curie point isotherms estimates. The calculated temperature-depth profile shows temperature of 80–120 °C at the basement, which is equivalent to oil window temperature in Mesozoic sediments and around 570–635 °C at Moho depth of 38–43 km and the thermal lithosphere is about 110 km thick, which is comparatively higher than those of adjoining regions. The present study reveals the brittle–ductile transition zone at 14–41 km depth (temperature around 250–600 °C) where earthquake nucleation takes place.     

Ore petrology,hydrothermal alteration,fluid inclusions,and sulfur stable isotopes of the Milin Kamak intermediate sulfidation epithermal Au-Ag deposit in Western Srednogorie,Bulgaria

The Milin Kamak gold-silver deposit is located in Western Srednogorie zone, 50 km west of Sofia, Bulgaria. This zone belongs to the Late Cretaceous Apuseni-Banat-Timok-Srednogorie magmatic and metallogenic belt. The deposit is hosted by altered trachybasalt to andesitic trachybasalt volcanic and volcanoclastic rocks with Upper Cretaceous age, which are considered to be products of the Breznik paleovolcano. Milin Kamak is the first gold-silver intermediate sulfidation type epithermal deposit recognized in Srednogorie zone in Bulgaria. It consists of eight ore zones with lengths ranging from 400 to 1000 m, widths from several cm to 3–4 m, rarely to 10–15 m, an average of 80–90 m depth (a maximum of 200 m) and dip steeply to the south. The average content of gold is 5.04 g/t and silver – 13.01 g/t. The styles of alteration are propylitic, sericite, argillic, and advanced argillic. Ore mineralization consists of three stages. Quartz-pyrite stage I is dominated by quartz, euhedral to subhedral pyrite, trace pyrrhotite and hematite in the upper levels of the deposit. Quartz-polymetallic stage II is represented by major anhedral pyrite, galena, Fe-poor sphalerite; minor chalcopyrite, tennantite, bournonite, tellurides and electrum; and trace pyrrhotite, arsenopyrite, marcasite. Gangue minerals are quartz and carbonates. The carbonate-gold stage III is defined by deposition of carbonate minerals and barite with native gold and stibnite.Fluid inclusions in quartz are liquid H₂O-rich with homogenization temperature (T_h) ranging from 238 to 345 °C as the majority of the measurements are in the range 238–273 °C. Ice-melting temperatures (T_m) range from −2.2 to −4.1 °C, salinity – from 3.7 to 6.6 wt.% NaCl equiv. These measurements imply an epithermal environment and low- to moderate salinity of the ore-forming fluids.δ³⁴S values of pyrite range from −0.49 to +2.44‰. The average calculated δ³⁴S values are 1.35‰. The total range of δ³⁴S values for pyrite are close to zero suggesting a magmatic source for the sulfur.     

First documented deep submarine explosive eruptions at the Marsili Seamount (Tyrrhenian Sea,Italy): A case of historical volcanism in the Mediterranean Sea

The Marsili Seamount (MS) is an about 3200 m high volcanic complex measuring 70 × 30 km with the top at ~ 500 m b.s.l. MS is interpreted as the ridge of the 2 Ma old Marsili back-arc basin belonging to the Calabrian Arc–Ionian Sea subduction system (Southern Tyrrhenian Sea, Italy). Previous studies indicate that the MS activity developed between 1 and 0.1 Ma through effusions of lava flows. Here, new stratigraphic, textural, geochemical, and ¹⁴C geochronological data from a 95 cm long gravity core (COR02) recovered at 839 m bsl in the MS central sector are presented. COR02 contains mud and two tephras consisting of 98 to 100 area% of volcanic ash. The thickness of the upper tephra (TEPH01) is 15 cm, and that of the lower tephra (TEPH02) is 60 cm. The tephras have poor to moderate sorting, loose to partly welded levels, and erosive contacts, which imply a short distance source of the pyroclastics. ¹⁴C dating on fossils above and below TEPH01 gives an age of 3 ka BP. Calculations of the sedimentation rates from the mud sediments above and between the tephras suggest that a formation of TEPH02 at 5 ka BP MS ashes has a high-K calcalkaline affinity with 53 wt.% < SiO₂ < 68 wt.%, and their composition overlaps that of the MS lava flows. The trace element pattern is consistent with fractional crystallization from a common, OIB-like basalt. The source area of ashes is the central sector of MS and not a subaerial volcano of the Campanian and/or Aeolian Quaternary volcanic districts. Submarine, explosive eruptions occurred at MS in historical times: this is the first evidence of explosive volcanic activity at a significant (500–800 m bsl) water depth in the Mediterranean Sea. MS is still active, the monitoring and an evaluation of the different types of hazards are highly recommended.     

Evolution of the Mo-rich scheelite skarn mineralization at Kozbudaklar,Western Anatolia,Turkey: Evidence from mineral chemistry and fluid inclusions

The Kozbudaklar scheelite skarn deposit in the Tavşanlı Zone, located approximately 22 km southeast of Bursa, is hosted by the Triassic calcic İnönü Marble and Eocene Topuk Pluton. At least four stages have been recognized through skarn evolution. Scheelite skarn distributed close to the Topuk Pluton occurred during the early (stage 1) and late (stage 2) prograde substages. The early prograde endo and exoskarn are composed of hedenbergite (Hd₉₆Joh₄)–plagioclase (An_55–64) and hedenbergite (Hd_61–94Joh_4–7), accompanied by calcic garnet (Grs_38–94Sps_1–5Alm₀) and scheelite (Pow_1–6). The second stage represents a relatively oxidized mineralogy dominated by diopside (Hd_16–48Joh_0–9), subcalcic garnet (Grs_24–92Sps_0–11Alm_0–31) and scheelite (Pow_7–32). The stage 3 and 4 mineral assemblages are characterized by few hydrous minerals in the retrograde stage and intense fracturing.Fluid inclusions from skarn rocks are indicative of multiple fluid events: (1) low-moderate salinity (5–16 wt.%NaCl equiv.) inclusions homogenized dominantly by a high-temperature (308 °C to > 600 °C) liquid phase in stage 1. Fluid inclusions in an early garnet homogenized over a similar temperature range (440 °C and 459 °C) into both liquid and vapor phases. Eutectic temperatures ranging from − 61.7 °C to − 35.0 °C that indicate the presence of H₂O–NaCl–(± MgCl₂ ± CaCl₂)–CO₂ solutions; (2) coexisting daughter mineral-bearing high salinity (29.5 ≥ 70 wt.%NaCl equiv.) and vapor-rich moderate salinity (11.5–16.7 wt.%NaCl equiv.) inclusions that homogenized in the liquid phase by the disappearance of the vapor phase at a similar temperature range (308 °C to > 600 °C) in stage 2. Eutectic temperatures range from − 67.9°C to − 51.8°C that shows the presence of H₂O–NaCl–CO₂–(± CH₄/N₂) solutions; (3) low-moderate salinity (12.5–7.6 wt.%NaCl equiv.) and temperature (320 °C to 215 °C) inclusions homogenized by the liquid phase in stage 3. Eutectic temperatures range from − 59.5 °C to − 44.2 °C indicating the presence of H₂O–NaCl–(± MgCl₂ ± CaCl₂)–CO₂ solutions; (4) inclusions of low salinity (9.9–0.9 wt.%NaCl equiv.) and homogenization temperature (183 °C to 101 °C) in stage 4.These data show that the Kozbudaklar skarn deposit was formed in a magmatic–hydrothermal system. In this model, carbonaceous fluids may have been exsolved from the plutonic rock during its emplacement and crystallization. Fluid inclusion data indicate that fluid boiling and immiscibility occurred at temperatures between 440 °C and 459 °C and pressures ranging from 50 MPa to 60 MPa based on hydrostatic considerations. Early scheelite was precipitated with relatively reduced mineral compositions. As a result of depressurization, Mo-rich scheelite with oxidized minerals formed via high salinity and vapor-rich inclusions. The second scheelite mineralization occurred in a normal hydrothermal system by an infiltration mechanism at pressures between approximately 40 and 1.5 MPa. At shallow depths (< 1.5 MPa) with increasing permeability, sulfide and oxide minerals were deposited in the retrograde stage, greatly assisted by meteoric water. Finally, as a result of the diminishing of ore-forming fluids, post-depositional barren quartz and calcite veins were formed.     

The exhumation history of collision-related mineralizing systems in Tibet: Insights from thermal studies of the Sharang and Yaguila deposits,central Lhasa

The large, newly discovered Sharang porphyry Mo deposit and nearby Yaguila skarn Pb–Zn–Ag (–Mo) deposit reside in the central Lhasa terrane, northern Gangdese metallogenic belt, Tibet. Multiple mineral chronometers (zircon U–Pb, sericite ⁴⁰Ar–³⁹Ar, and zircon and apatite (U–Th)/He) reveal that ore-forming porphyritic intrusions experienced rapid cooling (> 100 °C/Ma) during a monotonic magmatic–hydrothermal evolution. The magmatic–hydrothermal ore-forming event at Sharang lasted ~ 6.0 Myr (~ 1.8 Myr for cooling from > 900 to 350 °C and ~ 4.0 Myr for cooling from 350 to 200 °C) whereas cooling was more prolonged during ore formation at Yaguila (~ 1.8 Myr from > 900 to 500 °C and a maximum of ~ 16 Myr from > 900 to 350 °C). All porphyritic intrusions in the ore district experienced exhumation at a rate of 0.07–0.09 mm/yr (apatite He ages between ~ 37 and 30 Ma). Combined with previous studies, this work implies that uplift of the eastern section of the Lhasa terrane expanded from central Lhasa (37–30 Ma) to southern Lhasa (15–12 Ma) at an increasing exhumation rate. All available geochronologic data reveal that magmatic–hydrothermal–exhumation activities in the Sharang–Yaguila ore district occurred within four periods of magmatism with related mineralization. Significant porphyry-type Mo mineralization was associated with Late Cretaceous–Eocene felsic porphyritic intrusions in the central Lhasa terrane, resulting from Neotethyan oceanic subduction and India–Asia continental collision.     

3D thermal structure of the continental lithosphere beneath China and adjacent regions

Based on the Crust2.0 model and the topography data of Chinese continent and its adjacent regions, a three-dimensional finite element model is constructed in terms of the spherical coordinate system. In our numerical model, the average annual ground temperature from 195 meteorological stations and temperature of upper mantle derived from the seismic velocities are adopted as the top and bottom boundary conditions, respectively. The observed thermal conductivity and heat production from P wave velocity based on empirical formula are employed in our numerical model as well. The comparison between the calculated and observed surface heat flow proved that our results are reliable. The temperature beneath the Precambrian cratons is lower than that of other areas for 100–300 °C also. The typical temperature rang at the Moho is estimated to be 800–1000 °C beneath the Tibetan plateau and 500–700 °C beneath the Precambrian cratons (such as Indian plate, Sichuan basin, South China, North China and Tarim), respectively. The thermal state in the eastern part of Sino-Korean craton at the depth deeper than 60 km indicates that it was destructed. The thermal structure in center of Tibetan plateau (especially beneath Qiangtang area) supports the proposed flow of lower crustal or upper mantle material to the east. Generally, the distribution of volcanoes in Chinese continent is consistent with the high temperature areas in the crust or upper mantle. There are many obvious thermal transition zones across the orogenic belts. The thermal transition zone between eastern and western parts in the crust of Chinese continent is consistent with the north–south seismic zone.